Mercurial > repos > public > sbplib

comparison +scheme/LaplaceCurvilinearVirtaMin.m @ 1141:91058813b6e7 feature/laplace_curvilinear_test

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Multiply by 2 in VirtaMin so that scheme is actually provably stable.
author Martin Almquist <malmquist@stanford.edu>
date Fri, 21 Jun 2019 17:28:33 +0200
parents afd06a84b69c
children
comparison
equal deleted inserted replaced
1140:738de3a4058b 1141:91058813b6e7
293 case {'s', 'n'} 293 case {'s', 'n'}
294 b1 = h11 *lambda./obj.a11.^2; 294 b1 = h11 *lambda./obj.a11.^2;
295 b2 = gamm*lambda_min./obj.a22.^2; 295 b2 = gamm*lambda_min./obj.a22.^2;
296 end 296 end
297 297
298 tau1 = tuning * spdiag(-1./b1 - 1./b2); 298 tau1 = 2*tuning * spdiag(-1./b1 - 1./b2);
299 tau2 = 1; 299 tau2 = 1;
300 300
301 tau = (tau1*e + tau2*d)*H_b; 301 tau = (tau1*e + tau2*d)*H_b;
302 302
303 closure = obj.a*obj.Hi*tau*e'; 303 closure = obj.a*obj.Hi*tau*e';
421 b1_v = h11_v *lambda_v(I_v)./v.a11(I_v).^2; 421 b1_v = h11_v *lambda_v(I_v)./v.a11(I_v).^2;
422 b2_v = gamm_v*lambda_min_v(I_v)./v.a22(I_v).^2; 422 b2_v = gamm_v*lambda_min_v(I_v)./v.a22(I_v).^2;
423 end 423 end
424 424
425 tau1 = -1./(4*b1_u) -1./(4*b1_v) -1./(4*b2_u) -1./(4*b2_v); 425 tau1 = -1./(4*b1_u) -1./(4*b1_v) -1./(4*b2_u) -1./(4*b2_v);
426 tau1 = tuning * spdiag(tau1); 426 tau1 = 2*tuning * spdiag(tau1);
427 tau2 = 1/2; 427 tau2 = 1/2;
428 428
429 sig1 = -1/2; 429 sig1 = -1/2;
430 sig2 = 0; 430 sig2 = 0;
431 431