sbplib: +scheme/Beam2d.m annotate

annotate +scheme/Beam2d.m @ 998:2b1b944deae1 feature/getBoundaryOp

Add getBoundaryOperator to all 1d schemes. Did not add getBoundaryQuadrature because it doesnt make sense in 1d (?)

author	Martin Almquist <malmquist@stanford.edu>
date	Sat, 12 Jan 2019 13:35:19 -0800
parents	78db023a7fe3
children

rev	line source
141 cb2b12246b7e Fixed path to some superclasses. Jonatan Werpers <jonatan@werpers.com> parents: 125 diff changeset	1 classdef Beam2d < scheme.Scheme
0 48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	2 properties
175 8f22829b69d0 Added and upgraded schemes for the beam equation in 1d and 2d. Jonatan Werpers <jonatan@werpers.com> parents: 141 diff changeset	3 grid
0 48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	4 order % Order accuracy for the approximation
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	5
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	6 D % non-stabalized scheme operator
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	7 M % Derivative norm
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	8 alpha
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	9
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	10 H % Discrete norm
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	11 Hi
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	12 H_x, H_y % Norms in the x and y directions
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	13 Hx,Hy % Kroneckerd norms. 1'Hxv corresponds to integration in the x dir.
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	14 Hi_x, Hi_y
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	15 Hix, Hiy
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	16 e_w, e_e, e_s, e_n
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	17 d1_w, d1_e, d1_s, d1_n
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	18 d2_w, d2_e, d2_s, d2_n
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	19 d3_w, d3_e, d3_s, d3_n
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	20 gamm_x, gamm_y
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	21 delt_x, delt_y
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	22 end
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	23
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	24 methods
125 d52e5cdb6eff Fixed some class name, file name errors. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	25 function obj = Beam2d(m,lim,order,alpha,opsGen)
175 8f22829b69d0 Added and upgraded schemes for the beam equation in 1d and 2d. Jonatan Werpers <jonatan@werpers.com> parents: 141 diff changeset	26 default_arg('alpha',1);
0 48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	27 default_arg('opsGen',@sbp.Higher);
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	28
176 d095b5396103 Fixed some bugs in Beam schemes. Jonatan Werpers <jonatan@werpers.com> parents: 175 diff changeset	29 if ~isa(grid, 'grid.Cartesian') \|\| grid.D() ~= 2
175 8f22829b69d0 Added and upgraded schemes for the beam equation in 1d and 2d. Jonatan Werpers <jonatan@werpers.com> parents: 141 diff changeset	30 error('Grid must be 2d cartesian');
0 48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	31 end
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	32
175 8f22829b69d0 Added and upgraded schemes for the beam equation in 1d and 2d. Jonatan Werpers <jonatan@werpers.com> parents: 141 diff changeset	33 obj.grid = grid;
8f22829b69d0 Added and upgraded schemes for the beam equation in 1d and 2d. Jonatan Werpers <jonatan@werpers.com> parents: 141 diff changeset	34 obj.alpha = alpha;
8f22829b69d0 Added and upgraded schemes for the beam equation in 1d and 2d. Jonatan Werpers <jonatan@werpers.com> parents: 141 diff changeset	35 obj.order = order;
0 48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	36
175 8f22829b69d0 Added and upgraded schemes for the beam equation in 1d and 2d. Jonatan Werpers <jonatan@werpers.com> parents: 141 diff changeset	37 m_x = grid.m(1);
8f22829b69d0 Added and upgraded schemes for the beam equation in 1d and 2d. Jonatan Werpers <jonatan@werpers.com> parents: 141 diff changeset	38 m_y = grid.m(2);
0 48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	39
175 8f22829b69d0 Added and upgraded schemes for the beam equation in 1d and 2d. Jonatan Werpers <jonatan@werpers.com> parents: 141 diff changeset	40 h = grid.scaling();
8f22829b69d0 Added and upgraded schemes for the beam equation in 1d and 2d. Jonatan Werpers <jonatan@werpers.com> parents: 141 diff changeset	41 h_x = h(1);
8f22829b69d0 Added and upgraded schemes for the beam equation in 1d and 2d. Jonatan Werpers <jonatan@werpers.com> parents: 141 diff changeset	42 h_y = h(2);
0 48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	43
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	44 ops_x = opsGen(m_x,h_x,order);
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	45 ops_y = opsGen(m_y,h_y,order);
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	46
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	47 I_x = speye(m_x);
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	48 I_y = speye(m_y);
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	49
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	50 D4_x = sparse(ops_x.derivatives.D4);
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	51 H_x = sparse(ops_x.norms.H);
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	52 Hi_x = sparse(ops_x.norms.HI);
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	53 e_l_x = sparse(ops_x.boundary.e_1);
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	54 e_r_x = sparse(ops_x.boundary.e_m);
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	55 d1_l_x = sparse(ops_x.boundary.S_1);
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	56 d1_r_x = sparse(ops_x.boundary.S_m);
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	57 d2_l_x = sparse(ops_x.boundary.S2_1);
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	58 d2_r_x = sparse(ops_x.boundary.S2_m);
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	59 d3_l_x = sparse(ops_x.boundary.S3_1);
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	60 d3_r_x = sparse(ops_x.boundary.S3_m);
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	61
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	62 D4_y = sparse(ops_y.derivatives.D4);
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	63 H_y = sparse(ops_y.norms.H);
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	64 Hi_y = sparse(ops_y.norms.HI);
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	65 e_l_y = sparse(ops_y.boundary.e_1);
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	66 e_r_y = sparse(ops_y.boundary.e_m);
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	67 d1_l_y = sparse(ops_y.boundary.S_1);
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	68 d1_r_y = sparse(ops_y.boundary.S_m);
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	69 d2_l_y = sparse(ops_y.boundary.S2_1);
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	70 d2_r_y = sparse(ops_y.boundary.S2_m);
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	71 d3_l_y = sparse(ops_y.boundary.S3_1);
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	72 d3_r_y = sparse(ops_y.boundary.S3_m);
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	73
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	74
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	75 D4 = kr(D4_x, I_y) + kr(I_x, D4_y);
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	76
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	77 % Norms
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	78 obj.H = kr(H_x,H_y);
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	79 obj.Hx = kr(H_x,I_x);
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	80 obj.Hy = kr(I_x,H_y);
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	81 obj.Hix = kr(Hi_x,I_y);
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	82 obj.Hiy = kr(I_x,Hi_y);
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	83 obj.Hi = kr(Hi_x,Hi_y);
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	84
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	85 % Boundary operators
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	86 obj.e_w = kr(e_l_x,I_y);
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	87 obj.e_e = kr(e_r_x,I_y);
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	88 obj.e_s = kr(I_x,e_l_y);
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	89 obj.e_n = kr(I_x,e_r_y);
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	90 obj.d1_w = kr(d1_l_x,I_y);
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	91 obj.d1_e = kr(d1_r_x,I_y);
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	92 obj.d1_s = kr(I_x,d1_l_y);
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	93 obj.d1_n = kr(I_x,d1_r_y);
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	94 obj.d2_w = kr(d2_l_x,I_y);
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	95 obj.d2_e = kr(d2_r_x,I_y);
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	96 obj.d2_s = kr(I_x,d2_l_y);
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	97 obj.d2_n = kr(I_x,d2_r_y);
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	98 obj.d3_w = kr(d3_l_x,I_y);
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	99 obj.d3_e = kr(d3_r_x,I_y);
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	100 obj.d3_s = kr(I_x,d3_l_y);
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	101 obj.d3_n = kr(I_x,d3_r_y);
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	102
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	103 obj.D = alpha*D4;
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	104
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	105 obj.gamm_x = h_x*ops_x.borrowing.N.S2/2;
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	106 obj.delt_x = h_x^3*ops_x.borrowing.N.S3/2;
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	107
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	108 obj.gamm_y = h_y*ops_y.borrowing.N.S2/2;
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	109 obj.delt_y = h_y^3*ops_y.borrowing.N.S3/2;
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	110 end
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	111
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	112
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	113 % Closure functions return the opertors applied to the own doamin to close the boundary
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	114 % Penalty functions return the opertors to force the solution. In the case of an interface it returns the operator applied to the other doamin.
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	115 % boundary is a string specifying the boundary e.g. 'l','r' or 'e','w','n','s'.
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	116 % type is a string specifying the type of boundary condition if there are several.
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	117 % data is a function returning the data that should be applied at the boundary.
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	118 % neighbour_scheme is an instance of Scheme that should be interfaced to.
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	119 % neighbour_boundary is a string specifying which boundary to interface to.
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	120 function [closure, penalty_e,penalty_d] = boundary_condition(obj,boundary,type,data)
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	121 default_arg('type','dn');
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	122 default_arg('data',0);
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	123
997 78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	124 [e, d1, d2, d3] = obj.getBoundaryOperator({'e', 'd1', 'd2', 'd3'}, boundary);
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	125 s = obj.getBoundarySign(boundary);
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	126 [gamm, delt] = obj.getBoundaryBorrowing(boundary);
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	127 halfnorm_inv = obj.getHalfnormInv(boundary);
0 48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	128
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	129 switch type
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	130 % Dirichlet-neumann boundary condition
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	131 case {'dn'}
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	132 alpha = obj.alpha;
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	133
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	134 % tau1 < -alpha^2/gamma
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	135 tuning = 1.1;
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	136
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	137 tau1 = tuning * alpha/delt;
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	138 tau4 = s*alpha;
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	139
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	140 sig2 = tuning * alpha/gamm;
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	141 sig3 = -s*alpha;
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	142
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	143 tau = tau1e+tau4d3;
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	144 sig = sig2d1+sig3d2;
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	145
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	146 closure = halfnorm_inv(taue' + sig*d1');
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	147
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	148 pp_e = halfnorm_inv*tau;
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	149 pp_d = halfnorm_inv*sig;
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	150 switch class(data)
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	151 case 'double'
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	152 penalty_e = pp_e*data;
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	153 penalty_d = pp_d*data;
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	154 case 'function_handle'
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	155 penalty_e = @(t)pp_e*data(t);
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	156 penalty_d = @(t)pp_d*data(t);
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	157 otherwise
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	158 error('Wierd data argument!')
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	159 end
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	160
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	161 % Unknown, boundary condition
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	162 otherwise
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	163 error('No such boundary condition: type = %s',type);
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	164 end
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	165 end
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	166
944 a35ed1d124d3 Change from opts to type in a few schemes Jonatan Werpers <jonatan@werpers.com> parents: 910 diff changeset	167 function [closure, penalty] = interface(obj,boundary,neighbour_scheme,neighbour_boundary, type)
0 48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	168 % u denotes the solution in the own domain
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	169 % v denotes the solution in the neighbour domain
997 78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	170 [e_u, d1_u, d2_u, d3_u] = obj.getBoundaryOperator({'e', 'd1', 'd2', 'd3'}, boundary);
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	171 s_u = obj.getBoundarySign(boundary);
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	172 [gamm_u, delt_u] = obj.getBoundaryBorrowing(boundary);
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	173 halfnorm_inv = obj.getHalfnormInv(boundary);
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	174
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	175 [e_v, d1_v, d2_v, d3_v] = neighbour_scheme.getBoundaryOperator({'e', 'd1', 'd2', 'd3'}, neighbour_boundary);
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	176 s_v = neighbour_scheme.getBoundarySign(neighbour_boundary);
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	177 [gamm_v, delt_v] = neighbour_scheme.getBoundaryBorrowing(neighbour_boundary);
0 48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	178
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	179 tuning = 2;
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	180
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	181 alpha_u = obj.alpha;
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	182 alpha_v = neighbour_scheme.alpha;
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	183
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	184 tau1 = ((alpha_u/2)/delt_u + (alpha_v/2)/delt_v)/2*tuning;
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	185 % tau1 = (alpha_u/2 + alpha_v/2)/(2delt_u)tuning;
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	186 tau4 = s_u*alpha_u/2;
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	187
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	188 sig2 = ((alpha_u/2)/gamm_u + (alpha_v/2)/gamm_v)/2*tuning;
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	189 sig3 = -s_u*alpha_u/2;
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	190
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	191 phi2 = s_u*1/2;
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	192
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	193 psi1 = -s_u*1/2;
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	194
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	195 tau = tau1e_u + tau4d3_u;
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	196 sig = sig2d1_u + sig3d2_u ;
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	197 phi = phi2*d1_u ;
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	198 psi = psi1*e_u ;
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	199
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	200 closure = halfnorm_inv(taue_u' + sigd1_u' + phialpha_ud2_u' + psialpha_u*d3_u');
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	201 penalty = -halfnorm_inv(taue_v' + sigd1_v' + phialpha_vd2_v' + psialpha_v*d3_v');
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	202 end
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	203
997 78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	204 % Returns the boundary operator op for the boundary specified by the string boundary.
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	205 % op -- string or a cell array of strings
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	206 % boundary -- string
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	207 function varargout = getBoundaryOperator(obj, op, boundary)
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	208
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	209 if ~iscell(op)
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	210 op = {op};
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	211 end
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	212
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	213 for i = 1:numel(op)
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	214 switch op{i}
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	215 case 'e'
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	216 switch boundary
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	217 case 'w'
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	218 e = obj.e_w;
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	219 case 'e'
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	220 e = obj.e_e;
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	221 case 's'
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	222 e = obj.e_s;
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	223 case 'n'
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	224 e = obj.e_n;
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	225 otherwise
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	226 error('No such boundary: boundary = %s',boundary);
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	227 end
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	228 varargout{i} = e;
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	229
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	230 case 'd1'
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	231 switch boundary
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	232 case 'w'
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	233 d1 = obj.d1_w;
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	234 case 'e'
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	235 d1 = obj.d1_e;
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	236 case 's'
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	237 d1 = obj.d1_s;
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	238 case 'n'
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	239 d1 = obj.d1_n;
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	240 otherwise
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	241 error('No such boundary: boundary = %s',boundary);
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	242 end
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	243 varargout{i} = d1;
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	244 end
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	245
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	246 case 'd2'
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	247 switch boundary
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	248 case 'w'
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	249 d2 = obj.d2_w;
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	250 case 'e'
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	251 d2 = obj.d2_e;
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	252 case 's'
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	253 d2 = obj.d2_s;
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	254 case 'n'
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	255 d2 = obj.d2_n;
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	256 otherwise
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	257 error('No such boundary: boundary = %s',boundary);
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	258 end
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	259 varargout{i} = d2;
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	260 end
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	261
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	262 case 'd3'
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	263 switch boundary
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	264 case 'w'
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	265 d3 = obj.d3_w;
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	266 case 'e'
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	267 d3 = obj.d3_e;
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	268 case 's'
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	269 d3 = obj.d3_s;
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	270 case 'n'
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	271 d3 = obj.d3_n;
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	272 otherwise
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	273 error('No such boundary: boundary = %s',boundary);
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	274 end
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	275 varargout{i} = d3;
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	276 end
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	277 end
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	278 end
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	279
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	280 % Returns square boundary quadrature matrix, of dimension
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	281 % corresponding to the number of boundary points
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	282 %
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	283 % boundary -- string
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	284 function H_b = getBoundaryQuadrature(obj, boundary)
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	285
0 48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	286 switch boundary
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	287 case 'w'
997 78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	288 H_b = obj.H_y;
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	289 case 'e'
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	290 H_b = obj.H_y;
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	291 case 's'
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	292 H_b = obj.H_x;
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	293 case 'n'
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	294 H_b = obj.H_x;
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	295 otherwise
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	296 error('No such boundary: boundary = %s',boundary);
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	297 end
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	298 end
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	299
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	300 % Returns the boundary sign. The right boundary is considered the positive boundary
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	301 % boundary -- string
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	302 function s = getBoundarySign(obj, boundary)
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	303 switch boundary
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	304 case {'e','n'}
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	305 s = 1;
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	306 case {'w','s'}
0 48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	307 s = -1;
997 78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	308 otherwise
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	309 error('No such boundary: boundary = %s',boundary);
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	310 end
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	311 end
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	312
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	313 % Returns the halfnorm_inv used in SATs. TODO: better notation
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	314 function Hinv = getHalfnormInv(obj, boundary)
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	315 switch boundary
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	316 case 'w'
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	317 Hinv = obj.Hix;
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	318 case 'e'
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	319 Hinv = obj.Hix;
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	320 case 's'
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	321 Hinv = obj.Hiy;
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	322 case 'n'
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	323 Hinv = obj.Hiy;
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	324 otherwise
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	325 error('No such boundary: boundary = %s',boundary);
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	326 end
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	327 end
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	328
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	329 % Returns borrowing constant gamma
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	330 % boundary -- string
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	331 function [gamm, delt] = getBoundaryBorrowing(obj, boundary)
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	332 switch boundary
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	333 case {'w','e'}
0 48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	334 gamm = obj.gamm_x;
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	335 delt = obj.delt_x;
997 78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 944 diff changeset	336 case {'s','n'}
0 48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	337 gamm = obj.gamm_y;
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	338 delt = obj.delt_y;
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	339 otherwise
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	340 error('No such boundary: boundary = %s',boundary);
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	341 end
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	342 end
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	343
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	344 function N = size(obj)
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	345 N = prod(obj.m);
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	346 end
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	347
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	348 end
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	349 end

Mercurial > repos > public > sbplib

annotate +scheme/Beam2d.m @ 998:2b1b944deae1 feature/getBoundaryOp