Changeset 16751

Timestamp:

11/13/13 20:50:52 (11 years ago)

Author:

seroussi

Message:

NEW: added InputUpdateFromSolution for AdjointHoriz

File:

• issm/trunk-jpl/src/c/analyses/AdjointHorizAnalysis.cpp (modified) (1 diff)

issm/trunk-jpl/src/c/analyses/AdjointHorizAnalysis.cpp ¶

@@ -30 +30 @@
 }/*}}}*/
 void AdjointHorizAnalysis::InputUpdateFromSolution(IssmDouble* solution,Element* element){/*{{{*/
-        _error_("not implemented yet");
+        int  i;
+        int* doflist=NULL;
+
+        /*Fetch number of nodes and dof for this finite element*/
+        int numnodes = element->GetNumberOfNodes();
+        int numdof   = numnodes*2;
+
+        /*Fetch dof list and allocate solution vectors*/
+        element->GetDofList(&doflist,NoneApproximationEnum,GsetEnum);
+        IssmDouble* values  = xNew<IssmDouble>(numdof);
+        IssmDouble* lambdax = xNew<IssmDouble>(numnodes);
+        IssmDouble* lambday = xNew<IssmDouble>(numnodes);
+
+        /*Use the dof list to index into the solution vector: */
+        for(i=0;i<numdof;i++) values[i]=solution[doflist[i]];
+
+        /*Transform solution in Cartesian Space*/
+        element->TransformSolutionCoord(&values[0],XYEnum);
+
+        /*Ok, we have vx and vy in values, fill in vx and vy arrays: */
+        for(i=0;i<numnodes;i++){
+                lambdax[i]=values[i*NDOF2+0];
+                lambday[i]=values[i*NDOF2+1];
+
+                /*Check solution*/
+                if(xIsNan<IssmDouble>(lambdax[i])) _error_("NaN found in solution vector");
+                if(xIsNan<IssmDouble>(lambday[i])) _error_("NaN found in solution vector");
+        }
+
+        /*Add vx and vy as inputs to the tria element: */
+        element->AddInput(AdjointxEnum,lambdax,P1Enum);
+        element->AddInput(AdjointyEnum,lambday,P1Enum);
+
+        /*Free ressources:*/
+        xDelete<IssmDouble>(values);
+        xDelete<IssmDouble>(lambdax);
+        xDelete<IssmDouble>(lambday);
+        xDelete<int>(doflist);
 }/*}}}*/