Changeset 17671


Ignore:
Timestamp:
04/08/14 13:15:59 (11 years ago)
Author:
Mathieu Morlighem
Message:

CHG: moving some stuff from Tria to Element

Location:
issm/trunk-jpl/src/c/classes/Elements
Files:
9 edited

Legend:

Unmodified
Added
Removed

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

    r17636 r17671  
    526526                lidlist[i]=nodes[i]->Lid();
    527527        }
     528}
     529/*}}}*/
     530void  Element::GetVectorFromInputs(Vector<IssmDouble>* vector,int input_enum){/*{{{*/
     531
     532        /*Fetch number vertices for this element and allocate arrays*/
     533        int numvertices = this->GetNumberOfVertices();
     534        int*        vertexpidlist = xNew<int>(numvertices);
     535        IssmDouble* values        = xNew<IssmDouble>(numvertices);
     536
     537        /*Fill in values*/
     538        this->GetVertexPidList(vertexpidlist);
     539        this->GetInputListOnNodesVelocity(values,input_enum);
     540        vector->SetValues(numvertices,vertexpidlist,values,INS_VAL);
     541
     542        /*Clean up*/
     543        xDelete<int>(vertexpidlist);
     544        xDelete<IssmDouble>(values);
     545
     546}
     547/*}}}*/
     548void  Element::GetVertexPidList(int* pidlist){/*{{{*/
     549
     550        int numvertices = this->GetNumberOfVertices();
     551        for(int i=0;i<numvertices;i++) pidlist[i]=vertices[i]->Pid();
     552
    528553}
    529554/*}}}*/

  • issm/trunk-jpl/src/c/classes/Elements/Element.h

    r17585 r17671  
    8686                void       GetNodesLidList(int* lidlist);
    8787                void       GetPhi(IssmDouble* phi, IssmDouble*  epsilon, IssmDouble viscosity);
     88                void       GetVectorFromInputs(Vector<IssmDouble>* vector, int name_enum);
     89                void         GetVertexPidList(int* pidlist);
    8890                void       GetVerticesCoordinates(IssmDouble** xyz_list);
    8991                void       GetVerticesSidList(int* sidlist);
    … …  
    220222
    221223                virtual void   InputScale(int enum_type,IssmDouble scale_factor)=0;
    222                 virtual void   GetVectorFromInputs(Vector<IssmDouble>* vector, int name_enum)=0;
    223224                virtual IssmDouble TimeAdapt()=0;
    224225                virtual void   PositiveDegreeDay(IssmDouble* pdds,IssmDouble* pds,IssmDouble signorm)=0;

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

    r17641 r17671  
    825825}
    826826/*}}}*/
    827 /*FUNCTION Penta::GetVertexPidList {{{*/
    828 void  Penta::GetVertexPidList(int* doflist){
    829 
    830         for(int i=0;i<6;i++) doflist[i]=vertices[i]->Pid();
    831 
    832 }
    833 /*}}}*/
    834827/*FUNCTION Penta::GetElementType {{{*/
    835828int Penta::GetElementType(){
    … …  
    966959
    967960        return tau_parameter;
    968 }
    969 /*}}}*/
    970 /*FUNCTION Penta::GetVectorFromInputs{{{*/
    971 void  Penta::GetVectorFromInputs(Vector<IssmDouble>* vector,int input_enum){
    972 
    973         int vertexpidlist[NUMVERTICES];
    974 
    975         /*Get out if this is not an element input*/
    976         if (!IsInput(input_enum)) return;
    977 
    978         /*Prepare index list*/
    979         this->GetVertexPidList(&vertexpidlist[0]);
    980 
    981         /*Get input (either in element or material)*/
    982         Input* input=inputs->GetInput(input_enum);
    983         if(!input) _error_("Input " << EnumToStringx(input_enum) << " not found in element");
    984 
    985         /*We found the enum.  Use its values to fill into the vector, using the vertices ids: */
    986         input->GetVectorFromInputs(vector,&vertexpidlist[0]);
    987961}
    988962/*}}}*/

  • issm/trunk-jpl/src/c/classes/Elements/Penta.h

    r17641 r17671  
    8080                IssmDouble GetYcoord(Gauss* gauss);
    8181                IssmDouble GetZcoord(Gauss* gauss);
    82                 void   GetVectorFromInputs(Vector<IssmDouble>* vector,int name_enum);
    8382                void   GetVerticesCoordinatesBase(IssmDouble** pxyz_list);
    8483                void   GetVerticesCoordinatesTop(IssmDouble** pxyz_list);
    … …  
    178177                ElementMatrix* CreateBasalMassMatrix(void);
    179178                void           GetAreaCoordinates(IssmDouble *area_coordinates,IssmDouble* xyz_zero,IssmDouble* xyz_list,int numpoints);
    180 
    181                 void             GetVertexPidList(int* doflist);
    182179                int            GetElementType(void);
    183180                void           GetInputValue(IssmDouble* pvalue,Node* node,int enumtype);

  • issm/trunk-jpl/src/c/classes/Elements/Seg.h

    r17585 r17671  
    128128
    129129                void        GetSolutionFromInputsOneDof(Vector<IssmDouble>* solution,int enum_type){_error_("not implemented yet");};
    130                 void        GetVectorFromInputs(Vector<IssmDouble>* vector, int name_enum){_error_("not implemented yet");};
    131130                void        InputDepthAverageAtBase(int enum_type,int average_enum_type){_error_("not implemented yet");};
    132131                void        InputDuplicate(int original_enum,int new_enum){_error_("not implemented yet");};

  • issm/trunk-jpl/src/c/classes/Elements/Tetra.cpp

    r17610 r17671  
    844844                        tetra_node_ids[8]=iomodel->nodecounter+iomodel->numberofvertices+iomodel->numberofelements+iomodel->elements[4*index+3];
    845845                        break;
     846                case TaylorHoodEnum:
    846847                        numnodes        = 14;
    847                         tetra_node_ids   = xNew<int>(numnodes);
     848                        tetra_node_ids  = xNew<int>(numnodes);
    848849                        tetra_node_ids[0]=iomodel->nodecounter+iomodel->elements[4*index+0];
    849850                        tetra_node_ids[1]=iomodel->nodecounter+iomodel->elements[4*index+1];

  • issm/trunk-jpl/src/c/classes/Elements/Tetra.h

    r17585 r17671  
    133133
    134134                void        GetSolutionFromInputsOneDof(Vector<IssmDouble>* solution,int enum_type){_error_("not implemented yet");};
    135                 void        GetVectorFromInputs(Vector<IssmDouble>* vector, int name_enum){_error_("not implemented yet");};
    136135                void        InputDepthAverageAtBase(int enum_type,int average_enum_type){_error_("not implemented yet");};
    137136                void        InputDuplicate(int original_enum,int new_enum){_error_("not implemented yet");};

  • issm/trunk-jpl/src/c/classes/Elements/Tria.cpp

    r17670 r17671  
    834834
    835835}/*}}}*/
    836 /*FUNCTION Tria::GetVertexPidList {{{*/
    837 void  Tria::GetVertexPidList(int* doflist){
    838 
    839         for(int i=0;i<3;i++) doflist[i]=vertices[i]->Pid();
    840 
    841 }
    842 /*}}}*/
    843 /*FUNCTION Tria::GetVectorFromInputs{{{*/
    844 void  Tria::GetVectorFromInputs(Vector<IssmDouble>* vector,int input_enum){
    845 
    846         int vertexpidlist[NUMVERTICES];
    847 
    848         /*Get out if this is not an element input*/
    849         if(!IsInput(input_enum)) return;
    850 
    851         /*Prepare index list*/
    852         this->GetVertexPidList(&vertexpidlist[0]);
    853 
    854         /*Get input (either in element or material)*/
    855         Input* input=inputs->GetInput(input_enum);
    856         if(!input) _error_("Input " << EnumToStringx(input_enum) << " not found in element");
    857 
    858         /*We found the enum.  Use its values to fill into the vector, using the vertices ids: */
    859         input->GetVectorFromInputs(vector,&vertexpidlist[0]);
    860 }
    861 /*}}}*/
    862836/*FUNCTION Tria::GetXcoord {{{*/
    863837IssmDouble Tria::GetXcoord(Gauss* gauss){

  • issm/trunk-jpl/src/c/classes/Elements/Tria.h

    r17585 r17671  
    8282                int         NumberofNodesPressure(void);
    8383                void        GetSolutionFromInputsOneDof(Vector<IssmDouble>* solution,int enum_type);
    84                 void        GetVectorFromInputs(Vector<IssmDouble>* vector, int name_enum);
    8584                void        GetVerticesCoordinatesBase(IssmDouble** pxyz_list);
    8685                void        GetVerticesCoordinatesTop(IssmDouble** pxyz_list);
    … …  
    180179                void           GetAreaCoordinates(IssmDouble *area_coordinates,IssmDouble* xyz_zero,IssmDouble* xyz_list,int numpoints);
    181180                int            GetElementType(void);
    182 
    183                 void             GetVertexPidList(int* doflist);
    184181                IssmDouble     GetXcoord(Gauss* gauss);
    185182                IssmDouble     GetYcoord(Gauss* gauss);
Note: See TracChangeset for help on using the changeset viewer.