Changeset 4478

Timestamp:

07/08/10 15:34:15 (15 years ago)

Author:

Mathieu Morlighem

Message:

Added UpdateFromSolutionAdjointStokes

File:

• issm/trunk/src/c/objects/Elements/Penta.cpp (modified) (6 diffs)

issm/trunk/src/c/objects/Elements/Penta.cpp

@@ -374 +374 @@
                 InputUpdateFromSolutionAdjointHoriz( solution);
         }
+        else if (analysis_type==AdjointStokesAnalysisEnum){
+                InputUpdateFromSolutionAdjointStokes( solution);
+        }
         else if (analysis_type==BedSlopeXAnalysisEnum){
                 InputUpdateFromBedSlopeX( solution);
@@ -979 +982 @@
         bool onbed;
         bool shelf;
-        int analysis_type,sub_analysis_type;
+        int analysis_type;
 
         /*retrieve parameters: */
         parameters->FindParam(&analysis_type,AnalysisTypeEnum);
-        parameters->FindParam(&sub_analysis_type,AnalysisTypeEnum);
 
         /*retrieve inputs :*/
@@ -999 +1001 @@
         if (!onbed) return;
 
-        if (sub_analysis_type==HorizAnalysisEnum){
+        if (analysis_type==AdjointHorizAnalysisEnum){
 
                 /*MacAyeal or Pattyn*/
@@ -1006 +1008 @@
                 delete tria;
         }
-        else if (sub_analysis_type==StokesAnalysisEnum){
+        else if (analysis_type==AdjointStokesAnalysisEnum){
 
                 /*Stokes*/
@@ -1013 +1015 @@
                 delete tria;
         }
-        else ISSMERROR("%s%i%s\n","sub_analysis: ",sub_analysis_type," not supported yet");
+        else ISSMERROR("analysis %s not supported yet",EnumAsString(analysis_type));
 
 
@@ -2441 +2443 @@
         /*Add surfaceslopey as inputs to the tria element: */
         this->inputs->AddInput(new PentaVertexInput(SurfaceSlopeYEnum,values));
+}
+/*}}}*/
+/*FUNCTION Penta::InputUpdateFromSolutionAdjointStokes {{{1*/
+void  Penta::InputUpdateFromSolutionAdjointStokes(double* solution){
+
+        int i;
+
+        const int    numvertices=6;
+        const int    numdofpervertex=4;
+        const int    numdof=numdofpervertex*numvertices;
+
+        int          doflist[numdof];
+        double       values[numdof];
+        double       lambdax[numvertices];
+        double       lambday[numvertices];
+        double       lambdaz[numvertices];
+        double       lambdap[numvertices];
+
+        int          dummy;
+
+        /*Get dof list: */
+        GetDofList(&doflist[0],&dummy);
+
+        /*Use the dof list to index into the solution vector: */
+        for(i=0;i<numdof;i++){
+                values[i]=solution[doflist[i]];
+        }
+
+        /*Ok, we have vx and vy in values, fill in vx and vy arrays: */
+        for(i=0;i<numvertices;i++){
+                lambdax[i]=values[i*numdofpervertex+0];
+                lambday[i]=values[i*numdofpervertex+1];
+                lambdaz[i]=values[i*numdofpervertex+2];
+                lambdap[i]=values[i*numdofpervertex+3];
+        }
+
+        /*Add vx and vy as inputs to the tria element: */
+        this->inputs->AddInput(new PentaVertexInput(AdjointxEnum,lambdax));
+        this->inputs->AddInput(new PentaVertexInput(AdjointyEnum,lambday));
+        this->inputs->AddInput(new PentaVertexInput(AdjointzEnum,lambdaz));
+        this->inputs->AddInput(new PentaVertexInput(AdjointpEnum,lambdap));
+
 }
 /*}}}*/