Mercurial > repos > public > sbplib

comparison +scheme/Euler1d.m @ 1033:037f203b9bf5 feature/burgers1d

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Merge with branch feature/advectioRV to utilize the +rv package
author Vidar Stiernström <vidar.stiernstrom@it.uu.se>
date Thu, 17 Jan 2019 10:44:12 +0100
parents a35ed1d124d3
children 2b1b944deae1
comparison
equal deleted inserted replaced
854:18162a0a5bb5 1033:037f203b9bf5
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);