Mercurial > repos > public > sbplib

changeset 1085:49c0b8c7330a feature/laplace_curvilinear_test

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Clean up Dirichlet BC code in LaplaceCurvilinearNew.
author Martin Almquist <malmquist@stanford.edu>
date Fri, 29 Mar 2019 14:24:39 -0700
parents d5290a056049
children 90c23fd08f59
files +scheme/LaplaceCurvilinearNew.m
diffstat 1 files changed, 1 insertions(+), 6 deletions(-) [+]
line wrap: on
line diff
diff -r d5290a056049 -r 49c0b8c7330a +scheme/LaplaceCurvilinearNew.m
--- a/+scheme/LaplaceCurvilinearNew.m	Wed Mar 27 15:56:42 2019 -0700
+++ b/+scheme/LaplaceCurvilinearNew.m	Fri Mar 29 14:24:39 2019 -0700
@@ -301,12 +301,7 @@
                         sigma = sigma + e'*J*K{i,m}*K{i,m}*e;
                     end
                     sigma = sigma/s_b;
-                    % tau = tuning*(1/th_R + obj.dim/th_H)*sigma;
-
-                    tau_R = 1/th_R*sigma;
-                    tau_H = obj.dim*1/th_H*sigma;
-
-                    tau = tuning*(tau_R + tau_H);
+                    tau = tuning*(1/th_R + obj.dim/th_H)*sigma;
 
                     closure = a*Hi*d*b_b*H_b*e' ...
                              -a*Hi*e*tau*b_b*H_b*e';