Mercurial > repos > public > sbplib

changeset 1064:1341c6cea9c1 feature/laplace_curvilinear_test

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Set tuning to 1.0 in both old and new Dirichlet.
author Martin Almquist <malmquist@stanford.edu>
date Thu, 17 Jan 2019 18:18:00 -0800
parents e913cdb34dcb
children c2bd7f15da48
files +scheme/LaplaceCurvilinear.m +scheme/LaplaceCurvilinearNew.m
diffstat 2 files changed, 2 insertions(+), 15 deletions(-) [+]
line wrap: on
line diff
--- a/+scheme/LaplaceCurvilinear.m	Tue Jan 15 15:08:29 2019 -0800
+++ b/+scheme/LaplaceCurvilinear.m	Thu Jan 17 18:18:00 2019 -0800
@@ -244,8 +244,7 @@
             switch type
                 % Dirichlet boundary condition
                 case {'D','d','dirichlet'}
-                    tuning = 1.2;
-                    % tuning = 20.2;
+                    tuning = 1.0;
 
                     b1 = gamm*obj.lambda./obj.a11.^2;
                     b2 = gamm*obj.lambda./obj.a22.^2;
--- a/+scheme/LaplaceCurvilinearNew.m	Tue Jan 15 15:08:29 2019 -0800
+++ b/+scheme/LaplaceCurvilinearNew.m	Thu Jan 17 18:18:00 2019 -0800
@@ -276,19 +276,7 @@
             switch type
                 % Dirichlet boundary condition
                 case {'D','d','dirichlet'}
-                    tuning = 1.2;
-                    % tuning = 20.2;
-
-                    % b1 = gamm*obj.lambda./obj.a11.^2;
-                    % b2 = gamm*obj.lambda./obj.a22.^2;
-
-                    % tau1 = tuning * spdiag(-1./b1 - 1./b2);
-                    % tau2 = 1;
-
-                    % tau = (tau1*e + tau2*d)*H_b;
-
-                    % closure =  obj.a*Hi*tau*e';
-                    % penalty = -obj.a*Hi*tau;
+                    tuning = 1.0;
 
                     sigma = 0;
                     for i = 1:obj.dim