sbplib: +scheme/LaplaceCurvilinearNewCorner.m annotate

annotate +scheme/LaplaceCurvilinearNewCorner.m @ 1060:e40899094f20 feature/poroelastic

Elastic2dVariable: Clean up alpha in getBoundaryOperator. Add alpha1 and alpha2 as boundary operators.

author	Martin Almquist <malmquist@stanford.edu>
date	Sat, 26 Jan 2019 16:56:35 -0800
parents	b828e589d540
children	753de514ae77

rev	line source
1055 b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	1 classdef LaplaceCurvilinearNewCorner < scheme.Scheme
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	2 properties
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	3 m % Number of points in each direction, possibly a vector
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	4 h % Grid spacing
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	5 dim % Number of spatial dimensions
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	6
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	7 grid
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	8
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	9 order % Order of accuracy for the approximation
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	10
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	11 a,b % Parameters of the operator
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	12
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	13
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	14 % Inner products and operators for physical coordinates
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	15 D % Laplace operator
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	16 H, Hi % Inner product
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	17 e_w, e_e, e_s, e_n
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	18 d_w, d_e, d_s, d_n % Normal derivatives at the boundary
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	19 H_w, H_e, H_s, H_n % Boundary inner products
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	20 Dx, Dy % Physical derivatives
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	21 M % Gradient inner product
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	22
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	23 % Metric coefficients
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	24 J, Ji
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	25 a11, a12, a22
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	26 K
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	27 x_u
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	28 x_v
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	29 y_u
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	30 y_v
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	31 s_w, s_e, s_s, s_n % Boundary integral scale factors
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	32
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	33 % Inner product and operators for logical coordinates
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	34 H_u, H_v % Norms in the x and y directions
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	35 Hi_u, Hi_v
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	36 Hu,Hv % Kroneckerd norms. 1'Hxv corresponds to integration in the x dir.
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	37 Hiu, Hiv
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	38 du_w, dv_w
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	39 du_e, dv_e
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	40 du_s, dv_s
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	41 du_n, dv_n
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	42
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	43 % Borrowing constants
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	44 theta_M_u, theta_M_v
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	45 theta_R_u, theta_R_v
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	46 theta_H_u, theta_H_v
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	47
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	48 % Temporary
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	49 lambda
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	50 gamm_u, gamm_v
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	51
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	52 end
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	53
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	54 methods
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	55 % Implements adiv(bgrad(u)) as a SBP scheme
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	56 % TODO: Implement proper H, it should be the real physical quadrature, the logic quadrature may be but in a separate variable (H_logic?)
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	57
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	58 function obj = LaplaceCurvilinearNewCorner(g, order, a, b, opSet)
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	59 default_arg('opSet',@sbp.D2Variable);
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	60 default_arg('a', 1);
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	61 default_arg('b', @(x,y) 0*x + 1);
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	62
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	63 % assert(isa(g, 'grid.Curvilinear'))
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	64 if isa(a, 'function_handle')
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	65 a = grid.evalOn(g, a);
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	66 a = spdiag(a);
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	67 end
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	68
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	69 if isa(b, 'function_handle')
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	70 b = grid.evalOn(g, b);
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	71 b = spdiag(b);
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	72 end
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	73
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	74 % If b is scalar
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	75 if length(b) == 1
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	76 b = b*speye(g.N(), g.N());
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	77 end
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	78
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	79 dim = 2;
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	80 m = g.size();
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	81 m_u = m(1);
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	82 m_v = m(2);
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	83 m_tot = g.N();
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	84
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	85 if isa(g, 'grid.Curvilinear')
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	86 h = g.scaling();
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	87 h_u = h(1);
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	88 h_v = h(2);
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	89 else
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	90 h_u = 1/(m_u - 1);
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	91 h_v = 1/(m_v - 1);
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	92 h = [h_u, h_v];
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	93 end
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	94
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	95 % 1D operators
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	96 ops_u = opSet(m_u, {0, 1}, order);
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	97 ops_v = opSet(m_v, {0, 1}, order);
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	98
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	99 I_u = speye(m_u);
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	100 I_v = speye(m_v);
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	101
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	102 D1_u = ops_u.D1;
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	103 D2_u = ops_u.D2;
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	104 H_u = ops_u.H;
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	105 Hi_u = ops_u.HI;
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	106 e_l_u = ops_u.e_l;
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	107 e_r_u = ops_u.e_r;
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	108 d1_l_u = ops_u.d1_l;
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	109 d1_r_u = ops_u.d1_r;
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	110
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	111 D1_v = ops_v.D1;
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	112 D2_v = ops_v.D2;
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	113 H_v = ops_v.H;
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	114 Hi_v = ops_v.HI;
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	115 e_l_v = ops_v.e_l;
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	116 e_r_v = ops_v.e_r;
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	117 d1_l_v = ops_v.d1_l;
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	118 d1_r_v = ops_v.d1_r;
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	119
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	120
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	121 % Logical operators
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	122 Du = kr(D1_u,I_v);
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	123 Dv = kr(I_u,D1_v);
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	124 obj.Hu = kr(H_u,I_v);
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	125 obj.Hv = kr(I_u,H_v);
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	126 obj.Hiu = kr(Hi_u,I_v);
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	127 obj.Hiv = kr(I_u,Hi_v);
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	128
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	129 e_w = kr(e_l_u,I_v);
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	130 e_e = kr(e_r_u,I_v);
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	131 e_s = kr(I_u,e_l_v);
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	132 e_n = kr(I_u,e_r_v);
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	133 obj.du_w = kr(d1_l_u,I_v);
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	134 obj.dv_w = (e_w'*Dv)';
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	135 obj.du_e = kr(d1_r_u,I_v);
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	136 obj.dv_e = (e_e'*Dv)';
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	137 obj.du_s = (e_s'*Du)';
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	138 obj.dv_s = kr(I_u,d1_l_v);
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	139 obj.du_n = (e_n'*Du)';
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	140 obj.dv_n = kr(I_u,d1_r_v);
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	141
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	142
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	143 % Metric coefficients
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	144 coords = g.points();
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	145 x = coords(:,1);
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	146 y = coords(:,2);
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	147
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	148 x_u = Du*x;
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	149 x_v = Dv*x;
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	150 y_u = Du*y;
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	151 y_v = Dv*y;
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	152
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	153 J = x_u.y_v - x_v.y_u;
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	154 a11 = 1./J .* (x_v.^2 + y_v.^2);
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	155 a12 = -1./J .* (x_u.x_v + y_u.y_v);
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	156 a22 = 1./J .* (x_u.^2 + y_u.^2);
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	157 lambda = 1/2 * (a11 + a22 - sqrt((a11-a22).^2 + 4*a12.^2));
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	158
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	159 K = cell(dim, dim);
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	160 K{1,1} = spdiag(y_v./J);
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	161 K{1,2} = spdiag(-y_u./J);
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	162 K{2,1} = spdiag(-x_v./J);
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	163 K{2,2} = spdiag(x_u./J);
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	164 obj.K = K;
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	165
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	166 obj.x_u = x_u;
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	167 obj.x_v = x_v;
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	168 obj.y_u = y_u;
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	169 obj.y_v = y_v;
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	170
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	171 % Assemble full operators
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	172 L_12 = spdiag(a12);
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	173 Duv = DubL_12*Dv;
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	174 Dvu = DvbL_12*Du;
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	175
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	176 Duu = sparse(m_tot);
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	177 Dvv = sparse(m_tot);
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	178 ind = grid.funcToMatrix(g, 1:m_tot);
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	179
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	180 for i = 1:m_v
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	181 b_a11 = b*a11;
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	182 D = D2_u(b_a11(ind(:,i)));
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	183 p = ind(:,i);
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	184 Duu(p,p) = D;
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	185 end
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	186
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	187 for i = 1:m_u
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	188 b_a22 = b*a22;
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	189 D = D2_v(b_a22(ind(i,:)));
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	190 p = ind(i,:);
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	191 Dvv(p,p) = D;
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	192 end
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	193
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	194
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	195 % Physical operators
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	196 obj.J = spdiag(J);
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	197 obj.Ji = spdiag(1./J);
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	198
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	199 obj.D = obj.Jia(Duu + Duv + Dvu + Dvv);
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	200 obj.H = obj.J*kr(H_u,H_v);
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	201 obj.Hi = obj.Ji*kr(Hi_u,Hi_v);
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	202
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	203 obj.e_w = e_w;
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	204 obj.e_e = e_e;
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	205 obj.e_s = e_s;
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	206 obj.e_n = e_n;
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	207
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	208 %% normal derivatives
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	209 I_w = ind(1,:);
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	210 I_e = ind(end,:);
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	211 I_s = ind(:,1);
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	212 I_n = ind(:,end);
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	213
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	214 a11_w = spdiag(a11(I_w));
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	215 a12_w = spdiag(a12(I_w));
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	216 a11_e = spdiag(a11(I_e));
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	217 a12_e = spdiag(a12(I_e));
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	218 a22_s = spdiag(a22(I_s));
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	219 a12_s = spdiag(a12(I_s));
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	220 a22_n = spdiag(a22(I_n));
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	221 a12_n = spdiag(a12(I_n));
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	222
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	223 s_w = sqrt((e_w'x_v).^2 + (e_w'y_v).^2);
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	224 s_e = sqrt((e_e'x_v).^2 + (e_e'y_v).^2);
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	225 s_s = sqrt((e_s'x_u).^2 + (e_s'y_u).^2);
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	226 s_n = sqrt((e_n'x_u).^2 + (e_n'y_u).^2);
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	227
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	228 obj.d_w = -1(spdiag(1./s_w)(a11_wobj.du_w' + a12_wobj.dv_w'))';
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	229 obj.d_e = (spdiag(1./s_e)(a11_eobj.du_e' + a12_e*obj.dv_e'))';
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	230 obj.d_s = -1(spdiag(1./s_s)(a22_sobj.dv_s' + a12_sobj.du_s'))';
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	231 obj.d_n = (spdiag(1./s_n)(a22_nobj.dv_n' + a12_n*obj.du_n'))';
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	232
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	233 obj.Dx = spdiag( y_v./J)Du + spdiag(-y_u./J)Dv;
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	234 obj.Dy = spdiag(-x_v./J)Du + spdiag( x_u./J)Dv;
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	235
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	236 %% Boundary inner products
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	237 obj.H_w = H_v*spdiag(s_w);
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	238 obj.H_e = H_v*spdiag(s_e);
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	239 obj.H_s = H_u*spdiag(s_s);
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	240 obj.H_n = H_u*spdiag(s_n);
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	241
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	242 % Misc.
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	243 obj.m = m;
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	244 obj.h = [h_u h_v];
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	245 obj.order = order;
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	246 obj.grid = g;
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	247 obj.dim = dim;
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	248
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	249 obj.a = a;
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	250 obj.b = b;
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	251 obj.a11 = a11;
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	252 obj.a12 = a12;
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	253 obj.a22 = a22;
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	254 obj.s_w = spdiag(s_w);
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	255 obj.s_e = spdiag(s_e);
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	256 obj.s_s = spdiag(s_s);
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	257 obj.s_n = spdiag(s_n);
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	258
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	259 obj.theta_M_u = h_u*ops_u.borrowing.M.d1;
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	260 obj.theta_M_v = h_v*ops_v.borrowing.M.d1;
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	261
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	262 obj.theta_R_u = h_u*ops_u.borrowing.R.delta_D;
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	263 obj.theta_R_v = h_v*ops_v.borrowing.R.delta_D;
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	264
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	265 obj.theta_H_u = h_u*ops_u.borrowing.H11;
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	266 obj.theta_H_v = h_v*ops_v.borrowing.H11;
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	267
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	268 % Temporary
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	269 obj.lambda = lambda;
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	270 obj.gamm_u = h_u*ops_u.borrowing.M.d1;
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	271 obj.gamm_v = h_v*ops_v.borrowing.M.d1;
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	272 end
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	273
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	274
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	275 % Closure functions return the opertors applied to the own doamin to close the boundary
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	276 % Penalty functions return the opertors to force the solution. In the case of an interface it returns the operator applied to the other doamin.
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	277 % boundary is a string specifying the boundary e.g. 'l','r' or 'e','w','n','s'.
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	278 % type is a string specifying the type of boundary condition if there are several.
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	279 % data is a function returning the data that should be applied at the boundary.
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	280 % neighbour_scheme is an instance of Scheme that should be interfaced to.
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	281 % neighbour_boundary is a string specifying which boundary to interface to.
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	282 function [closure, penalty] = boundary_condition(obj, boundary, type, parameter)
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	283 default_arg('type','neumann');
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	284 default_arg('parameter', []);
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	285
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	286 e = obj.getBoundaryOperator('e', boundary);
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	287 d = obj.getBoundaryOperator('d', boundary);
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	288 H_b = obj.getBoundaryQuadrature(boundary);
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	289 s_b = obj.getBoundaryScaling(boundary);
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	290 [th_H, ~, th_R] = obj.getBoundaryBorrowing(boundary);
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	291 m = obj.getBoundaryNumber(boundary);
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	292
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	293 K = obj.K;
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	294 J = obj.J;
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	295 Hi = obj.Hi;
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	296 a = obj.a;
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	297 b_b = e'obj.be;
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	298
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	299 switch type
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	300 % Dirichlet boundary condition
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	301 case {'D','d','dirichlet'}
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	302 tuning = 1.0;
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	303
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	304 sigma = 0;
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	305 for i = 1:obj.dim
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	306 sigma = sigma + e'JK{i,m}K{i,m}e;
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	307 end
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	308 sigma = sigma/s_b;
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	309 % tau = tuning(1/th_R + obj.dim/th_H)sigma;
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	310
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	311 tau_R = 1/th_R*sigma;
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	312
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	313 tau_H = 1/th_H*sigma;
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	314 tau_H(1,1) = obj.dim*tau_H(1,1);
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	315 tau_H(end,end) = obj.dim*tau_H(end,end);
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	316
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	317 tau = tuning*(tau_R + tau_H);
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	318
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	319 closure = aHidb_bH_b*e' ...
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	320 -aHietaub_bH_be';
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	321
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	322 penalty = -aHidb_bH_b ...
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	323 +aHietaub_b*H_b;
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	324
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	325
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	326 % Neumann boundary condition. Note that the penalty is for du/dn and not b*du/dn.
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	327 case {'N','n','neumann'}
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	328 tau1 = -1;
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	329 tau2 = 0;
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	330 tau = (tau1e + tau2d)*H_b;
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	331
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	332 closure = aHitaub_bd';
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	333 penalty = -aHitau*b_b;
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	334
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	335
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	336 % Unknown boundary condition
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	337 otherwise
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	338 error('No such boundary condition: type = %s',type);
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	339 end
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	340 end
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	341
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	342 % type Struct that specifies the interface coupling.
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	343 % Fields:
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	344 % -- tuning: penalty strength, defaults to 1.2
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	345 % -- interpolation: type of interpolation, default 'none'
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	346 function [closure, penalty] = interface(obj,boundary,neighbour_scheme,neighbour_boundary,type)
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	347
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	348 % error('Not implemented')
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	349
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	350 defaultType.tuning = 1.0;
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	351 defaultType.interpolation = 'none';
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	352 default_struct('type', defaultType);
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	353
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	354 switch type.interpolation
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	355 case {'none', ''}
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	356 [closure, penalty] = interfaceStandard(obj,boundary,neighbour_scheme,neighbour_boundary,type);
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	357 case {'op','OP'}
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	358 [closure, penalty] = interfaceNonConforming(obj,boundary,neighbour_scheme,neighbour_boundary,type);
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	359 otherwise
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	360 error('Unknown type of interpolation: %s ', type.interpolation);
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	361 end
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	362 end
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	363
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	364 function [closure, penalty] = interfaceStandard(obj,boundary,neighbour_scheme,neighbour_boundary,type)
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	365 tuning = type.tuning;
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	366
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	367 dim = obj.dim;
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	368 % u denotes the solution in the own domain
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	369 % v denotes the solution in the neighbour domain
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	370 u = obj;
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	371 v = neighbour_scheme;
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	372
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	373 % Boundary operators, u
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	374 e_u = u.getBoundaryOperator('e', boundary);
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	375 d_u = u.getBoundaryOperator('d', boundary);
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	376 gamm_u = u.getBoundaryBorrowing(boundary);
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	377 s_b_u = u.getBoundaryScaling(boundary);
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	378 [th_H_u, ~, th_R_u] = u.getBoundaryBorrowing(boundary);
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	379 m_u = u.getBoundaryNumber(boundary);
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	380
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	381 % Coefficients, u
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	382 K_u = u.K;
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	383 J_u = u.J;
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	384 b_b_u = e_u'u.be_u;
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	385
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	386 % Boundary operators, v
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	387 e_v = v.getBoundaryOperator('e', neighbour_boundary);
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	388 d_v = v.getBoundaryOperator('d', neighbour_boundary);
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	389 gamm_v = v.getBoundaryBorrowing(neighbour_boundary);
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	390 s_b_v = v.getBoundaryScaling(neighbour_boundary);
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	391 [th_H_v, ~, th_R_v] = v.getBoundaryBorrowing(neighbour_boundary);
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	392 m_v = v.getBoundaryNumber(neighbour_boundary);
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	393
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	394 % Coefficients, v
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	395 K_v = v.K;
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	396 J_v = v.J;
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	397 b_b_v = e_v'v.be_v;
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	398
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	399 %--- Penalty strength tau -------------
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	400 sigma_u = 0;
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	401 sigma_v = 0;
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	402 for i = 1:obj.dim
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	403 sigma_u = sigma_u + e_u'J_uK_u{i,m_u}K_u{i,m_u}e_u;
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	404 sigma_v = sigma_v + e_v'J_vK_v{i,m_v}K_v{i,m_v}e_v;
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	405 end
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	406 sigma_u = sigma_u/s_b_u;
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	407 sigma_v = sigma_v/s_b_v;
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	408
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	409 tau_R_u = 1/th_R_u*sigma_u;
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	410 tau_R_v = 1/th_R_v*sigma_v;
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	411
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	412 tau_H_u = 1/th_H_u*sigma_u;
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	413 tau_H_u(1,1) = dim*tau_H_u(1,1);
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	414 tau_H_u(end,end) = dim*tau_H_u(end,end);
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	415
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	416 tau_H_v = 1/th_H_v*sigma_v;
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	417 tau_H_v(1,1) = dim*tau_H_v(1,1);
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	418 tau_H_v(end,end) = dim*tau_H_v(end,end);
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	419
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	420 tau = 1/4tuning(b_b_u(tau_R_u + tau_H_u) + b_b_v(tau_R_v + tau_H_v));
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	421 %--------------------------------------
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	422
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	423 % Operators/coefficients that are only required from this side
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	424 Hi = u.Hi;
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	425 H_b = u.getBoundaryQuadrature(boundary);
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	426 a = u.a;
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	427
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	428 closure = 1/2aHid_ub_b_uH_be_u' ...
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	429 -1/2aHie_uH_bb_b_ud_u' ...
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	430 -aHie_utauH_b*e_u';
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	431
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	432 penalty = -1/2aHid_ub_b_uH_be_v' ...
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	433 -1/2aHie_uH_bb_b_vd_v' ...
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	434 +aHie_utauH_b*e_v';
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	435 end
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	436
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	437 function [closure, penalty] = interfaceNonConforming(obj,boundary,neighbour_scheme,neighbour_boundary,type)
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	438
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	439 % TODO: Make this work for curvilinear grids
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	440 warning('LaplaceCurvilinear: Non-conforming grid interpolation only works for Cartesian grids.');
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	441 warning('LaplaceCurvilinear: Non-conforming interface uses Virtas penalty strength');
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	442 warning('LaplaceCurvilinear: Non-conforming interface assumes that b is constant');
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	443
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	444 % User can request special interpolation operators by specifying type.interpOpSet
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	445 default_field(type, 'interpOpSet', @sbp.InterpOpsOP);
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	446 interpOpSet = type.interpOpSet;
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	447 tuning = type.tuning;
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	448
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	449
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	450 % u denotes the solution in the own domain
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	451 % v denotes the solution in the neighbour domain
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	452 e_u = obj.getBoundaryOperator('e', boundary);
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	453 d_u = obj.getBoundaryOperator('d', boundary);
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	454 H_b_u = obj.getBoundaryQuadrature(boundary);
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	455 I_u = obj.getBoundaryIndices(boundary);
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	456 gamm_u = obj.getBoundaryBorrowing(boundary);
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	457
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	458 e_v = neighbour_scheme.getBoundaryOperator('e', neighbour_boundary);
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	459 d_v = neighbour_scheme.getBoundaryOperator('d', neighbour_boundary);
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	460 H_b_v = neighbour_scheme.getBoundaryQuadrature(neighbour_boundary);
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	461 I_v = neighbour_scheme.getBoundaryIndices(neighbour_boundary);
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	462 gamm_v = neighbour_scheme.getBoundaryBorrowing(neighbour_boundary);
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	463
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	464
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	465 % Find the number of grid points along the interface
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	466 m_u = size(e_u, 2);
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	467 m_v = size(e_v, 2);
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	468
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	469 Hi = obj.Hi;
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	470 a = obj.a;
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	471
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	472 u = obj;
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	473 v = neighbour_scheme;
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	474
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	475 b1_u = gamm_u*u.lambda(I_u)./u.a11(I_u).^2;
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	476 b2_u = gamm_u*u.lambda(I_u)./u.a22(I_u).^2;
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	477 b1_v = gamm_v*v.lambda(I_v)./v.a11(I_v).^2;
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	478 b2_v = gamm_v*v.lambda(I_v)./v.a22(I_v).^2;
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	479
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	480 tau_u = -1./(4b1_u) -1./(4b2_u);
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	481 tau_v = -1./(4b1_v) -1./(4b2_v);
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	482
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	483 tau_u = tuning * spdiag(tau_u);
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	484 tau_v = tuning * spdiag(tau_v);
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	485 beta_u = tau_v;
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	486
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	487 % Build interpolation operators
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	488 intOps = interpOpSet(m_u, m_v, obj.order, neighbour_scheme.order);
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	489 Iu2v = intOps.Iu2v;
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	490 Iv2u = intOps.Iv2u;
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	491
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	492 closure = aHie_utau_uH_b_u*e_u' + ...
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	493 aHie_uH_b_uIv2u.badbeta_uIu2v.good*e_u' + ...
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	494 a1/2Hid_uH_b_u*e_u' + ...
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	495 -a1/2Hie_uH_b_u*d_u';
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	496
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	497 penalty = -aHie_utau_uH_b_uIv2u.goode_v' + ...
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	498 -aHie_uH_b_uIv2u.badbeta_ue_v' + ...
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	499 -a1/2Hid_uH_b_uIv2u.goode_v' + ...
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	500 -a1/2Hie_uH_b_uIv2u.badd_v';
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	501
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	502 end
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	503
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	504 % Returns the boundary operator op for the boundary specified by the string boundary.
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	505 % op -- string
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	506 % boundary -- string
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	507 function o = getBoundaryOperator(obj, op, boundary)
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	508 assertIsMember(op, {'e', 'd'})
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	509 assertIsMember(boundary, {'w', 'e', 's', 'n'})
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	510
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	511 o = obj.([op, '_', boundary]);
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	512 end
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	513
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	514 % Returns square boundary quadrature matrix, of dimension
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	515 % corresponding to the number of boundary points
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	516 %
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	517 % boundary -- string
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	518 function H_b = getBoundaryQuadrature(obj, boundary)
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	519 assertIsMember(boundary, {'w', 'e', 's', 'n'})
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	520
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	521 H_b = obj.(['H_', boundary]);
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	522 end
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	523
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	524 % Returns square boundary quadrature scaling matrix, of dimension
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	525 % corresponding to the number of boundary points
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	526 %
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	527 % boundary -- string
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	528 function s_b = getBoundaryScaling(obj, boundary)
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	529 assertIsMember(boundary, {'w', 'e', 's', 'n'})
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	530
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	531 s_b = obj.(['s_', boundary]);
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	532 end
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	533
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	534 % Returns the coordinate number corresponding to the boundary
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	535 %
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	536 % boundary -- string
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	537 function m = getBoundaryNumber(obj, boundary)
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	538 assertIsMember(boundary, {'w', 'e', 's', 'n'})
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	539
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	540 switch boundary
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	541 case {'w', 'e'}
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	542 m = 1;
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	543 case {'s', 'n'}
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	544 m = 2;
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	545 end
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	546 end
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	547
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	548 % Returns the indices of the boundary points in the grid matrix
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	549 % boundary -- string
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	550 function I = getBoundaryIndices(obj, boundary)
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	551 ind = grid.funcToMatrix(obj.grid, 1:prod(obj.m));
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	552 switch boundary
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	553 case 'w'
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	554 I = ind(1,:);
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	555 case 'e'
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	556 I = ind(end,:);
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	557 case 's'
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	558 I = ind(:,1)';
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	559 case 'n'
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	560 I = ind(:,end)';
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	561 otherwise
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	562 error('No such boundary: boundary = %s',boundary);
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	563 end
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	564 end
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	565
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	566 % Returns borrowing constant gamma
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	567 % boundary -- string
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	568 function [theta_H, theta_M, theta_R] = getBoundaryBorrowing(obj, boundary)
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	569 switch boundary
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	570 case {'w','e'}
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	571 theta_H = obj.theta_H_u;
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	572 theta_M = obj.theta_M_u;
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	573 theta_R = obj.theta_R_u;
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	574 case {'s','n'}
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	575 theta_H = obj.theta_H_v;
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	576 theta_M = obj.theta_M_v;
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	577 theta_R = obj.theta_R_v;
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	578 otherwise
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	579 error('No such boundary: boundary = %s',boundary);
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	580 end
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	581 end
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	582
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	583 function N = size(obj)
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	584 N = prod(obj.m);
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	585 end
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	586 end
b828e589d540 Copy-paste LaplCurviNewCorner from laplace branch. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	587 end

Mercurial > repos > public > sbplib

annotate +scheme/LaplaceCurvilinearNewCorner.m @ 1060:e40899094f20 feature/poroelastic