Changeset 17538

Timestamp:

03/25/14 11:47:44 (11 years ago)

Author:

Mathieu Morlighem

Message:

CHG: done with create P Vector

File:

• issm/trunk-jpl/src/c/analyses/StressbalanceAnalysis.cpp (modified) (3 diffs)

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

@@ -3424 +3424 @@
         IssmDouble*    vdbasis = xNew<IssmDouble>(dim*vnumnodes);
         IssmDouble*    tbasis  = xNew<IssmDouble>(tnumnodes);
+        IssmDouble*    D       = xNewZeroInit<IssmDouble>(tausize*tnumnodes*tausize*tnumnodes);
 
         /*Retrieve all inputs and parameters*/
@@ -3530 +3531 @@
                         }
                 }
+
+                /*contribution -Dstar tau*/
+                for(i=0;i<tausize*tnumnodes;i++) D[i*(tausize*tnumnodes)+i] = -gauss->weight*Jdet;
+                TripleMultiply(Dstar,dim*tnumnodes,tausize*tnumnodes,0,
+                                        D,tausize*tnumnodes,tausize*tnumnodes,0,
+                                        tau,tausize*tnumnodes,1,0,
+                                        &pe->values[0],1);
+
+                /*contribution + r Dstar d*/
+                for(i=0;i<tausize*tnumnodes;i++) D[i*(tausize*tnumnodes)+i] = +r*gauss->weight*Jdet;
+                TripleMultiply(Dstar,dim*tnumnodes,tausize*tnumnodes,0,
+                                        D,tausize*tnumnodes,tausize*tnumnodes,0,
+                                        d,tausize*tnumnodes,1,0,
+                                        &pe->values[0],1);
         }
 
@@ -3541 +3556 @@
         xDelete<IssmDouble>(Dstar);
         xDelete<IssmDouble>(d);
+        xDelete<IssmDouble>(D);
         xDelete<IssmDouble>(tau);
         return pe;