Changeset 4353 for issm/trunk/src/c/objects/Elements/Tria.cpp

Timestamp:

06/30/10 14:17:09 (15 years ago)

Author:

Mathieu Morlighem

Message:

fixing control methods

File:

• issm/trunk/src/c/objects/Elements/Tria.cpp (modified) (7 diffs)

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

@@ -302 +302 @@
 /*FUNCTION Tria::InputUpdateFromConstant(int value, int name);{{{1*/
 void  Tria::InputUpdateFromConstant(int constant, int name){
-        /*Nothing updated for now*/
+        /*Check that name is an element input*/
+        if (!IsInput(name)) return;
+
+        /*update input*/
+        this->inputs->AddInput(new IntInput(name,constant));
 }
 /*}}}*/
 /*FUNCTION Tria::InputUpdateFromConstant(double value, int name);{{{1*/
 void  Tria::InputUpdateFromConstant(double constant, int name){
-        /*Nothing updated for now*/
+        /*Check that name is an element input*/
+        if (!IsInput(name)) return;
+
+        /*update input*/
+        this->inputs->AddInput(new DoubleInput(name,constant));
 }
 /*}}}*/
 /*FUNCTION Tria::InputUpdateFromConstant(bool value, int name);{{{1*/
 void  Tria::InputUpdateFromConstant(bool constant, int name){
-        /*Nothing updated for now*/
+        /*Check that name is an element input*/
+        if (!IsInput(name)) return;
+
+        /*update input*/
+        this->inputs->AddInput(new BoolInput(name,constant));
 }
 /*}}}*/
@@ -329 +341 @@
         else if (analysis_type==DiagnosticHorizAnalysisEnum){
                 InputUpdateFromSolutionDiagnosticHoriz( solution);
+        }
+        else if (analysis_type==AdjointAnalysisEnum){
+                InputUpdateFromSolutionAdjoint( solution);
         }
         else if (analysis_type==BedSlopeAnalysisEnum || analysis_type==SurfaceSlopeAnalysisEnum){
@@ -721 +736 @@
         double scale=0;
         double S=0;
-        double fit=-1;
+        int    fit=-1;
 
         /*retrieve some parameters: */
@@ -1887 +1902 @@
         double scaley=1;
         double S=0;
-        double fit=-1;
+        int    fit=-1;
 
         /*inputs: */
@@ -2120 +2135 @@
         /*inputs: */
         bool onwater;
-        int  fit;
+        int fit;
 
         /*retrieve inputs :*/
@@ -5319 +5334 @@
 }
 /*}}}*/
+/*FUNCTION Tria::InputUpdateFromSolutionAdjoint {{{1*/
+void  Tria::InputUpdateFromSolutionAdjoint(double* solution){
+
+        int i;
+
+        const int    numvertices=3;
+        const int    numdofpervertex=2;
+        const int    numdof=numdofpervertex*numvertices;
+
+        int          doflist[numdof];
+        double       values[numdof];
+        double       lambdax[numvertices];
+        double       lambday[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];
+        }
+
+        /*Add vx and vy as inputs to the tria element: */
+        this->inputs->AddInput(new TriaVertexInput(AdjointxEnum,lambdax));
+        this->inputs->AddInput(new TriaVertexInput(AdjointyEnum,lambday));
+
+}
+/*}}}*/
 /*FUNCTION Tria::InputUpdateFromSolutionDiagnosticHoriz {{{1*/
 void  Tria::InputUpdateFromSolutionDiagnosticHoriz(double* solution){
@@ -5473 +5524 @@
 /*FUNCTION Tria::IsInput{{{1*/
 bool Tria::IsInput(int name){
-        if (name==SurfaceSlopeXEnum ||
-                                name==SurfaceSlopeYEnum){
+        if (
+                                name==SurfaceSlopeXEnum ||
+                                name==SurfaceSlopeYEnum ||
+                                name==FitEnum
+                ){
                 return true;
         }