Mercurial > repos > public > sbplib
changeset 96:19d0c9325a3e
Wave2dCurve: Removed some comments.
| author | Jonatan Werpers <jonatan@werpers.com> |
|---|---|
| date | Wed, 02 Dec 2015 16:21:16 +0100 |
| parents | 480ef676742f |
| children | 33dba20b4b9d ce4eecbcb915 |
| files | +scheme/Wave2dCurve.m |
| diffstat | 1 files changed, 2 insertions(+), 3 deletions(-) [+] |
line wrap: on
line diff
diff -r 480ef676742f -r 19d0c9325a3e +scheme/Wave2dCurve.m --- a/+scheme/Wave2dCurve.m Wed Dec 02 15:51:28 2015 +0100 +++ b/+scheme/Wave2dCurve.m Wed Dec 02 16:21:16 2015 +0100 @@ -264,12 +264,11 @@ b2_v = gamm_v*v.lambda(I_v)./v.a22(I_v).^2; tau = -1./(4*b1_u) -1./(4*b1_v) -1./(4*b2_u) -1./(4*b2_v); - tau = tuning * spdiag(tau(:)); % Probably correct until here, see eq 27 + tau = tuning * spdiag(tau(:)); sig1 = 1/2; sig2 = -1/2; - % penalty_parameter_1 = halfnorm_inv_u_n*(tau + sig1*halfnorm_inv_u_t*F_u*e_u'*halfnorm_u_t)*e_u; %% This is what is in the paper, but there is an error in dimensions. - penalty_parameter_1 = halfnorm_inv_u_n*(e_u*tau + sig1*halfnorm_inv_u_t*F_u*e_u'*halfnorm_u_t*e_u); %% Random guess at a fix, should check theory for this. + penalty_parameter_1 = halfnorm_inv_u_n*(e_u*tau + sig1*halfnorm_inv_u_t*F_u*e_u'*halfnorm_u_t*e_u); penalty_parameter_2 = halfnorm_inv_u_n * sig2 * e_u;