source: issm/oecreview/Archive/23390-24306/ISSM-23546-23547.diff

../trunk-jpl/src/c/classes/FemModel.cpp

@@ -4829 +4829 @@
                 else{
                         for(int j=0;j<elements->Size();j++){
                                 /*Intermediaries*/
-                                Element*   element            =xDynamicCast<Element*>(elements->GetObjectByOffset(j));
-                                int        numvertices        =element->GetNumberOfVertices();
-                                IssmDouble zeros[numvertices] ={0.0};
+                                Element*    element     = xDynamicCast<Element*>(elements->GetObjectByOffset(j));
+                                int         numvertices = element->GetNumberOfVertices();
+                                IssmDouble* zeros       = xNewZeroInit<IssmDouble>(numvertices);
                                 switch(element->ObjectEnum()){
-                                case TriaEnum:
-                                        element->inputs->AddInput(new TriaInput(stackedinput_enum[i],&zeros[0],P1Enum));
-                                        break;
-                                case PentaEnum:
-                                        element->inputs->AddInput(new PentaInput(stackedinput_enum[i],&zeros[0],P1Enum));
-                                        break;
-                                case TetraEnum:
-                                        element->inputs->AddInput(new TetraInput(stackedinput_enum[i],&zeros[0],P1Enum));
-                                        break;
-                                default: _error_("Not implemented yet");
+                                        case TriaEnum:
+                                                element->inputs->AddInput(new TriaInput(stackedinput_enum[i],&zeros[0],P1Enum));
+                                                break;
+                                        case PentaEnum:
+                                                element->inputs->AddInput(new PentaInput(stackedinput_enum[i],&zeros[0],P1Enum));
+                                                break;
+                                        case TetraEnum:
+                                                element->inputs->AddInput(new TetraInput(stackedinput_enum[i],&zeros[0],P1Enum));
+                                                break;
+                                        default: _error_("Not implemented yet");
                                 }
+                                xDelete<IssmDouble>(zeros);
                         }
                 }
         }