comparison +scheme/LaplaceCurvilinear.m @ 921:19c05eefc8c6 feature/utux2D

Add warning that interfaceNonConforming in LaplaceCurve only works for Cartesian grids.
author Martin Almquist <malmquist@stanford.edu>
date Sun, 02 Dec 2018 16:27:49 -0800
parents b9c98661ff5d
children 1d91c2a8aada
comparison
equal deleted inserted replaced
915:af9a63a4a5d0 921:19c05eefc8c6
322 closure = obj.a*obj.Hi*( tau*e_u' + sig*d_u'); 322 closure = obj.a*obj.Hi*( tau*e_u' + sig*d_u');
323 penalty = obj.a*obj.Hi*(-tau*e_v' + sig*d_v'); 323 penalty = obj.a*obj.Hi*(-tau*e_v' + sig*d_v');
324 end 324 end
325 325
326 function [closure, penalty] = interfaceNonConforming(obj,boundary,neighbour_scheme,neighbour_boundary,opts) 326 function [closure, penalty] = interfaceNonConforming(obj,boundary,neighbour_scheme,neighbour_boundary,opts)
327
328 % TODO: Make this work for curvilinear grids
329 warning('LaplaceCurvilinear: Non-conforming grid interpolation only works for Cartesian grids.');
327 330
328 default_field(opts, 'tuning', 1.2); 331 default_field(opts, 'tuning', 1.2);
329 tuning = opts.tuning; 332 tuning = opts.tuning;
330 333
331 % u denotes the solution in the own domain 334 % u denotes the solution in the own domain