Changeset 4493


Ignore:
Timestamp:
07/09/10 07:38:28 (15 years ago)
Author:
Mathieu Morlighem
Message:

Now SUCCESS in macayeal sheet (TO BE REMOVED AT SOME POINT)

File:
1 edited

Legend:

Unmodified
Added
Removed
  • issm/trunk/src/c/objects/Elements/Tria.cpp

    r4492 r4493  
    29632963void  Tria::CreateKMatrixDiagnosticHorizFriction(Mat Kgg){
    29642964
    2965 
    29662965        /* local declarations */
    2967         int             i,j;
    2968         int analysis_type;
     2966        int       i,j;
     2967        int       analysis_type;
    29692968
    29702969        /* node data: */
    2971         const int    numgrids=3;
    2972         const int    numdof=2*numgrids;
    2973         double       xyz_list[numgrids][3];
    2974         int          doflist[numdof];
    2975         int          numberofdofspernode;
    2976        
     2970        const int numvertices = 3;
     2971        const int numdof   = 2 *numvertices;
     2972        double    xyz_list[numvertices][3];
     2973        int       doflist[numdof];
     2974        int       numberofdofspernode;
     2975
    29772976        /* gaussian points: */
    29782977        int     num_gauss,ig;
    2979         double* first_gauss_area_coord  = NULL;
    2980         double* second_gauss_area_coord = NULL;
    2981         double* third_gauss_area_coord  = NULL;
    2982         double* gauss_weights           = NULL;
     2978        double *first_gauss_area_coord  = NULL;
     2979        double *second_gauss_area_coord = NULL;
     2980        double *third_gauss_area_coord  = NULL;
     2981        double *gauss_weights           = NULL;
    29832982        double  gauss_weight;
    29842983        double  gauss_l1l2l3[3];
     
    30002999
    30013000        /*friction: */
    3002         Friction* friction=NULL;
    3003         double alpha2;
     3001        Friction *friction = NULL;
     3002        double    alpha2;
     3003        double    alpha2_list[numvertices];                                       //TO BE DELETED
     3004        double    gauss[numvertices][numvertices] = {{1,0,0},{0,1,0},{0,0,1}}; //TO BE DELETED
     3005
    30043006
    30053007        double MAXSLOPE=.06; // 6 %
     
    30183020       
    30193021        /* Get node coordinates and dof list: */
    3020         GetVerticesCoordinates(&xyz_list[0][0], nodes, numgrids);
     3022        GetVerticesCoordinates(&xyz_list[0][0], nodes, numvertices);
    30213023        GetDofList(&doflist[0],&numberofdofspernode);
    30223024
     
    30293031        if (drag_type!=2)ISSMERROR(" non-viscous friction not supported yet!");
    30303032        friction=new Friction("2d",inputs,matpar,analysis_type);
     3033
     3034        /*COMPUT alpha2_list (TO BE DELETED)*/
     3035        for(i=0;i<numvertices;i++){
     3036                friction->GetAlpha2(&alpha2_list[i],&gauss[i][0],VxEnum,VyEnum,VzEnum); // TO UNCOMMENT
     3037        }
    30313038
    30323039        /* Get gaussian points and weights (make this a statically initialized list of points? fstd): */
     
    30423049
    30433050                /*Friction: */
    3044                 friction->GetAlpha2(&alpha2, gauss_l1l2l3,VxEnum,VyEnum,VzEnum);
     3051                // friction->GetAlpha2(&alpha2, gauss_l1l2l3,VxEnum,VyEnum,VzEnum); // TO UNCOMMENT
     3052                GetParameterValue(&alpha2,&alpha2_list[0],gauss_l1l2l3); // TO BE DELETED
     3053
    30453054
    30463055                // If we have a slope > 6% for this element,  it means  we are on a mountain. In this particular case,
Note: See TracChangeset for help on using the changeset viewer.