Changeset 16698


Ignore:
Timestamp:
11/09/13 21:20:38 (11 years ago)
Author:
seroussi
Message:

CHG: working on UpdateInputFromSolution

Location:
issm/trunk-jpl/src/c
Files:
8 edited

Legend:

Unmodified
Added
Removed

  • TabularUnified issm/trunk-jpl/src/c/analyses/StressbalanceAnalysis.cpp

    r16697 r16698  
    952952void StressbalanceAnalysis::InputUpdateFromSolutionHoriz(IssmDouble* solution,Element* element){/*{{{*/
    953953
    954         int        i;
    955         IssmDouble rho_ice,g;
    956         int*       doflist=NULL;
     954        int         i,meshtype;
     955        IssmDouble  rho_ice,g;
     956        int*        doflist=NULL;
     957        IssmDouble* xyz_list=NULL;
    957958
    958959        /*Fetch number of nodes and dof for this finite element*/
    … …  
    969970        IssmDouble* pressure  = xNew<IssmDouble>(numnodes);
    970971        IssmDouble* thickness = xNew<IssmDouble>(numnodes);
     972        IssmDouble* surface   = xNew<IssmDouble>(numnodes);
    971973
    972974        /*Use the dof list to index into the solution vector: */
    … …  
    975977        /*Transform solution in Cartesian Space*/
    976978        element->TransformSolutionCoord(&values[0],XYEnum);
     979        element->FindParam(&meshtype,MeshTypeEnum);
    977980
    978981        /*Ok, we have vx and vy in values, fill in vx and vy arrays: */
    … …  
    994997        rho_ice = element->GetMaterialParameter(MaterialsRhoIceEnum);
    995998        g       = element->GetMaterialParameter(ConstantsGEnum);
    996         element->GetInputListOnNodes(&thickness[0],ThicknessEnum);
    997         for(i=0;i<numnodes;i++) pressure[i]=rho_ice*g*thickness[i];
     999        switch(meshtype){
     1000                case Mesh2DhorizontalEnum:
     1001                        element->GetInputListOnNodes(&thickness[0],ThicknessEnum);
     1002                        for(i=0;i<numnodes;i++) pressure[i]=rho_ice*g*thickness[i];
     1003                        break;
     1004                case Mesh3DEnum:   
     1005                        element->GetVerticesCoordinates(&xyz_list);
     1006                        element->GetInputListOnNodes(&surface[0],SurfaceEnum);
     1007                        for(i=0;i<numnodes;i++) pressure[i]=rho_ice*g*(surface[i]-xyz_list[i*3+2]);
     1008                        break;
     1009                default: _error_("mesh "<<EnumToStringx(meshtype)<<" not supported yet");
     1010        }
    9981011
    9991012        /*Now, we have to move the previous Vx and Vy inputs  to old
    … …  
    10111024        /*Free ressources:*/
    10121025        xDelete<IssmDouble>(thickness);
     1026        xDelete<IssmDouble>(surface);
    10131027        xDelete<IssmDouble>(pressure);
    10141028        xDelete<IssmDouble>(vel);
    … …  
    10171031        xDelete<IssmDouble>(vx);
    10181032        xDelete<IssmDouble>(values);
     1033        xDelete<IssmDouble>(xyz_list);
    10191034        xDelete<int>(doflist);
    10201035

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

    r16697 r16698  
    7171                virtual void   GetInputValue(int* pvalue,int enum_type)=0;
    7272                virtual void   GetInputValue(IssmDouble* pvalue,int enum_type)=0;
     73                virtual void   GetVerticesCoordinates(IssmDouble** xyz_list)=0;
    7374                virtual void   GetMaterialInputValue(IssmDouble* pvalue,Node* node,int enumtype)=0;
    7475

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

    r16697 r16698  
    121121
    122122/*Other*/
     123/*FUNCTION Penta::AddInput{{{*/
     124void  Penta::AddInput(int input_enum,IssmDouble* values, int interpolation_enum){
     125
     126        /*Call inputs method*/
     127        this->inputs->AddInput(new PentaInput(input_enum,values,interpolation_enum));
     128}
     129/*}}}*/
    123130/*FUNCTION Penta::BedNormal {{{*/
    124131void Penta::BedNormal(IssmDouble* bed_normal, IssmDouble xyz_list[3][3]){
    … …  
    234241
    235242        /* Get node coordinates and dof list: */
    236         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     243        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    237244        for(i=0;i<3;i++) for(j=0;j<3;j++) xyz_list_tria[i][j]=xyz_list[i][j];
    238245
    … …  
    302309
    303310        /* Get node coordinates and dof list: */
    304         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     311        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    305312
    306313        /*Retrieve all inputs we will be needing: */
    … …  
    12731280
    12741281        /*Get xyz list: */
    1275         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     1282        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    12761283        xmin=xyz_list[0][0]; xmax=xyz_list[0][0];
    12771284        ymin=xyz_list[0][1]; ymax=xyz_list[0][1];
    … …  
    14681475        if(!input) _error_("Input " << EnumToStringx(inputenum) << " not found in element");
    14691476        input->GetInputValue(pvalue);
     1477
     1478}/*}}}*/
     1479/*FUNCTION Penta::GetVerticesCoordinates(IssmDouble** pxyz_list){{{*/
     1480void Penta::GetVerticesCoordinates(IssmDouble** pxyz_list){
     1481
     1482        IssmDouble* xyz_list = xNew<IssmDouble>(NUMVERTICES*3);
     1483        ::GetVerticesCoordinates(xyz_list,this->vertices,NUMVERTICES);
     1484
     1485        /*Assign output pointer*/
     1486        *pxyz_list = xyz_list;
    14701487
    14711488}/*}}}*/
    … …  
    16511668        IssmDouble z_list[NUMVERTICES];
    16521669
    1653         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     1670        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    16541671        for(i=0;i<NUMVERTICES;i++) z_list[i]=xyz_list[i][2];
    16551672        PentaRef::GetInputValue(&z,z_list,gauss);
    … …  
    19551972
    19561973                /*Step2: Create element thickness input*/
    1957                 GetVerticesCoordinates(&xyz_list[0][0],penta->vertices,NUMVERTICES);
     1974                ::GetVerticesCoordinates(&xyz_list[0][0],penta->vertices,NUMVERTICES);
    19581975                for(i=0;i<3;i++){
    19591976                        Helem_list[i]=xyz_list[i+3][2]-xyz_list[i][2];
    … …  
    32153232
    32163233        /* Get node coordinates and dof list: */
    3217         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     3234        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    32183235
    32193236        minx=xyz_list[0][0];
    … …  
    34703487
    34713488        /*Retrieve all inputs and parameters*/
    3472         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     3489        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    34733490        Input* vx_input=inputs->GetInput(VxEnum); _assert_(vx_input);
    34743491        Input* vy_input=inputs->GetInput(VyEnum); _assert_(vy_input);
    … …  
    35233540        if(NoIceInElement())return 0;
    35243541
    3525         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     3542        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    35263543
    35273544        /*First calculate the area of the base (cross section triangle)
    … …  
    35493566        rho_ice=matpar->GetRhoIce();
    35503567        rho_water=matpar->GetRhoWater();
    3551         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     3568        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    35523569
    35533570        /*First calculate the area of the base (cross section triangle)
    … …  
    37703787        if(NoIceInElement() || !IsOnSurface()) return 0.;
    37713788
    3772         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     3789        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    37733790
    37743791        /*First calculate the area of the base (cross section triangle)
    … …  
    38453862
    38463863        /*Retrieve all inputs and parameters*/
    3847         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     3864        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    38483865        rho_water=matpar->GetRhoWater();
    38493866        rho_ice=matpar->GetRhoIce();
    … …  
    39944011        rho_ice=matpar->GetRhoIce();
    39954012        heatcapacity=matpar->GetHeatCapacity();
    3996         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     4013        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    39974014        for(i=0;i<NUMVERTICES2D;i++) for(j=0;j<3;j++) xyz_list_tria[i][j]=xyz_list[i][j];
    39984015
    … …  
    40764093
    40774094        /*Retrieve all inputs and parameters*/
    4078         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     4095        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    40794096        rho_water=matpar->GetRhoWater();
    40804097        rho_ice=matpar->GetRhoIce();
    … …  
    42234240        rho_ice=matpar->GetRhoIce();
    42244241        heatcapacity=matpar->GetHeatCapacity();
    4225         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     4242        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    42264243        for(i=0;i<NUMVERTICES2D;i++) for(j=0;j<3;j++) xyz_list_tria[i][j]=xyz_list[i][j];
    42274244
    … …  
    42924309
    42934310        /*Retrieve all inputs and parameters*/
    4294         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     4311        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    42954312        rho_ice=matpar->GetRhoIce();
    42964313        heatcapacity=matpar->GetHeatCapacity();
    … …  
    43834400
    43844401        /*Retrieve all inputs and parameters*/
    4385         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     4402        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    43864403        for(i=0;i<NUMVERTICES2D;i++) for(j=0;j<3;j++) xyz_list_tria[i][j]=xyz_list[i][j];
    43874404        mixed_layer_capacity=matpar->GetMixedLayerCapacity();
    … …  
    44464463
    44474464        /*Retrieve all inputs and parameters*/
    4448         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     4465        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    44494466        for(i=0;i<NUMVERTICES2D;i++) for(j=0;j<3;j++) xyz_list_tria[i][j]=xyz_list[i][j];
    44504467        parameters->FindParam(&analysis_type,AnalysisTypeEnum);
    … …  
    45654582
    45664583        /*Retrieve all inputs and parameters*/
    4567         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     4584        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    45684585        rho_ice=matpar->GetRhoIce();
    45694586        heatcapacity=matpar->GetHeatCapacity();
    … …  
    46494666
    46504667        /*Retrieve all inputs and parameters*/
    4651         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     4668        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    46524669        for(i=0;i<NUMVERTICES2D;i++) for(j=0;j<3;j++) xyz_list_tria[i][j]=xyz_list[i][j];
    46534670        mixed_layer_capacity=matpar->GetMixedLayerCapacity();
    … …  
    47084725
    47094726        /*Retrieve all inputs and parameters*/
    4710         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     4727        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    47114728        for(i=0;i<NUMVERTICES2D;i++) for(j=0;j<3;j++) xyz_list_tria[i][j]=xyz_list[i][j];
    47124729        parameters->FindParam(&analysis_type,AnalysisTypeEnum);
    … …  
    47884805
    47894806        /*Get all inputs and parameters*/
    4790         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     4807        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    47914808        Input* surface_input=inputs->GetInput(SurfaceEnum); _assert_(surface_input);
    47924809
    … …  
    48564873
    48574874        /*Get all inputs and parameters*/
    4858         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     4875        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    48594876        GetInputListOnVertices(&pressure[0],PressureEnum);
    48604877        Input* surface_input=inputs->GetInput(SurfaceEnum); _assert_(surface_input);
    … …  
    50295046        Input* enthalpy_input=inputs->GetInput(EnthalpyEnum); _assert_(enthalpy_input);
    50305047
    5031         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     5048        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    50325049        for(i=0;i<NUMVERTICES2D;i++) for(j=0;j<3;j++) xyz_list_tria[i][j]=xyz_list[i][j];
    50335050
    … …  
    53465363
    53475364        /*Retrieve all inputs and parameters*/
    5348         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     5365        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    53495366        Input* vx_input=inputs->GetInput(VxEnum);       _assert_(vx_input);
    53505367        Input* vy_input=inputs->GetInput(VyEnum);       _assert_(vy_input);
    … …  
    54165433
    54175434        /*Retrieve all inputs and parameters*/
    5418         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     5435        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    54195436        Input* vx_input=inputs->GetInput(VxEnum);       _assert_(vx_input);
    54205437        Input* vy_input=inputs->GetInput(VyEnum);       _assert_(vy_input);
    … …  
    55305547
    55315548        /*Retrieve all inputs and parameters*/
    5532         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     5549        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    55335550        this->parameters->FindParam(&num_responses,InversionNumCostFunctionsEnum);
    55345551        this->parameters->FindParam(&responses,NULL,InversionCostFunctionsEnum);
    … …  
    57265743
    57275744        /*Retrieve all inputs and parameters*/
    5728         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     5745        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    57295746        this->parameters->FindParam(&num_responses,InversionNumCostFunctionsEnum);
    57305747        this->parameters->FindParam(&responses,NULL,InversionCostFunctionsEnum);
    … …  
    60336050        parameters->FindParam(&analysis_type,AnalysisTypeEnum);
    60346051        GradientIndexing(&vertexpidlist[0],control_index);
    6035         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     6052        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    60366053        for(i=0;i<NUMVERTICES2D;i++) for(j=0;j<3;j++) xyz_list_tria[i][j]=xyz_list[i][j];
    60376054        Input* adjointx_input=inputs->GetInput(AdjointxEnum);               _assert_(adjointx_input);
    … …  
    61056122        /*Retrieve all inputs and parameters*/
    61066123        parameters->FindParam(&analysis_type,AnalysisTypeEnum);
    6107         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     6124        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    61086125        for(i=0;i<NUMVERTICES2D;i++) for(j=0;j<3;j++) xyz_list_tria[i][j]=xyz_list[i][j];
    61096126        GradientIndexing(&vertexpidlist[0],control_index);
    … …  
    69236940
    69246941        /* Get node coordinates and dof list: */
    6925         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     6942        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    69266943        this->parameters->FindParam(&viscosity_overshoot,StressbalanceViscosityOvershootEnum);
    69276944        Input* vx_input=inputs->GetInput(VxEnum);       _assert_(vx_input);
    … …  
    70107027
    70117028        /*retrieve inputs :*/
    7012         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     7029        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    70137030        for(i=0;i<NUMVERTICES2D;i++) for(j=0;j<3;j++) xyz_list_tria[i][j]=xyz_list[i][j];
    70147031        parameters->FindParam(&analysis_type,AnalysisTypeEnum);
    … …  
    71297146
    71307147        /* Get node coordinates and dof list: */
    7131         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     7148        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    71327149        parameters->FindParam(&FSreconditioning,StressbalanceFSreconditioningEnum);
    71337150        Input* vx_input=inputs->GetInput(VxEnum);       _assert_(vx_input);
    … …  
    72437260
    72447261        /*Retrieve all inputs and parameters*/
    7245         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     7262        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    72467263        parameters->FindParam(&analysis_type,AnalysisTypeEnum);
    72477264        parameters->FindParam(&FSreconditioning,StressbalanceFSreconditioningEnum);
    … …  
    75437560
    75447561        /*Retrieve all inputs and parameters*/
    7545         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     7562        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    75467563        this->parameters->FindParam(&viscosity_overshoot,StressbalanceViscosityOvershootEnum);
    75477564        Input* vx_input=inputs->GetInput(VxEnum);       _assert_(vx_input);
    … …  
    76947711
    76957712        /*Retrieve all inputs and parameters*/
    7696         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     7713        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    76977714        Input* vx_input=inputs->GetInput(VxEnum);        _assert_(vx_input);
    76987715        Input* vy_input=inputs->GetInput(VyEnum);        _assert_(vy_input);
    … …  
    77907807        /*Retrieve all inputs and parameters*/
    77917808        inputs->GetInputValue(&approximation,ApproximationEnum);
    7792         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     7809        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    77937810        this->parameters->FindParam(&viscosity_overshoot,StressbalanceViscosityOvershootEnum);
    77947811        Input* vx_input=inputs->GetInput(VxEnum);       _assert_(vx_input);
    … …  
    78647881
    78657882        /*Retrieve all inputs and parameters*/
    7866         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     7883        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    78677884        for(i=0;i<NUMVERTICES2D;i++) for(j=0;j<3;j++) xyz_list_tria[i][j]=xyz_list[i][j];
    78687885        parameters->FindParam(&migration_style,GroundinglineMigrationEnum);
    … …  
    79888005
    79898006        /*Retrieve all inputs and parameters*/
    7990         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     8007        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    79918008        parameters->FindParam(&FSreconditioning,StressbalanceFSreconditioningEnum);
    79928009        Input* vx_input=inputs->GetInput(VxEnum); _assert_(vx_input);
    … …  
    81008117
    81018118        /*Retrieve all inputs and parameters*/
    8102         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     8119        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    81038120        parameters->FindParam(&FSreconditioning,StressbalanceFSreconditioningEnum);
    81048121        Input* vx_input=inputs->GetInput(VxEnum); _assert_(vx_input);
    … …  
    81718188
    81728189        /*Retrieve all inputs and parameters*/
    8173         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     8190        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    81748191        parameters->FindParam(&analysis_type,AnalysisTypeEnum);
    81758192        Input* vx_input=inputs->GetInput(VxEnum); _assert_(vx_input);
    … …  
    82428259
    82438260        /*Retrieve all inputs and parameters*/
    8244         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     8261        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    82458262
    82468263        /* Start  looping on the number of gaussian points: */
    … …  
    82888305
    82898306        /*Retrieve all inputs and parameters*/
    8290         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     8307        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    82918308        for(i=0;i<NUMVERTICES2D;i++) for(j=0;j<3;j++) xyz_list_tria[i][j]=xyz_list[i+3][j];
    82928309        SurfaceNormal(&surface_normal[0],xyz_list_tria);
    … …  
    83558372
    83568373        /*Retrieve all inputs and parameters*/
    8357         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     8374        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    83588375        this->parameters->FindParam(&FSreconditioning,StressbalanceFSreconditioningEnum);
    83598376        Input* vx_input=inputs->GetInput(VxEnum);               _assert_(vx_input);
    … …  
    84268443
    84278444        /*Retrieve all inputs and parameters*/
    8428         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     8445        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    84298446        parameters->FindParam(&analysis_type,AnalysisTypeEnum);
    84308447        this->parameters->FindParam(&FSreconditioning,StressbalanceFSreconditioningEnum);
    … …  
    85178534
    85188535        /*Retrieve all inputs and parameters*/
    8519         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     8536        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    85208537        this->parameters->FindParam(&FSreconditioning,StressbalanceFSreconditioningEnum);
    85218538        Input* vx_input=inputs->GetInput(VxEnum);       _assert_(vx_input);
    … …  
    85918608
    85928609        /*Retrieve all inputs and parameters*/
    8593         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     8610        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    85948611        parameters->FindParam(&analysis_type,AnalysisTypeEnum);
    85958612        this->parameters->FindParam(&FSreconditioning,StressbalanceFSreconditioningEnum);
    … …  
    87568773
    87578774        /*Retrieve all inputs and parameters*/
    8758         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     8775        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    87598776        rho_ice=matpar->GetRhoIce();
    87608777        gravity=matpar->GetG();
    … …  
    88828899
    88838900        /*Retrieve all inputs and parameters*/
    8884         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     8901        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    88858902        Input* thickness_input=inputs->GetInput(ThicknessEnum); _assert_(thickness_input);
    88868903        Input* surface_input=inputs->GetInput(SurfaceEnum);     _assert_(surface_input);
    … …  
    89508967
    89518968        /*Initialize Element vector and other vectors*/
    8952         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     8969        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    89538970        ElementVector* pe    = new ElementVector(nodes,numnodes,this->parameters,HOApproximationEnum);
    89548971        IssmDouble*    basis = xNew<IssmDouble>(numnodes);
    … …  
    90659082
    90669083        /*Retrieve all inputs and parameters*/
    9067         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     9084        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    90689085        rho_water=matpar->GetRhoWater();
    90699086        rho_ice  =matpar->GetRhoIce();
    … …  
    91319148        gravity=matpar->GetG();
    91329149        B=material->GetB();
    9133         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     9150        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    91349151        Input* vx_input=inputs->GetInput(VxEnum); _assert_(vx_input);
    91359152        Input* vy_input=inputs->GetInput(VyEnum); _assert_(vy_input);
    … …  
    92239240
    92249241        /*Retrieve all inputs and parameters*/
    9225         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     9242        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    92269243        Input* loadingforcex_input=inputs->GetInput(LoadingforceXEnum);  _assert_(loadingforcex_input);
    92279244        Input* loadingforcey_input=inputs->GetInput(LoadingforceYEnum);  _assert_(loadingforcey_input);
    … …  
    92989315        rho_water=matpar->GetRhoWater();
    92999316        gravity=matpar->GetG();
    9300         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     9317        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    93019318        Input* bed_input=inputs->GetInput(BedEnum); _assert_(bed_input);
    93029319        Input* vx_input=inputs->GetInput(VxEnum);   _assert_(vx_input);
    … …  
    93799396
    93809397        /*Retrieve all inputs and parameters*/
    9381         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     9398        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    93829399        inputs->GetInputValue(&approximation,ApproximationEnum);
    93839400        Input* vx_input=inputs->GetInput(VxEnum); _assert_(vx_input);
    … …  
    94389455
    94399456        /*Retrieve all inputs and parameters*/
    9440         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     9457        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    94419458        for(i=0;i<NUMVERTICES2D;i++) for(j=0;j<3;j++) xyz_list_tria[i][j]=xyz_list[i][j];
    94429459        inputs->GetInputValue(&approximation,ApproximationEnum);
    … …  
    95499566
    95509567        /*Retrieve all inputs and parameters*/
    9551         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     9568        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    95529569        Input* vx_input=inputs->GetInput(VxEnum);       _assert_(vx_input);
    95539570        Input* vy_input=inputs->GetInput(VyEnum);       _assert_(vy_input);
    … …  
    96259642
    96269643        /*Retrieve all inputs and parameters*/
    9627         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     9644        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    96289645        Input* vx_input=inputs->GetInput(VxEnum); _assert_(vx_input);
    96299646        Input* vy_input=inputs->GetInput(VyEnum); _assert_(vy_input);
    … …  
    98439860
    98449861                /*Get node data: */
    9845                 GetVerticesCoordinates(&xyz_list[0][0],penta->vertices,NUMVERTICES);
     9862                ::GetVerticesCoordinates(&xyz_list[0][0],penta->vertices,NUMVERTICES);
    98469863
    98479864                /*Now Compute vel*/
    … …  
    99069923
    99079924        /*Get node data: */
    9908         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     9925        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    99099926
    99109927        /*Use the dof list to index into the solution vector: */
    … …  
    1000110018
    1000210019        /*Get node data: */
    10003         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     10020        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    1000410021
    1000510022        /*Use the dof list to index into the solution vector: */
    … …  
    1011710134
    1011810135                /*Get node data: */
    10119                 GetVerticesCoordinates(&xyz_list[0][0],penta->vertices,NUMVERTICES);
     10136                ::GetVerticesCoordinates(&xyz_list[0][0],penta->vertices,NUMVERTICES);
    1012010137
    1012110138                /*Now Compute vel*/
    … …  
    1017710194        /*Transform solution in Cartesian Space*/
    1017810195        ::TransformSolutionCoord(&values[0],nodes,NUMVERTICES,XYEnum);
    10179         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     10196        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    1018010197
    1018110198        /*Ok, we have vx and vy in values, fill in vx and vy arrays: */
    … …  
    1025810275
    1025910276        /*Get node data: */
    10260         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     10277        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    1026110278
    1026210279        /*Use the dof list to index into the solution vector: */
    … …  
    1034510362
    1034610363        /*Get node data: */
    10347         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     10364        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    1034810365
    1034910366        /*Use the dof list to index into the solution vector: */
    … …  
    1041710434        /*Get dof list and vertices coordinates: */
    1041810435        GetDofList(&doflist,NoneApproximationEnum,GsetEnum);
    10419         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     10436        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    1042010437
    1042110438        /*Use the dof list to index into the solution vector vz: */

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

    r16697 r16698  
    9595                IssmDouble GetZcoord(GaussPenta* gauss);
    9696                void   GetVectorFromInputs(Vector<IssmDouble>* vector,int name_enum);
     97                void   GetVerticesCoordinates(IssmDouble** pxyz_list);
    9798
    9899                int    Sid();
    … …  
    180181                /*}}}*/
    181182                /*Penta specific routines:{{{*/
    182                 void           AddInput(int input_enum, IssmDouble* values, int interpolation_enum){_error_("not implemented yet");};
     183                void           AddInput(int input_enum, IssmDouble* values, int interpolation_enum);
    183184                void             BedNormal(IssmDouble* bed_normal, IssmDouble xyz_list[3][3]);
    184185                ElementMatrix* CreateBasalMassMatrix(void);

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

    r16697 r16698  
    160160
    161161        /*Get xyz list: */
    162         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     162        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    163163        x1=xyz_list[0][0]; y1=xyz_list[0][1];
    164164        x2=xyz_list[1][0]; y2=xyz_list[1][1];
    … …  
    186186
    187187        /*Retrieve all inputs and parameters*/
    188         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     188        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    189189        this->parameters->FindParam(&dt,TimesteppingTimeStepEnum);
    190190        this->parameters->FindParam(&stabilization,MasstransportStabilizationEnum);
    … …  
    264264
    265265        /*Retrieve all inputs and parameters*/
    266         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     266        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    267267        this->parameters->FindParam(&dt,TimesteppingTimeStepEnum);
    268268        this->parameters->FindParam(&stabilization,MasstransportStabilizationEnum);
    … …  
    338338
    339339        /*Retrieve all inputs and parameters*/
    340         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     340        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    341341
    342342        /* Start looping on the number of gaussian points: */
    … …  
    381381
    382382        /*Retrieve all inputs and parameters*/
    383         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     383        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    384384        this->parameters->FindParam(&input_enum,InputToL2ProjectEnum);
    385385        switch(input_enum){
    … …  
    430430
    431431        /*Retrieve all inputs and parameters*/
    432         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     432        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    433433        this->parameters->FindParam(&dt,TimesteppingTimeStepEnum);
    434434        Input* vy_input     = inputs->GetInput(VyEnum);                         _assert_(vy_input);
    … …  
    475475
    476476        /*Retrieve all inputs and parameters*/
    477         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     477        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    478478        this->parameters->FindParam(&dt,TimesteppingTimeStepEnum);
    479479        Input* vy_input            = inputs->GetInput(VyEnum);                         _assert_(vy_input);

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

    r16697 r16698  
    9292                void        GetNodesLidList(int* lidlist){_error_("not implemented yet");};
    9393                int         GetNumberOfNodes(void){_error_("not implemented yet");};
     94                void        GetVerticesCoordinates(IssmDouble** pxyz_list){_error_("not implemented yet");};
    9495                int         Sid(){_error_("not implemented yet");};
    9596                void        InputChangeName(int input_enum, int enum_type_old){_error_("not implemented yet");};

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

    r16697 r16698  
    342342
    343343        /*Retrieve all inputs and parameters*/
    344         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     344        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    345345
    346346        /* Start looping on the number of gaussian points: */
    … …  
    592592
    593593        /*Retrieve all inputs and parameters*/
    594         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     594        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    595595        this->parameters->FindParam(&input_enum,InputToL2ProjectEnum);
    596596        switch(input_enum){
    … …  
    685685
    686686        /* Get node coordinates and dof list: */
    687         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     687        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    688688
    689689        /*Retrieve all inputs we will be needing: */
    … …  
    907907
    908908        /*Get xyz list: */
    909         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     909        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    910910        x1=xyz_list[0][0]; y1=xyz_list[0][1];
    911911        x2=xyz_list[1][0]; y2=xyz_list[1][1];
    … …  
    27002700        if(NoIceInElement())return 0;
    27012701
    2702         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     2702        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    27032703
    27042704        for(int i=0;i<3;i++){
    … …  
    27642764
    27652765        /* Get node coordinates and dof list: */
    2766         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     2766        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    27672767
    27682768        minx=xyz_list[0][0];
    … …  
    30293029        rho_ice=matpar->GetRhoIce();
    30303030        rho_water=matpar->GetRhoWater();
    3031         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     3031        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    30323032
    30333033        /*First calculate the area of the base (cross section triangle)
    … …  
    30683068
    30693069        /*Get xyz list: */
    3070         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     3070        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    30713071
    30723072        /*get area coordinates of 0 and 1 locations: */
    … …  
    31353135
    31363136        /*Get xyz list: */
    3137         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     3137        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    31383138
    31393139        /*get area coordinates of 0 and 1 locations: */
    … …  
    33293329   if(NoIceInElement())return 0;
    33303330
    3331         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     3331        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    33323332
    33333333        /*First calculate the area of the base (cross section triangle)
    … …  
    34293429
    34303430        /*figure out gravity center of our element: */
    3431         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     3431        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    34323432        x0=(xyz_list[0][0]+xyz_list[1][0]+xyz_list[2][0])/3.0;
    34333433        y0=(xyz_list[0][1]+xyz_list[1][1]+xyz_list[2][1])/3.0;
    … …  
    35233523
    35243524        /*Retrieve all inputs and parameters*/
    3525         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     3525        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    35263526        parameters->FindParam(&FSreconditioning,StressbalanceFSreconditioningEnum);
    35273527        Input* vx_input=inputs->GetInput(VxEnum); _assert_(vx_input);
    … …  
    35943594
    35953595        /*Retrieve all inputs and parameters*/
    3596         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     3596        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    35973597        parameters->FindParam(&analysis_type,AnalysisTypeEnum);
    35983598        Input* vx_input=inputs->GetInput(VxEnum); _assert_(vx_input);
    … …  
    36703670
    36713671        /*Retrieve all inputs and parameters*/
    3672         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     3672        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    36733673        Input* thickness_input=inputs->GetInput(ThicknessEnum); _assert_(thickness_input);
    36743674        Input* vx_input=inputs->GetInput(VxEnum);               _assert_(vx_input);
    … …  
    37443744
    37453745        /*Retrieve all inputs and parameters*/
    3746         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     3746        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    37473747        parameters->FindParam(&migration_style,GroundinglineMigrationEnum);
    37483748        Input* surface_input=inputs->GetInput(SurfaceEnum);       _assert_(surface_input);
    … …  
    38863886
    38873887        /*Initialize Element matrix and vectors*/
    3888         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     3888        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    38893889        ElementVector* pe     = new ElementVector(nodes,vnumnodes+pnumnodes,this->parameters,FSvelocityEnum);
    38903890        IssmDouble*    vbasis = xNew<IssmDouble>(vnumnodes);
    … …  
    39583958
    39593959        /*Retrieve all inputs and parameters*/
    3960         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     3960        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    39613961        Input* loadingforcex_input=inputs->GetInput(LoadingforceXEnum);  _assert_(loadingforcex_input);
    39623962        Input* loadingforcey_input=inputs->GetInput(LoadingforceYEnum);  _assert_(loadingforcey_input);
    … …  
    40254025        rho_water=matpar->GetRhoWater();
    40264026        gravity=matpar->GetG();
    4027         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     4027        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    40284028        Input* bed_input=inputs->GetInput(BedEnum); _assert_(bed_input);
    40294029
    … …  
    40964096
    40974097        /*Retrieve all inputs and parameters*/
    4098         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     4098        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    40994099        Input* thickness_input=inputs->GetInput(ThicknessEnum); _assert_(thickness_input);
    41004100        Input* surface_input=inputs->GetInput(SurfaceEnum);     _assert_(surface_input);
    … …  
    41654165
    41664166        /*Initialize Element vector and other vectors*/
    4167         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     4167        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    41684168        ElementVector* pe    = new ElementVector(nodes,numnodes,this->parameters,SSAApproximationEnum);
    41694169        IssmDouble*    basis = xNew<IssmDouble>(numnodes);
    … …  
    42844284
    42854285        /*Retrieve all inputs and parameters*/
    4286         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     4286        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    42874287        Input* vx_input=inputs->GetInput(VxEnum);       _assert_(vx_input);
    42884288        Input* vy_input=inputs->GetInput(VyEnum);       _assert_(vy_input);
    … …  
    43354335        IssmDouble y_list[NUMVERTICES];
    43364336
    4337         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     4337        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    43384338        for(int i=0;i<NUMVERTICES;i++) y_list[i]=xyz_list[i][1];
    43394339        TriaRef::GetInputValue(&y,&y_list[0],gauss,P1Enum);
    … …  
    45834583
    45844584        /*Retrieve all inputs we will be needing: */
    4585         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     4585        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    45864586        Input* weights_input              = inputs->GetInput(InversionCostFunctionsCoefficientsEnum);   _assert_(weights_input);
    45874587        Input* thickness_input            = inputs->GetInput(ThicknessEnum); _assert_(thickness_input);
    … …  
    48364836
    48374837        /*Retrieve all inputs we will be needing: */
    4838         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     4838        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    48394839        GradientIndexing(&vertexpidlist[0],control_index);
    48404840        Input* rheologyb_input=material->inputs->GetInput(MaterialsRheologyBbarEnum); _assert_(rheologyb_input);
    … …  
    48784878
    48794879        /* Get node coordinates and dof list: */
    4880         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     4880        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    48814881        GradientIndexing(&doflist[0],control_index);
    48824882
    … …  
    49354935
    49364936        /* Get node coordinates and dof list: */
    4937         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     4937        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    49384938        GradientIndexing(&doflist[0],control_index);
    49394939
    … …  
    49974997        /*retrive parameters: */
    49984998        parameters->FindParam(&analysis_type,AnalysisTypeEnum);
    4999         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     4999        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    50005000        GradientIndexing(&vertexpidlist[0],control_index);
    50015001        this->GetConnectivityList(&connectivity[0]);
    … …  
    50775077        /*Retrieve all inputs we will be needing: */
    50785078        if(IsFloating())return;
    5079         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     5079        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    50805080        GradientIndexing(&vertexpidlist[0],control_index);
    50815081        Input* dragcoefficient_input=inputs->GetInput(FrictionCoefficientEnum); _assert_(dragcoefficient_input);
    … …  
    51375137
    51385138        /* Get node coordinates and dof list: */
    5139         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     5139        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    51405140        GradientIndexing(&vertexpidlist[0],control_index);
    51415141
    … …  
    51805180
    51815181        /* Get node coordinates and dof list: */
    5182         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     5182        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    51835183        GradientIndexing(&vertexpidlist[0],control_index);
    51845184
    … …  
    52295229
    52305230        /* Get node coordinates and dof list: */
    5231         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     5231        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    52325232        GradientIndexing(&vertexpidlist[0],control_index);
    52335233
    52345234        /*Retrieve all inputs and parameters*/
    5235         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     5235        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    52365236        this->parameters->FindParam(&num_responses,InversionNumCostFunctionsEnum);
    52375237        this->parameters->FindParam(&responses,NULL,InversionCostFunctionsEnum);
    … …  
    53485348
    53495349        /*Retrieve all inputs we will be needing: */
    5350         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     5350        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    53515351        Input* weights_input  =inputs->GetInput(InversionCostFunctionsCoefficientsEnum);              _assert_(weights_input);
    53525352        Input* rheologyb_input=material->inputs->GetInput(MaterialsRheologyBbarEnum); _assert_(rheologyb_input);
    … …  
    53875387
    53885388        /* Get node coordinates and dof list: */
    5389         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     5389        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    53905390
    53915391        /*Retrieve all inputs we will be needing: */
    … …  
    54465446
    54475447        /* Get node coordinates and dof list: */
    5448         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     5448        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    54495449
    54505450        /*Retrieve all inputs we will be needing: */
    … …  
    55055505
    55065506        /* Get node coordinates and dof list: */
    5507         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     5507        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    55085508
    55095509        /*Retrieve all inputs we will be needing: */
    … …  
    55635563
    55645564        /* Get node coordinates and dof list: */
    5565         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     5565        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    55665566
    55675567        /*Retrieve all inputs we will be needing: */
    … …  
    56225622
    56235623        /* Get node coordinates and dof list: */
    5624         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     5624        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    56255625
    56265626        /*Retrieve all inputs we will be needing: */
    … …  
    56845684
    56855685        /*Retrieve all inputs we will be needing: */
    5686         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     5686        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    56875687        Input* weights_input  =inputs->GetInput(InversionCostFunctionsCoefficientsEnum);   _assert_(weights_input);
    56885688        Input* thickness_input=inputs->GetInput(ThicknessEnum); _assert_(thickness_input);
    … …  
    57285728
    57295729        /*Retrieve all inputs we will be needing: */
    5730         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     5730        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    57315731        Input* weights_input  = inputs->GetInput(InversionCostFunctionsCoefficientsEnum); _assert_(weights_input);
    57325732        Input* thickness_input= inputs->GetInput(ThicknessEnum);                          _assert_(thickness_input);
    … …  
    57795779
    57805780        /*Retrieve all inputs we will be needing: */
    5781         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     5781        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    57825782        Input* weights_input  = inputs->GetInput(InversionCostFunctionsCoefficientsEnum); _assert_(weights_input);
    57835783        Input* thickness_input= inputs->GetInput(ThicknessEnum);                          _assert_(thickness_input);
    … …  
    58275827
    58285828        /*Retrieve all inputs we will be needing: */
    5829         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     5829        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    58305830        Input* thickness_input   =inputs->GetInput(ThicknessEnum);   _assert_(thickness_input);
    58315831        Input* thicknessobs_input=inputs->GetInput(InversionThicknessObsEnum);_assert_(thicknessobs_input);
    … …  
    58785878
    58795879        /*Retrieve all inputs and parameters*/
    5880         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     5880        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    58815881        this->parameters->FindParam(&num_responses,InversionNumCostFunctionsEnum);
    58825882        this->parameters->FindParam(&responses,NULL,InversionCostFunctionsEnum);
    … …  
    59685968
    59695969        /*Retrieve all inputs and parameters*/
    5970         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     5970        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    59715971        this->parameters->FindParam(&num_responses,InversionNumCostFunctionsEnum);
    59725972        this->parameters->FindParam(&responses,NULL,InversionCostFunctionsEnum);
    … …  
    61466146
    61476147        /*Retrieve all inputs we will be needing: */
    6148         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     6148        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    61496149        Input* weights_input=inputs->GetInput(InversionCostFunctionsCoefficientsEnum);         _assert_(weights_input);
    61506150        Input* drag_input   =inputs->GetInput(FrictionCoefficientEnum); _assert_(drag_input);
    … …  
    62186218
    62196219        /*Retrieve all inputs and parameters*/
    6220         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     6220        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    62216221        Input* vx_input=inputs->GetInput(VxEnum);       _assert_(vx_input);
    62226222        Input* vy_input=inputs->GetInput(VyEnum);       _assert_(vy_input);
    … …  
    63906390
    63916391        /*Retrieve all inputs and parameters*/
    6392         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     6392        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    63936393        latentheat   = matpar->GetLatentHeat();
    63946394        heatcapacity = matpar->GetHeatCapacity();
    … …  
    65266526
    65276527        /*Retrieve all inputs and parameters*/
    6528         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     6528        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    65296529        this->parameters->FindParam(&dt,TimesteppingTimeStepEnum);
    65306530        this->parameters->FindParam(&diffusivity,HydrologyshreveStabilizationEnum);
    … …  
    66166616
    66176617        /*Retrieve all inputs and parameters*/
    6618         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     6618        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    66196619        this->parameters->FindParam(&dt,TimesteppingTimeStepEnum);
    66206620        sediment_storing       = matpar->GetSedimentStoring();
    … …  
    66836683
    66846684        /*Retrieve all inputs and parameters*/
    6685         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     6685        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    66866686        Input* thickness_input=inputs->GetInput(HydrologydcEplThicknessEnum); _assert_(thickness_input);
    66876687        this->parameters->FindParam(&dt,TimesteppingTimeStepEnum);
    … …  
    67456745
    67466746        /*Retrieve all inputs and parameters*/
    6747         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     6747        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    67486748
    67496749        /* Start looping on the number of gaussian points: */
    … …  
    67906790
    67916791        /*Retrieve all inputs and parameters*/
    6792         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     6792        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    67936793        this->parameters->FindParam(&dt,TimesteppingTimeStepEnum);
    67946794        Input* basal_melting_input=inputs->GetInput(BasalforcingsMeltingRateEnum); _assert_(basal_melting_input);
    … …  
    68406840        /*Retrieve all inputs and parameters*/
    68416841        sediment_storing = matpar->GetSedimentStoring();
    6842         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     6842        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    68436843        this->parameters->FindParam(&dt,TimesteppingTimeStepEnum);
    68446844        Input* water_input=inputs->GetInput(BasalforcingsMeltingRateEnum);  _assert_(water_input);
    … …  
    69066906        /*Retrieve all inputs and parameters*/
    69076907        epl_specificstoring = matpar->GetEplSpecificStoring();
    6908         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     6908        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    69096909        this->parameters->FindParam(&dt,TimesteppingTimeStepEnum);
    69106910        Input* residual_input=inputs->GetInput(SedimentHeadResidualEnum);     _assert_(residual_input);
    … …  
    69766976
    69776977        /*Retrieve all inputs and parameters*/
    6978         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     6978        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    69796979        this->parameters->FindParam(&input_enum,InputToL2ProjectEnum);
    69806980        switch(input_enum){
    … …  
    74077407
    74087408        /*Retrieve all inputs and parameters*/
    7409         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     7409        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    74107410
    74117411        /* Start  looping on the number of gaussian points: */
    … …  
    74537453
    74547454        /*Retrieve all inputs and parameters*/
    7455         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     7455        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    74567456
    74577457        /*Get vertex indices that lie on bed*/
    … …  
    75037503
    75047504        /*Retrieve all inputs and parameters*/
    7505         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     7505        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    75067506
    75077507        /*Get vertex indices that lie on bed*/
    … …  
    75677567
    75687568        /*Retrieve all inputs and parameters*/
    7569         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     7569        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    75707570        this->parameters->FindParam(&dt,TimesteppingTimeStepEnum);
    75717571        this->parameters->FindParam(&meshtype,MeshTypeEnum);
    … …  
    76837683
    76847684        /*Retrieve all inputs and parameters*/
    7685         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     7685        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    76867686        this->parameters->FindParam(&dt,TimesteppingTimeStepEnum);
    76877687        this->parameters->FindParam(&meshtype,MeshTypeEnum);
    … …  
    77607760
    77617761        /*Retrieve all inputs and parameters*/
    7762         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     7762        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    77637763        this->parameters->FindParam(&dt,TimesteppingTimeStepEnum);
    77647764        this->parameters->FindParam(&stabilization,MasstransportStabilizationEnum);
    … …  
    78737873
    78747874        /*Retrieve all inputs and parameters*/
    7875         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     7875        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    78767876        this->parameters->FindParam(&dt,TimesteppingTimeStepEnum);
    78777877        this->parameters->FindParam(&stabilization,MasstransportStabilizationEnum);
    … …  
    79957995
    79967996        /*Retrieve all inputs and parameters*/
    7997         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     7997        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    79987998        this->parameters->FindParam(&dt,TimesteppingTimeStepEnum);
    79997999        Input* ms_input     = inputs->GetInput(SurfaceforcingsMassBalanceEnum);      _assert_(ms_input);
    … …  
    80468046        /*Retrieve all inputs and parameters*/
    80478047        this->parameters->FindParam(&dt,TimesteppingTimeStepEnum);
    8048         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     8048        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    80498049        Input* ms_input=inputs->GetInput(SurfaceforcingsMassBalanceEnum); _assert_(ms_input);
    80508050        Input* mb_input=inputs->GetInput(BasalforcingsMeltingRateEnum);   _assert_(mb_input);
    … …  
    80898089
    80908090        /*Retrieve all inputs and parameters*/
    8091         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     8091        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    80928092        this->parameters->FindParam(&dt,TimesteppingTimeStepEnum);
    80938093        Input* vz_input     = inputs->GetInput(VzEnum);                         _assert_(vz_input);
    … …  
    81508150
    81518151        /*Retrieve all inputs and parameters*/
    8152         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     8152        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    81538153        this->parameters->FindParam(&dt,TimesteppingTimeStepEnum);
    81548154        Input* vz_input  = inputs->GetInput(VzEnum);                         _assert_(vz_input);
    … …  
    82388238
    82398239        /*Retrieve all inputs and parameters*/
    8240         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     8240        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    82418241        this->parameters->FindParam(&dt,TimesteppingTimeStepEnum);
    82428242        this->parameters->FindParam(&meshtype,MeshTypeEnum);
    … …  
    83658365
    83668366        /*Retrieve all inputs and parameters*/
    8367         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     8367        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    83688368        this->parameters->FindParam(&dt,TimesteppingTimeStepEnum);
    83698369        damage_input  = this->material->inputs->GetInput(DamageDbarEnum);     _assert_(damage_input);
    … …  
    86628662
    86638663        /*Retrieve all Inputs and parameters: */
    8664         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     8664        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    86658665        this->parameters->FindParam(&stabilization,BalancethicknessStabilizationEnum);
    86668666        this->parameters->FindParam(&meshtype,MeshTypeEnum);
    … …  
    87678767
    87688768        /*Retrieve all inputs and parameters*/
    8769         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     8769        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    87708770        Input* vx_input=inputs->GetInput(VxEnum); _assert_(vx_input);
    87718771        Input* vy_input=inputs->GetInput(VyEnum); _assert_(vy_input);
    … …  
    88298829
    88308830        /*Retrieve all Inputs and parameters: */
    8831         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     8831        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    88328832        Input* H_input =inputs->GetInput(ThicknessEnum); _assert_(H_input);
    88338833        h=sqrt(2.*this->GetArea());
    … …  
    89028902
    89038903        /*Retrieve all inputs and parameters*/
    8904         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     8904        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    89058905        Input* thickness_input=inputs->GetInput(ThicknessEnum); _assert_(thickness_input);
    89068906
    … …  
    89648964
    89658965        /*Retrieve all inputs and parameters*/
    8966         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     8966        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    89678967        Input* ms_input=inputs->GetInput(SurfaceforcingsMassBalanceEnum); _assert_(ms_input);
    89688968        Input* mb_input=inputs->GetInput(BasalforcingsMeltingRateEnum);   _assert_(mb_input);
    … …  
    90069006
    90079007        /*Retrieve all inputs and parameters*/
    9008         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     9008        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    90099009        Input* ms_input=inputs->GetInput(SurfaceforcingsMassBalanceEnum); _assert_(ms_input);
    90109010        Input* mb_input=inputs->GetInput(BasalforcingsMeltingRateEnum);   _assert_(mb_input);
    … …  
    90569056
    90579057        /*Retrieve all inputs and parameters*/
    9058         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     9058        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    90599059        Input* ms_input=inputs->GetInput(SurfaceforcingsMassBalanceEnum); _assert_(ms_input);
    90609060        Input* mb_input=inputs->GetInput(BasalforcingsMeltingRateEnum);   _assert_(mb_input);
    … …  
    91329132
    91339133        /*Retrieve all inputs and parameters*/
    9134         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     9134        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    91359135        Input* H_input       = inputs->GetInput(ThicknessEnum); _assert_(H_input);
    91369136        Input* surface_input = inputs->GetInput(SurfaceEnum);   _assert_(surface_input);
    … …  
    91849184
    91859185        /*Retrieve all inputs and parameters*/
    9186         GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
     9186        ::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
    91879187        Input* H_input       = inputs->GetInput(ThicknessEnum); _assert_(H_input);
    91889188        Input* surface_input = inputs->GetInput(SurfaceEnum);   _assert_(surface_input);

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

    r16697 r16698  
    106106                void        GetSolutionFromInputsOneDof(Vector<IssmDouble>* solution,int enum_type);
    107107                void        GetVectorFromInputs(Vector<IssmDouble>* vector, int name_enum);
     108                void        GetVerticesCoordinates(IssmDouble** pxyz_list){_error_("not implemented yet");};
    108109                void        InputCreate(IssmDouble* vector,IoModel* iomodel,int M,int N,int vector_type,int vector_enum,int code);
    109110                void        InputDepthAverageAtBase(int enum_type,int average_enum_type,int object_enum=MeshElementsEnum);
Note: See TracChangeset for help on using the changeset viewer.