Mercurial > repos > public > sbplib

changeset 1084:d5290a056049 feature/laplace_curvilinear_test

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Fix bug in new implementations that made matrices non-sparse.
author Martin Almquist <malmquist@stanford.edu>
date Wed, 27 Mar 2019 15:56:42 -0700
parents 7a55a72729e6
children 49c0b8c7330a
files +scheme/LaplaceCurvilinearNew.m +scheme/LaplaceCurvilinearNewCorner.m
diffstat 2 files changed, 6 insertions(+), 6 deletions(-) [+]
line wrap: on
line diff
--- a/+scheme/LaplaceCurvilinearNew.m	Wed Jan 23 17:11:48 2019 -0800
+++ b/+scheme/LaplaceCurvilinearNew.m	Wed Mar 27 15:56:42 2019 -0700
@@ -296,7 +296,7 @@
                 case {'D','d','dirichlet'}
                     tuning = 1.0;
 
-                    sigma = 0;
+                    sigma = 0*b_b;
                     for i = 1:obj.dim
                         sigma = sigma + e'*J*K{i,m}*K{i,m}*e;
                     end
@@ -389,8 +389,8 @@
             b_b_v = e_v'*v.b*e_v;
 
             %--- Penalty strength tau -------------
-            sigma_u = 0;
-            sigma_v = 0;
+            sigma_u = 0*b_b_u;
+            sigma_v = 0*b_b_v;
             for i = 1:obj.dim
                 sigma_u = sigma_u + e_u'*J_u*K_u{i,m_u}*K_u{i,m_u}*e_u;
                 sigma_v = sigma_v + e_v'*J_v*K_v{i,m_v}*K_v{i,m_v}*e_v;
--- a/+scheme/LaplaceCurvilinearNewCorner.m	Wed Jan 23 17:11:48 2019 -0800
+++ b/+scheme/LaplaceCurvilinearNewCorner.m	Wed Mar 27 15:56:42 2019 -0700
@@ -296,7 +296,7 @@
                 case {'D','d','dirichlet'}
                     tuning = 1.0;
 
-                    sigma = 0;
+                    sigma = 0*b_b;
                     for i = 1:obj.dim
                         sigma = sigma + e'*J*K{i,m}*K{i,m}*e;
                     end
@@ -392,8 +392,8 @@
             b_b_v = e_v'*v.b*e_v;
 
             %--- Penalty strength tau -------------
-            sigma_u = 0;
-            sigma_v = 0;
+            sigma_u = 0*b_b_u;
+            sigma_v = 0*b_b_v;
             for i = 1:obj.dim
                 sigma_u = sigma_u + e_u'*J_u*K_u{i,m_u}*K_u{i,m_u}*e_u;
                 sigma_v = sigma_v + e_v'*J_v*K_v{i,m_v}*K_v{i,m_v}*e_v;