Changeset 25211

Timestamp:

07/06/20 10:57:47 (5 years ago)

Author:

Mathieu Morlighem

Message:

BUG: problems for 1D simulations

File:

• issm/trunk-jpl/src/c/analyses/MasstransportAnalysis.cpp (modified) (2 diffs)

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

@@ -365 +365 @@
                         for(int i=0;i<numnodes;i++){
                                 for(int j=0;j<numnodes;j++){
-                                        Ke->values[i*numnodes+j] += D_scalar*dvxdx*dbasis[0*numnodes+i]*dbasis[0*numnodes+j];
+                                        Ke->values[i*numnodes+j] += D_scalar*dvxdx*basis[i]*basis[j];
                                         Ke->values[i*numnodes+j] += D_scalar*vx*dbasis[0*numnodes+j]*basis[i];
                                 }
@@ -431 +431 @@
                                 for(int i=0;i<numnodes;i++){
                                         for(int j=0;j<numnodes;j++){
-                                                Ke->values[i*numnodes+j] += D_scalar*h/(2.*vel)*dbasis[0*numnodes+i] *D[0]*dbasis[0*numnodes+j];
+                                                Ke->values[i*numnodes+j] += dbasis[0*numnodes+i]*D[0]*dbasis[0*numnodes+j];
                                         }
                                 }