Changeset 5663


Ignore:
Timestamp:
09/03/10 09:47:15 (15 years ago)
Author:
Mathieu Morlighem
Message:

more macros

Location:
issm/trunk/src/c/objects/Elements
Files:
2 edited

Legend:

Unmodified
Added
Removed

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

    r5662 r5663  
    2020
    2121/*Element macros*/
    22 #define NUMVERTICES 6
     22#define NUMVERTICES   6
     23#define NUMVERTICES2D 3
     24#define NDOF1 1
     25#define NDOF2 2
     26#define NDOF3 3
     27#define NDOF4 4
    2328
    2429/*Penta constructors and destructor*/
    … …  
    23022307        /*Collapsed formulation: */
    23032308        int       i;
    2304         const int NDOF2=2;
    23052309        const int numdofs=NDOF2*NUMVERTICES;
    23062310        int*      doflist=NULL;
    … …  
    23862390
    23872391        /* node data: */
    2388         const int    NUMVERTICES3d=6;
    2389         const int    NUMVERTICES2d=3;
    2390         const int    numdof2d=2*NUMVERTICES2d;
    2391         double       xyz_list[NUMVERTICES3d][3];
     2392        const int    numdof2d=2*NUMVERTICES2D;
     2393        double       xyz_list[NUMVERTICES][3];
    23922394        int*         doflist=NULL;
    23932395
    … …  
    24972499
    24982500                /* Get node coordinates and dof list: */
    2499                 GetVerticesCoordinates(&xyz_list[0][0], nodes, NUMVERTICES3d);
     2501                GetVerticesCoordinates(&xyz_list[0][0], nodes,NUMVERTICES);
    25002502                tria->GetDofList(&doflist,MacAyealApproximationEnum);  //Pattyn dof list
    25012503
    … …  
    27152717
    27162718        int i,j;
    2717         const int DOFPERGRID=4;
    2718         const int numdof=NUMVERTICES*DOFPERGRID;
     2719        const int numdof=NUMVERTICES*NDOF4;
    27192720        int*      doflist=NULL;
    27202721
    2721         const int NUMVERTICES2d=3;
    2722         const int numdof2d=NUMVERTICES2d*DOFPERGRID;
     2722        const int numdof2d=NUMVERTICES2D*NDOF4;
    27232723
    27242724        /*Collapsed formulation: */
    … …  
    27272727        /*Grid data: */
    27282728        double     xyz_list[NUMVERTICES][3];
    2729         double    xyz_list_tria[NUMVERTICES2d][3];
     2729        double    xyz_list_tria[NUMVERTICES2D][3];
    27302730        double    bed_normal[3];
    27312731
    … …  
    28722872                friction=new Friction("3d",inputs,matpar,analysis_type);
    28732873
    2874                 for(i=0;i<NUMVERTICES2d;i++){
     2874                for(i=0;i<NUMVERTICES2D;i++){
    28752875                        for(j=0;j<3;j++){
    28762876                                xyz_list_tria[i][j]=xyz_list[i][j];
    … …  
    29702970
    29712971        /* node data: */
    2972         const int    NDOF1=1;
    29732972        const int    numdof=NDOF1*NUMVERTICES;
    29742973        double       xyz_list[NUMVERTICES][3];
    … …  
    31483147
    31493148        /* node data: */
    3150         const int    NDOF1=1;
    31513149        const int    numdof=NDOF1*NUMVERTICES;
    31523150        double       xyz_list[NUMVERTICES][3];
    … …  
    34803478        int i,j,k;
    34813479       
    3482         const int NDOF2=2;
    34833480        const int numdofs=NDOF2*NUMVERTICES;
    34843481        int*      doflist=NULL;
    … …  
    36673664
    36683665        /* node data: */
    3669         const int    NDOF2=2;
    36703666        const int    numdof=NDOF2*NUMVERTICES;
    36713667        double       xyz_list[NUMVERTICES][3];
    … …  
    37563752        /*indexing: */
    37573753        int i,j;
    3758 
    3759         const int DOFPERGRID=4;
    3760         const int numdof=NUMVERTICES*DOFPERGRID;
    3761         const int NUMVERTICES2d=3;
    3762         int       numdof2d=NUMVERTICES2d*DOFPERGRID;
     3754        const int numdof=NUMVERTICES*NDOF4;
     3755        int       numdof2d=NUMVERTICES2D*NDOF4;
    37633756        int*      doflist=NULL;
    37643757
    … …  
    37673760
    37683761        /*parameters: */
    3769         double             xyz_list_tria[NUMVERTICES2d][3];
     3762        double             xyz_list_tria[NUMVERTICES2D][3];
    37703763        double         xyz_list[NUMVERTICES][3];
    37713764        double             bed_normal[3];
    … …  
    38893882                        /* Build gaussian vector */
    38903883                        for(i=0;i<NUMVERTICES+1;i++){
    3891                                 Pe_gaussian[i*DOFPERGRID+2]=-rho_ice*gravity*Jdet*gauss_weight*l1l7[i];
     3884                                Pe_gaussian[i*NDOF4+2]=-rho_ice*gravity*Jdet*gauss_weight*l1l7[i];
    38923885                        }
    38933886
    … …  
    39343927        if ( (onbed==1) && (shelf==1)){
    39353928
    3936                 for(i=0;i<NUMVERTICES2d;i++){
     3929                for(i=0;i<NUMVERTICES2D;i++){
    39373930                        for(j=0;j<3;j++){
    39383931                                xyz_list_tria[i][j]=xyz_list[i][j];
    … …  
    39703963                        BedNormal(&bed_normal[0],xyz_list_tria);
    39713964
    3972                         for(i=0;i<NUMVERTICES2d;i++){
     3965                        for(i=0;i<NUMVERTICES2D;i++){
    39733966                                for(j=0;j<3;j++){
    3974                                         Pe_temp[i*DOFPERGRID+j]+=water_pressure*gauss_weight*Jdet2d*L[i]*bed_normal[j];
     3967                                        Pe_temp[i*NDOF4+j]+=water_pressure*gauss_weight*Jdet2d*L[i]*bed_normal[j];
    39753968                                }
    39763969                        }
    … …  
    40254018
    40264019        /* node data: */
    4027         const int    NDOF1=1;
    40284020        const int    numdof=NDOF1*NUMVERTICES;
    40294021        double       xyz_list[NUMVERTICES][3];
    … …  
    41894181        int found=0;
    41904182
    4191         const int  NDOF1=1;
    41924183        const int  numdof=NUMVERTICES*NDOF1;
    41934184        int*       doflist=NULL;
    … …  
    52675258        int i;
    52685259
    5269         const int    NUMVERTICES2d=3;
    52705260        const int    numdofpervertex=2;
    52715261        const int    numdof=numdofpervertex*NUMVERTICES;
    5272         const int    numdof2d=numdofpervertex*NUMVERTICES2d;
     5262        const int    numdof2d=numdofpervertex*NUMVERTICES2D;
    52735263        int*         doflistp=NULL;
    52745264        int*         doflistm=NULL;

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

    r5662 r5663  
    2121/*Element macros*/
    2222#define NUMVERTICES 3
     23#define NDOF1 1
     24#define NDOF2 2
     25#define NDOF3 3
     26#define NDOF4 4
    2327
    2428/*Tria constructors and destructor*/
    … …  
    628632
    629633        /*constants: */
    630         const int    NDOF2=2;
    631634        const int    numdof=NDOF2*NUMVERTICES;
    632635
    … …  
    899902void  Tria::GradjB(Vec gradient){
    900903
    901         /* constants*/
    902         const int    NDOF2=2;
    903 
    904904        /*Intermediaries*/
    905905        int        i,ig;
    … …  
    978978
    979979        /* node data: */
    980         const int    NDOF2=2;
    981980        double       xyz_list[NUMVERTICES][3];
    982981        int          doflist1[NUMVERTICES];
    … …  
    16031602        /* node data: */
    16041603        const int    numdof=2*NUMVERTICES;
    1605         const int    NDOF2=2;
    16061604        double       xyz_list[NUMVERTICES][3];
    16071605
    … …  
    17051703        /* node data: */
    17061704        const int    numdof=2*NUMVERTICES;
    1707         const int    NDOF2=2;
    17081705        double       xyz_list[NUMVERTICES][3];
    17091706
    … …  
    18161813        /* node data: */
    18171814        const int    numdof=2*NUMVERTICES;
    1818         const int    NDOF2=2;
    18191815        double       xyz_list[NUMVERTICES][3];
    18201816
    … …  
    19241920        /* node data: */
    19251921        const int    numdof=2*NUMVERTICES;
    1926         const int    NDOF2=2;
    19271922        double       xyz_list[NUMVERTICES][3];
    19281923
    … …  
    20352030        /* node data: */
    20362031        const int    numdof=2*NUMVERTICES;
    2037         const int    NDOF2=2;
    20382032        double       xyz_list[NUMVERTICES][3];
    20392033
    … …  
    24752469
    24762470        /* node data: */
    2477         const int    NDOF1=1;
    24782471        const int    numdof=NDOF1*NUMVERTICES;
    24792472        double       xyz_list[NUMVERTICES][3];
    … …  
    26302623
    26312624        /* node data: */
    2632         const int    NDOF1=1;
    26332625        const int    numdof=NDOF1*NUMVERTICES;
    26342626        double       xyz_list[NUMVERTICES][3];
    … …  
    27182710
    27192711        /* node data: */
    2720         const int    NDOF1=1;
    27212712        const int    numdof=NDOF1*NUMVERTICES;
    27222713        double       xyz_list[NUMVERTICES][3];
    … …  
    33783369        int    i;
    33793370        int    connectivity;
    3380         const int NDOF2=2;
    33813371        const int numdofs=NUMVERTICES*NDOF2;
    33823372        int*         doflist=NULL;
    … …  
    34143404
    34153405        /* node data: */
    3416         const int    NDOF1=1;
    34173406        const int    numdof=NDOF1*NUMVERTICES;
    34183407        double       xyz_list[NUMVERTICES][3];
    … …  
    35163505        int i,j;
    35173506
    3518         const int  NDOF1=1;
    35193507        const int  numdof=NUMVERTICES*NDOF1;
    35203508        int*       doflist=NULL;
    … …  
    35873575
    35883576        /* node data: */
    3589         const int    NDOF1=1;
    35903577        const int    numdof=NDOF1*NUMVERTICES;
    35913578        double       xyz_list[NUMVERTICES][3];
    … …  
    37563743
    37573744        /* node data: */
    3758         const int    NDOF1=1;
    37593745        const int    numdof=NDOF1*NUMVERTICES;
    37603746        double       xyz_list[NUMVERTICES][3];
    … …  
    38653851        /*constants: */
    38663852        const int    numnodes=3;
    3867         const int    NDOF1=1;
    38683853        const int    numdof=NDOF1*numnodes;
    38693854
    … …  
    39153900       
    39163901        /* node data: */
    3917         const int    NDOF1=1;
    39183902        const int    numdof=NDOF1*NUMVERTICES;
    39193903        double       xyz_list[NUMVERTICES][3];
    … …  
    39993983
    40003984        /* node data: */
    4001         const int    NDOF1=1;
    40023985        const int    numdof=NDOF1*NUMVERTICES;
    40033986        double       xyz_list[NUMVERTICES][3];
    … …  
    40704053
    40714054        /* node data: */
    4072         const int    NDOF1=1;
    40734055        const int    numdof=NDOF1*NUMVERTICES;
    40744056        double       xyz_list[NUMVERTICES][3];
    … …  
    41414123
    41424124        /* node data: */
    4143         const int    NDOF1=1;
    41444125        const int    numdof=NDOF1*NUMVERTICES;
    41454126        double       xyz_list[NUMVERTICES][3];
    … …  
    42074188
    42084189        /* node data: */
    4209         const int    NDOF1=1;
    42104190        const int    numdof=NDOF1*NUMVERTICES;
    42114191        double       xyz_list[NUMVERTICES][3];
    … …  
    43014281        /* node data: */
    43024282        const int    numdof=2*NUMVERTICES;
    4303         const int    NDOF2=2;
    43044283        double       xyz_list[NUMVERTICES][3];
    43054284        int*         doflist=NULL;
    … …  
    44994478        /* node data: */
    45004479        const int    numdof=2*NUMVERTICES;
    4501         const int    NDOF2=2;
    45024480        double       xyz_list[NUMVERTICES][3];
    45034481        int*         doflist=NULL;
    … …  
    47204698
    47214699        /* node data: */
    4722         const int    NDOF4=4;
    47234700        const int    numdof=NDOF4*NUMVERTICES;
    47244701        double       xyz_list[NUMVERTICES][3];
    … …  
    49364913        /*Collapsed formulation: */
    49374914        int       i;
    4938         const int NDOF2=2;
    49394915        const int numdofs=NDOF2*NUMVERTICES;
    49404916        int*         doflist=NULL;
    … …  
    50124988
    50134989        /* node data: */
    5014         const int    NDOF1=1;
    50154990        const int    numdof=NDOF1*NUMVERTICES;
    50164991        double       xyz_list[NUMVERTICES][3];
    … …  
    50885063
    50895064        /* node data: */
    5090         const int    NDOF1=1;
    50915065        const int    numdof=NDOF1*NUMVERTICES;
    50925066        double       xyz_list[NUMVERTICES][3];
    … …  
    51605134
    51615135        /* node data: */
    5162         const int    NDOF1=1;
    51635136        const int    numdof=NDOF1*NUMVERTICES;
    51645137        double       xyz_list[NUMVERTICES][3];
    … …  
    52395212        int i,found;
    52405213       
    5241         const int  NDOF1=1;
    52425214        const int  numdof=NUMVERTICES*NDOF1;
    52435215        int*         doflist=NULL;
    … …  
    53325304        int i,found;
    53335305       
    5334         const int  NDOF1=1;
    53355306        const int  numdof=NUMVERTICES*NDOF1;
    53365307        int*       doflist=NULL;
    … …  
    58805851
    58815852        /* node data: */
    5882         const int    NDOF2=2;
    58835853        double       xyz_list[NUMVERTICES][3];
    58845854        int          doflist1[NUMVERTICES];
Note: See TracChangeset for help on using the changeset viewer.