Changeset 21746 for issm/trunk-jpl/src/c/modules/GetVectorFromInputsx/GetVectorFromInputsx.cpp

Timestamp:

05/26/17 19:47:25 (8 years ago)

Author:

Eric.Larour

Message:

CHG: now possible to retrieve without an vector type.

File:

• issm/trunk-jpl/src/c/modules/GetVectorFromInputsx/GetVectorFromInputsx.cpp (modified) (1 diff)

issm/trunk-jpl/src/c/modules/GetVectorFromInputsx/GetVectorFromInputsx.cpp

@@ -54 +54 @@
         *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));
+        Input* input=element->GetInput(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!");
+}