Changeset 25439


Ignore:
Timestamp:
08/21/20 10:49:02 (5 years ago)
Author:
Mathieu Morlighem
Message:

CHG: removing gauss loop in analyses

Location:
issm/trunk-jpl/src/c/analyses
Files:
26 edited

Legend:

Unmodified
Added
Removed

  • issm/trunk-jpl/src/c/analyses/AdjointBalancethickness2Analysis.cpp

    r25379 r25439  
    7474        /* Start  looping on the number of gaussian points: */
    7575        Gauss* gauss=element->NewGauss(2);
    76         for(int ig=gauss->begin();ig<gauss->end();ig++){
    77                 gauss->GaussPoint(ig);
     76        while(gauss->next()){
    7877
    7978                element->JacobianDeterminant(&Jdet,xyz_list,gauss);
    … …  
    175174
    176175        Gauss* gauss=element->NewGauss(2);
    177         for(int ig=gauss->begin();ig<gauss->end();ig++){
    178                 gauss->GaussPoint(ig);
     176        while(gauss->next()){
    179177
    180178                element->JacobianDeterminant(&Jdet,xyz_list,gauss);
    … …  
    224222
    225223        Gauss* gauss=element->NewGauss(2);
    226         for(int ig=gauss->begin();ig<gauss->end();ig++){
    227                 gauss->GaussPoint(ig);
     224        while(gauss->next()){
    228225
    229226                element->JacobianDeterminant(&Jdet,xyz_list,gauss);
    … …  
    278275        /* Start  looping on the number of gaussian points: */
    279276        Gauss* gauss=element->NewGauss(2);
    280         for(int ig=gauss->begin();ig<gauss->end();ig++){
    281                 gauss->GaussPoint(ig);
     277        while(gauss->next()){
    282278
    283279                element->JacobianDeterminant(&Jdet,xyz_list,gauss);
    … …  
    328324        /* Start  looping on the number of gaussian points: */
    329325        Gauss* gauss=element->NewGauss(2);
    330         for(int ig=gauss->begin();ig<gauss->end();ig++){
    331                 gauss->GaussPoint(ig);
     326        while(gauss->next()){
    332327
    333328                element->JacobianDeterminant(&Jdet,xyz_list,gauss);

  • issm/trunk-jpl/src/c/analyses/AdjointBalancethicknessAnalysis.cpp

    r25379 r25439  
    9494        /* Start  looping on the number of gaussian points: */
    9595        Gauss* gauss=basalelement->NewGauss(2);
    96         for(int ig=gauss->begin();ig<gauss->end();ig++){
    97                 gauss->GaussPoint(ig);
     96        while(gauss->next()){
    9897
    9998                basalelement->JacobianDeterminant(&Jdet,xyz_list,gauss);
    … …  
    288287        /* Start  looping on the number of gaussian points: */
    289288        Gauss* gauss=element->NewGauss(4);
    290         for(int ig=gauss->begin();ig<gauss->end();ig++){
    291                 gauss->GaussPoint(ig);
     289        while(gauss->next()){
    292290
    293291                adjoint_input->GetInputDerivativeValue(&Dlambda[0],xyz_list,gauss);
    … …  
    335333        /* Start  looping on the number of gaussian points: */
    336334        Gauss* gauss=element->NewGauss(4);
    337         for(int ig=gauss->begin();ig<gauss->end();ig++){
    338                 gauss->GaussPoint(ig);
     335        while(gauss->next()){
    339336
    340337                adjoint_input->GetInputDerivativeValue(&Dlambda[0],xyz_list,gauss);

  • issm/trunk-jpl/src/c/analyses/AdjointHorizAnalysis.cpp

    r25379 r25439  
    102102        /* Start  looping on the number of gaussian points: */
    103103        Gauss* gauss=element->NewGauss(5);
    104         for(int ig=gauss->begin();ig<gauss->end();ig++){
    105                 gauss->GaussPoint(ig);
     104        while(gauss->next()){
    106105
    107106                element->JacobianDeterminant(&Jdet,xyz_list,gauss);
    … …  
    179178        /* Start  looping on the number of gaussian points: */
    180179        Gauss* gauss=element->NewGauss(5);
    181         for(int ig=gauss->begin();ig<gauss->end();ig++){
    182                 gauss->GaussPoint(ig);
     180        while(gauss->next()){
    183181
    184182                element->JacobianDeterminant(&Jdet,xyz_list,gauss);
    … …  
    285283        /* Start  looping on the number of gaussian points: */
    286284        Gauss* gauss=basalelement->NewGauss(2);
    287         for(int ig=gauss->begin();ig<gauss->end();ig++){
    288                 gauss->GaussPoint(ig);
     285        while(gauss->next()){
    289286
    290287                basalelement->JacobianDeterminant(&Jdet,xyz_list,gauss);
    … …  
    397394        /* Start  looping on the number of gaussian points: */
    398395        Gauss* gauss=element->NewGaussTop(4);
    399         for(int ig=gauss->begin();ig<gauss->end();ig++){
    400                 gauss->GaussPoint(ig);
     396        while(gauss->next()){
    401397
    402398                element->JacobianDeterminantTop(&Jdet,xyz_list_top,gauss);
    … …  
    636632        /* Start  looping on the number of gaussian points: */
    637633        Gauss* gauss=element->NewGaussTop(4);
    638         for(int ig=gauss->begin();ig<gauss->end();ig++){
    639                 gauss->GaussPoint(ig);
     634        while(gauss->next()){
    640635
    641636                element->JacobianDeterminantTop(&Jdet,xyz_list_top,gauss);
    … …  
    889884        /* Start  looping on the number of gaussian points: */
    890885        Gauss* gauss=basalelement->NewGauss(4);
    891         for(int ig=gauss->begin();ig<gauss->end();ig++){
    892                 gauss->GaussPoint(ig);
     886        while(gauss->next()){
    893887
    894888                basalelement->JacobianDeterminant(&Jdet,xyz_list,gauss);
    … …  
    12281222        /* Start  looping on the number of gaussian points: */
    12291223        Gauss* gauss=basalelement->NewGauss(2);
    1230         for(int ig=gauss->begin();ig<gauss->end();ig++){
    1231                 gauss->GaussPoint(ig);
     1224        while(gauss->next()){
    12321225
    12331226                basalelement->JacobianDeterminant(&Jdet,xyz_list,gauss);
    … …  
    13211314        /* Start  looping on the number of gaussian points: */
    13221315        Gauss* gauss=basalelement->NewGauss(4);
    1323         for(int ig=gauss->begin();ig<gauss->end();ig++){
    1324                 gauss->GaussPoint(ig);
     1316        while(gauss->next()){
    13251317
    13261318                thickness_input->GetInputValue(&thickness,gauss);
    … …  
    14131405        /* Start  looping on the number of gaussian points: */
    14141406        Gauss* gauss=element->NewGauss(2);
    1415         for(int ig=gauss->begin();ig<gauss->end();ig++){
    1416                 gauss->GaussPoint(ig);
     1407        while(gauss->next()){
    14171408                element->JacobianDeterminant(&Jdet,xyz_list,gauss);
    14181409                element->NodalFunctionsP1Derivatives(dbasis,xyz_list,gauss);
    … …  
    14811472        /* Start  looping on the number of gaussian points: */
    14821473        Gauss* gauss=element->NewGauss(2);
    1483         for(int ig=gauss->begin();ig<gauss->end();ig++){
    1484                 gauss->GaussPoint(ig);
     1474        while(gauss->next()){
    14851475                element->JacobianDeterminant(&Jdet,xyz_list,gauss);
    14861476                element->NodalFunctionsP1(basis,gauss);
    … …  
    15421532        /* Start  looping on the number of gaussian points: */
    15431533        Gauss* gauss=element->NewGauss(4);
    1544         for(int ig=gauss->begin();ig<gauss->end();ig++){
    1545                 gauss->GaussPoint(ig);
     1534        while(gauss->next()){
    15461535
    15471536                thickness_input->GetInputValue(&thickness,gauss);
    … …  
    16341623        /* Start  looping on the number of gaussian points: */
    16351624        Gauss* gauss=basalelement->NewGauss(4);
    1636         for(int ig=gauss->begin();ig<gauss->end();ig++){
    1637                 gauss->GaussPoint(ig);
     1625        while(gauss->next()){
    16381626
    16391627                thickness_input->GetInputValue(&thickness,gauss);
    … …  
    17291717        /* Start  looping on the number of gaussian points: */
    17301718        Gauss* gauss=basalelement->NewGauss(2);
    1731         for(int ig=gauss->begin();ig<gauss->end();ig++){
    1732                 gauss->GaussPoint(ig);
     1719        while(gauss->next()){
    17331720
    17341721                basalelement->JacobianDeterminant(&Jdet,xyz_list,gauss);
    … …  
    18071794        /* Start  looping on the number of gaussian points: */
    18081795        Gauss* gauss=element->NewGaussBase(4);
    1809         for(int ig=gauss->begin();ig<gauss->end();ig++){
    1810                 gauss->GaussPoint(ig);
     1796        while(gauss->next()){
    18111797
    18121798                adjointx_input->GetInputValue(&lambda, gauss);
    … …  
    19031889        /* Start  looping on the number of gaussian points: */
    19041890        Gauss* gauss=element->NewGaussBase(4);
    1905         for(int ig=gauss->begin();ig<gauss->end();ig++){
    1906                 gauss->GaussPoint(ig);
     1891        while(gauss->next()){
    19071892
    19081893                adjointx_input->GetInputValue(&lambda, gauss);
    … …  
    20081993        /* Start  looping on the number of gaussian points: */
    20091994        Gauss* gauss=basalelement->NewGauss(4);
    2010         for(int ig=gauss->begin();ig<gauss->end();ig++){
    2011                 gauss->GaussPoint(ig);
     1995        while(gauss->next()){
    20121996
    20131997                adjointx_input->GetInputValue(&lambda, gauss);
    … …  
    20842068        /* Start  looping on the number of gaussian points: */
    20852069        Gauss* gauss=element->NewGaussBase(4);
    2086         for(int ig=gauss->begin();ig<gauss->end();ig++){
    2087                 gauss->GaussPoint(ig);
     2070        while(gauss->next()){
    20882071
    20892072                adjointx_input->GetInputValue(&lambda, gauss);
    … …  
    21782161        /* Start  looping on the number of gaussian points: */
    21792162        Gauss* gauss=element->NewGaussBase(4);
    2180         for(int ig=gauss->begin();ig<gauss->end();ig++){
    2181                 gauss->GaussPoint(ig);
     2163        while(gauss->next()){
    21822164
    21832165                adjointx_input->GetInputValue(&lambda, gauss);
    … …  
    22832265        /* Start  looping on the number of gaussian points: */
    22842266        Gauss* gauss=basalelement->NewGauss(4);
    2285         for(int ig=gauss->begin();ig<gauss->end();ig++){
    2286                 gauss->GaussPoint(ig);
     2267        while(gauss->next()){
    22872268
    22882269                adjointx_input->GetInputValue(&lambda, gauss);
    … …  
    23872368        /* Start  looping on the number of gaussian points: */
    23882369        Gauss* gauss=basalelement->NewGauss(4);
    2389         for(int ig=gauss->begin();ig<gauss->end();ig++){
    2390                 gauss->GaussPoint(ig);
     2370        while(gauss->next()){
    23912371
    23922372                thickness_input->GetInputValue(&thickness,gauss);

  • issm/trunk-jpl/src/c/analyses/Balancethickness2Analysis.cpp

    r25379 r25439  
    9090        /* Start  looping on the number of gaussian points: */
    9191        Gauss* gauss=element->NewGauss(2);
    92         for(int ig=gauss->begin();ig<gauss->end();ig++){
    93                 gauss->GaussPoint(ig);
     92        while(gauss->next()){
    9493                element->NodalFunctionsDerivatives(dbasis,xyz_list,gauss);
    9594                element->JacobianDeterminant(&Jdet,xyz_list,gauss);
    … …  
    148147        /* Start  looping on the number of gaussian points: */
    149148        Gauss* gauss=element->NewGauss(2);
    150         for(int ig=gauss->begin();ig<gauss->end();ig++){
    151                 gauss->GaussPoint(ig);
     149        while(gauss->next()){
    152150
    153151                element->JacobianDeterminant(&Jdet,xyz_list,gauss);

  • issm/trunk-jpl/src/c/analyses/BalancethicknessAnalysis.cpp

    r25379 r25439  
    186186        /* Start  looping on the number of gaussian points: */
    187187        Gauss* gauss=element->NewGauss(2);
    188         for(int ig=gauss->begin();ig<gauss->end();ig++){
    189                 gauss->GaussPoint(ig);
     188        while(gauss->next()){
    190189
    191190                element->JacobianDeterminant(&Jdet,xyz_list,gauss);
    … …  
    284283        /* Start  looping on the number of gaussian points: */
    285284        Gauss* gauss=element->NewGauss(2);
    286         for(int ig=gauss->begin();ig<gauss->end();ig++){
    287                 gauss->GaussPoint(ig);
     285        while(gauss->next()){
    288286
    289287                element->JacobianDeterminant(&Jdet,xyz_list,gauss);
    … …  
    361359        /* Start  looping on the number of gaussian points: */
    362360        Gauss* gauss=element->NewGauss(2);
    363         for(int ig=gauss->begin();ig<gauss->end();ig++){
    364                 gauss->GaussPoint(ig);
     361        while(gauss->next()){
    365362
    366363                element->JacobianDeterminant(&Jdet,xyz_list,gauss);
    … …  
    402399        /* Start  looping on the number of gaussian points: */
    403400        Gauss* gauss=element->NewGauss(2);
    404         for(int ig=gauss->begin();ig<gauss->end();ig++){
    405                 gauss->GaussPoint(ig);
     401        while(gauss->next()){
    406402
    407403                element->JacobianDeterminant(&Jdet,xyz_list,gauss);
    … …  
    468464        /* Start  looping on the number of gaussian points: */
    469465        Gauss* gauss=element->NewGauss(4);
    470         for(int ig=gauss->begin();ig<gauss->end();ig++){
    471                 gauss->GaussPoint(ig);
     466        while(gauss->next()){
    472467
    473468                thickness_input->GetInputValue(&thickness, gauss);

  • issm/trunk-jpl/src/c/analyses/BalancevelocityAnalysis.cpp

    r25379 r25439  
    106106        /*Start looping on the number of gaussian points:*/
    107107        Gauss* gauss=element->NewGauss(2);
    108         for(int ig=gauss->begin();ig<gauss->end();ig++){
    109                 gauss->GaussPoint(ig);
     108        while(gauss->next()){
    110109
    111110                H_input->GetInputValue(&thickness,gauss);
    … …  
    200199        /* Start  looping on the number of gaussian points: */
    201200        Gauss* gauss=basalelement->NewGauss(2);
    202         for(int ig=gauss->begin();ig<gauss->end();ig++){
    203                 gauss->GaussPoint(ig);
     201        while(gauss->next()){
    204202
    205203                basalelement->JacobianDeterminant(&Jdet,xyz_list,gauss);

  • issm/trunk-jpl/src/c/analyses/DamageEvolutionAnalysis.cpp

    r25379 r25439  
    472472        /* Start  looping on the number of gaussian points: */
    473473        Gauss* gauss=element->NewGauss(2);
    474         for(int ig=gauss->begin();ig<gauss->end();ig++){
    475                 gauss->GaussPoint(ig);
     474        while(gauss->next()){
    476475
    477476                element->JacobianDeterminant(&Jdet,xyz_list,gauss);
    … …  
    655654        /* Start  looping on the number of gaussian points: */
    656655        Gauss* gauss=element->NewGauss(2);
    657         for(int ig=gauss->begin();ig<gauss->end();ig++){
    658                 gauss->GaussPoint(ig);
     656        while(gauss->next()){
    659657
    660658                element->JacobianDeterminant(&Jdet,xyz_list,gauss);
    … …  
    749747        /* Start  looping on the number of gaussian points: */
    750748        Gauss* gauss=element->NewGauss(2);
    751         for(int ig=gauss->begin();ig<gauss->end();ig++){
    752                 gauss->GaussPoint(ig);
     749        while(gauss->next()){
    753750
    754751                element->JacobianDeterminant(&Jdet,xyz_list,gauss);
    … …  
    795792        /* Start  looping on the number of gaussian points: */
    796793        Gauss* gauss=element->NewGauss(2);
    797         for(int ig=gauss->begin();ig<gauss->end();ig++){
    798                 gauss->GaussPoint(ig);
    799 
     794        while(gauss->next()){
    800795                element->JacobianDeterminant(&Jdet,xyz_list,gauss);
    801796                element->NodalFunctions(basis,gauss);

  • issm/trunk-jpl/src/c/analyses/DepthAverageAnalysis.cpp

    r25379 r25439  
    7070        /* Start  looping on the number of gaussian points: */
    7171        Gauss* gauss=element->NewGauss(2);
    72         for(int ig=gauss->begin();ig<gauss->end();ig++){
    73                 gauss->GaussPoint(ig);
     72        while(gauss->next()){
    7473
    7574                element->JacobianDeterminant(&Jdet,xyz_list,gauss);
    … …  
    118117        /* Start  looping on the number of gaussian points: */
    119118        Gauss* gauss=element->NewGauss(3);
    120         for(int ig=gauss->begin();ig<gauss->end();ig++){
    121                 gauss->GaussPoint(ig);
     119        while(gauss->next()){
    122120
    123121                element->JacobianDeterminant(&Jdet,xyz_list,gauss);

  • issm/trunk-jpl/src/c/analyses/ExtrapolationAnalysis.cpp

    r25379 r25439  
    138138        /* Start  looping on the number of gaussian points: */
    139139        Gauss* gauss=workelement->NewGauss(2);
    140         for(int ig=gauss->begin();ig<gauss->end();ig++){
    141                 gauss->GaussPoint(ig);
     140        while(gauss->next()){
    142141
    143142                workelement->JacobianDeterminant(&Jdet,xyz_list,gauss);

  • issm/trunk-jpl/src/c/analyses/ExtrudeFromBaseAnalysis.cpp

    r25379 r25439  
    6969        /* Start  looping on the number of gaussian points: */
    7070        Gauss* gauss=element->NewGauss(2);
    71         for(int ig=gauss->begin();ig<gauss->end();ig++){
    72                 gauss->GaussPoint(ig);
     71        while(gauss->next()){
    7372
    7473                element->JacobianDeterminant(&Jdet,xyz_list,gauss);

  • issm/trunk-jpl/src/c/analyses/ExtrudeFromTopAnalysis.cpp

    r25379 r25439  
    6969        /* Start  looping on the number of gaussian points: */
    7070        Gauss* gauss=element->NewGauss(2);
    71         for(int ig=gauss->begin();ig<gauss->end();ig++){
    72                 gauss->GaussPoint(ig);
     71        while(gauss->next()){
    7372
    7473                element->JacobianDeterminant(&Jdet,xyz_list,gauss);

  • issm/trunk-jpl/src/c/analyses/FreeSurfaceBaseAnalysis.cpp

    r25379 r25439  
    183183        /* Start  looping on the number of gaussian points: */
    184184        Gauss* gauss=basalelement->NewGauss(2);
    185         for(int ig=gauss->begin();ig<gauss->end();ig++){
    186                 gauss->GaussPoint(ig);
     185        while(gauss->next()){
    187186
    188187                basalelement->JacobianDeterminant(&Jdet,xyz_list,gauss);
    … …  
    320319        /* Start  looping on the number of gaussian points: */
    321320        Gauss* gauss=basalelement->NewGauss(2);
    322         for(int ig=gauss->begin();ig<gauss->end();ig++){
    323                 gauss->GaussPoint(ig);
     321        while(gauss->next()){
    324322
    325323                basalelement->JacobianDeterminant(&Jdet,xyz_list,gauss);

  • issm/trunk-jpl/src/c/analyses/FreeSurfaceTopAnalysis.cpp

    r25379 r25439  
    159159        /* Start  looping on the number of gaussian points: */
    160160        Gauss* gauss=topelement->NewGauss(2);
    161         for(int ig=gauss->begin();ig<gauss->end();ig++){
    162                 gauss->GaussPoint(ig);
     161        while(gauss->next()){
    163162
    164163                topelement->JacobianDeterminant(&Jdet,xyz_list,gauss);
    … …  
    298297        /* Start  looping on the number of gaussian points: */
    299298        Gauss* gauss=topelement->NewGauss(2);
    300         for(int ig=gauss->begin();ig<gauss->end();ig++){
    301                 gauss->GaussPoint(ig);
     299        while(gauss->next()){
    302300
    303301                topelement->JacobianDeterminant(&Jdet,xyz_list,gauss);

  • issm/trunk-jpl/src/c/analyses/GLheightadvectionAnalysis.cpp

    r25379 r25439  
    103103        /* Start  looping on the number of gaussian points: */
    104104        Gauss* gauss=basalelement->NewGauss(4);
    105         for(int ig=gauss->begin();ig<gauss->end();ig++){
    106                 gauss->GaussPoint(ig);
     105        while(gauss->next()){
    107106
    108107                basalelement->JacobianDeterminant(&Jdet,xyz_list,gauss);

  • issm/trunk-jpl/src/c/analyses/HydrologyDCEfficientAnalysis.cpp

    r25379 r25439  
    215215        /* Start  looping on the number of gaussian points: */
    216216        Gauss* gauss = basalelement->NewGauss(2);
    217         for(int ig=gauss->begin();ig<gauss->end();ig++){
    218                 gauss       ->GaussPoint(ig);
     217        while(gauss->next()){
    219218                basalelement->JacobianDeterminant(&Jdet,xyz_list,gauss);
    220219                basalelement->NodalFunctionsDerivatives(dbasis,xyz_list,gauss);
    … …  
    339338        /* Start  looping on the number of gaussian points: */
    340339        Gauss* gauss = basalelement->NewGauss(2);
    341         for(int ig=gauss->begin();ig<gauss->end();ig++){
    342                 gauss->GaussPoint(ig);
     340        while(gauss->next()){
    343341                basalelement ->JacobianDeterminant(&Jdet,xyz_list,gauss);
    344342                basalelement ->NodalFunctions(basis,gauss);

  • issm/trunk-jpl/src/c/analyses/HydrologyDCInefficientAnalysis.cpp

    r25379 r25439  
    245245        /* Start  looping on the number of gaussian points: */
    246246        Gauss* gauss=basalelement->NewGauss(2);
    247 
    248         for(int ig=gauss -> begin();ig<gauss->end();ig++){
    249                 gauss       ->GaussPoint(ig);
     247        while(gauss->next()){
    250248                basalelement->JacobianDeterminant(&Jdet,xyz_list,gauss);
    251249                basalelement->NodalFunctionsDerivatives(dbasis,xyz_list,gauss);
    … …  
    384382        /* Start  looping on the number of gaussian points: */
    385383        Gauss* gauss=basalelement->NewGauss(2);
    386 
    387         for(int ig=gauss->begin();ig<gauss->end();ig++){
    388                 gauss->GaussPoint(ig);
     384        while(gauss->next()){
    389385                basalelement->JacobianDeterminant(&Jdet,xyz_list,gauss);
    390386                basalelement->NodalFunctions(basis,gauss);

  • issm/trunk-jpl/src/c/analyses/HydrologyShreveAnalysis.cpp

    r25379 r25439  
    177177        /* Start  looping on the number of gaussian points: */
    178178        Gauss* gauss=element->NewGauss(2);
    179         for(int ig=gauss->begin();ig<gauss->end();ig++){
    180                 gauss->GaussPoint(ig);
     179        while(gauss->next()){
    181180
    182181                element->JacobianDeterminant(&Jdet,xyz_list,gauss);
    … …  
    255254        /* Start  looping on the number of gaussian points: */
    256255        Gauss* gauss=element->NewGauss(2);
    257         for(int ig=gauss->begin();ig<gauss->end();ig++){
    258                 gauss->GaussPoint(ig);
     256        while(gauss->next()){
    259257
    260258                element->JacobianDeterminant(&Jdet,xyz_list,gauss);

  • issm/trunk-jpl/src/c/analyses/L2ProjectionBaseAnalysis.cpp

    r25379 r25439  
    102102        /* Start  looping on the number of gaussian points: */
    103103        Gauss* gauss=basalelement->NewGauss(2);
    104         for(int ig=gauss->begin();ig<gauss->end();ig++){
    105                 gauss->GaussPoint(ig);
     104        while(gauss->next()){
    106105
    107106                basalelement->JacobianDeterminant(&Jdet,xyz_list,gauss);
    … …  
    172171        /* Start  looping on the number of gaussian points: */
    173172        Gauss* gauss=basalelement->NewGauss(2);
    174         for(int ig=gauss->begin();ig<gauss->end();ig++){
    175                 gauss->GaussPoint(ig);
     173        while(gauss->next()){
    176174
    177175                basalelement->JacobianDeterminant(&Jdet,xyz_list,gauss);

  • issm/trunk-jpl/src/c/analyses/L2ProjectionEPLAnalysis.cpp

    r25379 r25439  
    132132        /* Start  looping on the number of gaussian points: */
    133133        Gauss* gauss=basalelement->NewGauss(2);
    134         for(int ig=gauss->begin();ig<gauss->end();ig++){
    135                 gauss->GaussPoint(ig);
     134        while(gauss->next()){
    136135
    137136                basalelement->JacobianDeterminant(&Jdet,xyz_list,gauss);
    … …  
    203202        /* Start  looping on the number of gaussian points: */
    204203        Gauss* gauss=basalelement->NewGauss(2);
    205         for(int ig=gauss->begin();ig<gauss->end();ig++){
    206                 gauss->GaussPoint(ig);
     204        while(gauss->next()){
    207205
    208206                basalelement->JacobianDeterminant(&Jdet,xyz_list,gauss);

  • issm/trunk-jpl/src/c/analyses/LevelsetAnalysis.cpp

    r25379 r25439  
    308308        /* Start  looping on the number of gaussian points: */
    309309        Gauss* gauss=basalelement->NewGauss(2);
    310         for(int ig=gauss->begin();ig<gauss->end();ig++){
    311                 gauss->GaussPoint(ig);
     310        while(gauss->next()){
    312311
    313312                basalelement->JacobianDeterminant(&Jdet,xyz_list,gauss);
    … …  
    551550
    552551        /*Intermediaries */
    553         int i, ig, domaintype;
     552        int domaintype;
    554553        IssmDouble  Jdet,dt;
    555554        IssmDouble  lsf;
    … …  
    573572                /* Start  looping on the number of gaussian points: */
    574573                Gauss* gauss=basalelement->NewGauss(2);
    575                 for(ig=gauss->begin();ig<gauss->end();ig++){
    576                         gauss->GaussPoint(ig);
    577 
     574                while(gauss->next()){
    578575                        basalelement->JacobianDeterminant(&Jdet,xyz_list,gauss);
    579576                        basalelement->NodalFunctions(basis,gauss);
    … …  
    581578                        /* old function value */
    582579                        levelset_input->GetInputValue(&lsf,gauss);
    583                         for(i=0;i<numnodes;i++) pe->values[i]+=Jdet*gauss->weight*lsf*basis[i];
     580                        for(int i=0;i<numnodes;i++) pe->values[i]+=Jdet*gauss->weight*lsf*basis[i];
    584581                }
    585582

  • issm/trunk-jpl/src/c/analyses/MasstransportAnalysis.cpp

    r25379 r25439  
    328328        /* Start  looping on the number of gaussian points: */
    329329        Gauss* gauss=element->NewGauss(2);
    330         for(int ig=gauss->begin();ig<gauss->end();ig++){
    331                 gauss->GaussPoint(ig);
     330        while(gauss->next()){
    332331
    333332                element->JacobianDeterminant(&Jdet,xyz_list,gauss);
    … …  
    517516        /* Start  looping on the number of gaussian points: */
    518517        Gauss* gauss=element->NewGauss(2);
    519         for(int ig=gauss->begin();ig<gauss->end();ig++){
    520                 gauss->GaussPoint(ig);
     518        while(gauss->next()){
    521519
    522520                element->JacobianDeterminant(&Jdet,xyz_list,gauss);
    … …  
    642640
    643641        /* Start  looping on the number of gaussian points: */
    644         for(int ig=gauss->begin();ig<gauss->end();ig++){
    645                 gauss->GaussPoint(ig);
     642        while(gauss->next()){
    646643
    647644                element->JacobianDeterminant(&Jdet,xyz_list,gauss);
    … …  
    748745
    749746        /* Start  looping on the number of gaussian points: */
    750         for(int ig=gauss->begin();ig<gauss->end();ig++){
    751                 gauss->GaussPoint(ig);
     747        while(gauss->next()){
    752748
    753749                element->JacobianDeterminant(&Jdet,xyz_list,gauss);
    … …  
    953949        /* Start  looping on the number of gaussian points: */
    954950        Gauss* gauss=element->NewGauss(2);
    955         for(int ig=gauss->begin();ig<gauss->end();ig++){
    956                 gauss->GaussPoint(ig);
     951        while(gauss->next()){
    957952
    958953                element->JacobianDeterminant(&Jdet,xyz_list,gauss);
    … …  
    10061001        /* Start  looping on the number of gaussian points: */
    10071002        Gauss* gauss=element->NewGauss(2);
    1008         for(int ig=gauss->begin();ig<gauss->end();ig++){
    1009                 gauss->GaussPoint(ig);
     1003        while(gauss->next()){
    10101004
    10111005                element->JacobianDeterminant(&Jdet,xyz_list,gauss);
    … …  
    10791073
    10801074        /* Start  looping on the number of gaussian points: */
    1081         for(int ig=gauss->begin();ig<gauss->end();ig++){
    1082                 gauss->GaussPoint(ig);
     1075        while(gauss->next()){
    10831076
    10841077                element->JacobianDeterminant(&Jdet,xyz_list,gauss);

  • issm/trunk-jpl/src/c/analyses/MeltingAnalysis.cpp

    r25379 r25439  
    102102        /* Start  looping on the number of gaussian points: */
    103103        Gauss* gauss=basalelement->NewGauss(2);
    104         for(int ig=gauss->begin();ig<gauss->end();ig++){
    105                 gauss->GaussPoint(ig);
     104        while(gauss->next()){
    106105
    107106                basalelement->JacobianDeterminant(&Jdet,xyz_list,gauss);

  • issm/trunk-jpl/src/c/analyses/SmoothAnalysis.cpp

    r25379 r25439  
    7474        /* Start looping on the number of gaussian points: */
    7575        Gauss* gauss=element->NewGauss(2);
    76         for(int ig=gauss->begin();ig<gauss->end();ig++){
    77                 gauss->GaussPoint(ig);
     76        while(gauss->next()){
    7877
    7978                element->JacobianDeterminant(&Jdet,xyz_list,gauss);
    … …  
    155154        /* Start  looping on the number of gaussian points: */
    156155        Gauss* gauss=element->NewGauss(2);
    157         for(int ig=gauss->begin();ig<gauss->end();ig++){
    158                 gauss->GaussPoint(ig);
     156        while(gauss->next()){
    159157
    160158                element->JacobianDeterminant(&Jdet,xyz_list,gauss);

  • issm/trunk-jpl/src/c/analyses/StressbalanceSIAAnalysis.cpp

    r25379 r25439  
    444444
    445445                Gauss* gauss=element->NewGaussLine(nodedown,nodeup,3);
    446                 for(int ig=gauss->begin();ig<gauss->end();ig++){
    447                         gauss->GaussPoint(ig);
     446                while(gauss->next()){
    448447
    449448                        B_input->GetInputValue(&B,gauss);

  • issm/trunk-jpl/src/c/analyses/ThermalAnalysis.cpp

    r25379 r25439  
    354354        /* Start  looping on the number of gaussian points: */
    355355        Gauss* gauss=element->NewGaussBase(4);
    356         for(int ig=gauss->begin();ig<gauss->end();ig++){
    357                 gauss->GaussPoint(ig);
     356        while(gauss->next()){
    358357
    359358                element->JacobianDeterminantBase(&Jdet,xyz_list_base,gauss);
    … …  
    412411        /* Start  looping on the number of gaussian points: */
    413412        Gauss* gauss=element->NewGauss(4);
    414         for(int ig=gauss->begin();ig<gauss->end();ig++){
    415                 gauss->GaussPoint(ig);
     413        while(gauss->next()){
    416414
    417415                element->JacobianDeterminant(&Jdet,xyz_list,gauss);
    … …  
    564562        /* Start  looping on the number of gaussian points: */
    565563        Gauss* gauss   = element->NewGaussBase(4);
    566         for(int ig=gauss->begin();ig<gauss->end();ig++){
    567                 gauss->GaussPoint(ig);
     564        while(gauss->next()){
    568565
    569566                element->JacobianDeterminantBase(&Jdet,xyz_list_base,gauss);
    … …  
    623620        /* Start  looping on the number of gaussian points: */
    624621        Gauss* gauss=element->NewGaussBase(4);
    625         for(int ig=gauss->begin();ig<gauss->end();ig++){
    626                 gauss->GaussPoint(ig);
     622        while(gauss->next()){
    627623
    628624                element->JacobianDeterminantBase(&Jdet,xyz_list_base,gauss);
    … …  
    683679        /* Start  looping on the number of gaussian points: */
    684680        Gauss* gauss=element->NewGauss(4);
    685         for(int ig=gauss->begin();ig<gauss->end();ig++){
    686                 gauss->GaussPoint(ig);
     681        while(gauss->next()){
    687682
    688683                element->JacobianDeterminant(&Jdet,xyz_list,gauss);

  • issm/trunk-jpl/src/c/analyses/UzawaPressureAnalysis.cpp

    r25379 r25439  
    8787
    8888        //Gauss* gauss = element->NewGauss(5);
    89         //for(int ig=gauss->begin();ig<gauss->end();ig++){
    90         //      gauss->GaussPoint(ig);
    91 
     89        //while(gauss->next()){
    9290        //      element->JacobianDeterminant(&Jdet,xyz_list,gauss);
    9391        //      this->GetM(M,element,gauss);
    … …  
    136134
    137135        Gauss* gauss = element->NewGauss(5);
    138         for(int ig=gauss->begin();ig<gauss->end();ig++){
    139                 gauss->GaussPoint(ig);
     136        while(gauss->next()){
    140137
    141138                element->JacobianDeterminant(&Jdet,xyz_list,gauss);
    … …  
    251248                delete gauss;
    252249                gauss = element->NewGaussBase(3);
    253                 for(int ig=gauss->begin();ig<gauss->end();ig++){
    254                         gauss->GaussPoint(ig);
     250                while(gauss->next()){
    255251
    256252                        /*Compute Jlambda*/
Note: See TracChangeset for help on using the changeset viewer.