Changeset 26887

Timestamp:

02/16/22 07:27:08 (3 years ago)

Author:

Mathieu Morlighem

Message:

NEW: added method to update input from passive double, needed for checkpointing!

Location:

issm/trunk-jpl/src/c/modules/SetControlInputsFromVectorx

Files:

• SetControlInputsFromVectorx.cpp (modified) (1 diff)
• SetControlInputsFromVectorx.h (modified) (1 diff)

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

@@ -40 +40 @@
         xDelete<IssmDouble>(serial_vector);
 }
+void SetControlInputsFromVectorx(FemModel* femmodel,IssmPDouble* vector){
+
+        /*Get total size and recast*/
+        int  num_controls;
+        int* M = NULL;
+        int* N = NULL;
+        femmodel->parameters->FindParam(&num_controls,InversionNumControlParametersEnum);
+        femmodel->parameters->FindParam(&M,NULL,ControlInputSizeMEnum);
+        femmodel->parameters->FindParam(&N,NULL,ControlInputSizeNEnum);
+
+        int size = 0;
+        for(int i=0;i<num_controls;i++) size += M[i]*N[i]; 
+
+        IssmDouble* serial_vector=xNew<IssmDouble>(size);
+        for(int i=0;i<size;i++) serial_vector[i] = reCast<IssmDouble>(vector[i]);
+
+        SetControlInputsFromVectorx(femmodel,serial_vector);
+
+        xDelete<IssmDouble>(serial_vector);
+        xDelete<int>(M);
+        xDelete<int>(N);
+}

issm/trunk-jpl/src/c/modules/SetControlInputsFromVectorx/SetControlInputsFromVectorx.h

@@ -10 +10 @@
 void SetControlInputsFromVectorx(FemModel* femmodel,Vector<IssmDouble>* vector);
 void SetControlInputsFromVectorx(FemModel* femmodel,IssmDouble* vector);
+void SetControlInputsFromVectorx(FemModel* femmodel,IssmPDouble* vector);
 
 #endif