Changeset 23547

Timestamp:

12/13/18 08:14:18 (6 years ago)

Author:

Mathieu Morlighem

Message:

BUG: has to allocate memory for vector zeros

File:

• issm/trunk-jpl/src/c/classes/FemModel.cpp (modified) (1 diff)

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

@@ -4830 +4830 @@
                         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);
                         }
                 }