Changeset 10573

Timestamp:

11/10/11 09:35:02 (13 years ago)

Author:

Eric.Larour

Message:

3D Dakota transient update

File:

• issm/trunk/src/c/objects/Elements/Penta.cpp (modified) (1 diff)

issm/trunk/src/c/objects/Elements/Penta.cpp

@@ -5081 +5081 @@
 }
 /*}}}*/
+/*FUNCTION Penta::InputUpdateFromVectorMatrix(double* matrix, int nrows, int ncols, int name, int type);{{{1*/
+void  Penta::InputUpdateFromVectorMatrix(double* matrix, int nrows, int ncols, int name, int type){
+        
+        int i,j,t;
+
+        /*Check that name is an element input*/
+        if (!IsInput(name)) return;
+
+        switch(type){
+
+                case VertexEnum:
+                        
+                        double values[6];
+                        double time;
+
+                        /*Ok, create transient input: */
+                        for(t=0;t<ncols;t++){ //ncols is the number of times
+
+                                /*create input values: */
+                                for(i=0;i<6;i++){
+                                        row=this->nodes[i]->GetSidList();
+                                        values[i]=(double)matrix[ncols*row+t];
+                                }
+
+                                /*time? :*/
+                                time=(double)matrix[(nrows-1)*ncols+t]*yts;
+
+                                if(t==0) transientinput=new TransientInput(vector_enum);
+                                transientinput->AddTimeInput(new PentaVertexInput(vector_enum,nodeinputs),time);
+                        }
+                        this->inputs->AddInput(transientinput);
+                        break;
+
+                default:
+                        _error_("type %i (%s) not implemented yet",type,EnumToStringx(type));
+        }
+
+}
+/*}}}*/
 #endif