Mercurial > repos > public > sbplib

comparison +scheme/Euler1d.m @ 905:459eeb99130f feature/utux2D

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Include type as (optional) input parameter in the interface method of all schemes.
author Martin Almquist <malmquist@stanford.edu>
date Thu, 22 Nov 2018 22:03:44 -0800
parents f39f98b59f61
children b9c98661ff5d
comparison
equal deleted inserted replaced
904:14b093a344eb 905:459eeb99130f
444 end 444 end
445 445
446 closure = @closure_fun; 446 closure = @closure_fun;
447 end 447 end
448 448
449 function [closure, penalty] = interface(obj,boundary,neighbour_scheme,neighbour_boundary) 449 function [closure, penalty] = interface(obj,boundary,neighbour_scheme,neighbour_boundary,type)
450 error('NOT DONE') 450 error('NOT DONE')
451 % u denotes the solution in the own domain 451 % u denotes the solution in the own domain
452 % v denotes the solution in the neighbour domain 452 % v denotes the solution in the neighbour domain
453 [e_u,d1_u,d2_u,d3_u,s_u,gamm_u,delt_u, halfnorm_inv] = obj.get_boundary_ops(boundary); 453 [e_u,d1_u,d2_u,d3_u,s_u,gamm_u,delt_u, halfnorm_inv] = obj.get_boundary_ops(boundary);
454 [e_v,d1_v,d2_v,d3_v,s_v,gamm_v,delt_v] = neighbour_scheme.get_boundary_ops(neighbour_boundary); 454 [e_v,d1_v,d2_v,d3_v,s_v,gamm_v,delt_v] = neighbour_scheme.get_boundary_ops(neighbour_boundary);