Changeset 17543


Ignore:
Timestamp:
03/26/14 09:25:12 (11 years ago)
Author:
seroussi
Message:

NEW: added gauss points for segments up to 7 order

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

Legend:

Unmodified
Added
Removed

  • issm/trunk-jpl/src/c/shared/Numerics/GaussPoints.cpp

    r17479 r17543  
    4242        static IssmPDouble xi4[]={-0.861136311594053,-0.339981043584856, 0.339981043584856, 0.861136311594053};
    4343
    44         static IssmPDouble* wgtp[MAX_LINE_GAUS_PTS]={wgt1 ,wgt2 ,wgt3 ,wgt4 };
    45         static IssmPDouble* xip [MAX_LINE_GAUS_PTS]={xi1  ,xi2  ,xi3  ,xi4  };
     44        /*npoint= 5*/
     45        static IssmPDouble wgt5[]={0.5688888888888889,0.4786286704993665,0.4786286704993665,0.2369268850561891,0.2369268850561891};
     46        static IssmPDouble xi5[]={0.0000000000000000,-0.5384693101056831,0.5384693101056831,-0.9061798459386640,0.9061798459386640};
     47
     48        /*npoint= 6*/
     49        static IssmPDouble wgt6[]={0.3607615730481386,0.3607615730481386,0.4679139345726910,0.4679139345726910,0.1713244923791704,0.1713244923791704};
     50        static IssmPDouble xi6[]={-0.6612093864662645,0.6612093864662645,-0.2386191860831969,0.2386191860831969,-0.9324695142031521,0.9324695142031521};
     51
     52        /*npoint= 7*/
     53        static IssmPDouble wgt7[]={0.4179591836734694,0.3818300505051189,0.3818300505051189,0.2797053914892766,0.2797053914892766,0.1294849661688697,0.1294849661688697};
     54        static IssmPDouble xi7[]={0.0000000000000000,-0.4058451513773972,0.4058451513773972,-0.7415311855993945,0.7415311855993945,-0.9491079123427585,0.9491079123427585};
     55
     56        static IssmPDouble* wgtp[MAX_LINE_GAUS_PTS]={wgt1 ,wgt2 ,wgt3 ,wgt4, wgt5, wgt6, wgt7};
     57        static IssmPDouble* xip [MAX_LINE_GAUS_PTS]={xi1  ,xi2  ,xi3  ,xi4, xi5, xi6, xi7};
    4658
    4759        *pxgaus =xNew<IssmPDouble>(ngaus);

  • issm/trunk-jpl/src/c/shared/Numerics/GaussPoints.h

    r17479 r17543  
    77#define _GAUSSPOINTS_H
    88
    9 #define MAX_LINE_GAUS_PTS    4
     9#define MAX_LINE_GAUS_PTS    7
    1010void GaussLegendreLinear(IssmPDouble** pxgaus, IssmPDouble** pxwgt, int ngaus);
    1111#define MAX_TRIA_SYM_ORD    20
Note: See TracChangeset for help on using the changeset viewer.