Changeset 24426

Timestamp:

12/03/19 21:05:47 (5 years ago)

Author:

Mathieu Morlighem

Message:

CHG: making InputDepthAverage more efficient

File:

• issm/trunk-jpl/src/c/classes/Elements/Penta.cpp (modified) (3 diffs)

issm/trunk-jpl/src/c/classes/Elements/Penta.cpp

@@ -2214 +2214 @@
         Penta* penta = this;
         int    step  = 0;
+        Gauss* gauss[3];
+        for(int iv=0;iv<3;iv++) gauss[iv] = penta->NewGaussLine(iv,iv+3,3);
+
         for(;;){
 
@@ -2229 +2232 @@
                         }
                         /*Integrate over edge*/
-                        Gauss* gauss=penta->NewGaussLine(iv,iv+3,3);
-                        for(int ig=gauss->begin();ig<gauss->end();ig++){
-                                gauss->GaussPoint(ig);
-                                penta->JacobianDeterminantLine(&Jdet,&xyz_list_line[0][0],gauss);
-                                original_input->GetInputValue(&value,gauss);
-                                total[iv] += value*Jdet*gauss->weight;
-                                intz[iv]  += Jdet*gauss->weight;
+                        for(int ig=gauss[iv]->begin();ig<gauss[iv]->end();ig++){
+                                gauss[iv]->GaussPoint(ig);
+                                penta->JacobianDeterminantLine(&Jdet,&xyz_list_line[0][0],gauss[iv]);
+                                original_input->GetInputValue(&value,gauss[iv]);
+                                total[iv] += value*Jdet*gauss[iv]->weight;
+                                intz[iv]  += Jdet*gauss[iv]->weight;
                         }
-                        delete gauss;
                 }
 
@@ -2247 +2248 @@
                 step++;
         }
+        for(int iv=0;iv<3;iv++) delete gauss[iv];
 
         /*Now we only need to divide the depth integrated input by the total thickness!*/