Changeset 3028

Timestamp:

02/12/10 19:19:23 (15 years ago)

Author:

Eric.Larour

Message:

Added misfti response.

File:

• issm/trunk/src/c/parallel/DakotaResponses.cpp (modified) (1 diff)

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

@@ -182 +182 @@
                         if(my_rank==0)responses[i]=max_abs_vz;
                 }
+                else if(strcmp(response_descriptor,"misfit")==0){
+                        
+                        ParameterInputs* inputs=NULL; //needed to plug velocity into elements.
+                        double J=0;
+                        int analysis_type,sub_analysis_type;
+                        int numberofdofspernode,numberofnodes;
+                        Vec u_g=NULL;
+                
+                        /*Recover analysis_type: */
+                        fem->parameters->FindParam(&analysis_type,"analysis_type");
+                        fem->parameters->FindParam(&sub_analysis_type,"sub_analysis_type");
+                        fem->parameters->FindParam(&numberofdofspernode,"numberofdofspernode");
+                        fem->parameters->FindParam(&numberofnodes,"numberofnodes");
+        
+                        /*Recover velocity: */
+                        found=results->FindResult(&u_g,"u_g");
+                        if(!found)throw ErrorException(__FUNCT__," could not find velocity to compute misfit");
+
+                        /*Add to inputs: */
+                        inputs->Add("velocity",u_g,numberofdofspernode,numberofnodes);
+
+                        /*Compute misfit: */
+                        Misfitx( &J, fem->elements,fem->nodes, fem->loads, fem->materials, fem->parameters,inputs,analysis_type,sub_analysis_type);
+        
+                        if(my_rank==0)responses[i]=J;
+
+                        /*Some cleanup: */
+                        VecFree(&u_g);
+
+                }
                 else{
                         if(my_rank==0)printf("%s%s%s"," response descriptor : ",response_descriptor," not supported yet!");