SegRef Class Reference

#include <SegRef.h>

Inheritance diagram for SegRef:

Public Member Functions
	SegRef ()
	~SegRef ()
void	GetInputDerivativeValue (IssmDouble *p, IssmDouble *plist, IssmDouble *xyz_list, GaussSeg *gauss, int finiteelement)
void	GetInputValue (IssmDouble *p, IssmDouble *plist, GaussSeg *gauss, int finiteelement)
void	GetJacobian (IssmDouble *J, IssmDouble *xyz_list, GaussSeg *gauss)
void	GetJacobianDeterminant (IssmDouble *Jdet, IssmDouble *xyz_list, GaussSeg *gauss)
void	GetJacobianInvert (IssmDouble *Jinv, IssmDouble *xyz_list, GaussSeg *gauss)
void	GetNodalFunctions (IssmDouble *basis, GaussSeg *gauss, int finiteelement)
void	GetNodalFunctionsDerivatives (IssmDouble *dbasis, IssmDouble *xyz_list, GaussSeg *gauss, int finiteelement)
void	GetNodalFunctionsDerivativesReference (IssmDouble *dbasis, GaussSeg *gauss, int finiteelement)
void	Marshall (char **pmarshalled_data, int *pmarshalled_data_size, int marshall_direction)
int	NumberofNodes (int finiteelement)

Detailed Description

Definition at line 12 of file SegRef.h.

Constructor & Destructor Documentation

SegRef()

SegRef::SegRef

(

)

Definition at line 23 of file SegRef.cpp.

              {/*{{{*/
}

~SegRef()

SegRef::~SegRef

(

)

Definition at line 26 of file SegRef.cpp.

               {/*{{{*/
}

Member Function Documentation

GetInputDerivativeValue()

void SegRef::GetInputDerivativeValue	(	IssmDouble *	p,
		IssmDouble *	plist,
		IssmDouble *	xyz_list,
		GaussSeg *	gauss,
		int	finiteelement
	)

Definition at line 31 of file SegRef.cpp.

                                                                                                                            {/*{{{*/
 
   /*From node values of parameter p (plist[0],plist[1]), return parameter derivative value at gaussian 
    * point specified by gauss_basis:
    *   dp/dx=plist[0]*dh1/dx+plist[1]*dh2/dx
    *
    * p is a vector already allocated.
    *
    * WARNING: For a significant gain in performance, it is better to use
    * static memory allocation instead of dynamic.
    */
 
   /*Allocate derivatives of basis functions*/
   IssmDouble  dbasis[1*NUMNODESMAX];
 
   /*Fetch number of nodes for this finite element*/
   int numnodes = this->NumberofNodes(finiteelement);
   _assert_(numnodes<=NUMNODESMAX);
 
   /*Get basis functions derivatives at this point*/
   GetNodalFunctionsDerivatives(&dbasis[0],xyz_list,gauss,finiteelement);
 
   /*Calculate parameter for this Gauss point*/
   IssmDouble dpx=0.;
   for(int i=0;i<numnodes;i++) dpx += dbasis[0*numnodes+i]*plist[i];
 
   /*Assign values*/
   p[0]=dpx;
}

GetInputValue()

void SegRef::GetInputValue	(	IssmDouble *	p,
		IssmDouble *	plist,
		GaussSeg *	gauss,
		int	finiteelement
	)

Definition at line 61 of file SegRef.cpp.

                                                                                             {/*{{{*/
   /* WARNING: For a significant gain in performance, it is better to use
    * static memory allocation instead of dynamic.*/
 
   /*Allocate basis functions*/
   IssmDouble  basis[NUMNODESMAX];
 
   /*Fetch number of nodes for this finite element*/
   int numnodes = this->NumberofNodes(finiteelement);
   _assert_(numnodes<=NUMNODESMAX);
 
   /*Get basis functions at this point*/
   GetNodalFunctions(&basis[0],gauss,finiteelement);
 
   /*Calculate parameter for this Gauss point*/
   IssmDouble value =0.;
   for(int i=0;i<numnodes;i++) value += basis[i]*plist[i];
 
   /*Assign output pointer*/
   *p = value;
}

GetJacobian()

void SegRef::GetJacobian	(	IssmDouble *	J,
		IssmDouble *	xyz_list,
		GaussSeg *	gauss
	)

Definition at line 83 of file SegRef.cpp.

                                                                           {/*{{{*/
   /*The Jacobian is constant over the element, discard the gaussian points. 
    * J is assumed to have been allocated of size 1*/
 
   IssmDouble x1=xyz_list[3*0+0];
   IssmDouble x2=xyz_list[3*1+0];
 
   *J=.5*fabs(x2-x1);
}

GetJacobianDeterminant()

void SegRef::GetJacobianDeterminant	(	IssmDouble *	Jdet,
		IssmDouble *	xyz_list,
		GaussSeg *	gauss
	)

