Mercurial > repos > public > sbplib
changeset 516:afff85574ddb feature/quantumTriangles
Changed the penalty on the first derivative penalty terms, did not solve the problem
| author | Ylva Rydin <ylva.rydin@telia.com> |
|---|---|
| date | Fri, 07 Jul 2017 14:52:47 +0200 |
| parents | ba92b27da5a0 |
| children | 7a091a3527df |
| files | +scheme/Schrodinger2dCurve.m |
| diffstat | 1 files changed, 2 insertions(+), 1 deletions(-) [+] |
line wrap: on
line diff
diff -r ba92b27da5a0 -r afff85574ddb +scheme/Schrodinger2dCurve.m --- a/+scheme/Schrodinger2dCurve.m Fri Jul 07 14:22:02 2017 +0200 +++ b/+scheme/Schrodinger2dCurve.m Fri Jul 07 14:52:47 2017 +0200 @@ -197,6 +197,7 @@ end % Closure functions return the opertors applied to the own doamin to close the boundary + % Penalty functions return the opertors to force the solution. In the case of an interface it returns the operator applied to the other doamin. % boundary is a string specifying the boundary e.g. 'l','r' or 'e','w','n','s'. % type is a string specifying the type of boundary condition if there are several. @@ -213,7 +214,7 @@ F = @(t)(s * a_n(t)*d_n' + s * a_t(t) *d_t')'; tau1 = 1; a = @(t)spdiag(g(t)); - tau2 = @(t) (-1*s*a(t) - abs(a(t)))/4; + tau2 = @(t) (-1*s*a(t))/2; penalty_parameter_1 = @(t) 1*1i*halfnorm_inv_n*halfnorm_inv_t*F(t)*e'*halfnorm_t*e; penalty_parameter_2 = @(t) halfnorm_inv_n*e*tau2(t);