Ignore:
Timestamp:
03/14/12 14:15:28 (13 years ago)
Author:
cborstad
Message:

merged trunk-jpl through revision 11707 into branches/trunk-jpl-damage

Location:
issm/branches/trunk-jpl-damage
Files:
2 edited

Legend:

Unmodified
Added
Removed

  • issm/branches/trunk-jpl-damage

  • issm/branches/trunk-jpl-damage/src/c/objects/Elements/Tria.cpp

    r11684 r11708  
    279279/*Other*/
    280280/*FUNCTION Tria::AverageOntoPartition {{{1*/
    281 void  Tria::AverageOntoPartition(Vec partition_contributions,Vec partition_areas,double* vertex_response,double* qmu_part){
     281void  Tria::AverageOntoPartition(Vector* partition_contributions,Vector* partition_areas,double* vertex_response,double* qmu_part){
    282282
    283283        bool      already=false;
    … …  
    312312                }
    313313                if(!already){
    314                         VecSetValue(partition_contributions,partition[i],mean*area,ADD_VALUES);
    315                         VecSetValue(partition_areas,partition[i],area,ADD_VALUES);
     314                        partition_contributions->SetValue(partition[i],mean*area,ADD_VAL);
     315                        partition_areas->SetValue(partition[i],area,ADD_VAL);
    316316                };
    317317        }
    … …  
    925925/*}}}*/
    926926/*FUNCTION Tria::ComputeBasalStress {{{1*/
    927 void  Tria::ComputeBasalStress(Vec eps){
     927void  Tria::ComputeBasalStress(Vector* eps){
    928928        _error_("Not Implemented yet");
    929929}
    930930/*}}}*/
    931931/*FUNCTION Tria::ComputeStrainRate {{{1*/
    932 void  Tria::ComputeStrainRate(Vec eps){
     932void  Tria::ComputeStrainRate(Vector* eps){
    933933        _error_("Not Implemented yet");
    934934}
    … …  
    13241324/*}}}*/
    13251325/*FUNCTION Tria::GetVectorFromInputs{{{1*/
    1326 void  Tria::GetVectorFromInputs(Vec vector,int input_enum){
     1326void  Tria::GetVectorFromInputs(Vector* vector,int input_enum){
    13271327
    13281328        int doflist1[NUMVERTICES];
    … …  
    13431343/*}}}*/
    13441344/*FUNCTION Tria::GetVectorFromResults{{{1*/
    1345 void  Tria::GetVectorFromResults(Vec vector,int offset,int interp){
     1345void  Tria::GetVectorFromResults(Vector* vector,int offset,int interp){
    13461346
    13471347        /*Get result*/
    … …  
    22072207/*}}}*/
    22082208/*FUNCTION Tria::PotentialSheetUngrounding{{{1*/
    2209 void  Tria::PotentialSheetUngrounding(Vec potential_sheet_ungrounding){
     2209void  Tria::PotentialSheetUngrounding(Vector* potential_sheet_ungrounding){
    22102210
    22112211        int     i;
    … …  
    22292229                        if (bed_hydro>ba[i]){
    22302230                                /*Vertex that could potentially unground, flag it*/
    2231                                 VecSetValue(potential_sheet_ungrounding,nodes[i]->Sid(),1,INSERT_VALUES);
     2231                                potential_sheet_ungrounding->SetValue(nodes[i]->Sid(),1,INS_VAL);
    22322232                        }
    22332233                }
    … …  
    22872287/*}}}1*/
    22882288/*FUNCTION Tria::SmearFunction {{{1*/
    2289 void  Tria::SmearFunction(Vec smearedvector,double (*WeightFunction)(double distance,double radius),double radius){
     2289void  Tria::SmearFunction(Vector* smearedvector,double (*WeightFunction)(double distance,double radius),double radius){
    22902290        _error_("not implemented yet");
    22912291
    … …  
    25082508/*}}}*/
    25092509/*FUNCTION Tria::UpdatePotentialSheetUngrounding{{{1*/
    2510 int Tria::UpdatePotentialSheetUngrounding(double* vertices_potentially_ungrounding,Vec vec_nodes_on_iceshelf,double* nodes_on_iceshelf){
     2510int Tria::UpdatePotentialSheetUngrounding(double* vertices_potentially_ungrounding,Vector* vec_nodes_on_iceshelf,double* nodes_on_iceshelf){
    25112511
    25122512        int i;
    … …  
    25162516        for(i=0;i<3;i++){
    25172517                if (vertices_potentially_ungrounding[nodes[i]->Sid()]){
    2518                         VecSetValue(vec_nodes_on_iceshelf,nodes[i]->Sid(),1,INSERT_VALUES);
     2518                        vec_nodes_on_iceshelf->SetValue(nodes[i]->Sid(),1,INS_VAL);
    25192519               
    25202520                        /*If node was not on ice shelf, we flipped*/
    … …  
    31743174        }
    31753175
    3176         solution->SetValues(numdof,doflist,values,INSERT_VALUES);
     3176        solution->SetValues(numdof,doflist,values,INS_VAL);
    31773177
    31783178        /*Free ressources:*/
    … …  
    32133213        }
    32143214
    3215         solution->SetValues(numdof,doflist,values,INSERT_VALUES);
     3215        solution->SetValues(numdof,doflist,values,INS_VAL);
    32163216
    32173217        /*Free ressources:*/
    … …  
    33813381/*}}}*/
    33823382/*FUNCTION Tria::ControlInputGetGradient{{{1*/
    3383 void Tria::ControlInputGetGradient(Vec gradient,int enum_type,int control_index){
     3383void Tria::ControlInputGetGradient(Vector* gradient,int enum_type,int control_index){
    33843384
    33853385        int doflist1[NUMVERTICES];
    … …  
    34403440}/*}}}*/
    34413441/*FUNCTION Tria::Gradj {{{1*/
    3442 void  Tria::Gradj(Vec gradient,int control_type,int control_index){
     3442void  Tria::Gradj(Vector* gradient,int control_type,int control_index){
    34433443        /*dJ/dalpha = ∂L/∂alpha = ∂J/∂alpha + ∂/∂alpha(KU-F)*/
    34443444
    … …  
    35023502/*}}}*/
    35033503/*FUNCTION Tria::GradjBGradient{{{1*/
    3504 void  Tria::GradjBGradient(Vec gradient,int weight_index,int control_index){
     3504void  Tria::GradjBGradient(Vector* gradient,int weight_index,int control_index){
    35053505
    35063506        int        i,ig;
    … …  
    35353535                for (i=0;i<NUMVERTICES;i++) grade_g[i]+=-weight*Jdet*gauss->weight*(dbasis[0][i]*dk[0]+dbasis[1][i]*dk[1]);
    35363536        }
    3537         VecSetValues(gradient,NUMVERTICES,doflist1,(const double*)grade_g,ADD_VALUES);
     3537        gradient->SetValues(NUMVERTICES,doflist1,grade_g,ADD_VAL);
    35383538
    35393539        /*Clean up and return*/
    … …  
    35423542/*}}}*/
    35433543/*FUNCTION Tria::GradjZGradient{{{1*/
    3544 void  Tria::GradjZGradient(Vec gradient,int weight_index,int control_index){
     3544void  Tria::GradjZGradient(Vector* gradient,int weight_index,int control_index){
    35453545
    35463546        int        i,ig;
    … …  
    35753575                for (i=0;i<NUMVERTICES;i++) grade_g[i]+=-weight*Jdet*gauss->weight*(dbasis[0][i]*dk[0]+dbasis[1][i]*dk[1]);
    35763576        }
    3577         VecSetValues(gradient,NUMVERTICES,doflist1,(const double*)grade_g,ADD_VALUES);
     3577        gradient->SetValues(NUMVERTICES,doflist1,grade_g,ADD_VAL);
    35783578
    35793579        /*Clean up and return*/
    … …  
    35823582/*}}}*/
    35833583/*FUNCTION Tria::GradjBMacAyeal{{{1*/
    3584 void  Tria::GradjBMacAyeal(Vec gradient,int control_index){
     3584void  Tria::GradjBMacAyeal(Vector* gradient,int control_index){
    35853585
    35863586        /*Intermediaries*/
    … …  
    36323632        }
    36333633
    3634         VecSetValues(gradient,NUMVERTICES,doflist,(const double*)grad,ADD_VALUES);
     3634        gradient->SetValues(NUMVERTICES,doflist,grad,ADD_VAL);
    36353635
    36363636        /*clean-up*/
    … …  
    36393639/*}}}*/
    36403640/*FUNCTION Tria::GradjZMacAyeal{{{1*/
    3641 void  Tria::GradjZMacAyeal(Vec gradient,int control_index){
     3641void  Tria::GradjZMacAyeal(Vector* gradient,int control_index){
    36423642
    36433643        /*Intermediaries*/
    … …  
    36893689        }
    36903690
    3691         VecSetValues(gradient,NUMVERTICES,doflist,(const double*)grad,ADD_VALUES);
     3691        gradient->SetValues(NUMVERTICES,doflist,grad,ADD_VAL);
    36923692
    36933693        /*clean-up*/
    … …  
    36963696/*}}}*/
    36973697/*FUNCTION Tria::GradjDragMacAyeal {{{1*/
    3698 void  Tria::GradjDragMacAyeal(Vec gradient,int control_index){
     3698void  Tria::GradjDragMacAyeal(Vector* gradient,int control_index){
    36993699
    37003700        int        i,ig;
    … …  
    37763776        /*End Analytical gradient*/
    37773777
    3778         VecSetValues(gradient,NUMVERTICES,doflist1,(const double*)grade_g,ADD_VALUES);
     3778        gradient->SetValues(NUMVERTICES,doflist1,grade_g,ADD_VAL);
    37793779
    37803780        /*Clean up and return*/
    … …  
    37843784/*}}}*/
    37853785/*FUNCTION Tria::GradjDragGradient{{{1*/
    3786 void  Tria::GradjDragGradient(Vec gradient, int weight_index,int control_index){
     3786void  Tria::GradjDragGradient(Vector* gradient, int weight_index,int control_index){
    37873787
    37883788        int        i,ig;
    … …  
    38213821                }
    38223822        }
    3823         VecSetValues(gradient,NUMVERTICES,doflist1,(const double*)grade_g,ADD_VALUES);
     3823        gradient->SetValues(NUMVERTICES,doflist1,grade_g,ADD_VAL);
    38243824
    38253825        /*Clean up and return*/
    … …  
    38283828/*}}}*/
    38293829/*FUNCTION Tria::GradjDhDtBalancedthickness{{{1*/
    3830 void  Tria::GradjDhDtBalancedthickness(Vec gradient,int control_index){
     3830void  Tria::GradjDhDtBalancedthickness(Vector* gradient,int control_index){
    38313831
    38323832        /*Intermediaries*/
    … …  
    38403840        for(int i=0;i<NUMVERTICES;i++) gradient_g[i]=-lambda[i];
    38413841
    3842         VecSetValues(gradient,NUMVERTICES,doflist1,(const double*)gradient_g,INSERT_VALUES);
     3842        gradient->SetValues(NUMVERTICES,doflist1,gradient_g,INS_VAL);
    38433843}
    38443844/*}}}*/
    38453845/*FUNCTION Tria::GradjVxBalancedthickness{{{1*/
    3846 void  Tria::GradjVxBalancedthickness(Vec gradient,int control_index){
     3846void  Tria::GradjVxBalancedthickness(Vector* gradient,int control_index){
    38473847
    38483848        /*Intermediaries*/
    … …  
    38803880        }
    38813881
    3882         VecSetValues(gradient,NUMVERTICES,doflist1,(const double*)grade_g,ADD_VALUES);
     3882        gradient->SetValues(NUMVERTICES,doflist1,grade_g,ADD_VAL);
    38833883
    38843884        /*Clean up and return*/
    … …  
    38873887/*}}}*/
    38883888/*FUNCTION Tria::GradjVyBalancedthickness{{{1*/
    3889 void  Tria::GradjVyBalancedthickness(Vec gradient,int control_index){
     3889void  Tria::GradjVyBalancedthickness(Vector* gradient,int control_index){
    38903890
    38913891        /*Intermediaries*/
    … …  
    39223922                for(i=0;i<NUMVERTICES;i++) grade_g[i]+=thickness*Dlambda[1]*Jdet*gauss->weight*basis[i];
    39233923        }
    3924         VecSetValues(gradient,NUMVERTICES,doflist1,(const double*)grade_g,ADD_VALUES);
     3924        gradient->SetValues(NUMVERTICES,doflist1,grade_g,ADD_VAL);
    39253925
    39263926        /*Clean up and return*/
    … …  
    50165016/*}}}*/
    50175017/*FUNCTION Tria::GetVectorFromControlInputs{{{1*/
    5018 void  Tria::GetVectorFromControlInputs(Vec vector,int control_enum,int control_index,const char* data){
     5018void  Tria::GetVectorFromControlInputs(Vector* vector,int control_enum,int control_index,const char* data){
    50195019
    50205020        int doflist1[NUMVERTICES];
    … …  
    53185318        }
    53195319
    5320         solution->SetValues(numdof,doflist,values,INSERT_VALUES);
     5320        solution->SetValues(numdof,doflist,values,INS_VAL);
    53215321
    53225322        /*Free ressources:*/
Note: See TracChangeset for help on using the changeset viewer.