Mercurial > repos > public > sbplib

comparison +scheme/Laplace1D.m @ 946:706d1c2b4199 feature/utux2D

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Raname opts to type in a bunch of interface methods
author Jonatan Werpers <jonatan@werpers.com>
date Wed, 05 Dec 2018 15:57:41 +0100
parents 21394c78c72e
children
comparison
equal deleted inserted replaced
945:af397cc35239 946:706d1c2b4199
81 otherwise 81 otherwise
82 error('No such boundary condition: type = %s',type); 82 error('No such boundary condition: type = %s',type);
83 end 83 end
84 end 84 end
85 85
86 function [closure, penalty] = interface(obj,boundary,neighbour_scheme,neighbour_boundary,opts) 86 function [closure, penalty] = interface(obj, boundary, neighbour_scheme, neighbour_boundary, type)
87 % u denotes the solution in the own domain 87 % u denotes the solution in the own domain
88 % v denotes the solution in the neighbour domain 88 % v denotes the solution in the neighbour domain
89 89
90 [e_u,d_u,s_u] = obj.get_boundary_ops(boundary); 90 [e_u,d_u,s_u] = obj.get_boundary_ops(boundary);
91 [e_v,d_v,s_v] = neighbour_scheme.get_boundary_ops(neighbour_boundary); 91 [e_v,d_v,s_v] = neighbour_scheme.get_boundary_ops(neighbour_boundary);
96 96
97 gamm_u = obj.gamm; 97 gamm_u = obj.gamm;
98 gamm_v = neighbour_scheme.gamm; 98 gamm_v = neighbour_scheme.gamm;
99 99
100 tuning = 1.1; 100 tuning = 1.1;
101 101
102 tau1 = -(a_u/gamm_u + a_v/gamm_v) * tuning; 102 tau1 = -(a_u/gamm_u + a_v/gamm_v) * tuning;
103 tau2 = 1/2*a_u; 103 tau2 = 1/2*a_u;
104 sig1 = -1/2; 104 sig1 = -1/2;
105 sig2 = 0; 105 sig2 = 0;
106 106