Mercurial > repos > public > sbplib

comparison +scheme/LaplaceCurvilinear.m @ 563:8f780429dab2 feature/grids/laplace_refactor

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Fix sign error in Neumann BC
author Jonatan Werpers <jonatan@werpers.com>
date Wed, 30 Aug 2017 17:29:33 +0200
parents 11d8d6ccbcd7
children 9bf49338f8e6
comparison
equal deleted inserted replaced
562:11d8d6ccbcd7 563:8f780429dab2
258 258
259 % Neumann boundary condition 259 % Neumann boundary condition
260 case {'N','n','neumann'} 260 case {'N','n','neumann'}
261 tau1 = -1; 261 tau1 = -1;
262 tau2 = 0; 262 tau2 = 0;
263 tau = s*obj.a*obj.Ji*(tau1*e + tau2*d); 263 tau = obj.a*obj.Ji*(tau1*e + tau2*d);
264 264
265 closure = halfnorm_inv*tau*d'; 265 closure = halfnorm_inv*tau*d';
266 penalty = -halfnorm_inv*tau; 266 penalty = -halfnorm_inv*tau;
267 267
268 % Characteristic boundary condition 268 % Characteristic boundary condition