Mercurial > repos > public > sbplib

comparison +multiblock/LaplaceSquared.m @ 1073:95113a592421 feature/grids/LaplaceSquared

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Fix bug in d3 boundary condition
author Jonatan Werpers <jonatan@werpers.com>
date Wed, 13 Feb 2019 10:24:06 +0100
parents f6b3af6febf3
children 2aaced07d1e5
comparison
equal deleted inserted replaced
1072:6468a5f6ec79 1073:95113a592421
85 tau = obj.Hi*A*d1*B_b*H_b; 85 tau = obj.Hi*A*d1*B_b*H_b;
86 closure = tau*d2'; 86 closure = tau*d2';
87 penalty = -tau; 87 penalty = -tau;
88 case 'd3' 88 case 'd3'
89 e = obj.getBoundaryOperator('e', boundary); 89 e = obj.getBoundaryOperator('e', boundary);
90 d3 = obj.getBoundaryOperator('d1', boundary); 90 d3 = obj.getBoundaryOperator('d3', boundary);
91 H_b = obj.getBoundaryQuadrature(boundary); 91 H_b = obj.getBoundaryQuadrature(boundary);
92 92
93 A = spdiag(obj.a); 93 A = spdiag(obj.a);
94 94
95 tau = -obj.Hi*A*e*H_b; 95 tau = -obj.Hi*A*e*H_b;