Mercurial > repos > public > sbplib
comparison +scheme/Laplace1d.m @ 1116:33c378e508d2
Bugfix minus sign in Laplace1d Dirichlet penalty.
| author | Martin Almquist <malmquist@stanford.edu> |
|---|---|
| date | Thu, 02 May 2019 13:25:14 -0700 |
| parents | a22b23098021 |
| children |
comparison
equal
deleted
inserted
replaced
| 1107:3230e4cbdbb4 | 1116:33c378e508d2 |
|---|---|
| 68 tau2 = 1; | 68 tau2 = 1; |
| 69 | 69 |
| 70 tau = tau1*e + tau2*d; | 70 tau = tau1*e + tau2*d; |
| 71 | 71 |
| 72 closure = obj.a*obj.Hi*tau*e'; | 72 closure = obj.a*obj.Hi*tau*e'; |
| 73 penalty = obj.a*obj.Hi*tau; | 73 penalty = -obj.a*obj.Hi*tau; |
| 74 | 74 |
| 75 % Neumann boundary condition | 75 % Neumann boundary condition |
| 76 case {'N','n','neumann'} | 76 case {'N','n','neumann'} |
| 77 tau = -e; | 77 tau = -e; |
| 78 | 78 |