Changeset 25413

Timestamp:

08/16/20 21:35:14 (5 years ago)

Author:

Mathieu Morlighem

Message:

CHG: fixing a few things but transient P1 inversion not quite working yet still

File:

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

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

@@ -1121 +1121 @@
                 if(gradient_input->GetInputInterpolationType()==P1Enum){
                         _assert_(N==1);
-                        for(int i=0;i<NUMVERTICES;i++){
-                                values[i] = gradient[idlist[i]];
-                        }
+                        for(int i=0;i<NUMVERTICES;i++) values[i] = gradient[idlist[i]];
                         gradient_input->SetInput(P1Enum,NUMVERTICES,&lidlist[0],&values[0]);
                 }
@@ -1144 +1142 @@
                 for(int n=0;n<N;n++){
                         if(interp[control_index]==P1Enum){
-                                _error_("not implemented");
+                                for(int i=0;i<NUMVERTICES;i++) values[i] = gradient[idlist[i]];
+                                gradient_input->AddTriaTimeInput(n,NUMVERTICES,&lidlist[0],&values[0],P1Enum);
                         }
                         else if(interp[control_index]==P0Enum){
@@ -2474 +2473 @@
                         TriaInput* input_n = input->GetTriaInput(n); _assert_(input_n);
                         if(input_n->GetInputInterpolationType()==P1Enum){
-                                IssmDouble values[NUMVERTICES];
                                 input_n->Serve(NUMVERTICES,&lidlist[0]);
                                 for(int i=0;i<NUMVERTICES;i++) values[count+i] = input_n->element_values[i];