Mercurial > repos > public > sbplib
changeset 82:df642750e8f7
Added Dirichelt BC to Wave2dCurve.
| author | Jonatan Werpers <jonatan@werpers.com> |
|---|---|
| date | Thu, 26 Nov 2015 16:24:25 +0100 |
| parents | 9c0192cf099f |
| children | 08fbc284718f |
| files | +scheme/Wave2dCurve.m |
| diffstat | 1 files changed, 20 insertions(+), 11 deletions(-) [+] |
line wrap: on
line diff
--- a/+scheme/Wave2dCurve.m Thu Nov 26 16:24:02 2015 +0100 +++ b/+scheme/Wave2dCurve.m Thu Nov 26 16:24:25 2015 +0100 @@ -175,19 +175,29 @@ switch type % Dirichlet boundary condition case {'D','d','dirichlet'} - error('not implemented') - alpha = obj.alpha; + % v denotes the solution in the neighbour domain + tuning = 1.2; + % tuning = 20.2; + [e, d_n, d_t, coeff_n, coeff_t, s, gamm, halfnorm_inv_n, halfnorm_inv_t, halfnorm_t] = obj.get_boundary_ops(boundary); + + a_n = spdiag(coeff_n); + a_t = spdiag(coeff_t); + + F = (s * a_n * d_n' + s * a_t*d_t')'; + + u = obj; - % tau1 < -alpha^2/gamma - tuning = 1.1; - tau1 = -tuning*alpha/gamm; - tau2 = s*alpha; + b1 = gamm*u.lambda./u.a11.^2; + b2 = gamm*u.lambda./u.a22.^2; - p = tau1*e + tau2*d; + tau = -1./b1 - 1./b2; + tau = tuning * spdiag(tau(:)); + sig1 = 1/2; - closure = halfnorm_inv*p*e'; + penalty_parameter_1 = halfnorm_inv_n*(tau + sig1*halfnorm_inv_t*F*e'*halfnorm_t)*e; - pp = halfnorm_inv*p; + closure = obj.Ji*obj.c^2 * penalty_parameter_1*e'; + pp = -obj.Ji*obj.c^2 * penalty_parameter_1; switch class(data) case 'double' penalty = pp*data; @@ -254,8 +264,7 @@ b2_v = gamm_v*v.lambda./v.a22.^2; - tau = -1./(4*b1_u) -1./(4*b1_v) -1./(4*b2_u) -1./(4*b2_v); - m_tot = obj.m(1)*obj.m(2); + tau = -1./(4*b1_u) -1./(4*b1_v) -1./(4*b2_u) -1./(4*b2_v); tau = tuning * spdiag(tau(:)); sig1 = 1/2; sig2 = -1/2*s_u;