Changeset 1206

Timestamp:

07/01/09 16:33:06 (16 years ago)

Author:

Mathieu Morlighem

Message:

use numberofdofspernode

File:

• issm/trunk/src/c/parallel/GradJCompute.cpp (modified) (3 diffs)

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

@@ -49 +49 @@
         femmodel->parameters->FindParam((void*)&sub_analysis_type,"sub_analysis_type");
         femmodel->parameters->FindParam((void*)&numberofnodes,"numberofnodes");
+        femmodel->parameters->FindParam((void*)&numberofdofspernode,"numberofdofspernode");
         femmodel->parameters->FindParam((void*)&solverstring,"solverstring");
         femmodel->parameters->FindParam((void*)&control_type,"control_type");
@@ -57 +58 @@
         diagnostic_core_nonlinear(&u_g,&K_ff0,&K_fs0,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);
 
         _printf_("%s\n","      buid Du, difference between observed velocity and model velocity:");
@@ -76 +77 @@
         _printf_("%s\n","      compute gradJ:");
         VecToMPISerial(&lambda_g_double,lambda_g);VecFree(&lambda_g);
-        inputs->Add("adjoint",lambda_g_double,2,numberofnodes);
+        inputs->Add("adjoint",lambda_g_double,numberofdofspernode,numberofnodes);
 
         Gradjx( &grad_g, numberofnodes,femmodel->elements,femmodel->nodes, femmodel->loads, femmodel->materials,