Changeset 669

Timestamp:

06/01/09 11:32:42 (16 years ago)

Author:

Eric.Larour

Message:

Process result for control methods

Location:

issm/trunk/src/c/parallel

Files:

• OutputControl.cpp (modified) (1 diff)
• ProcessResults.cpp (modified) (2 diffs)

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

@@ -32 +32 @@
         results->AddObject(result);
 
-        result=new Result(results->Size()+1,0,1,"p_g",p_g,numberofnodes);
+        result=new Result(results->Size()+1,0,1,"param_g",p_g,numberofnodes);
         results->AddObject(result);
 

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

@@ -57 +57 @@
         double* p_g_serial=NULL;
         double* pressure=NULL;
+        double* parameter=NULL;
         double* partition=NULL;
         double  yts;
@@ -216 +217 @@
                         xfree((void**)&partition);
                 }
+                else if(strcmp(result->GetFieldName(),"param_g")==0){
+                        /*easy, param_g is of size numberofnodes, on 1 dof, just repartition: */
+                        result->GetField(&p_g);
+                        VecToMPISerial(&p_g_serial,p_g);
+                        fem_p->parameters->FindParam((void*)&numberofnodes,"numberofnodes");
+                        VecToMPISerial(&partition,fem_p->partition);
+
+                        parameter=(double*)xmalloc(numberofnodes*sizeof(double));
+
+                        for(i=0;i<numberofnodes;i++){
+                                parameter[i]=p_g_serial[(int)partition[i]];
+                        }
+                        
+                        /*Ok, add parameter to newresults: */
+                        newresult=new Result(newresults->Size()+1,result->GetTime(),result->GetStep(),"parameter",parameter,numberofnodes);
+                        newresults->AddObject(newresult);
+
+                        /*do some cleanup: */
+                        xfree((void**)&p_g_serial);
+                        xfree((void**)&partition);
+                }
+                else{
+                        /*Just copy the result into the new results dataset: */
+                        newresults->AddObject(result);
+                }
         }