Changeset 25437


Ignore:
Timestamp:
08/21/20 08:57:26 (5 years ago)
Author:
Mathieu Morlighem
Message:

CHG: fixing AD compilation

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

Legend:

Unmodified
Added
Removed

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

    r25379 r25437  
    614614        /* Start  looping on the number of gaussian points: */
    615615        Gauss* gauss=element->NewGauss(4);
    616         for(int ig=gauss->begin();ig<gauss->end();ig++){
    617                 gauss->GaussPoint(ig);
     616        while(gauss->next()){
    618617
    619618                element->JacobianDeterminant(&Jdet,xyz_list,gauss);
    … …  
    752751        /* Start  looping on the number of gaussian points: */
    753752        Gauss* gauss=element->NewGaussBase(4);
    754         for(int ig=gauss->begin();ig<gauss->end();ig++){
    755                 gauss->GaussPoint(ig);
     753        while(gauss->next()){
    756754
    757755                element->JacobianDeterminantBase(&Jdet,xyz_list_base,gauss);
    … …  
    836834        /* Start  looping on the number of gaussian points: */
    837835        Gauss* gauss=element->NewGauss(4);
    838         for(int ig=gauss->begin();ig<gauss->end();ig++){
    839                 gauss->GaussPoint(ig);
     836        while(gauss->next()){
    840837
    841838                element->JacobianDeterminant(&Jdet,xyz_list,gauss);
    … …  
    964961        Gauss* gauss=element->NewGaussBase(4);
    965962        Gauss* gaussup=element->NewGaussTop(4);
    966         for(int ig=gauss->begin();ig<gauss->end();ig++){
    967                 gauss->GaussPoint(ig);
    968                 gaussup->GaussPoint(ig);
     963        while(gauss->next() && gaussup->next()){
    969964
    970965                element->JacobianDeterminantBase(&Jdet,xyz_list_base,gauss);
    … …  
    10511046        /* Start  looping on the number of gaussian points: */
    10521047        Gauss* gauss=element->NewGaussBase(4);
    1053         for(int ig=gauss->begin();ig<gauss->end();ig++){
    1054                 gauss->GaussPoint(ig);
     1048        while(gauss->next()){
    10551049
    10561050                element->JacobianDeterminantBase(&Jdet,xyz_list_base,gauss);

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

    r25379 r25437  
    259259        /* Start  looping on the number of gaussian points: */
    260260        Gauss* gauss=element->NewGauss(2);
    261         for(int ig=gauss->begin();ig<gauss->end();ig++){
    262                 gauss->GaussPoint(ig);
     261        while(gauss->next()){
    263262
    264263                element->JacobianDeterminant(&Jdet,xyz_list,gauss);
    … …  
    367366        /* Start  looping on the number of gaussian points: */
    368367        Gauss* gauss=element->NewGauss(2);
    369         for(int ig=gauss->begin();ig<gauss->end();ig++){
    370                 gauss->GaussPoint(ig);
     368        while(gauss->next()){
    371369
    372370                element->JacobianDeterminant(&Jdet,xyz_list,gauss);

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

    r25392 r25437  
    196196        /* Start  looping on the number of gaussian points: */
    197197        Gauss* gauss=element->NewGauss(1);
    198         for(int ig=gauss->begin();ig<gauss->end();ig++){
    199                 gauss->GaussPoint(ig);
     198        while(gauss->next()){
    200199
    201200                element->JacobianDeterminant(&Jdet,xyz_list,gauss);
    … …  
    271270        /* Start  looping on the number of gaussian points: */
    272271        Gauss* gauss=element->NewGauss(2);
    273         for(int ig=gauss->begin();ig<gauss->end();ig++){
    274                 gauss->GaussPoint(ig);
     272        while(gauss->next()){
    275273
    276274                element->JacobianDeterminant(&Jdet,xyz_list,gauss);
    … …  
    509507        /* Start  looping on the number of gaussian points: */
    510508        Gauss* gauss=element->NewGauss(2);
    511         for(int ig=gauss->begin();ig<gauss->end();ig++){
    512                 gauss->GaussPoint(ig);
     509        while(gauss->next()){
    513510
    514511                element->JacobianDeterminant(&Jdet,xyz_list,gauss);

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

    r25436 r25437  
    13201320        /* Start  looping on the number of gaussian points: */
    13211321        Gauss* gauss = basalelement->NewGauss(2);
    1322         for(int ig=gauss->begin();ig<gauss->end();ig++){
    1323                 gauss->GaussPoint(ig);
     1322        while(gauss->next()){
    13241323
    13251324                basalelement->JacobianDeterminant(&Jdet,xyz_list,gauss);
    … …  
    20892088        Gauss* gauss      = element->NewGauss(5);
    20902089        Gauss* gauss_base = basalelement->NewGauss();
    2091         for(int ig=gauss->begin();ig<gauss->end();ig++){
    2092 
    2093                 gauss->GaussPoint(ig);
     2090        while(gauss->next()){
    20942091                gauss->SynchronizeGaussBase(gauss_base);
    20952092
    … …  
    21752172        /* Start  looping on the number of gaussian points: */
    21762173        Gauss* gauss=element->NewGauss(2);
    2177         for(int ig=gauss->begin();ig<gauss->end();ig++){
    2178 
    2179                 gauss->GaussPoint(ig);
    2180 
     2174        while(gauss->next()){
    21812175                element->JacobianDeterminant(&Jdet,xyz_list,gauss);
    21822176                element->NodalFunctions(basis, gauss);
    … …  
    22322226        /*Start looping on Gaussian points*/
    22332227        Gauss* gauss=element->NewGauss(xyz_list,xyz_list_front,3);
    2234         for(int ig=gauss->begin();ig<gauss->end();ig++){
    2235 
    2236                 gauss->GaussPoint(ig);
     2228        while(gauss->next()){
    22372229                thickness_input->GetInputValue(&thickness,gauss);
    22382230                base_input->GetInputValue(&bed,gauss);
    … …  
    23842376        /* Start  looping on the number of gaussian points: */
    23852377        Gauss* gauss = element->NewGauss(5);
    2386         for(int ig=gauss->begin();ig<gauss->end();ig++){
    2387                 gauss->GaussPoint(ig);
     2378        while(gauss->next()){
    23882379
    23892380                element->JacobianDeterminant(&Jdet,xyz_list,gauss);
    … …  
    24822473
    24832474        /* Start  looping on the number of gaussian points: */
    2484         for(int ig=gauss->begin();ig<gauss->end();ig++){
    2485                 gauss->GaussPoint(ig);
     2475        while(gauss->next()){
    24862476
    24872477                friction->GetAlpha2(&alpha2,gauss);
    … …  
    25572547        /* Start  looping on the number of gaussian points: */
    25582548        Gauss* gauss = element->NewGauss(5);
    2559         for(int ig=gauss->begin();ig<gauss->end();ig++){
    2560                 gauss->GaussPoint(ig);
     2549        while(gauss->next()){
    25612550
    25622551                element->JacobianDeterminant(&Jdet,xyz_list,gauss);
    … …  
    26292618        /* Start  looping on the number of gaussian points: */
    26302619        Gauss* gauss=element->NewGauss(3);
    2631         for(int ig=gauss->begin();ig<gauss->end();ig++){
    2632                 gauss->GaussPoint(ig);
     2620        while(gauss->next()){
    26332621
    26342622                element->JacobianDeterminant(&Jdet,xyz_list,gauss);
    … …  
    27022690        /* Start  looping on the number of gaussian points: */
    27032691        Gauss* gauss=element->NewGauss(3);
    2704         for(int ig=gauss->begin();ig<gauss->end();ig++){
    2705                 gauss->GaussPoint(ig);
     2692        while(gauss->next()){
    27062693
    27072694                element->JacobianDeterminant(&Jdet,xyz_list,gauss);
    … …  
    27692756
    27702757        /* Start  looping on the number of gaussian points: */
    2771         for(int ig=gauss->begin();ig<gauss->end();ig++){
    2772 
    2773                 gauss->GaussPoint(ig);
     2758        while(gauss->next()){
    27742759                surface_input->GetInputValue(&surface,gauss);
    27752760                sealevel_input->GetInputValue(&sealevel,gauss);
    … …  
    30773062        /* Start  looping on the number of gaussian points: */
    30783063        Gauss* gauss = element->NewGauss(5);
    3079         for(int ig=gauss->begin();ig<gauss->end();ig++){
    3080                 gauss->GaussPoint(ig);
     3064        while(gauss->next()){
    30813065
    30823066                element->JacobianDeterminant(&Jdet,xyz_list,gauss);
    … …  
    31983182        /* Start  looping on the number of gaussian points: */
    31993183        gauss=element->NewGaussBase(3);
    3200         for(int ig=gauss->begin();ig<gauss->end();ig++){
    3201                 gauss->GaussPoint(ig);
     3184        while(gauss->next()){
    32023185
    32033186                base_input->GetInputDerivativeValue(&slope[0],xyz_list,gauss);
    … …  
    32593242        /* Start  looping on the number of gaussian points: */
    32603243        Gauss* gauss = element->NewGauss(5);
    3261         for(int ig=gauss->begin();ig<gauss->end();ig++){
    3262                 gauss->GaussPoint(ig);
     3244        while(gauss->next()){
    32633245
    32643246                element->JacobianDeterminant(&Jdet,xyz_list,gauss);
    … …  
    33563338        /* Start  looping on the number of gaussian points: */
    33573339        Gauss* gauss = element->NewGauss(5);
    3358         for(int ig=gauss->begin();ig<gauss->end();ig++){
    3359                 gauss->GaussPoint(ig);
     3340        while(gauss->next()){
    33603341
    33613342                element->JacobianDeterminant(&Jdet,xyz_list,gauss);
    … …  
    34103391        /* Start  looping on the number of gaussian points: */
    34113392        Gauss* gauss = element->NewGauss(5);
    3412         for(int ig=gauss->begin();ig<gauss->end();ig++){
    3413                 gauss->GaussPoint(ig);
     3393        while(gauss->next()){
    34143394
    34153395                element->JacobianDeterminant(&Jdet,xyz_list,gauss);
    … …  
    34833463        /* Start  looping on the number of gaussian points: */
    34843464        Gauss* gauss = element->NewGauss(5);
    3485         for(int ig=gauss->begin();ig<gauss->end();ig++){
    3486                 gauss->GaussPoint(ig);
     3465        while(gauss->next()){
    34873466
    34883467                element->JacobianDeterminant(&Jdet,xyz_list,gauss);
    … …  
    35263505                delete gauss;
    35273506                gauss = element->NewGaussBase(5);
    3528                 for(int ig=gauss->begin();ig<gauss->end();ig++){
    3529                         gauss->GaussPoint(ig);
     3507                while(gauss->next()){
    35303508
    35313509                        element->JacobianDeterminantBase(&Jdet,xyz_list_base,gauss);
    … …  
    36083586        /* Start  looping on the number of gaussian points: */
    36093587        Gauss* gauss = element->NewGauss(5);
    3610         for(int ig=gauss->begin();ig<gauss->end();ig++){
    3611                 gauss->GaussPoint(ig);
     3588        while(gauss->next()){
    36123589
    36133590                element->JacobianDeterminant(&Jdet,xyz_list,gauss);
    … …  
    36773654        /* Start  looping on the number of gaussian points: */
    36783655        gauss=element->NewGaussBase(10);
    3679         for(int ig=gauss->begin();ig<gauss->end();ig++){
    3680                 gauss->GaussPoint(ig);
     3656        while(gauss->next()){
    36813657
    36823658                x_coord=element->GetXcoord(xyz_list,gauss);
    … …  
    37523728        /* Start  looping on the number of gaussian points: */
    37533729        gauss=element->NewGaussBase(3);
    3754         for(int ig=gauss->begin();ig<gauss->end();ig++){
    3755                 gauss->GaussPoint(ig);
     3730        while(gauss->next()){
    37563731
    37573732                alpha2_input->GetInputValue(&alpha2, gauss);
    … …  
    38083783        /* Start  looping on the number of gaussian points: */
    38093784        gauss=element->NewGaussBase(3);
    3810         for(int ig=gauss->begin();ig<gauss->end();ig++){
    3811                 gauss->GaussPoint(ig);
     3785        while(gauss->next()){
    38123786
    38133787                sigmann_input->GetInputValue(&sigmann, gauss);
    … …  
    38663840        /* Start  looping on the number of gaussian points: */
    38673841        Gauss* gauss=element->NewGauss(5);
    3868         for(int ig=gauss->begin();ig<gauss->end();ig++){
    3869                 gauss->GaussPoint(ig);
     3842        while(gauss->next()){
    38703843
    38713844                element->JacobianDeterminant(&Jdet,xyz_list,gauss);
    … …  
    39663939
    39673940        /* Start  looping on the number of gaussian points: */
    3968         for(int ig=gauss->begin();ig<gauss->end();ig++){
    3969                 gauss->GaussPoint(ig);
     3941        while(gauss->next()){
    39703942
    39713943                friction->GetAlpha2(&alpha2,gauss);
    … …  
    40974069        /* Start  looping on the number of gaussian points: */
    40984070        Gauss* gauss=element->NewGauss(5);
    4099         for(int ig=gauss->begin();ig<gauss->end();ig++){
    4100                 gauss->GaussPoint(ig);
     4071        while(gauss->next()){
    41014072
    41024073                element->JacobianDeterminant(&Jdet,xyz_list,gauss);
    … …  
    41784149
    41794150        /* Start  looping on the number of gaussian points: */
    4180         for(int ig=gauss->begin();ig<gauss->end();ig++){
    4181                 gauss->GaussPoint(ig);
     4151        while(gauss->next()){
    41824152
    41834153                element->JacobianDeterminantSurface(&Jdet,xyz_list_front,gauss);
    … …  
    42414211        /* Start  looping on the number of gaussian points: */
    42424212        Gauss* gauss=element->NewGaussBase(5);
    4243         for(int ig=gauss->begin();ig<gauss->end();ig++){
    4244                 gauss->GaussPoint(ig);
     4213        while(gauss->next()){
    42454214
    42464215                element->JacobianDeterminantBase(&Jdet,xyz_list_base,gauss);
    … …  
    42624231                IssmDouble  dt,mb;
    42634232                element->FindParam(&dt,TimesteppingTimeStepEnum);
    4264                 for(int ig=gauss->begin();ig<gauss->end();ig++){
    4265                         gauss->GaussPoint(ig);
     4233                while(gauss->next()){
    42664234                        element->JacobianDeterminantBase(&Jdet,xyz_list_base,gauss);
    42674235                        element->NodalFunctionsVelocity(vbasis,gauss);
    … …  
    43144282
    43154283        gauss=element->NewGauss(5);
    4316         for(int ig=gauss->begin();ig<gauss->end();ig++){
    4317                 gauss->GaussPoint(ig);
     4284        while(gauss->next()){
    43184285                element->JacobianDeterminant(&Jdet,xyz_list,gauss);
    43194286
    … …  
    43374304                delete gauss;
    43384305                gauss=element->NewGaussBase(5);
    4339                 for(int ig=gauss->begin();ig<gauss->end();ig++){
    4340                         gauss->GaussPoint(ig);
     4306                while(gauss->next()){
    43414307
    43424308                        element->JacobianDeterminantBase(&Jdet,xyz_list_base,gauss);
    … …  
    44484414        delete gauss;
    44494415        gauss=element->NewGauss(5);
    4450         for(int ig=gauss->begin();ig<gauss->end();ig++){
    4451                 gauss->GaussPoint(ig);
     4416        while(gauss->next()){
    44524417                element->JacobianDeterminant(&Jdet,xyz_list,gauss);
    44534418
    … …  
    53805345
    53815346                Gauss* gauss=element->NewGauss(5);
    5382                 for(int ig=gauss->begin();ig<gauss->end();ig++){
    5383                         gauss->GaussPoint(ig);
     5347                while(gauss->next()){
    53845348                        element->JacobianDeterminant(&Jdet,xyz_list,gauss);
    53855349                        element->NodalFunctionsTensor(tbasis,gauss);
    … …  
    58135777        /* Start  looping on the number of gaussian points: */
    58145778        Gauss* gauss=element->NewGaussBase(2);
    5815         for(int ig=gauss->begin();ig<gauss->end();ig++){
    5816 
    5817                 gauss->GaussPoint(ig);
     5779        while(gauss->next()){
    58185780
    58195781                element->JacobianDeterminantBase(&Jdet2d,xyz_list_tria,gauss);
    … …  
    59325894        Gauss* gauss=element->NewGauss(5);
    59335895        Gauss* gauss_tria=new GaussTria();
    5934         for(int ig=gauss->begin();ig<gauss->end();ig++){
    5935 
    5936                 gauss->GaussPoint(ig);
     5896        while(gauss->next()){
    59375897                gauss->SynchronizeGaussBase(gauss_tria);
    59385898
    … …  
    60345994        /* Start  looping on the number of gaussian points: */
    60355995        Gauss* gauss=element->NewGaussBase(2);
    6036         for(int ig=gauss->begin();ig<gauss->end();ig++){
    6037 
    6038                 gauss->GaussPoint(ig);
     5996        while(gauss->next()){
    60395997
    60405998                /*Friction: */
    … …  
    61186076        Gauss* gauss=element->NewGauss(5);
    61196077        Gauss* gauss_tria=new GaussTria();
    6120         for(int ig=gauss->begin();ig<gauss->end();ig++){
    6121 
    6122                 gauss->GaussPoint(ig);
     6078        while(gauss->next()){
    61236079                gauss->SynchronizeGaussBase(gauss_tria);
    61246080
    … …  
    62766232        Gauss* gauss=element->NewGauss(5);
    62776233        Gauss* gauss_tria=new GaussTria();
    6278         for(int ig=gauss->begin();ig<gauss->end();ig++){
    6279 
    6280                 gauss->GaussPoint(ig);
     6234        while(gauss->next()){
    62816235                gauss->SynchronizeGaussBase(gauss_tria);
    62826236
    … …  
    63816335        /* Start looping on the number of gauss 2d (nodes on the bedrock) */
    63826336        Gauss* gauss=element->NewGaussBase(2);
    6383         for(int ig=gauss->begin();ig<gauss->end();ig++){
    6384 
    6385                 gauss->GaussPoint(ig);
     6337        while(gauss->next()){
    63866338
    63876339                element->JacobianDeterminantBase(&Jdet2d,xyz_list_tria,gauss);
    … …  
    64576409        /* Start  looping on the number of gaussian points: */
    64586410        Gauss* gauss=element->NewGauss(5);
    6459         for(int ig=gauss->begin();ig<gauss->end();ig++){
    6460 
    6461                 gauss->GaussPoint(ig);
     6411        while(gauss->next()){
    64626412
    64636413                element->JacobianDeterminant(&Jdet, xyz_list,gauss);
    … …  
    65486498        /* Start looping on the number of gauss 2d (nodes on the bedrock) */
    65496499        Gauss* gauss=element->NewGaussBase(2);
    6550         for(int ig=gauss->begin();ig<gauss->end();ig++){
    6551 
    6552                 gauss->GaussPoint(ig);
     6500        while(gauss->next()){
    65536501
    65546502                element->JacobianDeterminantBase(&Jdet2d,xyz_list_tria,gauss);
    … …  
    66236571        /* Start  looping on the number of gaussian points: */
    66246572        Gauss* gauss=element->NewGauss(5);
    6625         for(int ig=gauss->begin();ig<gauss->end();ig++){
    6626 
    6627                 gauss->GaussPoint(ig);
     6573        while(gauss->next()){
    66286574                element->JacobianDeterminant(&Jdet,xyz_list,gauss);
    66296575                element->NodalFunctionsP1(&basis[0], gauss);

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

    r25379 r25437  
    209209        Gauss* gauss = element->NewGaussBase(2);
    210210        element->NormalBase(&normal[0],xyz_list);
    211         for(int ig=gauss->begin();ig<gauss->end();ig++){
    212                 gauss->GaussPoint(ig);
     211        while(gauss->next()){
    213212
    214213                element->JacobianDeterminantBase(&Jdet,xyz_list,gauss);
    … …  
    246245        Gauss* gauss = element->NewGaussTop(2);
    247246        element->NormalTop(&normal[0],xyz_list);
    248         for(int ig=gauss->begin();ig<gauss->end();ig++){
    249                 gauss->GaussPoint(ig);
     247        while(gauss->next()){
    250248
    251249                element->JacobianDeterminantTop(&Jdet,xyz_list,gauss);
    … …  
    281279        /* Start  looping on the number of gaussian points: */
    282280        Gauss* gauss = element->NewGauss(2);
    283         for(int ig=gauss->begin();ig<gauss->end();ig++){
    284                 gauss->GaussPoint(ig);
     281        while(gauss->next()){
    285282
    286283                element->JacobianDeterminant(&Jdet,xyz_list,gauss);
    … …  
    360357        /* Start  looping on the number of gaussian points: */
    361358        Gauss* gauss=element->NewGaussBase(2);
    362         for(int ig=gauss->begin();ig<gauss->end();ig++){
    363                 gauss->GaussPoint(ig);
     359        while(gauss->next()){
    364360
    365361                groundedice_melting_input->GetInputValue(&gmb,gauss);
    … …  
    424420        /* Start  looping on the number of gaussian points: */
    425421        Gauss* gauss=element->NewGaussTop(2);
    426         for(int ig=gauss->begin();ig<gauss->end();ig++){
    427                 gauss->GaussPoint(ig);
     422        while(gauss->next()){
    428423
    429424                smb_input->GetInputValue(&smb,gauss);
    … …  
    477472        /* Start  looping on the number of gaussian points: */
    478473        Gauss* gauss=element->NewGauss(2);
    479         for(int ig=gauss->begin();ig<gauss->end();ig++){
    480                 gauss->GaussPoint(ig);
     474        while(gauss->next()){
    481475
    482476                element->JacobianDeterminant(&Jdet,xyz_list,gauss);

  • issm/trunk-jpl/src/c/classes/Cfdragcoeffabsgrad.cpp

    r25379 r25437  
    156156        /* Start  looping on the number of gaussian points: */
    157157        Gauss* gauss=basalelement->NewGauss(2);
    158         for(int ig=gauss->begin();ig<gauss->end();ig++){
    159 
    160                 gauss->GaussPoint(ig);
     158        while(gauss->next()){
    161159
    162160                /* Get Jacobian determinant: */

  • issm/trunk-jpl/src/c/modules/SurfaceMassBalancex/Gembx.cpp

    r25435 r25437  
    21012101                                // ERA5 new aIdx=1, swIdx=0
    21022102                                if (aIdx==1 && swIdx==0){
    2103                                         if (abs(adThresh - 820) < Dtol){
     2103                                        if (fabs(adThresh - 820) < Dtol){
    21042104                                                // ERA5 new aIdx=1, swIdx=0, MODIS 820
    21052105                                                M0 = max(1.8230 - (0.1645 * log(C)),0.25);
Note: See TracChangeset for help on using the changeset viewer.