GetVectorFromInputsx.cpp File Reference

#include "./GetVectorFromInputsx.h"
#include "../../shared/shared.h"
#include "../../toolkits/toolkits.h"

Go to the source code of this file.

Functions
void	GetVectorFromInputsx (Vector< IssmDouble > **pvector, FemModel *femmodel, int name, int type)
void	GetVectorFromInputsx (Vector< IssmDouble > **pvector, FemModel *femmodel, int name, int type, IssmDouble time)
void	GetVectorFromInputsx (IssmDouble **pvector, FemModel *femmodel, int name, int type)
void	GetVectorFromInputsx (IssmDouble **pvector, int *pvector_size, FemModel *femmodel, int name)

Function Documentation

GetVectorFromInputsx() [1/4]

void GetVectorFromInputsx	(	Vector< IssmDouble > **	pvector,
		FemModel *	femmodel,
		int	name,
		int	type
	)

Definition at line 9 of file GetVectorFromInputsx.cpp.

                                                                                            { /*{{{*/
 
   int i;
   Vector<IssmDouble>* vector=NULL;
 
   switch(type){
      case ElementSIdEnum:
         vector=new Vector<IssmDouble>(femmodel->elements->NumberOfElements());
         break;
      case VertexPIdEnum: case VertexSIdEnum:
         vector=new Vector<IssmDouble>(femmodel->vertices->NumberOfVertices());
         break;
      case NodesEnum:case NodeSIdEnum:
         vector=new Vector<IssmDouble>(femmodel->nodes->NumberOfNodes());
         break;
      default:
         _error_("vector type: " << EnumToStringx(type) << " not supported yet!");
   }
   /*Look up in elements*/
   for(i=0;i<femmodel->elements->Size();i++){
      Element* element=xDynamicCast<Element*>(femmodel->elements->GetObjectByOffset(i));
      element->GetVectorFromInputs(vector,name,type);
   }
 
   vector->Assemble();
 
   /*Assign output pointers:*/
   *pvector=vector;
} /*}}}*/

GetVectorFromInputsx() [2/4]

void GetVectorFromInputsx	(	Vector< IssmDouble > **	pvector,
		FemModel *	femmodel,
		int	name,
		int	type,
		IssmDouble	time
	)

Definition at line 38 of file GetVectorFromInputsx.cpp.

                                                                                                            {/*{{{*/
 
   int i;
   Vector<IssmDouble>* vector=NULL;
 
   switch(type){
   case VertexPIdEnum: case VertexSIdEnum:
      vector=new Vector<IssmDouble>(femmodel->vertices->NumberOfVertices());
      break;
   case NodesEnum:case NodeSIdEnum:
      vector=new Vector<IssmDouble>(femmodel->nodes->NumberOfNodes());
      break;
   default:
         _error_("vector type: " << EnumToStringx(type) << " not supported yet!");
   }
   /*Look up in elements*/
   for(i=0;i<femmodel->elements->Size();i++){
      Element* element=xDynamicCast<Element*>(femmodel->elements->GetObjectByOffset(i));
      element->GetVectorFromInputs(vector,name,type,time);
   }
 
   vector->Assemble();
 
   /*Assign output pointers:*/
   *pvector=vector;
}/*}}}*/

GetVectorFromInputsx() [3/4]

void GetVectorFromInputsx	(	IssmDouble **	pvector,
		FemModel *	femmodel,
		int	name,
		int	type
	)

Definition at line 64 of file GetVectorFromInputsx.cpp.

                                                                                     {/*{{{*/
 
   /*output: */
   IssmDouble* vector=NULL;
 
   /*intermediary: */
   Vector<IssmDouble>* vec_vector=NULL;
 
   GetVectorFromInputsx(&vec_vector,femmodel,name,type);
   vector=vec_vector->ToMPISerial();
 
   /*Free ressources:*/
   delete vec_vector;
 
   /*Assign output pointers:*/
   *pvector=vector;
}/*}}}*/

GetVectorFromInputsx() [4/4]

void GetVectorFromInputsx	(	IssmDouble **	pvector,
		int *	pvector_size,
		FemModel *	femmodel,
		int	name
	)

Definition at line 81 of file GetVectorFromInputsx.cpp.

                                                                                              { /*{{{*/
 
   int interpolation_type;
   /*this one is special: we don't specify the type, but let the nature of the inputs dictace. 
    * P0 -> ElementSIdEnum, P1 ->VertexSIdEnum: */
 
   /*We go find the input of the first element, and query its interpolation type: */
   Element* element=xDynamicCast<Element*>(femmodel->elements->GetObjectByOffset(0));
   Input2* input=element->GetInput2(name); 
   if (!input) _error_("could not find input: " << name);
 
   interpolation_type=input->GetInputInterpolationType(); 
   if(interpolation_type==P0Enum){
      *pvector_size=femmodel->elements->NumberOfElements();
      GetVectorFromInputsx(pvector,femmodel,name, ElementSIdEnum);
   }
   else if(interpolation_type==P1Enum){
      *pvector_size=femmodel->vertices->NumberOfVertices();
      GetVectorFromInputsx(pvector,femmodel,name, VertexSIdEnum);
   }
   else _error_("interpolation type : " << interpolation_type << " not supported yet!");
}/*}}}*/