Mercurial > repos > public > sbplib

changeset 1117:27019aca2f17 feature/poroelastic

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Fix bug in LaplCurveNewCorner that makes bc matrix non-sparse.
author Martin Almquist <malmquist@stanford.edu>
date Sat, 04 May 2019 14:44:34 -0700
parents 753de514ae77
children 07d0caf915e4
files +scheme/LaplaceCurvilinearNewCorner.m
diffstat 1 files changed, 5 insertions(+), 5 deletions(-) [+]
line wrap: on
line diff
--- a/+scheme/LaplaceCurvilinearNewCorner.m	Sun Mar 03 15:00:42 2019 -0800
+++ b/+scheme/LaplaceCurvilinearNewCorner.m	Sat May 04 14:44:34 2019 -0700
@@ -63,13 +63,13 @@
             % assert(isa(g, 'grid.Curvilinear'))
             if isa(a, 'function_handle')
                 a = grid.evalOn(g, a);
-                a = spdiag(a);
             end
+            a = spdiag(a);
 
             if isa(b, 'function_handle')
                 b = grid.evalOn(g, b);
-                b = spdiag(b);
             end
+            b = spdiag(b);
 
             % If b is scalar
             if length(b) == 1
@@ -305,7 +305,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
@@ -401,8 +401,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;