Changeset 19637


Ignore:
Timestamp:
10/14/15 14:36:29 (9 years ago)
Author:
seroussi
Message:

BUG: increased number of gauss points for thermal and enthalpy

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

Legend:

Unmodified
Added
Removed

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

    r19162 r19637  
    530530
    531531        /* Start  looping on the number of gaussian points: */
    532         Gauss* gauss=element->NewGauss(2);
     532        Gauss* gauss=element->NewGauss(4);
    533533        for(int ig=gauss->begin();ig<gauss->end();ig++){
    534534                gauss->GaussPoint(ig);
    … …  
    648648
    649649        /* Start  looping on the number of gaussian points: */
    650         Gauss* gauss=element->NewGaussBase(2);
     650        Gauss* gauss=element->NewGaussBase(4);
    651651        for(int ig=gauss->begin();ig<gauss->end();ig++){
    652652                gauss->GaussPoint(ig);
    … …  
    737737
    738738        /* Start  looping on the number of gaussian points: */
    739         Gauss* gauss=element->NewGauss(3);
     739        Gauss* gauss=element->NewGauss(4);
    740740        for(int ig=gauss->begin();ig<gauss->end();ig++){
    741741                gauss->GaussPoint(ig);
    … …  
    844844
    845845        /* Start  looping on the number of gaussian points: */
    846         Gauss* gauss=element->NewGaussBase(2);
    847         Gauss* gaussup=element->NewGaussTop(2);
     846        Gauss* gauss=element->NewGaussBase(4);
     847        Gauss* gaussup=element->NewGaussTop(4);
    848848        for(int ig=gauss->begin();ig<gauss->end();ig++){
    849849                gauss->GaussPoint(ig);
    … …  
    938938
    939939        /* Start  looping on the number of gaussian points: */
    940         Gauss* gauss=element->NewGaussBase(2);
     940        Gauss* gauss=element->NewGaussBase(4);
    941941        for(int ig=gauss->begin();ig<gauss->end();ig++){
    942942                gauss->GaussPoint(ig);

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

    r19527 r19637  
    247247
    248248        /* Start  looping on the number of gaussian points: */
    249         Gauss* gauss=element->NewGaussBase(2);
     249        Gauss* gauss=element->NewGaussBase(4);
    250250        for(int ig=gauss->begin();ig<gauss->end();ig++){
    251251                gauss->GaussPoint(ig);
    … …  
    313313
    314314        /* Start  looping on the number of gaussian points: */
    315         Gauss* gauss=element->NewGauss(2);
     315        Gauss* gauss=element->NewGauss(4);
    316316        for(int ig=gauss->begin();ig<gauss->end();ig++){
    317317                gauss->GaussPoint(ig);
    … …  
    452452
    453453        /* Start  looping on the number of gaussian points: */
    454         Gauss* gauss   = element->NewGaussBase(2);
     454        Gauss* gauss   = element->NewGaussBase(4);
    455455        for(int ig=gauss->begin();ig<gauss->end();ig++){
    456456                gauss->GaussPoint(ig);
    … …  
    511511
    512512        /* Start  looping on the number of gaussian points: */
    513         Gauss* gauss=element->NewGaussBase(2);
     513        Gauss* gauss=element->NewGaussBase(4);
    514514        for(int ig=gauss->begin();ig<gauss->end();ig++){
    515515                gauss->GaussPoint(ig);
    … …  
    571571
    572572        /* Start  looping on the number of gaussian points: */
    573         Gauss* gauss=element->NewGauss(3);
     573        Gauss* gauss=element->NewGauss(4);
    574574        for(int ig=gauss->begin();ig<gauss->end();ig++){
    575575                gauss->GaussPoint(ig);
    … …  
    714714        int        *doflist   = NULL;
    715715        IssmDouble *xyz_list  = NULL;
    716         bool        hack      = false;
     716        bool        hack      = true;
    717717
    718718        /*Fetch number of nodes and dof for this finite element*/
Note: See TracChangeset for help on using the changeset viewer.