Changeset 17926 for issm/trunk-jpl/src/c/classes/Elements/SegRef.cpp
- Timestamp:
- 05/04/14 21:04:34 (11 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
TabularUnified issm/trunk-jpl/src/c/classes/Elements/SegRef.cpp ΒΆ
r17874 r17926 49 49 50 50 /*Reference Element numerics*/ 51 /*FUNCTION SegRef::GetBprimeMasstransport{{{*/52 void SegRef::GetBprimeMasstransport(IssmDouble* Bprime, IssmDouble* xyz_list, GaussSeg* gauss){53 /*Compute B' matrix. B'=[B1' B2' B3'] where Bi' is of size 3*NDOF2.54 * For node i, Bi' can be expressed in the actual coordinate system55 * by:56 * Bi_prime=[ dN/dx ]57 * where N is the finiteelement function for node i.58 *59 * We assume B' has been allocated already, of size: 3x(NDOF2*numnodes)60 */61 62 /*Fetch number of nodes for this finite element*/63 int numnodes = this->NumberofNodes();64 65 /*Get nodal functions derivatives*/66 IssmDouble* dbasis=xNew<IssmDouble>(2*numnodes);67 GetNodalFunctionsDerivatives(dbasis,xyz_list,gauss);68 69 /*Build B': */70 for(int i=0;i<numnodes;i++){71 Bprime[i] = dbasis[i];72 }73 74 /*Clean-up*/75 xDelete<IssmDouble>(dbasis);76 }77 /*}}}*/78 51 /*FUNCTION SegRef::GetNodalFunctions(IssmDouble* basis,GaussSeg* gauss){{{*/ 79 52 void SegRef::GetNodalFunctions(IssmDouble* basis,GaussSeg* gauss){
Note:
See TracChangeset
for help on using the changeset viewer.