Mercurial > repos > public > sbplib

comparison +scheme/Euler1d.m @ 1018:eca4d9df9365 feature/advectionRV

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Merge with default for latest changes
author Vidar Stiernström <vidar.stiernstrom@it.uu.se>
date Tue, 11 Dec 2018 17:31:28 +0100
parents a35ed1d124d3
children 2b1b944deae1
comparison
equal deleted inserted replaced
1017:2d7c1333bd6c 1018:eca4d9df9365
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);