Mercurial > repos > public > sbplib

diff +scheme/Laplace1d.m @ 1083:a22b23098021 feature/laplace1dUpdate

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Divide interface penalty strength by 4, to make it sharp. Stability verified by eigenvalue computations.
author Martin Almquist <malmquist@stanford.edu>
date Tue, 05 Mar 2019 12:41:50 -0800
parents a20fc67e9ac0
children 33c378e508d2
line wrap: on
line diff
--- a/+scheme/Laplace1d.m	Tue Mar 05 11:26:42 2019 -0800
+++ b/+scheme/Laplace1d.m	Tue Mar 05 12:41:50 2019 -0800
@@ -104,7 +104,7 @@
 
             tuning = 1.1;
 
-            tau1 = -(a_u/gamm_u + a_v/gamm_v) * tuning;
+            tau1 = -1/4*(a_u/gamm_u + a_v/gamm_v) * tuning;
             tau2 = 1/2*a_u;
             sig1 = -1/2;
             sig2 = 0;