Mercurial > repos > public > sbplib

comparison +scheme/LaplaceCurvilinearNew.m @ 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
comparison
equal deleted inserted replaced
1084:d5290a056049 1085:49c0b8c7330a
299 sigma = 0*b_b; 299 sigma = 0*b_b;
300 for i = 1:obj.dim 300 for i = 1:obj.dim
301 sigma = sigma + e'*J*K{i,m}*K{i,m}*e; 301 sigma = sigma + e'*J*K{i,m}*K{i,m}*e;
302 end 302 end
303 sigma = sigma/s_b; 303 sigma = sigma/s_b;
304 % tau = tuning*(1/th_R + obj.dim/th_H)*sigma; 304 tau = tuning*(1/th_R + obj.dim/th_H)*sigma;
305
306 tau_R = 1/th_R*sigma;
307 tau_H = obj.dim*1/th_H*sigma;
308
309 tau = tuning*(tau_R + tau_H);
310 305
311 closure = a*Hi*d*b_b*H_b*e' ... 306 closure = a*Hi*d*b_b*H_b*e' ...
312 -a*Hi*e*tau*b_b*H_b*e'; 307 -a*Hi*e*tau*b_b*H_b*e';
313 308
314 penalty = -a*Hi*d*b_b*H_b ... 309 penalty = -a*Hi*d*b_b*H_b ...