Changeset 16841

Timestamp:

11/20/13 20:33:57 (11 years ago)

Author:

Mathieu Morlighem

Message:

BUG: fixing KMatrix Thermal

File:

• issm/trunk-jpl/src/c/analyses/ThermalAnalysis.cpp (modified) (1 diff)

issm/trunk-jpl/src/c/analyses/ThermalAnalysis.cpp

@@ -170 +170 @@
 
                 element->JacobianDeterminant(&Jdet,xyz_list,gauss);
-                D_scalar=gauss->weight*Jdet*kappa;
+                D_scalar=gauss->weight*Jdet;
                 if(dt!=0.) D_scalar=D_scalar*dt;
 
                 /*Conduction: */
                 GetBConduct(B,element,xyz_list,gauss); 
-                D[0][0]=D_scalar;
-                D[1][1]=D_scalar;
-                D[2][2]=D_scalar;
+                D[0][0]=D_scalar*kappa;
+                D[1][1]=D_scalar*kappa;
+                D[2][2]=D_scalar*kappa;
                 TripleMultiply(B,3,numnodes,1,
                                         &D[0][0],3,3,0,