Changeset 25445


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

CHG: moving some gauss to next function

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

Legend:

Unmodified
Added
Removed

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

    r25406 r25445  
    520520                /*Start looping on Gaussian points*/
    521521                Gauss* gauss=this->NewGaussBase(xyz_list,&xyz_front[0][0],3);
    522                 for(int ig=gauss->begin();ig<gauss->end();ig++){
    523 
    524                         gauss->GaussPoint(ig);
     522                while(gauss->next()){
    525523                        thickness_input->GetInputValue(&thickness,gauss);
    526524                        calvingratex_input->GetInputValue(&calvingratex,gauss);
    … …  
    648646                /*Start looping on Gaussian points*/
    649647                Gauss* gauss=this->NewGaussBase(xyz_list,&xyz_front[0][0],3);
    650                 for(int ig=gauss->begin();ig<gauss->end();ig++){
    651 
    652                         gauss->GaussPoint(ig);
     648                while(gauss->next()){
    653649                        thickness_input->GetInputValue(&thickness,gauss);
    654650                        calvingratex_input->GetInputValue(&calvingratex,gauss);
    … …  
    728724        /* Start  looping on the number of gaussian points: */
    729725        gauss=new GaussPenta(0,1,2,2);
    730         for(int ig=gauss->begin();ig<gauss->end();ig++){
    731 
    732                 gauss->GaussPoint(ig);
     726        while(gauss->next()){
    733727
    734728                /*Compute strain rate viscosity and pressure: */
    … …  
    20502044        /*Start looping on Gaussian points*/
    20512045        Gauss* gauss=this->NewGaussBase(xyz_list,&xyz_front[0][0],3);
    2052         for(int ig=gauss->begin();ig<gauss->end();ig++){
    2053 
    2054                 gauss->GaussPoint(ig);
     2046        while(gauss->next()){
    20552047                thickness_input->GetInputValue(&thickness,gauss);
    20562048                vx_input->GetInputValue(&vx,gauss);
    … …  
    21722164        /*Start looping on Gaussian points*/
    21732165        Gauss* gauss=this->NewGaussBase(xyz_list,&xyz_front[0][0],3);
    2174         for(int ig=gauss->begin();ig<gauss->end();ig++){
    2175 
    2176                 gauss->GaussPoint(ig);
     2166        while(gauss->next()){
    21772167                thickness_input->GetInputValue(&thickness,gauss);
    21782168                vx_input->GetInputValue(&vx,gauss);
    … …  
    22892279                        }
    22902280                        /*Integrate over edge*/
    2291                         for(int ig=gauss[iv]->begin();ig<gauss[iv]->end();ig++){
    2292                                 gauss[iv]->GaussPoint(ig);
     2281                        gauss[iv]->Reset();
     2282                        while(gauss[iv]->next()){
    22932283                                penta->JacobianDeterminantLine(&Jdet,&xyz_list_line[0][0],gauss[iv]);
    22942284                                original_input->GetInputValue(&value,gauss[iv]);
    … …  
    36193609                for(int ig=0;ig<3;ig++){
    36203610                        GaussPenta* gauss=new GaussPenta(ig,3+ig,11);
    3621                         for (int iv=gauss->begin();iv<gauss->end();iv++){
    3622                                 gauss->GaussPoint(iv);
     3611                        while(gauss->next()){
    36233612
    36243613                                /* Get the value we need*/
    … …  
    38393828        /*Start looping on Gaussian points*/
    38403829        Gauss* gauss=this->NewGaussBase(xyz_list,&xyz_front[0][0],3);
    3841         for(int ig=gauss->begin();ig<gauss->end();ig++){
    3842 
    3843                 gauss->GaussPoint(ig);
     3830        while(gauss->next()){
    38443831                thickness_input->GetInputValue(&thickness,gauss);
    38453832                calvingratex_input->GetInputValue(&calvingratex,gauss);
    … …  
    39663953        /*Start looping on Gaussian points*/
    39673954        Gauss* gauss=this->NewGaussBase(xyz_list,&xyz_front[0][0],3);
    3968         for(int ig=gauss->begin();ig<gauss->end();ig++){
    3969 
    3970                 gauss->GaussPoint(ig);
     3955        while(gauss->next()){
    39713956                thickness_input->GetInputValue(&thickness,gauss);
    39723957                calvingratex_input->GetInputValue(&calvingratex,gauss);
    … …  
    40154000        /* Start  looping on the number of gaussian points: */
    40164001        gauss = this->NewGauss(point1,fraction1,fraction2,1-mainlyfloating,3);
    4017         for(int ig=gauss->begin();ig<gauss->end();ig++){
    4018 
    4019                 gauss->GaussPoint(ig);
     4002        while(gauss->next()){
    40204003                this->JacobianDeterminantBase(&Jdet,&xyz_list[0][0],gauss);
    40214004                floatingmelt_input->GetInputValue(&floatingmelt,gauss);
    … …  
    40604043        /* Start  looping on the number of gaussian points: */
    40614044        gauss = this->NewGauss(point1,fraction1,fraction2,mainlyfloating,3);
    4062         for(int ig=gauss->begin();ig<gauss->end();ig++){
    4063 
    4064                 gauss->GaussPoint(ig);
     4045        while(gauss->next()){
    40654046                this->JacobianDeterminantBase(&Jdet,&xyz_list[0][0],gauss);
    40664047                groundedmelt_input->GetInputValue(&groundedmelt,gauss);

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

    r25442 r25445  
    672672                /*Start looping on Gaussian points*/
    673673                Gauss* gauss=this->NewGauss(xyz_list,&xyz_front[0][0],3);
    674                 for(int ig=gauss->begin();ig<gauss->end();ig++){
    675 
    676                         gauss->GaussPoint(ig);
     674                while(gauss->next()){
    677675                        thickness_input->GetInputValue(&thickness,gauss);
    678676                        calvingratex_input->GetInputValue(&calvingratex,gauss);
    … …  
    815813                /*Start looping on Gaussian points*/
    816814                Gauss* gauss=this->NewGauss(xyz_list,&xyz_front[0][0],3);
    817                 for(int ig=gauss->begin();ig<gauss->end();ig++){
    818 
    819                         gauss->GaussPoint(ig);
     815                while(gauss->next()){
    820816                        thickness_input->GetInputValue(&thickness,gauss);
    821817                        calvingratex_input->GetInputValue(&calvingratex,gauss);
    … …  
    26452641        /*Start looping on Gaussian points*/
    26462642        Gauss* gauss=this->NewGauss(xyz_list,xyz_front,3);
    2647         for(int ig=gauss->begin();ig<gauss->end();ig++){
    2648 
    2649                 gauss->GaussPoint(ig);
     2643        while(gauss->next()){
    26502644                thickness_input->GetInputValue(&thickness,gauss);
    26512645                vx_input->GetInputValue(&vx,gauss);
    … …  
    27742768        /*Start looping on Gaussian points*/
    27752769        Gauss* gauss=this->NewGauss(xyz_list,&xyz_front[0][0],3);
    2776         for(int ig=gauss->begin();ig<gauss->end();ig++){
    2777 
    2778                 gauss->GaussPoint(ig);
     2770        while(gauss->next()){
    27792771                thickness_input->GetInputValue(&thickness,gauss);
    27802772                vx_input->GetInputValue(&vx,gauss);
    … …  
    29012893        /*Start looping on Gaussian points*/
    29022894        Gauss* gauss=this->NewGauss(xyz_list,&xyz_front[0][0],3);
    2903         for(int ig=gauss->begin();ig<gauss->end();ig++){
    2904 
    2905                 gauss->GaussPoint(ig);
     2895        while(gauss->next()){
    29062896                thickness_input->GetInputValue(&thickness,gauss);
    29072897                vx_input->GetInputValue(&vx,gauss);
    … …  
    33673357        volume=0;
    33683358
    3369         for(int ig=gauss->begin();ig<gauss->end();ig++){
    3370                 gauss->GaussPoint(ig);
     3359        while(gauss->next()){
    33713360
    33723361                this->JacobianDeterminant(&Jdet,xyz_list,gauss);
    … …  
    34663455        /* Start  looping on the number of gaussian points: */
    34673456        gauss=new GaussTria(2);
    3468         for(int ig=gauss->begin();ig<gauss->end();ig++){
    3469 
    3470                 gauss->GaussPoint(ig);
     3457        while(gauss->next()){
    34713458
    34723459                /* Get Jacobian determinant: */
    … …  
    35053492        /* Start  looping on the number of gaussian points: */
    35063493        gauss=new GaussTria(2);
    3507         for(int ig=gauss->begin();ig<gauss->end();ig++){
    3508 
    3509                 gauss->GaussPoint(ig);
     3494        while(gauss->next()){
    35103495
    35113496                /* Get Jacobian determinant: */
    … …  
    43814366        /*Start looping on Gaussian points*/
    43824367        Gauss* gauss=this->NewGauss(xyz_list,&xyz_front[0][0],3);
    4383         for(int ig=gauss->begin();ig<gauss->end();ig++){
    4384 
    4385                 gauss->GaussPoint(ig);
     4368        while(gauss->next()){
    43864369                thickness_input->GetInputValue(&thickness,gauss);
    43874370                calvingratex_input->GetInputValue(&calvingratex,gauss);
    … …  
    45164499        /*Start looping on Gaussian points*/
    45174500        Gauss* gauss=this->NewGauss(xyz_list,&xyz_front[0][0],3);
    4518         for(int ig=gauss->begin();ig<gauss->end();ig++){
    4519 
    4520                 gauss->GaussPoint(ig);
     4501        while(gauss->next()){
    45214502                thickness_input->GetInputValue(&thickness,gauss);
    45224503                calvingratex_input->GetInputValue(&calvingratex,gauss);
    … …  
    45624543        /* Start  looping on the number of gaussian points: */
    45634544        gauss = this->NewGauss(point1,fraction1,fraction2,1-mainlyfloating,3);
    4564         for(int ig=gauss->begin();ig<gauss->end();ig++){
    4565 
    4566                 gauss->GaussPoint(ig);
     4545        while(gauss->next()){
    45674546                this->JacobianDeterminant(&Jdet,&xyz_list[0][0],gauss);
    45684547                floatingmelt_input->GetInputValue(&floatingmelt,gauss);
    … …  
    46074586        /* Start  looping on the number of gaussian points: */
    46084587        gauss = this->NewGauss(point1,fraction1,fraction2,mainlyfloating,2);
    4609         for(int ig=gauss->begin();ig<gauss->end();ig++){
    4610 
    4611                 gauss->GaussPoint(ig);
     4588        while(gauss->next()){
    46124589                this->JacobianDeterminant(&Jdet,&xyz_list[0][0],gauss);
    46134590                groundedmelt_input->GetInputValue(&groundedmelt,gauss);
    … …  
    57705747                total_weight=0;
    57715748                I=0;
    5772                 for(int ig=gauss->begin();ig<gauss->end();ig++){
     5749                while(gauss->next()){
    57735750                        IssmDouble Ig=0;
    5774                         gauss->GaussPoint(ig);
    57755751                        deltathickness_input->GetInputValue(&Ig,gauss);
    57765752                        I+=Ig*gauss->weight;
Note: See TracChangeset for help on using the changeset viewer.