Changeset 5250 for issm/trunk/src/c/modules/MassFluxx/MassFluxx.cpp

Timestamp:

08/13/10 14:54:02 (15 years ago)

Author:

Eric.Larour

Message:

New Response module, to centralize all response management.
Now used by DakotaResponsesx. Will be used by CostFunctionx.
Corrected parametercontrolB for new naming conventions.

File:

• issm/trunk/src/c/modules/MassFluxx/MassFluxx.cpp (modified) (3 diffs)

issm/trunk/src/c/modules/MassFluxx/MassFluxx.cpp

@@ -10 +10 @@
 #include "../../EnumDefinitions/EnumDefinitions.h"
 
-void MassFluxx(double* pmass_flux, Elements* elements,Nodes* nodes, Vertices* vertices,Loads* loads,Materials* materials,Parameters* parameters,  double* segments,int num_segments,bool process_units){
+void MassFluxx(double* pmass_flux, Elements* elements,Nodes* nodes, Vertices* vertices,Loads* loads,Materials* materials, Parameters* parameters,char* descriptor,bool process_units){
 
         int i,j;
@@ -23 +23 @@
         double all_mass_flux=0;
 
+        int  counter;
+
+        /*all segments: */
+        double** array=NULL;
+        int      M;
+        int*     mdims_array=NULL;
+        int*     ndims_array=NULL;
+
+        /*our segments of interest: */
+        double*  segments=NULL;
+        int      num_segments;
+
+        /*First, figure out which segment to compute our mass flux on. Start with retrieving qmu_mass_flux_segments: */
+        parameters->FindParam(&array,&M,&mdims_array,&ndims_array,QmuMassFluxSegmentsEnum);
+
+        /*figure out counter in the descriptor: */
+        sscanf(descriptor,"MassFlux%i",&counter);
+
+        /*retrieve segments from array: */
+        segments=array[counter-1]; //matlab to "C" indexing
+        num_segments=mdims_array[counter-1];
+
         /*Go through segments, and then elements, and figure out which elements belong to a segment. 
          * When we find one, use the element to compute the mass flux on the segment: */
-
         for(i=0;i<num_segments;i++){
                 element_id=(int)*(segments+5*i+4);
@@ -43 +64 @@
         #endif
 
+        /*Free ressources:*/
+        for(j=0;j<M;j++){
+                double* matrix=array[j];
+                xfree((void**)&matrix);
+        }
+        xfree((void**)&mdims_array);
+        xfree((void**)&ndims_array);
+        xfree((void**)&array);
+
         /*Assign output pointers: */
         *pmass_flux=mass_flux;