Mercurial > repos > public > sbplib

comparison +scheme/LaplaceCurvilinear.m @ 1115:0e5dc5dfa559 feature/laplace_curvilinear_test

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Temporary ugly fix for reversing direction of boundary in interface coupling
author Martin Almquist <malmquist@stanford.edu>
date Mon, 29 Apr 2019 12:55:51 -0700
parents 9ddd2e89d018
children 0aed89043ad6
comparison
equal deleted inserted replaced
1110:9ddd2e89d018 1115:0e5dc5dfa559
370 d_v = v.getBoundaryOperator('d', neighbour_boundary); 370 d_v = v.getBoundaryOperator('d', neighbour_boundary);
371 s_b_v = v.getBoundaryScaling(neighbour_boundary); 371 s_b_v = v.getBoundaryScaling(neighbour_boundary);
372 [th_H_v, ~, th_R_v] = v.getBoundaryBorrowing(neighbour_boundary); 372 [th_H_v, ~, th_R_v] = v.getBoundaryBorrowing(neighbour_boundary);
373 m_v = v.getBoundaryNumber(neighbour_boundary); 373 m_v = v.getBoundaryNumber(neighbour_boundary);
374 374
375 % BUGFIX?!?!?
376 if (strcmp(boundary,'s') && strcmp(neighbour_boundary,'e')) || (strcmp(boundary,'e') && strcmp(neighbour_boundary,'s'))
377 e_v = fliplr(e_v);
378 d_v = fliplr(d_v);
379 s_b_v = rot90(s_b_v,2);
380 end
381
375 % Coefficients, v 382 % Coefficients, v
376 K_v = v.K; 383 K_v = v.K;
377 J_v = v.J; 384 J_v = v.J;
378 b_b_v = e_v'*v.b*e_v; 385 b_b_v = e_v'*v.b*e_v;
379 386