Changeset 2112 for issm/trunk/src/c/Qmux/SpawnCoreParallel.cpp

Timestamp:

09/04/09 16:58:25 (16 years ago)

Author:

Eric.Larour

Message:

Heavy debugging of Mass flux capability

File:

• issm/trunk/src/c/Qmux/SpawnCoreParallel.cpp (modified) (3 diffs)

issm/trunk/src/c/Qmux/SpawnCoreParallel.cpp

@@ -42 +42 @@
         /*output from core solutions: */
         DataSet* results=NULL;
+        DataSet* processed_results=NULL;
 
         char** responses_descriptors=NULL;
@@ -110 +111 @@
         #endif
 
+        /*broadcast response descriptors: */
+        MPI_Bcast(&numresponses,1,MPI_INT,0,MPI_COMM_WORLD); 
+        if(my_rank!=0){
+                responses_descriptors=(char**)xmalloc(numresponses*sizeof(char*));
+        }
+        for(i=0;i<numresponses;i++){
+                if(my_rank==0){
+                        string=responses_descriptors[i];
+                        string_length=(strlen(string)+1)*sizeof(char);
+                }
+                MPI_Bcast(&string_length,1,MPI_INT,0,MPI_COMM_WORLD); 
+                if(my_rank!=0)string=(char*)xmalloc(string_length);
+                MPI_Bcast(string,string_length,MPI_CHAR,0,MPI_COMM_WORLD); 
+                if(my_rank!=0)responses_descriptors[i]=string;
+        }
+
+        #ifdef _ISSM_DEBUG_
+        for(i=0;i<numresponses;i++){
+                PetscSynchronizedPrintf(MPI_COMM_WORLD,"variable descriptor %i: %s value: %g\n",i,responses_descriptors[i],responses[i]);
+                PetscSynchronizedFlush(MPI_COMM_WORLD);
+        }
+        #endif
+
+
         _printf_("qmu iteration: %i\n",counter);
 
@@ -150 +175 @@
         /*Now process the outputs, before computing the dakota responses: */
         if(debug)_printf_("process results:\n");
-        ProcessResults(&results,model,analysis_type); 
-        
-
+        ProcessResults(&processed_results,results,model,analysis_type); 
+        
         /*compute responses on cpu 0: dummy for now! */
         if(debug)_printf_("compute dakota responses:\n");
-        DakotaResponses(responses,responses_descriptors,numresponses,model,results,analysis_type,sub_analysis_type);
+        DakotaResponses(responses,responses_descriptors,numresponses,model,results,processed_results,analysis_type,sub_analysis_type);
 
         /*Free ressources:*/
         delete results;
+        delete processed_results;
 
         //variables only on cpu != 0