Changeset 17472 for issm/trunk-jpl/src/c/classes/Elements/Element.cpp

Timestamp:

03/18/14 22:44:50 (11 years ago)

Author:

Mathieu Morlighem

Message:

NEW: working on Tetra elements

File:

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

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

@@ -429 +429 @@
 
 }/*}}}*/
+void Element::GetNodesSidList(int* sidlist){/*{{{*/
+
+        _assert_(sidlist);
+        _assert_(nodes);
+        int numnodes = this->GetNumberOfNodes();
+        for(int i=0;i<numnodes;i++){
+                sidlist[i]=nodes[i]->Sid();
+        }
+}
+/*}}}*/
+void Element::GetNodesLidList(int* lidlist){/*{{{*/
+
+        _assert_(lidlist);
+        _assert_(nodes);
+        int numnodes = this->GetNumberOfNodes();
+        for(int i=0;i<numnodes;i++){
+                lidlist[i]=nodes[i]->Lid();
+        }
+}
+/*}}}*/
+void Element::GetVerticesCoordinates(IssmDouble** pxyz_list){/*{{{*/
+
+        int         numvertices = this->GetNumberOfVertices();
+        IssmDouble* xyz_list    = xNew<IssmDouble>(numvertices*3);
+        ::GetVerticesCoordinates(xyz_list,this->vertices,numvertices);
+
+        *pxyz_list = xyz_list;
+
+}/*}}}*/
 void Element::GetVerticesSidList(int* sidlist){/*{{{*/
 
@@ -439 +468 @@
         int numvertices = this->GetNumberOfVertices();
         for(int i=0;i<numvertices;i++) connectivity[i]=this->vertices[i]->Connectivity();
+}
+/*}}}*/
+void Element::InputChangeName(int original_enum,int new_enum){/*{{{*/
+        this->inputs->ChangeEnum(original_enum,new_enum);
+}
+/*}}}*/
+void Element::InputCreate(IssmDouble* vector,IoModel* iomodel,int M,int N,int vector_type,int vector_enum,int code){/*{{{*/
+
+        /*Intermediaries*/
+        int        i,t,row;
+        IssmDouble time;
+
+        /*Branch on type of vector: nodal or elementary: */
+        if(vector_type==1){ //nodal vector
+
+                int         numvertices = this->GetNumberOfVertices();
+                int        *vertexids   = xNew<int>(numvertices);
+                IssmDouble *values      = xNew<IssmDouble>(numvertices);
+
+                /*Recover vertices ids needed to initialize inputs*/
+                _assert_(iomodel->elements);
+                for(i=0;i<numvertices;i++){ 
+                        vertexids[i]=reCast<int>(iomodel->elements[numvertices*this->Sid()+i]); //ids for vertices are in the elements array from Matlab
+                }
+
+                /*Are we in transient or static? */
+                if(M==iomodel->numberofvertices){
+                        for(i=0;i<numvertices;i++) values[i]=vector[vertexids[i]-1];
+                        this->AddInput(vector_enum,values,P1Enum);
+                }
+                else if(M==iomodel->numberofvertices+1){
+                        /*create transient input: */
+                        IssmDouble* times = xNew<IssmDouble>(N);
+                        for(t=0;t<N;t++) times[t] = vector[(M-1)*N+t];
+                        TransientInput* transientinput=new TransientInput(vector_enum,times,N);
+                        for(t=0;t<N;t++){
+                                for(i=0;i<numvertices;i++) values[i]=vector[N*(vertexids[i]-1)+t];
+                                switch(this->ObjectEnum()){
+                                        case TriaEnum:  transientinput->AddTimeInput(new TriaInput( vector_enum,values,P1Enum)); break;
+                                        case PentaEnum: transientinput->AddTimeInput(new PentaInput(vector_enum,values,P1Enum)); break;
+                                        case TetraEnum: transientinput->AddTimeInput(new TetraInput(vector_enum,values,P1Enum)); break;
+                                        default: _error_("Not implemented yet");
+                                }
+                        }
+                        this->inputs->AddInput(transientinput);
+                        xDelete<IssmDouble>(times);
+                }
+                else _error_("nodal vector is either numberofvertices or numberofvertices+1 long. Field provided (" << EnumToStringx(vector_enum) << ") is " << M << " long");
+
+                xDelete<IssmDouble>(values);
+                xDelete<int>(vertexids);
+        }
+        else if(vector_type==2){ //element vector
+                /*Are we in transient or static? */
+                if(M==iomodel->numberofelements){
+                        if (code==5){ //boolean
+                                this->inputs->AddInput(new BoolInput(vector_enum,reCast<bool>(vector[this->Sid()])));
+                        }
+                        else if (code==6){ //integer
+                                this->inputs->AddInput(new IntInput(vector_enum,reCast<int>(vector[this->Sid()])));
+                        }
+                        else if (code==7){ //IssmDouble
+                                this->inputs->AddInput(new DoubleInput(vector_enum,vector[this->Sid()]));
+                        }
+                        else _error_("could not recognize nature of vector from code " << code);
+                }
+                else {
+                        _error_("transient element inputs not supported yet!");
+                }
+        }
+        else{
+                _error_("Cannot add input for vector type " << vector_type << " (not supported)");
+        }
+}/*}}}*/
+void Element::InputUpdateFromConstant(int constant, int name){/*{{{*/
+
+        /*Check that name is an element input*/
+        if(!IsInput(name)) return;
+
+        /*update input*/
+        this->inputs->AddInput(new IntInput(name,constant));
+}
+/*}}}*/
+void Element::InputUpdateFromConstant(IssmDouble constant, int name){/*{{{*/
+
+        /*Check that name is an element input*/
+        if(!IsInput(name)) return;
+
+        /*update input*/
+        this->inputs->AddInput(new DoubleInput(name,constant));
+}
+/*}}}*/
+void Element::InputUpdateFromConstant(bool constant, int name){/*{{{*/
+
+        /*Check that name is an element input*/
+        if(!IsInput(name)) return;
+
+        /*update input*/
+        this->inputs->AddInput(new BoolInput(name,constant));
 }
 /*}}}*/
@@ -535 +663 @@
 }
 /*}}}*/
+ElementVector* Element::NewElementVector(int approximation_enum){/*{{{*/
+        return new ElementVector(nodes,this->GetNumberOfNodes(),this->parameters,approximation_enum);
+}
+/*}}}*/
+ElementMatrix* Element::NewElementMatrix(int approximation_enum){/*{{{*/
+        return new ElementMatrix(nodes,this->GetNumberOfNodes(),this->parameters,approximation_enum);
+}
+/*}}}*/
+ElementMatrix* Element::NewElementMatrixCoupling(int number_nodes,int approximation_enum){/*{{{*/
+        return new ElementMatrix(nodes,number_nodes,this->parameters,approximation_enum);
+}
+/*}}}*/
 void Element::ResultInterpolation(int* pinterpolation,int* pnodesperelement,int output_enum){/*{{{*/