Changeset 24430

Timestamp:

12/04/19 16:29:43 (6 years ago)

Author:

Mathieu Morlighem

Message:

CHG: fixing levelset

File:

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

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

@@ -523 +523 @@
                                         for(int i=0;i<numnodes;i++){
                                                 for(int j=0;j<numnodes;j++){
-                                                        Ke->values[i*numnodes+j] += D_scalar*h/(2.*vel)*w[i]*(w[j]*dbasis[0*numnodes+j]*dbasis[0*numnodes+i] + w[j]*dbasis[1*numnodes+j]*dbasis[1*numnodes+i]);
+                                                        Ke->values[i*numnodes+j] += D_scalar*h/(2.*vel)*(
+                                                                                dbasis[0*numnodes+i] *(w[0]*w[0]*dbasis[0*numnodes+j] + w[0]*w[1]*dbasis[1*numnodes+j]) +
+                                                                                dbasis[1*numnodes+i] *(w[1]*w[0]*dbasis[0*numnodes+j] + w[1]*w[1]*dbasis[1*numnodes+j]) 
+                                                                                );
                                                 }
                                         }