GetVectorFromInputsx.cpp

Go to the documentation of this file.

 
#include "./GetVectorFromInputsx.h"
#include "../../shared/shared.h"
#include "../../toolkits/toolkits.h"
 
void GetVectorFromInputsx(Vector<IssmDouble>** pvector,FemModel* femmodel,int name,int type){ /*{{{*/
 
   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;
} /*}}}*/
void GetVectorFromInputsx(Vector<IssmDouble>** pvector,FemModel* femmodel,int name,int type,IssmDouble time){/*{{{*/
 
   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;
}/*}}}*/
void GetVectorFromInputsx(IssmDouble** pvector,FemModel* femmodel,int name, int type){/*{{{*/
 
   /*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;
}/*}}}*/
void GetVectorFromInputsx(IssmDouble** pvector,int* pvector_size, FemModel* femmodel,int name){ /*{{{*/
 
   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!");
}/*}}}*/