Mercurial > repos > public > sbplib

comparison +scheme/Euler1d.m @ 968:a4ad90b37998 feature/poroelastic

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Merge with default.
author Martin Almquist <malmquist@stanford.edu>
date Sun, 23 Dec 2018 14:39:31 +0100
parents a35ed1d124d3
children 2b1b944deae1
comparison
equal deleted inserted replaced
967:368a2773f78b 968:a4ad90b37998
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);