Changeset 4493
- Timestamp:
- 07/09/10 07:38:28 (15 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk/src/c/objects/Elements/Tria.cpp
r4492 r4493 2963 2963 void Tria::CreateKMatrixDiagnosticHorizFriction(Mat Kgg){ 2964 2964 2965 2966 2965 /* local declarations */ 2967 int 2968 int analysis_type;2966 int i,j; 2967 int analysis_type; 2969 2968 2970 2969 /* node data: */ 2971 const int numgrids=3;2972 const int numdof=2*numgrids;2973 double xyz_list[numgrids][3];2974 int 2975 int 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 2977 2976 /* gaussian points: */ 2978 2977 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; 2983 2982 double gauss_weight; 2984 2983 double gauss_l1l2l3[3]; … … 3000 2999 3001 3000 /*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 3004 3006 3005 3007 double MAXSLOPE=.06; // 6 % … … 3018 3020 3019 3021 /* Get node coordinates and dof list: */ 3020 GetVerticesCoordinates(&xyz_list[0][0], nodes, num grids);3022 GetVerticesCoordinates(&xyz_list[0][0], nodes, numvertices); 3021 3023 GetDofList(&doflist[0],&numberofdofspernode); 3022 3024 … … 3029 3031 if (drag_type!=2)ISSMERROR(" non-viscous friction not supported yet!"); 3030 3032 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 } 3031 3038 3032 3039 /* Get gaussian points and weights (make this a statically initialized list of points? fstd): */ … … 3042 3049 3043 3050 /*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 3045 3054 3046 3055 // 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.