Changeset 1208

Timestamp:

07/01/09 17:14:38 (15 years ago)

Author:

Mathieu Morlighem

Message:

use numberofdofspernode needed for Stokes

Location:

issm/trunk/src/c/parallel

Files:

• GradJCompute.cpp (modified) (1 diff)
• objectivefunctionC.cpp (modified) (3 diffs)

issm/trunk/src/c/parallel/GradJCompute.cpp

@@ -22 +22 @@
         int sub_analysis_type;
         int numberofnodes;
+        int numberofdofspernode;
         char* solverstring=NULL;
         char* control_type=NULL;

issm/trunk/src/c/parallel/objectivefunctionC.cpp

@@ -27 +27 @@
         double* grad_g=NULL;
         ParameterInputs* inputs=NULL;
+        int numberofdofspernode;
         int n;
 
@@ -60 +61 @@
         femmodel->parameters->FindParam((void*)&sub_analysis_type,"sub_analysis_type");
         femmodel->parameters->FindParam((void*)&numberofnodes,"numberofnodes");
+        femmodel->parameters->FindParam((void*)&numberofdofspernode,"numberofdofspernode");
 
         /*First copy param_g so we don't modify it: */
@@ -77 +79 @@
         diagnostic_core_nonlinear(&u_g,NULL,NULL,femmodel,inputs,analysis_type,sub_analysis_type);
         VecToMPISerial(&u_g_double,u_g); VecFree(&u_g);
-        inputs->Add("velocity",u_g_double,2,numberofnodes);
+        inputs->Add("velocity",u_g_double,numberofdofspernode,numberofnodes);
 
         //Compute misfit for this velocity field.