Changeset 26073 for issm/trunk-jpl/src/c/classes/Elements/Element.cpp

Timestamp:

03/11/21 10:05:44 (4 years ago)

Author:

Eric.Larour

Message:

CHG: fixed test2001 by modifying the mass transport core. Now create inputs
directly in the core, not the MassTransportAnalysis::InputUpdateFromSolution.
Also need a new IoModel::InputToConstant routine.

File:

• issm/trunk-jpl/src/c/classes/Elements/Element.cpp (modified) (2 diffs)

issm/trunk-jpl/src/c/classes/Elements/Element.cpp

@@ -1709 +1709 @@
 }
 /*}}}*/
+void       Element::InputCreateP1FromConstant(Inputs* inputs,IoModel* iomodel,IssmDouble value_in,int vector_enum){/*{{{*/
+
+        const int NUM_VERTICES = this->GetNumberOfVertices();
+
+        int        vertexids[MAXVERTICES];
+        int        vertexlids[MAXVERTICES];
+        IssmDouble values[MAXVERTICES];
+
+        /*Recover vertices ids needed to initialize inputs*/
+        _assert_(iomodel->elements);
+        for(int i=0;i<NUM_VERTICES;i++){
+                vertexids[i] =reCast<int>(iomodel->elements[NUM_VERTICES*this->Sid()+i]); //ids for vertices are in the elements array from Matlab
+                vertexlids[i]=iomodel->my_vertices_lids[vertexids[i]-1];
+        }
+
+        for(int i=0;i<NUM_VERTICES;i++) values[i]=value_in;
+        this->SetElementInput(inputs,NUM_VERTICES,vertexlids,values,vector_enum);
+
+}
+/*}}}*/
 void       Element::ControlInputCreate(IssmDouble* vector,IssmDouble* min_vector,IssmDouble* max_vector,Inputs* inputs,IoModel* iomodel,int M,int N,IssmDouble scale,int input_enum,int id){/*{{{*/
 
@@ -1946 +1966 @@
         /*update input*/
         this->SetElementInput(name,constant);
+}
+/*}}}*/
+void       Element::InputUpdateFromConstant(IssmDouble constant, int name, int type){/*{{{*/
+
+        /*Check that name is an element input*/
+        if(!IsInputEnum(name)) return;
+
+        /*update input*/
+        this->SetElementInput(name,constant,type);
 }
 /*}}}*/