Mercurial > repos > public > sbplib

comparison +scheme/Staggered1DAcousticsVariable.m @ 1081:9c74d96fc9e2 feature/d1_staggered

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Bugfix in interface SATs, Staggered1DAcousticsVariable
author Martin Almquist <malmquist@stanford.edu>
date Tue, 05 Mar 2019 11:18:44 -0800
parents 18e10217dca9
children
comparison
equal deleted inserted replaced
920:386ef449df51 1081:9c74d96fc9e2
330 Z_u = sparse(Z_u); 330 Z_u = sparse(Z_u);
331 331
332 Z_v = TH_v*Ttilde_p*Sp_u; 332 Z_v = TH_v*Ttilde_p*Sp_u;
333 Z_v = sparse(Z_v); 333 Z_v = sparse(Z_v);
334 334
335 closure_u = Hi_u*e_u*Z_u*e_u'; 335 closure_u = Hi_u*e_u*(A_u\Z_u)*e_u';
336 penalty_u = -Hi_u*e_u*Z_u*e_v'; 336 penalty_u = -Hi_u*e_u*(A_u\Z_u)*e_v';
337 closure_v = Hi_v*e_v*Z_v*e_v'; 337 closure_v = Hi_v*e_v*(A_v\Z_v)*e_v';
338 penalty_v = -Hi_v*e_v*Z_v*e_u'; 338 penalty_v = -Hi_v*e_v*(A_v\Z_v)*e_u';
339 339
340 switch boundary 340 switch boundary
341 case 'l' 341 case 'l'
342 closure = closure_v; 342 closure = closure_v;
343 penalty = penalty_v; 343 penalty = penalty_v;