Changeset 25215

Timestamp:

07/06/20 18:02:00 (5 years ago)

Author:

Eric.Larour

Message:

CHG: handle distributed uncertain variables. These variables give rise to special
behaviour in terms of updating the inputs, using code specific to each enum coming in.

Location:

issm/trunk-jpl/src/c/modules/InputUpdateFromDakotax

Files:

• InputUpdateFromDakotax.cpp (modified) (3 diffs)
• InputUpdateFromDakotax.h (modified) (1 diff)

issm/trunk-jpl/src/c/modules/InputUpdateFromDakotax/InputUpdateFromDakotax.cpp

@@ -6 +6 @@
 #include "../../shared/shared.h"
 #include "../../toolkits/toolkits.h"
+#include "../../classes/Inputs2/TransientInput2.h"
+#include "../../classes/Inputs2/DatasetInput2.h"
+#include "../../classes/Inputs2/TriaInput2.h"
 #include "../InputUpdateFromMatrixDakotax/InputUpdateFromMatrixDakotax.h"
 #include "../InputUpdateFromConstantx/InputUpdateFromConstantx.h"
 #include "../InputUpdateFromVectorDakotax/InputUpdateFromVectorDakotax.h"
-
-void InputUpdateFromDakotax(FemModel* femmodel,double* variables,char* *variables_descriptors,int numdakotavariables){
+                        
+void InputUpdateFromDakotax(FemModel* femmodel,double* variables,char* *variables_descriptors,int numdakotavariables){ /*{{{*/
 
         int     i,j,k,l;
@@ -166 +169 @@
                         xDelete<double>(distributed_values);
                 }
+                else if (strncmp(descriptor,"distributed_",12)==0){
+                        
+                        /*recover partition vector: */
+                        variable_partition=variable_partitions[variablecount];
+                        npart=variable_partitions_npart[variablecount];
+
+                        /*Variable is distributed. Determine root name of variable (ex: distributed_DragCoefficient_1 -> DragCoefficient). 
+                         * Allocate distributed_values and fill the distributed_values with the next npart variables: */
+
+                        memcpy(root,strstr(descriptor,"_")+1,(strlen(strstr(descriptor,"_")+1)+1)*sizeof(char));
+                        *strstr(root,"_")='\0';
+
+                        distributed_values=xNew<double>(npart);
+                        for(j=0;j<npart;j++){
+                                distributed_values[j]=variables[i+j];
+                        }
+
+                        //for (int j=0;j<npart;j++)_printf_(j << ":" << distributed_values[j] << "\n");
+                        
+                        //Call specialty code:
+                        InputUpdateSpecialtyCode(femmodel,distributed_values,variable_partition,npart,root);
+                        
+                        /*increment i to skip the distributed values just collected: */
+                        i+=npart-1; //careful, the for loop will add 1.
+                        
+                        /*Free allocations: */
+                        xDelete<double>(parameter);
+                        xDelete<double>(distributed_values);
+                }
                 else if (strncmp(descriptor,"indexed_",8)==0){
                         _error_("indexed variables not supported yet!");
@@ -188 +220 @@
         xDelete<int>(variable_partitions_nt);
 
-}
+} /*}}}*/
+void  InputUpdateSpecialtyCode(FemModel* femmodel,IssmDouble* distributed_values,IssmDouble* variable_partition,int npart,char* root){ //{{{
+
+        /*Here, we put all the code that cannot be handled any other place: */
+        if (strncmp(root,"SurfaceloadModelid",18)==0){
+
+                if(VerboseQmu())_printf_("Updating SurfaceloadModelid MME, with ids: " << "\n");
+                for (int i=0;i<npart;i++)_printf_(distributed_values[i] << " ");
+
+
+                if (femmodel->inputs2->GetInputObjectEnum(SurfaceloadIceThicknessChangeEnum)==DatasetInput2Enum)
+                        MmeToInput(femmodel,distributed_values,variable_partition,npart,SurfaceloadIceThicknessChangeEnum, P0Enum);
+
+                if (femmodel->inputs2->GetInputObjectEnum(MaskIceLevelsetEnum)==DatasetInput2Enum)
+                        MmeToInput(femmodel,distributed_values,variable_partition,npart,MaskIceLevelsetEnum, P1Enum);
+
+                if (femmodel->inputs2->GetInputObjectEnum(MaskOceanLevelsetEnum)==DatasetInput2Enum)
+                        MmeToInput(femmodel,distributed_values,variable_partition,npart,MaskOceanLevelsetEnum, P1Enum);
+
+        }
+        else _error_("InputUpdateSpecialtyCode error message: " << root << " not supported yet!");
+
+}       //}}}
+void  MmeToInput(FemModel* femmodel,IssmDouble* distributed_values,IssmDouble* variable_partition,int npart,int rootenum, int interpolationenum){ /*{{{*/
+
+        TransientInput2* transientinput  = NULL;
+        TransientInput2* transientinput2 = NULL;
+        Tria* element                    = NULL;
+        IssmDouble value;
+        IssmDouble* values               = NULL;
+        IssmDouble* times                = NULL;
+        int N;
+        int id;
+
+        /*find thickness dataset: */
+        DatasetInput2* datasetinput = femmodel->inputs2->GetDatasetInput2(rootenum);
+
+        /*Initialize new transient input: */
+        transientinput = datasetinput->GetTransientInputByOffset(0); _assert_(transientinput);
+        transientinput->GetAllTimes(&times,&N);
+        femmodel->inputs2->SetTransientInput(DummyEnum,times,N);
+        transientinput2 = femmodel->inputs2->GetTransientInput(DummyEnum);
+                
+        for (int i=0;i<femmodel->elements->Size();i++){
+
+                Tria*   element=xDynamicCast<Tria*>(femmodel->elements->GetObjectByOffset(i));
+
+                if((int)variable_partition[element->Sid()]==-1)id=0; //grab background field
+                else id=distributed_values[(int)variable_partition[element->Sid()]]; //grab partition field
+
+                if (element->Sid()==577 && rootenum==MaskOceanLevelsetEnum)_printf_("577 -> " << id);
+
+                /*recover the right field from the mme: */
+                transientinput = datasetinput->GetTransientInputByOffset(id); _assert_(transientinput);
+
+                /*copy values from the transientinput to the final transientinput2: */
+                for (int j=0;j<N;j++){
+                        TriaInput2* tria_input=transientinput->GetTriaInput(j);
+                        element->InputServe(tria_input);
+                        if(interpolationenum==P0Enum){
+                                value=tria_input->element_values[0];
+                                transientinput2->AddTriaTimeInput( j,1,&(element->lid),&value,P0Enum); 
+                        }
+                        else if(interpolationenum==P1Enum){
+
+                                /*Get values and lid list*/
+                                const int   numvertices     = element->GetNumberOfVertices();
+                                int        *vertexlids      = xNew<int>(numvertices);
+                                int        *vertexsids      = xNew<int>(numvertices);
+
+                                /*Recover vertices ids needed to initialize inputs*/
+                                element->GetVerticesLidList(&vertexlids[0]);
+                                element->GetVerticesSidList(&vertexsids[0]);
+                                values=tria_input->element_values;
+                                if (element->Sid()==577 && rootenum==MaskOceanLevelsetEnum)_printf_(" sids: " << vertexsids[0] << " " << vertexsids[1] << " " << vertexsids[2] << " values: " << values[0] << " " << values[1] << " " << values[2] << "\n");
+                                transientinput2->AddTriaTimeInput( j,numvertices,vertexlids,values,P1Enum); 
+                        }
+                }
+        }
+
+        /*wipe out existing SurfaceloadIceThicknessChangeEnum dataset:*/
+        femmodel->inputs2->ChangeEnum(DummyEnum,rootenum);
+
+        //reconfigure:
+        transientinput2->Configure(femmodel->parameters);
+}       //}}}

issm/trunk-jpl/src/c/modules/InputUpdateFromDakotax/InputUpdateFromDakotax.h

@@ -9 +9 @@
 
 void  InputUpdateFromDakotax(FemModel* femmodel,double* variables,char* *variables_descriptors,int numvariables);
+void  InputUpdateSpecialtyCode(FemModel* femmodel,IssmDouble* distributed_values,IssmDouble* variable_partition,int npart,char* root);
+void  MmeToInput(FemModel* femmodel,IssmDouble* distributed_values,IssmDouble* variable_partition,int npart,int rootenum, int interpolationenum);
+
 
 #endif  /* _INPUTUPDATEFROMDAKOTAXX_H */