Changeset 16894
- Timestamp:
- 11/22/13 12:35:33 (11 years ago)
- Location:
- issm/trunk-jpl/src/c/analyses
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/src/c/analyses/EnthalpyAnalysis.cpp
r16888 r16894 287 287 vel=sqrt(vx*vx + vy*vy + vz*vz)+1.e-14; 288 288 h=sqrt( pow(hx*vx/vel,2) + pow(hy*vy/vel,2) + pow(hz*vz/vel,2)); 289 K[0][0]=h/(2.*vel)* fabs(vx*vx); K[0][1]=h/(2.*vel)*fabs(vx*vy); K[0][2]=h/(2.*vel)*fabs(vx*vz);290 K[1][0]=h/(2.*vel)* fabs(vy*vx); K[1][1]=h/(2.*vel)*fabs(vy*vy); K[1][2]=h/(2.*vel)*fabs(vy*vz);291 K[2][0]=h/(2.*vel)* fabs(vz*vx); K[2][1]=h/(2.*vel)*fabs(vz*vy); K[2][2]=h/(2.*vel)*fabs(vz*vz);289 K[0][0]=h/(2.*vel)*vx*vx; K[0][1]=h/(2.*vel)*vx*vy; K[0][2]=h/(2.*vel)*vx*vz; 290 K[1][0]=h/(2.*vel)*vy*vx; K[1][1]=h/(2.*vel)*vy*vy; K[1][2]=h/(2.*vel)*vy*vz; 291 K[2][0]=h/(2.*vel)*vz*vx; K[2][1]=h/(2.*vel)*vz*vy; K[2][2]=h/(2.*vel)*vz*vz; 292 292 for(int i=0;i<3;i++) for(int j=0;j<3;j++) K[i][j] = D_scalar*K[i][j]; 293 293 294 294 GetBAdvecprime(Bprime,element,xyz_list,gauss); 295 296 295 TripleMultiply(Bprime,3,numnodes,1, 297 296 &K[0][0],3,3,0, -
issm/trunk-jpl/src/c/analyses/StressbalanceAnalysis.cpp
r16893 r16894 1449 1449 ElementMatrix* Ke2=CreateKMatrixL1L2Friction(element); 1450 1450 ElementMatrix* Ke =new ElementMatrix(Ke1,Ke2); 1451 //if(element->Id()==62){1452 //Ke->Echo();1453 //Ke=element->CreateKMatrix();1454 //Ke->Echo();1455 //_error_("S");1456 //}1457 1451 1458 1452 /*clean-up and return*/
Note:
See TracChangeset
for help on using the changeset viewer.