OilfieldToolsNet/IFEM
4
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

changed: unify (EqualOrder|Compatible)Operators::Residual::Laplacian

Browse Source 
use positive sign in both
This commit is contained in:
Arne Morten Kvarving 2016-11-18 14:36:06 +01:00 committed by timovanopstal
parent 691565d575
commit 92b01e54a4
1 changed files with 3 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

6
Apps/Common/EqualOrderOperators.h
View File

@ -171,8 +171,8 @@ public:
diff += dUdX(k,l)*fe.grad(basis)(i,l); diff += dUdX(k,l)*fe.grad(basis)(i,l);
diff *= scale*fe.detJxW; diff *= scale*fe.detJxW;
// Add negative residual to rhs of momentum equation // Add residual to rhs of momentum equation
EV((i-1)*cmp + k) -= diff; EV((i-1)*cmp + k) += diff;
} }
} }
@ -182,7 +182,7 @@ public:
for (size_t k = 1;k <= cmp;k++) for (size_t k = 1;k <= cmp;k++)
for (size_t l = 1;l <= cmp;l++) for (size_t l = 1;l <= cmp;l++)
// Diffusion // Diffusion
EV((i-1)*cmp + k) -= scale*dUdX(l,k)*fe.grad(basis)(i,l)*fe.detJxW; EV((i-1)*cmp + k) += scale*dUdX(l,k)*fe.grad(basis)(i,l)*fe.detJxW;
} }
} }
}; };