Changeset 16842

Timestamp:

11/20/13 20:42:50 (11 years ago)

Author:

Mathieu Morlighem

Message:

BUG: fixing transient term in KMatrix

File:

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

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

@@ -199 +199 @@
                 /*Transient: */
                 if(dt!=0.){
+                        D_scalar=gauss->weight*Jdet;
                         element->NodalFunctions(basis,gauss);
-
                         TripleMultiply(basis,numnodes,1,0,
                                                 &D_scalar,1,1,0,
                                                 basis,1,numnodes,0,
                                                 &Ke->values[0],1);
+                        D_scalar=D_scalar*dt;
                 }