Changeset 3159


Ignore:
Timestamp:
03/02/10 14:18:24 (15 years ago)
Author:
Mathieu Morlighem
Message:

Added some macros and renamed GetElementNodeData -> GetVerticesCoordinates

Location:
issm/trunk/src/c
Files:
1 deleted
9 edited

Legend:

Unmodified
Added
Removed

  • issm/trunk/src/c/Makefile.am

    r3126 r3159  
    129129                                        ./shared/Elements/ResolvePointers.cpp\
    130130                                        ./shared/Elements/Paterson.cpp\
    131                                         ./shared/Elements/GetElementNodeData.cpp\
     131                                        ./shared/Elements/GetVerticesCoordinates.cpp\
    132132                                        ./shared/String/isdistributed.cpp\
    133133                                        ./shared/String/sharedstring.h\
    … …  
    456456                                        ./shared/Elements/ResolvePointers.cpp\
    457457                                        ./shared/Elements/Paterson.cpp\
    458                                         ./shared/Elements/GetElementNodeData.cpp\
     458                                        ./shared/Elements/GetVerticesCoordinates.cpp\
    459459                                        ./shared/String/isdistributed.cpp\
    460460                                        ./shared/String/sharedstring.h\

  • issm/trunk/src/c/include/typedefs.h

    r1 r3159  
    77
    88#define UNDEF -9999
     9#define ONETHIRD 0.333333333333333333333333333333333333333333333333333333333333
     10#define SQRT2 1.414213562373095048801688724209698078569671875376948073176679738
     11#define SQRT3 1.732050807568877293527446341505872366942805253810380628055806979
     12#define PI 3.141592653589793238462643383279502884197169399375105820974944592308
    913
    1014#endif //ifndef _ISSMTYPEDEFS_H_
    11 

  • issm/trunk/src/c/objects/Beam.cpp

    r3041 r3159  
    179179
    180180        /*Get node data: */
    181         GetElementNodeData( &xyz_list[0][0], nodes, numgrids);
     181        GetVerticesCoordinates(&xyz_list[0][0], nodes, numgrids);
    182182       
    183183        /*Get dof list on which we will plug the pressure values: */
    … …  
    369369
    370370        //Get all element grid data:
    371         GetElementNodeData( &xyz_list[0][0], nodes, numgrids);
     371        GetVerticesCoordinates(&xyz_list[0][0], nodes, numgrids);
    372372        for(i=0;i<numgrids;i++)z_list[i]=xyz_list[i][2];
    373373       

  • issm/trunk/src/c/objects/Icefront.cpp

    r3107 r3159  
    315315        /* Get dof list and node coordinates: */
    316316        GetDofList(&doflist[0],&numberofdofspernode);
    317         GetElementNodeData( &xyz_list[0][0], nodes, numgrids);
     317        GetVerticesCoordinates(&xyz_list[0][0], nodes, numgrids);
    318318       
    319319        /*Now build thickness_list_element and bed_list_element: */
    … …  
    423423        /* Get dof list and node coordinates: */
    424424        GetDofList(&doflist[0],&numberofdofspernode);
    425         GetElementNodeData( &xyz_list[0][0], nodes, numgrids);
     425        GetVerticesCoordinates(&xyz_list[0][0], nodes, numgrids);
    426426       
    427427        //Identify which grids are comprised in the quad:
    … …  
    582582        /* Get dof list and node coordinates: */
    583583        GetDofList(&doflist[0],&numberofdofspernode);
    584         GetElementNodeData( &xyz_list[0][0], nodes, numgrids);
     584        GetVerticesCoordinates(&xyz_list[0][0], nodes, numgrids);
    585585       
    586586        //Identify which grids are comprised in the quad:

  • issm/trunk/src/c/objects/Penta.cpp

    r3131 r3159  
    229229
    230230        /*Get node data: */
    231         GetElementNodeData( &xyz_list[0][0], nodes, numgrids);
     231        GetVerticesCoordinates(&xyz_list[0][0], nodes, numgrids);
    232232
    233233        /*pressure is lithostatic: */
    … …  
    510510
    511511                /* Get node coordinates and dof list: */
    512                 GetElementNodeData( &xyz_list[0][0], nodes, numgrids);
     512                GetVerticesCoordinates(&xyz_list[0][0], nodes, numgrids);
    513513                GetDofList(&doflist[0],&numberofdofspernode);
    514514
    … …  
    743743
    744744        /* Get node coordinates and dof list: */
    745         GetElementNodeData( &xyz_list[0][0], nodes, numgrids);
     745        GetVerticesCoordinates(&xyz_list[0][0], nodes, numgrids);
    746746        GetDofList(&doflist[0],&numberofdofspernode);
    747747
    … …  
    995995
    996996        /* Get node coordinates and dof list: */
    997         GetElementNodeData( &xyz_list[0][0], nodes, numgrids);
     997        GetVerticesCoordinates(&xyz_list[0][0], nodes, numgrids);
    998998        GetDofList(&doflist[0],&numberofdofspernode);
    999999
    … …  
    12281228
    12291229        /* Get node coordinates and dof list: */
    1230         GetElementNodeData( &xyz_list[0][0], nodes, numgrids);
     1230        GetVerticesCoordinates(&xyz_list[0][0], nodes, numgrids);
    12311231        GetDofList(&doflist[0],&numberofdofspernode);
    12321232
    … …  
    15761576
    15771577                /* Get node coordinates and dof list: */
    1578                 GetElementNodeData( &xyz_list[0][0], nodes, numgrids);
     1578                GetVerticesCoordinates(&xyz_list[0][0], nodes, numgrids);
    15791579                GetDofList(&doflist[0],&numberofdofspernode);
    15801580
    … …  
    17431743
    17441744        /* Get node coordinates and dof list: */
    1745         GetElementNodeData( &xyz_list[0][0], nodes, numgrids);
     1745        GetVerticesCoordinates(&xyz_list[0][0], nodes, numgrids);
    17461746        GetDofList(&doflist[0],&numberofdofspernode);
    17471747
    … …  
    19621962        /*Now, handle the standard penta element: */
    19631963        /* Get node coordinates and dof list: */
    1964         GetElementNodeData( &xyz_list[0][0], nodes, numgrids);
     1964        GetVerticesCoordinates(&xyz_list[0][0], nodes, numgrids);
    19651965        GetDofList(&doflist[0],&numberofdofspernode);
    19661966
    … …  
    21722172
    21732173        /* Get node coordinates and dof list: */
    2174         GetElementNodeData( &xyz_list[0][0], nodes, numgrids);
     2174        GetVerticesCoordinates(&xyz_list[0][0], nodes, numgrids);
    21752175        GetDofList(&doflist[0],&numberofdofspernode);
    21762176
    … …  
    30493049        double z1,z2,z3,z4,z5,z6;
    30503050
    3051         double sqrt3=sqrt(3.0);
     3051        double sqrt3=SQRT3;
    30523052
    30533053        /*Figure out xi,eta and zi (parametric coordinates), for this gaussian point: */
    … …  
    34613461        const int numgrids=6;
    34623462        double A1,A2,A3,z;
    3463         double sqrt3=sqrt(3.0);
    3464 
    3465         A1=gauss_coord[0]; //first area coordinate value. In term of xi and eta: A1=(1-xi)/2-eta/(2*sqrt(3));
    3466         A2=gauss_coord[1]; //second area coordinate value In term of xi and eta: A2=(1+xi)/2-eta/(2*sqrt(3));
    3467         A3=gauss_coord[2]; //third area coordinate value  In term of xi and eta: A3=y/sqrt(3);
     3463        double sqrt3=SQRT3;
     3464
     3465        A1=gauss_coord[0]; //first area coordinate value. In term of xi and eta: A1=(1-xi)/2-eta/(2*SQRT3);
     3466        A2=gauss_coord[1]; //second area coordinate value In term of xi and eta: A2=(1+xi)/2-eta/(2*SQRT3);
     3467        A3=gauss_coord[2]; //third area coordinate value  In term of xi and eta: A3=y/SQRT3;
    34683468        z=gauss_coord[3]; //fourth vertical coordinate value. Corresponding nodal function: (1-z)/2 and (1+z)/2
    34693469
    … …  
    35083508         * natural coordinate system) at the gaussian point. */
    35093509
    3510         double sqrt3=sqrt(3.0);
     3510        double sqrt3=SQRT3;
    35113511        int    numgrids=7; //six plus bubble grids
    35123512

  • issm/trunk/src/c/objects/Tria.cpp

    r3102 r3159  
    382382
    383383        /* Get node coordinates and dof list: */
    384         GetElementNodeData( &xyz_list[0][0], nodes, numgrids);
     384        GetVerticesCoordinates(&xyz_list[0][0], nodes, numgrids);
    385385        GetDofList(&doflist[0],&numberofdofspernode);
    386386
    … …  
    586586
    587587        /* Get node coordinates and dof list: */
    588         GetElementNodeData( &xyz_list[0][0], nodes, numgrids);
     588        GetVerticesCoordinates(&xyz_list[0][0], nodes, numgrids);
    589589        GetDofList(&doflist[0],&numberofdofspernode);
    590590
    … …  
    795795
    796796        /* Get node coordinates and dof list: */
    797         GetElementNodeData( &xyz_list[0][0], nodes, numgrids);
     797        GetVerticesCoordinates(&xyz_list[0][0], nodes, numgrids);
    798798        GetDofList(&doflist[0],&numberofdofspernode);
    799799
    … …  
    10051005       
    10061006        /* Get node coordinates and dof list: */
    1007         GetElementNodeData( &xyz_list[0][0], nodes, numgrids);
     1007        GetVerticesCoordinates(&xyz_list[0][0], nodes, numgrids);
    10081008        GetDofList(&doflist[0],&numberofdofspernode);
    10091009
    … …  
    11671167
    11681168        /* Get node coordinates and dof list: */
    1169         GetElementNodeData( &xyz_list[0][0], nodes, numgrids);
     1169        GetVerticesCoordinates(&xyz_list[0][0], nodes, numgrids);
    11701170        GetDofList(&doflist[0],&numberofdofspernode);
    11711171
    … …  
    12871287
    12881288        /* Get node coordinates and dof list: */
    1289         GetElementNodeData( &xyz_list[0][0], nodes, numgrids);
     1289        GetVerticesCoordinates(&xyz_list[0][0], nodes, numgrids);
    12901290        GetDofList(&doflist[0],&numberofdofspernode);
    12911291
    … …  
    13351335#define __FUNCT__ "Tria::CreateKMatrixPrognostic"
    13361336void  Tria::CreateKMatrixPrognostic(Mat Kgg,void* vinputs,int analysis_type,int sub_analysis_type){
    1337 
    13381337
    13391338        /* local declarations */
    … …  
    13641363        double DLprime[2][2]={0.0};
    13651364        double DL_scalar;
    1366         double Ke_gg[numdof][numdof]={0.0};//local element stiffness matrix
    1367         double Ke_gg_gaussian[numdof][numdof]={0.0}; //stiffness matrix evaluated at the gaussian point.
    1368         double Ke_gg_thickness1[numdof][numdof]={0.0}; //stiffness matrix evaluated at the gaussian point.
    1369         double Ke_gg_thickness2[numdof][numdof]={0.0}; //stiffness matrix evaluated at the gaussian point.
    1370 
     1365        double Ke_gg[numdof][numdof]={0.0};
     1366        double Ke_gg_gaussian[numdof][numdof]={0.0};
     1367        double Ke_gg_thickness1[numdof][numdof]={0.0};
     1368        double Ke_gg_thickness2[numdof][numdof]={0.0};
    13711369        double Jdettria;
    13721370
    … …  
    13841382        double  dt;
    13851383        int     dofs[2]={0,1};
    1386         int     found=0;
     1384        int     found;
    13871385
    13881386        ParameterInputs* inputs=NULL;
    … …  
    14041402
    14051403        /* Get node coordinates and dof list: */
    1406         GetElementNodeData( &xyz_list[0][0], nodes, numgrids);
     1404        GetVerticesCoordinates(&xyz_list[0][0], nodes, numgrids);
    14071405        GetDofList(&doflist[0],&numberofdofspernode);
    14081406
    14091407        //Create Artificial diffusivity once for all if requested
    1410         if(numpar->artdiff){
     1408        if(numpar->artdiff==1){
    14111409                //Get the Jacobian determinant
    14121410                gauss_l1l2l3[0]=1.0/3.0; gauss_l1l2l3[1]=1.0/3.0; gauss_l1l2l3[2]=1.0/3.0;
    … …  
    14261424        /* Start  looping on the number of gaussian points: */
    14271425        for (ig=0; ig<num_gauss; ig++){
     1426
    14281427                /*Pick up the gaussian point: */
    14291428                gauss_weight=*(gauss_weights+ig);
    … …  
    14871486                for( i=0; i<numdof; i++) for(j=0;j<numdof;j++) Ke_gg[i][j]+=Ke_gg_thickness2[i][j];
    14881487
    1489                 if(numpar->artdiff){
     1488                if(numpar->artdiff==1){
    14901489
    14911490                        /* Compute artificial diffusivity */
    … …  
    15881587       
    15891588        /* Get node coordinates and dof list: */
    1590         GetElementNodeData( &xyz_list[0][0], nodes, numgrids);
     1589        GetVerticesCoordinates(&xyz_list[0][0], nodes, numgrids);
    15911590        GetDofList(&doflist[0],&numberofdofspernode);
    15921591
    … …  
    16821681
    16831682        /* Get node coordinates and dof list: */
    1684         GetElementNodeData( &xyz_list[0][0], nodes, numgrids);
     1683        GetVerticesCoordinates(&xyz_list[0][0], nodes, numgrids);
    16851684        GetDofList(&doflist[0],&numberofdofspernode);
    16861685
    … …  
    18301829
    18311830        /* Get node coordinates and dof list: */
    1832         GetElementNodeData( &xyz_list[0][0], nodes, numgrids);
     1831        GetVerticesCoordinates(&xyz_list[0][0], nodes, numgrids);
    18331832        GetDofList(&doflist[0],&numberofdofspernode);
    18341833
    … …  
    19211920
    19221921        /* Get node coordinates and dof list: */
    1923         GetElementNodeData( &xyz_list[0][0], nodes, numgrids);
     1922        GetVerticesCoordinates(&xyz_list[0][0], nodes, numgrids);
    19241923        GetDofList(&doflist[0],&numberofdofspernode);
    19251924
    … …  
    20192018
    20202019        /* Get node coordinates and dof list: */
    2021         GetElementNodeData( &xyz_list[0][0], nodes, numgrids);
     2020        GetVerticesCoordinates(&xyz_list[0][0], nodes, numgrids);
    20222021        GetDofList(&doflist[0],&numberofdofspernode);
    20232022
    … …  
    21292128
    21302129        /* Get node coordinates and dof list: */
    2131         GetElementNodeData( &xyz_list[0][0], nodes, numgrids);
     2130        GetVerticesCoordinates(&xyz_list[0][0], nodes, numgrids);
    21322131        GetDofList(&doflist[0],&numberofdofspernode);
    21332132
    … …  
    23072306
    23082307        /* Get node coordinates and dof list: */
    2309         GetElementNodeData( &xyz_list[0][0], nodes, numgrids);
     2308        GetVerticesCoordinates(&xyz_list[0][0], nodes, numgrids);
    23102309        GetDofList(&doflist[0],&numberofdofspernode);
    23112310
    … …  
    23862385
    23872386        /* Get node coordinates and dof list: */
    2388         GetElementNodeData( &xyz_list[0][0], nodes, numgrids);
     2387        GetVerticesCoordinates(&xyz_list[0][0], nodes, numgrids);
    23892388        GetDofList(&doflist[0],&numberofdofspernode);
    23902389
    … …  
    24942493       
    24952494        /* Get node coordinates and dof list: */
    2496         GetElementNodeData( &xyz_list[0][0], nodes, numgrids);
     2495        GetVerticesCoordinates(&xyz_list[0][0], nodes, numgrids);
    24972496        GetDofList(&doflist[0],&numberofdofspernode);
    24982497
    … …  
    26092608       
    26102609        /* Get node coordinates and dof list: */
    2611         GetElementNodeData( &xyz_list[0][0], nodes, numgrids);
     2610        GetVerticesCoordinates(&xyz_list[0][0], nodes, numgrids);
    26122611        GetDofList(&doflist[0],&numberofdofspernode);
    26132612
    … …  
    28032802
    28042803        /* Get node coordinates and dof list: */
    2805         GetElementNodeData( &xyz_list[0][0], nodes, numgrids);
     2804        GetVerticesCoordinates(&xyz_list[0][0], nodes, numgrids);
    28062805        GetDofList(&doflist[0],&numberofdofspernode);
    28072806
    … …  
    30053004
    30063005        /*Get xyz list: */
    3007         GetElementNodeData( &xyz_list[0][0], nodes, numgrids);
     3006        GetVerticesCoordinates(&xyz_list[0][0], nodes, numgrids);
    30083007        x1=xyz_list[0][0]; y1=xyz_list[0][1];
    30093008        x2=xyz_list[1][0]; y2=xyz_list[1][1];
    … …  
    30273026
    30283027        /*Get xyz list: */
    3029         GetElementNodeData( &xyz_list[0][0], nodes, numgrids);
     3028        GetVerticesCoordinates(&xyz_list[0][0], nodes, numgrids);
    30303029        x1=xyz_list[0][0]; y1=xyz_list[0][1];
    30313030        x2=xyz_list[1][0]; y2=xyz_list[1][1];
    … …  
    32543253        const int numgrids=3;
    32553254        double x1,y1,x2,y2,x3,y3;
    3256         double sqrt3=sqrt(3.0);
    32573255       
    32583256        x1=*(xyz_list+numgrids*0+0);
    … …  
    32653263
    32663264        *(J+NDOF2*0+0)=1.0/2.0*(x2-x1);
    3267         *(J+NDOF2*1+0)=sqrt3/6.0*(2*x3-x1-x2);
     3265        *(J+NDOF2*1+0)=SQRT3/6.0*(2*x3-x1-x2);
    32683266        *(J+NDOF2*0+1)=1.0/2.0*(y2-y1);
    3269         *(J+NDOF2*1+1)=sqrt3/6.0*(2*y3-y1-y2);
     3267        *(J+NDOF2*1+1)=SQRT3/6.0*(2*y3-y1-y2);
    32703268}
    32713269/*}}}*/
    … …  
    32883286
    32893287
    3290         *Jdet=sqrt(3.0)/6.0*((x2-x1)*(y3-y1)-(y2-y1)*(x3-x1));
     3288        *Jdet=SQRT3/6.0*((x2-x1)*(y3-y1)-(y2-y1)*(x3-x1));
    32913289
    32923290
    … …  
    33183316
    33193317
    3320         *Jdet=sqrt(3.0)/6.0*pow(pow(((y2-y1)*(z3-z1)-(z2-z1)*(y3-y1)),2.0)+pow(((z2-z1)*(x3-x1)-(x2-x1)*(z3-z1)),2.0)+pow(((x2-x1)*(y3-y1)-(y2-y1)*(x3-x1)),2.0),0.5);
     3318        *Jdet=SQRT3/6.0*pow(pow(((y2-y1)*(z3-z1)-(z2-z1)*(y3-y1)),2.0)+pow(((z2-z1)*(x3-x1)-(x2-x1)*(z3-z1)),2.0)+pow(((x2-x1)*(y3-y1)-(y2-y1)*(x3-x1)),2.0),0.5);
    33213319
    33223320
    … …  
    34693467        const int numgrids=3;
    34703468
    3471         double sqrt3=sqrt(3.0);
    3472 
    34733469        /*First nodal function: */
    34743470        *(dl1dl3+numgrids*0+0)=-1.0/2.0;
    3475         *(dl1dl3+numgrids*1+0)=-1.0/(2.0*sqrt3);
     3471        *(dl1dl3+numgrids*1+0)=-1.0/(2.0*SQRT3);
    34763472
    34773473        /*Second nodal function: */
    34783474        *(dl1dl3+numgrids*0+1)=1.0/2.0;
    3479         *(dl1dl3+numgrids*1+1)=-1.0/(2.0*sqrt3);
     3475        *(dl1dl3+numgrids*1+1)=-1.0/(2.0*SQRT3);
    34803476
    34813477        /*Third nodal function: */
    34823478        *(dl1dl3+numgrids*0+2)=0;
    3483         *(dl1dl3+numgrids*1+2)=1.0/sqrt3;
     3479        *(dl1dl3+numgrids*1+2)=1.0/SQRT3;
    34843480
    34853481}
    … …  
    36663662
    36673663        /* Get node coordinates and dof list: */
    3668         GetElementNodeData( &xyz_list[0][0], nodes, numgrids);
     3664        GetVerticesCoordinates(&xyz_list[0][0], nodes, numgrids);
    36693665        GetDofList1(&doflist1[0]);
    36703666
    … …  
    38433839
    38443840        /* Get node coordinates and dof list: */
    3845         GetElementNodeData( &xyz_list[0][0], nodes, numgrids);
     3841        GetVerticesCoordinates(&xyz_list[0][0], nodes, numgrids);
    38463842        GetDofList1(&doflist1[0]);
    38473843
    … …  
    40654061
    40664062        /* Get node coordinates and dof list: */
    4067         GetElementNodeData( &xyz_list[0][0], nodes, numgrids);
     4063        GetVerticesCoordinates(&xyz_list[0][0], nodes, numgrids);
    40684064        GetDofList1(&doflist1[0]);
    40694065
    … …  
    42634259
    42644260        /*Get xyz list: */
    4265         GetElementNodeData( &xyz_list[0][0], nodes, numgrids);
     4261        GetVerticesCoordinates(&xyz_list[0][0], nodes, numgrids);
    42664262
    42674263        /*recover velocity at three element nodes: */
    … …  
    43834379
    43844380        /* Get node coordinates and dof list: */
    4385         GetElementNodeData( &xyz_list[0][0], nodes, numgrids);
     4381        GetVerticesCoordinates(&xyz_list[0][0], nodes, numgrids);
    43864382
    43874383        /* Recover input data: */

  • issm/trunk/src/c/shared/Elements/ResolvePointers.cpp

    r358 r3159  
    2525        Object* object=NULL;
    2626        int i;
    27        
    2827
    2928        for(i=0;i<num_objects;i++){
    … …  
    5756                }
    5857        }
    59        
    6058}
    61 

  • issm/trunk/src/c/shared/Elements/elements.h

    r1 r3159  
    1111void ResolvePointers(Object** objects,int* object_ids,int* object_offsets,int num_objects,DataSet* dataset);
    1212double Paterson(double temperature);
    13 int GetElementNodeData( double* xyz,  Node** nodes, int numgrids);
     13int GetVerticesCoordinates(double* xyz,  Node** nodes, int numgrids);
    1414
    1515#endif //ifndef _SHARED_ELEMENTS_H_

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

    r1439 r3159  
    66#include "../Alloc/alloc.h"
    77#include "../../include/macros.h"
     8#include "../../include/typedefs.h"
    89#include "../Exceptions/exceptions.h"
    910#include <math.h>
    1011#include <float.h>
    11 
    1212
    1313/*=======1=========2=========3=========4=========5=========6=========7=========8
    … …  
    16051605                        for (i=0; i<nigaus; i++) {
    16061606                                xi        =      1./2.*(1.-egaus[j])*xgaus[i];
    1607                                 eta       =sqrt(3.)/2.*(1.+egaus[j]);
    1608                                 (*pwgt)[ipt]=xwgt[i]*ewgt[j]*(sqrt(3.)/4.*(1.-egaus[j]));
    1609 
    1610                                 (*pl1 )[ipt]=(1.-xi-eta/sqrt(3.))/2.;
    1611                                 (*pl2 )[ipt]=(1.+xi-eta/sqrt(3.))/2.;
    1612                                 (*pl3 )[ipt]=       eta/sqrt(3.);
     1607                                eta       =SQRT3/2.*(1.+egaus[j]);
     1608                                (*pwgt)[ipt]=xwgt[i]*ewgt[j]*(SQRT3/4.*(1.-egaus[j]));
     1609
     1610                                (*pl1 )[ipt]=(1.-xi-eta/SQRT3)/2.;
     1611                                (*pl2 )[ipt]=(1.+xi-eta/SQRT3)/2.;
     1612                                (*pl3 )[ipt]=       eta/SQRT3;
    16131613
    16141614                                ipt++;
    … …  
    19551955                        (*pl3 )[i]=l3p [iord-1][i];
    19561956                        (*pl4 )[i]=l4p [iord-1][i];
    1957                         (*pwgt)[i]=wgtp[iord-1][i]*2./3.*sqrt(2.);
     1957                        (*pwgt)[i]=wgtp[iord-1][i]*2./3.*SQRT2;
    19581958                }
    19591959        }
    … …  
    19891989                                for (i=0; i<nigaus; i++) {
    19901990                                        xi        =1./4.*(1.-egaus[j])*(1.-zgaus[k])*xgaus[i];
    1991                                         eta       =1./4./sqrt(3.)
     1991                                        eta       =1./4./SQRT3
    19921992                                                           *(5.+3.*egaus[j]-zgaus[k]-3.*egaus[j]*zgaus[k]);
    19931993                                        zeta      =sqrt(2./3.)*(1.+zgaus[k]);
    19941994                                        (*pwgt)[ipt]=xwgt[i]*ewgt[j]*zwgt[k]
    1995                                                                  *(sqrt(2.)/16.
     1995                                                                 *(SQRT2/16.
    19961996                                                                   *(1.-egaus[j])*pow(1.-zgaus[k],2.));
    19971997
    1998                                         (*pl1 )[ipt]=(1.-xi-eta/sqrt(3.)-zeta/sqrt(6.))/2.;
    1999                                         (*pl2 )[ipt]=(1.+xi-eta/sqrt(3.)-zeta/sqrt(6.))/2.;
    2000                                         (*pl3 )[ipt]=(      eta         -zeta/sqrt(8.))/sqrt(3.);
    2001                                         (*pl4 )[ipt]=(                   zeta/sqrt(8.))*sqrt(3.);
     1998                                        (*pl1 )[ipt]=(1.-xi-eta/SQRT3-zeta/sqrt(6.))/2.;
     1999                                        (*pl2 )[ipt]=(1.+xi-eta/SQRT3-zeta/sqrt(6.))/2.;
     2000                                        (*pl3 )[ipt]=(      eta         -zeta/sqrt(8.))/SQRT3;
     2001                                        (*pl4 )[ipt]=(                   zeta/sqrt(8.))*SQRT3;
    20022002
    20032003                                        ipt++;
Note: See TracChangeset for help on using the changeset viewer.