Changeset 4654


Ignore:
Timestamp:
07/19/10 12:27:07 (15 years ago)
Author:
Mathieu Morlighem
Message:

fixed BIG bug in gauss points, to be hard coded by the way

File:
1 edited

Legend:

Unmodified
Added
Removed

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

    r4652 r4654  
    592592        double surface[numgrids];
    593593        double xyz_list[numgrids][3];
    594         double gauss[numgrids][numgrids]={{1,0,0,0},{0,1,0,0},{0,0,1,0},{1,0,0,1},{0,1,0,1},{0,0,1,1}};
     594        double gauss[numgrids][4]={{1,0,0,-1},{0,1,0,-1},{0,0,1,-1},{1,0,0,1},{0,1,0,1},{0,0,1,1}};
    595595
    596596        /*inputs: */
    … …  
    20952095        double       rho_ice,g;
    20962096        double       xyz_list[numvertices][3];
    2097         double       gauss[numvertices][numvertices]={{1,0,0,0},{0,1,0,0},{0,0,1,0},{1,0,0,1},{0,1,0,1},{0,0,1,1}};
     2097        double       gauss[numvertices][4]={{1,0,0,-1},{0,1,0,-1},{0,0,1,-1},{1,0,0,1},{0,1,0,1},{0,0,1,1}};
    20982098       
    20992099        Input*       VzInput=NULL;
    … …  
    21762176        double       rho_ice,g;
    21772177        double       xyz_list[numvertices][3];
    2178         double       gauss[numvertices][numvertices]={{1,0,0,0},{0,1,0,0},{0,0,1,0},{1,0,0,1},{0,1,0,1},{0,0,1,1}};
     2178        double       gauss[numvertices][4]={{1,0,0,-1},{0,1,0,-1},{0,0,1,-1},{1,0,0,1},{0,1,0,1},{0,0,1,1}};
    21792179       
    21802180        Input*       VzInput=NULL;
    … …  
    22572257        double       rho_ice,g;
    22582258        double       xyz_list[numvertices][3];
    2259         double       gauss[numvertices][numvertices]={{1,0,0,0},{0,1,0,0},{0,0,1,0},{1,0,0,1},{0,1,0,1},{0,0,1,1}};
     2259        double       gauss[numvertices][4]={{1,0,0,-1},{0,1,0,-1},{0,0,1,-1},{1,0,0,1},{0,1,0,1},{0,0,1,1}};
    22602260
    22612261        Input*       VxInput=NULL;
    … …  
    28142814                values[i*numdofpervertex+2]=vz;
    28152815                values[i*numdofpervertex+3]=p/stokesreconditioning;
     2816                if (
     2817                                        id==1089 ||
     2818                                        id==2761 ||
     2819                                        id==1121 ||
     2820                                        id==2762
     2821                        ) printf("element %i vertex %i p=%g  pstokes=%g\n",id,i,p,p/stokesreconditioning);
    28162822        }
    28172823
    … …  
    58695875        int       grid=-1;
    58705876        int       i;
    5871         double gauss[numnodes][numnodes]={{1,0,0,0},{0,1,0,0},{0,0,1,0},{1,0,0,1},{0,1,0,1},{0,0,1,1}};
     5877        double gauss[numnodes][4]={{1,0,0,-1},{0,1,0,-1},{0,0,1,-1},{1,0,0,1},{0,1,0,1},{0,0,1,1}};
    58725878
    58735879        /*go through 3 nodes (all nodes for tria) and identify 1st and 2nd nodes: */
Note: See TracChangeset for help on using the changeset viewer.