Mercurial > repos > public > sbplib

changeset 1141:91058813b6e7 feature/laplace_curvilinear_test

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Multiply by 2 in VirtaMin so that scheme is actually provably stable.
author Martin Almquist <malmquist@stanford.edu>
date Fri, 21 Jun 2019 17:28:33 +0200
parents 738de3a4058b
children 5e7692ed7c7c
files +scheme/LaplaceCurvilinearVirtaMin.m
diffstat 1 files changed, 2 insertions(+), 2 deletions(-) [+]
line wrap: on
line diff
--- a/+scheme/LaplaceCurvilinearVirtaMin.m	Fri Jun 21 17:28:02 2019 +0200
+++ b/+scheme/LaplaceCurvilinearVirtaMin.m	Fri Jun 21 17:28:33 2019 +0200
@@ -295,7 +295,7 @@
                         b2 = gamm*lambda_min./obj.a22.^2;
                     end
 
-                    tau1 = tuning * spdiag(-1./b1 - 1./b2);
+                    tau1 = 2*tuning * spdiag(-1./b1 - 1./b2);
                     tau2 = 1;
 
                     tau = (tau1*e + tau2*d)*H_b;
@@ -423,7 +423,7 @@
             end
 
             tau1 = -1./(4*b1_u) -1./(4*b1_v) -1./(4*b2_u) -1./(4*b2_v);
-            tau1 = tuning * spdiag(tau1);
+            tau1 = 2*tuning * spdiag(tau1);
             tau2 = 1/2;
 
             sig1 = -1/2;