Changeset 24431 for issm/trunk-jpl/src/c/analyses/ExtrapolationAnalysis.cpp

Timestamp:

12/04/19 21:11:16 (5 years ago)

Author:

Mathieu Morlighem

Message:

CHG: trying to speed up the code

File:

• issm/trunk-jpl/src/c/analyses/ExtrapolationAnalysis.cpp (modified) (2 diffs)

issm/trunk-jpl/src/c/analyses/ExtrapolationAnalysis.cpp

@@ -207 +207 @@
 ElementVector* ExtrapolationAnalysis::CreatePVector(Element* element){/*{{{*/
         return NULL;
-
-}/*}}}*/
-void           ExtrapolationAnalysis::GetB(IssmDouble* B,Element* element,IssmDouble* xyz_list,Gauss* gauss, int dim){/*{{{*/
-        /*Compute B  matrix. B=[B1 B2 B3] where Bi is of size 3*NDOF2. 
-         * For node i, Bi can be expressed in the actual coordinate system
-         * by: 
-         *       Bi=[ N ]
-         *          [ N ]
-         * where N is the finiteelement function for node i.
-         *
-         * We assume B_prog has been allocated already, of size: 2x(NDOF1*numnodes)
-         */
-
-        /*Fetch number of nodes for this finite element*/
-        int numnodes = element->GetNumberOfNodes();
-
-        /*Get nodal functions*/
-        IssmDouble* basis=xNew<IssmDouble>(numnodes);
-        element->NodalFunctions(basis,gauss);
-
-        /*Build B: */
-        for(int i=0;i<numnodes;i++)
-                for(int d=0;d<dim;d++)
-                        B[numnodes*d+i] = basis[i];
-
-        /*Clean-up*/
-        xDelete<IssmDouble>(basis);
-}/*}}}*/
-void           ExtrapolationAnalysis::GetBprime(IssmDouble* Bprime,Element* element,IssmDouble* xyz_list,Gauss* gauss, int dim){/*{{{*/
-        /*Compute B'  matrix. B'=[B1' B2' B3'] where Bi' is of size 3*NDOF2. 
-         * For node i, Bi' can be expressed in the actual coordinate system
-         * by: 
-         *       Bi_prime=[ dN/dx ]
-         *                [ dN/dy ]
-         * where N is the finiteelement function for node i.
-         *
-         * We assume B' has been allocated already, of size: 3x(NDOF2*numnodes)
-         */
-
-        /*Fetch number of nodes for this finite element*/
-        int numnodes = element->GetNumberOfNodes();
-
-        /*Get nodal functions derivatives*/
-        IssmDouble* dbasis=xNew<IssmDouble>(dim*numnodes);
-        element->NodalFunctionsDerivatives(dbasis,xyz_list,gauss);
-
-        /*Build B': */
-        for(int i=0;i<numnodes;i++)
-                for(int d=0;d<dim;d++)
-                        Bprime[numnodes*d+i] = dbasis[d*numnodes+i];
-
-        /*Clean-up*/
-        xDelete<IssmDouble>(dbasis);
-
 }/*}}}*/
 void           ExtrapolationAnalysis::GetSolutionFromInputs(Vector<IssmDouble>* solution,Element* element){/*{{{*/
@@ -288 +234 @@
         }
 }/*}}}*/
-int                             ExtrapolationAnalysis::GetExtrapolationCase(Element* element){/*{{{*/
+int            ExtrapolationAnalysis::GetExtrapolationCase(Element* element){/*{{{*/
 
         /* Get case of extrapolation, depending on domain quality, and extrapolation variable */