Changeset 16188 for issm/trunk-jpl/src/c/modules/ModelProcessorx/Damage/CreateLoadsDamage.cpp

Timestamp:

09/19/13 17:33:08 (12 years ago)

Author:

Eric.Larour

Message:

CHG: finished putting together damage evolution framework. Need validation, verification.

File:

• issm/trunk-jpl/src/c/modules/ModelProcessorx/Damage/CreateLoadsDamage.cpp (modified) (1 diff)

issm/trunk-jpl/src/c/modules/ModelProcessorx/Damage/CreateLoadsDamage.cpp

@@ -8 +8 @@
 
 void    CreateLoadsDamage(Loads** ploads, IoModel* iomodel){
+
+        /*Recover pointer: */
+        Loads* loads=*ploads;
+
+        /*create penalties for nodes: no node can have a damage > 1*/
+        iomodel->FetchData(1,DamageSpcdamageEnum);
+        CreateSingleNodeToElementConnectivity(iomodel);
+
+        for(int i=0;i<iomodel->numberofvertices;i++){
+
+                /*keep only this partition's nodes:*/
+                if((iomodel->my_vertices[i]==1)){
+                        if (xIsNan<IssmDouble>(iomodel->Data(DamageSpcdamageEnum)[i])){ //No penalty applied on spc nodes!
+                                loads->AddObject(new Pengrid(iomodel->loadcounter+i+1,i,iomodel,DamageEvolutionAnalysisEnum));
+                        }
+                }
+        }
+        iomodel->DeleteData(1,DamageSpcdamageEnum);
+
+        /*Assign output pointer: */
+        *ploads=loads;
         
-        /*No loads*/
-
 }