Changeset 22535

Timestamp:

03/13/18 07:06:32 (7 years ago)

Author:

Mathieu Morlighem

Message:

CHG: removing GetB and GetBprime

Location:

issm/trunk-jpl/src/c/analyses

Files:

• StressbalanceVerticalAnalysis.cpp (modified) (1 diff)
• StressbalanceVerticalAnalysis.h (modified) (1 diff)

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

@@ -479 +479 @@
         return pe;
 }/*}}}*/
-void           StressbalanceVerticalAnalysis::GetB(IssmDouble* B,Element* element,IssmDouble* xyz_list,Gauss* gauss){/*{{{*/
-        /*      Compute B  matrix. B=[dh1/dz dh2/dz dh3/dz dh4/dz dh5/dz dh6/dz];
-                where hi is the interpolation function for node i.*/
-
-        /*Fetch number of nodes for this finite element*/
-        int numnodes = element->GetNumberOfNodes();
-
-        /*Get nodal functions derivatives*/
-        IssmDouble* dbasis=xNew<IssmDouble>(3*numnodes);
-        element->NodalFunctionsDerivatives(dbasis,xyz_list,gauss);
-
-        /*Build B: */
-        for(int i=0;i<numnodes;i++){
-                B[i] = dbasis[2*numnodes+i];  
-        }
-
-        /*Clean-up*/
-        xDelete<IssmDouble>(dbasis);
-}/*}}}*/
-void           StressbalanceVerticalAnalysis::GetBprime(IssmDouble* Bprime,Element* element,IssmDouble* xyz_list,Gauss* gauss){/*{{{*/
-
-        element->NodalFunctions(Bprime,gauss);
-
-}/*}}}*/
 void           StressbalanceVerticalAnalysis::GetSolutionFromInputs(Vector<IssmDouble>* solution,Element* element){/*{{{*/
         element->GetSolutionFromInputsOneDof(solution,VzEnum);

issm/trunk-jpl/src/c/analyses/StressbalanceVerticalAnalysis.h

@@ -32 +32 @@
                 ElementVector* CreatePVectorSurface(Element* element);
                 ElementVector* CreatePVectorVolume(Element* element);
-                void           GetB(IssmDouble* B,Element* element,IssmDouble* xyz_list,Gauss* gauss);
-                void           GetBprime(IssmDouble* B,Element* element,IssmDouble* xyz_list,Gauss* gauss);
                 void           GetSolutionFromInputs(Vector<IssmDouble>* solution,Element* element);
                 void           GradientJ(Vector<IssmDouble>* gradient,Element* element,int control_type,int control_index);