Definition at line 93 of file SegRef.cpp.

                                                                                          {/*{{{*/
   /*The Jacobian determinant is constant over the element, discard the gaussian points. 
    * J is assumed to have been allocated of size 1.*/
 
   /*Call Jacobian routine to get the jacobian:*/
   GetJacobian(Jdet, xyz_list, gauss);
   if(*Jdet<0) _error_("negative jacobian determinant!");
 
}

GetJacobianInvert()

void SegRef::GetJacobianInvert	(	IssmDouble *	Jinv,
		IssmDouble *	xyz_list,
		GaussSeg *	gauss
	)

Definition at line 103 of file SegRef.cpp.

                                                                                    {/*{{{*/
 
   /*Jacobian*/
   IssmDouble J;
 
   /*Call Jacobian routine to get the jacobian:*/
   GetJacobian(&J, xyz_list, gauss);
 
   /*Invert Jacobian matrix: */
   *Jinv = 1./J;
}

GetNodalFunctions()

void SegRef::GetNodalFunctions	(	IssmDouble *	basis,
		GaussSeg *	gauss,
		int	finiteelement
	)

Definition at line 115 of file SegRef.cpp.

                                                                                 {/*{{{*/
   /*This routine returns the values of the nodal functions  at the gaussian point.*/
 
   _assert_(basis);
 
   switch(finiteelement){
      case P0Enum:
         basis[0]=1.;
         return;
      case P1Enum: case P1DGEnum:
         basis[0]=(1.-gauss->coord1)/2.;
         basis[1]=(1.+gauss->coord1)/2.;
         return;
      case P2Enum:
         basis[0]=(gauss->coord1-1.)*gauss->coord1/2.;
         basis[1]=gauss->coord1*(1.+gauss->coord1)/2.;
         basis[2]=(1.-gauss->coord1)*(1.+gauss->coord1);
         return;
      default:
         _error_("Element type "<<EnumToStringx(finiteelement)<<" not supported yet");
   }
}

GetNodalFunctionsDerivatives()

void SegRef::GetNodalFunctionsDerivatives	(	IssmDouble *	dbasis,
		IssmDouble *	xyz_list,
		GaussSeg *	gauss,
		int	finiteelement
	)

Definition at line 138 of file SegRef.cpp.

                                                                                                                   {/*{{{*/
 
   /*This routine returns the values of the nodal functions derivatives  (with respect to the 
    * actual coordinate system): */
   IssmDouble    Jinv;
 
   /*Fetch number of nodes for this finite element*/
   int numnodes = this->NumberofNodes(finiteelement);
 
   /*Get nodal functions derivatives in reference triangle*/
   IssmDouble dbasis_ref[1*NUMNODESMAX];
   GetNodalFunctionsDerivativesReference(dbasis_ref,gauss,finiteelement); 
 
   /*Get Jacobian invert: */
   GetJacobianInvert(&Jinv, xyz_list, gauss);
 
   /*Build dbasis: 
    * [dhi/dx]= Jinv*[dhi/dr]
    */
   for(int i=0;i<numnodes;i++){
      dbasis[i] = Jinv*dbasis_ref[i];
   }
}

GetNodalFunctionsDerivativesReference()

void SegRef::GetNodalFunctionsDerivativesReference	(	IssmDouble *	dbasis,
		GaussSeg *	gauss,
		int	finiteelement
	)

Definition at line 162 of file SegRef.cpp.

                                                                                                      {/*{{{*/
   /*This routine returns the values of the nodal functions derivatives  (with respect to the 
    * natural coordinate system) at the gaussian point. */
 
   _assert_(dbasis && gauss);
 
   switch(finiteelement){
      case P0Enum:
         /*Nodal function 1*/
         dbasis[0] = 0.;
         break;
      case P1Enum: case P1DGEnum:
         /*Nodal function 1*/
         dbasis[0] = -0.5;
         /*Nodal function 2*/
         dbasis[1] = 0.5;
         return;
      case P2Enum:
         /*Nodal function 1*/
         dbasis[0] = (gauss->coord1-1.)/2. + gauss->coord1/2.;
         dbasis[1] = (1.+gauss->coord1)/2. + gauss->coord1/2.;
         dbasis[2] = -2.*gauss->coord1;
         return;
      default:
         _error_("Element type "<<EnumToStringx(finiteelement)<<" not supported yet");
   }
 
}

Marshall()

void SegRef::Marshall ( char **  pmarshalled_data, int *  pmarshalled_data_size, int  marshall_direction  )

inline

Definition at line 26 of file SegRef.h.

{ /*do nothing */};

NumberofNodes()

int SegRef::NumberofNodes

(

int

finiteelement

)

Definition at line 191 of file SegRef.cpp.

                                           {/*{{{*/
 
   switch(finiteelement){
      case P0Enum:                return NUMNODESP0;
      case P1Enum:                return NUMNODESP1;
      case P1DGEnum:              return NUMNODESP1;
      default: _error_("Element type "<<EnumToStringx(finiteelement)<<" not supported yet");
   }
 
   return -1;
}

---

The documentation for this class was generated from the following files:

• src/c/classes/Elements/SegRef.h
• src/c/classes/Elements/SegRef.cpp