Mercurial > repos > public > sbplib

comparison +scheme/Wave2dCurve.m @ 451:4e266dfe9edc feature/grids

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Add warning for possible bug
author Jonatan Werpers <jonatan@werpers.com>
date Wed, 10 May 2017 13:13:06 +0200
parents 8d455e49364f
children 459eeb99130f
comparison
equal deleted inserted replaced
450:8d455e49364f 451:4e266dfe9edc
230 a_t = spdiags(coeff_t,0,length(coeff_t),length(coeff_t)); 230 a_t = spdiags(coeff_t,0,length(coeff_t),length(coeff_t));
231 d = s*(a_n * d_n' + a_t*d_t')'; % outward facing normal derivative 231 d = s*(a_n * d_n' + a_t*d_t')'; % outward facing normal derivative
232 232
233 tau = -c.^2 * 1/beta*obj.Ji*e; 233 tau = -c.^2 * 1/beta*obj.Ji*e;
234 234
235 warning('is this right?! /c?')
235 closure{1} = halfnorm_inv*tau/c*spdiag(scale_factor)*e'; 236 closure{1} = halfnorm_inv*tau/c*spdiag(scale_factor)*e';
236 closure{2} = halfnorm_inv*tau*beta*d'; 237 closure{2} = halfnorm_inv*tau*beta*d';
237 penalty = -halfnorm_inv*tau; 238 penalty = -halfnorm_inv*tau;
238 239
239 % Unknown, boundary condition 240 % Unknown, boundary condition