Mercurial > repos > public > sbplib
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 |