Changeset 16894


Ignore:
Timestamp:
11/22/13 12:35:33 (11 years ago)
Author:
Mathieu Morlighem
Message:

BUG: fixing enthalpy (fabs not used in stabilization)

Location:
issm/trunk-jpl/src/c/analyses
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • issm/trunk-jpl/src/c/analyses/EnthalpyAnalysis.cpp

    r16888 r16894  
    287287                        vel=sqrt(vx*vx + vy*vy + vz*vz)+1.e-14;
    288288                        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;
    292292                        for(int i=0;i<3;i++) for(int j=0;j<3;j++) K[i][j] = D_scalar*K[i][j];
    293293
    294294                        GetBAdvecprime(Bprime,element,xyz_list,gauss);
    295 
    296295                        TripleMultiply(Bprime,3,numnodes,1,
    297296                                                &K[0][0],3,3,0,
  • issm/trunk-jpl/src/c/analyses/StressbalanceAnalysis.cpp

    r16893 r16894  
    14491449        ElementMatrix* Ke2=CreateKMatrixL1L2Friction(element);
    14501450        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         //}
    14571451
    14581452        /*clean-up and return*/
Note: See TracChangeset for help on using the changeset viewer.