Changeset 13042
- Timestamp:
- 08/15/12 11:24:43 (13 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/src/c/classes/objects/Elements/Penta.cpp
r13036 r13042 6327 6327 /*Intermediaries */ 6328 6328 int i,j,ig,approximation; 6329 IssmDouble 6330 IssmDouble viscosity, oldviscosity, newviscosity, viscosity_overshoot;6331 IssmDouble epsilon[5],oldepsilon[5];/* epsilon=[exx,eyy,exy,exz,eyz];*/6332 IssmDouble epsilons[6];//6 for stokes6333 IssmDouble 6334 IssmDouble 6335 IssmDouble 6336 IssmDouble D[3][3]={0.0};// material matrix, simple scalar matrix.6337 IssmDouble 6338 IssmDouble Ke_gg_gaussian[numdof2d][numdof2d]; //stiffness matrix evaluated at the gaussian point.6339 Tria * tria=NULL;6340 Penta * pentabase=NULL;6341 GaussPenta *gauss =NULL;6342 GaussTria *gauss_tria =NULL;6329 IssmDouble Jdet; 6330 IssmDouble viscosity , oldviscosity, newviscosity, viscosity_overshoot; 6331 IssmDouble epsilon[5],oldepsilon[5]; /* epsilon=[exx,eyy,exy,exz,eyz];*/ 6332 IssmDouble epsilons[6]; //6 for stokes 6333 IssmDouble xyz_list[NUMVERTICES][3]; 6334 IssmDouble B[3][numdof2d]; 6335 IssmDouble Bprime[3][numdof2d]; 6336 IssmDouble D[3][3]= {0.0}; // material matrix, simple scalar matrix. 6337 IssmDouble D_scalar; 6338 IssmDouble Ke_gg_gaussian[numdof2d][numdof2d]; 6339 Tria *tria = NULL; 6340 Penta *pentabase = NULL; 6341 GaussPenta *gauss = NULL; 6342 GaussTria *gauss_tria = NULL; 6343 6343 6344 6344 /*Find penta on bed as this is a macayeal elements: */
Note:
See TracChangeset
for help on using the changeset viewer.