sbplib: +scheme/Elastic2dVariable.m annotate

annotate +scheme/Elastic2dVariable.m @ 972:104f0af001e0 feature/poroelastic

Clean up Elastic2dVariable.interfaceStandard

author	Martin Almquist <malmquist@stanford.edu>
date	Tue, 25 Dec 2018 08:33:35 +0100
parents	23d9ca6755be
children	f5a99cca4654

rev	line source
687 e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	1 classdef Elastic2dVariable < scheme.Scheme
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	2
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	3 % Discretizes the elastic wave equation:
813 b374a8aa9246 Correct interface penalty strength in Elastic2dVariable Martin Almquist <malmquist@stanford.edu> parents: 795 diff changeset	4 % rho u_{i,tt} = di lambda dj u_j + dj mu di u_j + dj mu dj u_i
687 e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	5 % opSet should be cell array of opSets, one per dimension. This
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	6 % is useful if we have periodic BC in one direction.
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	7
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	8 properties
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	9 m % Number of points in each direction, possibly a vector
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	10 h % Grid spacing
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	11
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	12 grid
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	13 dim
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	14
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	15 order % Order of accuracy for the approximation
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	16
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	17 % Diagonal matrices for varible coefficients
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	18 LAMBDA % Variable coefficient, related to dilation
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	19 MU % Shear modulus, variable coefficient
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	20 RHO, RHOi % Density, variable
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	21
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	22 D % Total operator
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	23 D1 % First derivatives
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	24
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	25 % Second derivatives
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	26 D2_lambda
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	27 D2_mu
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	28
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	29 % Traction operators used for BC
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	30 T_l, T_r
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	31 tau_l, tau_r
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	32
855 5751262b323b Add 1D quadrature matrices as property in Elastic2dVariable. Martin Almquist <malmquist@stanford.edu> parents: 795 diff changeset	33 H, Hi, H_1D % Inner products
687 e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	34 e_l, e_r
963 c75ddd568fcc Turn alpha into a boundary operator. Add properties H_w etc for getBoundaryQuadrature to work. Martin Almquist <malmquist@stanford.edu> parents: 961 diff changeset	35 e_w, e_e, e_s, e_n
813 b374a8aa9246 Correct interface penalty strength in Elastic2dVariable Martin Almquist <malmquist@stanford.edu> parents: 795 diff changeset	36
b374a8aa9246 Correct interface penalty strength in Elastic2dVariable Martin Almquist <malmquist@stanford.edu> parents: 795 diff changeset	37
687 e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	38 d1_l, d1_r % Normal derivatives at the boundary
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	39 E % E{i}^T picks out component i
813 b374a8aa9246 Correct interface penalty strength in Elastic2dVariable Martin Almquist <malmquist@stanford.edu> parents: 795 diff changeset	40
687 e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	41 H_boundary % Boundary inner products
963 c75ddd568fcc Turn alpha into a boundary operator. Add properties H_w etc for getBoundaryQuadrature to work. Martin Almquist <malmquist@stanford.edu> parents: 961 diff changeset	42 H_w, H_e, H_s, H_n
687 e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	43
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	44 % Kroneckered norms and coefficients
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	45 RHOi_kron
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	46 Hi_kron
861 607c631f175e Small changes to Elastic2dVariable to facilitate adjoing gradient computation. Martin Almquist <malmquist@stanford.edu> parents: 855 diff changeset	47
882 14fee299ada2 In Elastic2dVariable: Improve notation for borrowing constants. Update interface constant to match corrected derivation. Tests ok. Martin Almquist <malmquist@stanford.edu> parents: 861 diff changeset	48 % Borrowing constants of the form gamma*h, where gamma is a dimensionless constant.
14fee299ada2 In Elastic2dVariable: Improve notation for borrowing constants. Update interface constant to match corrected derivation. Tests ok. Martin Almquist <malmquist@stanford.edu> parents: 861 diff changeset	49 theta_R % Borrowing (d1- D1)^2 from R
14fee299ada2 In Elastic2dVariable: Improve notation for borrowing constants. Update interface constant to match corrected derivation. Tests ok. Martin Almquist <malmquist@stanford.edu> parents: 861 diff changeset	50 theta_H % First entry in norm matrix
14fee299ada2 In Elastic2dVariable: Improve notation for borrowing constants. Update interface constant to match corrected derivation. Tests ok. Martin Almquist <malmquist@stanford.edu> parents: 861 diff changeset	51 theta_M % Borrowing d1^2 from M.
14fee299ada2 In Elastic2dVariable: Improve notation for borrowing constants. Update interface constant to match corrected derivation. Tests ok. Martin Almquist <malmquist@stanford.edu> parents: 861 diff changeset	52
861 607c631f175e Small changes to Elastic2dVariable to facilitate adjoing gradient computation. Martin Almquist <malmquist@stanford.edu> parents: 855 diff changeset	53 % Structures used for adjoint optimization
607c631f175e Small changes to Elastic2dVariable to facilitate adjoing gradient computation. Martin Almquist <malmquist@stanford.edu> parents: 855 diff changeset	54 B
687 e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	55 end
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	56
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	57 methods
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	58
861 607c631f175e Small changes to Elastic2dVariable to facilitate adjoing gradient computation. Martin Almquist <malmquist@stanford.edu> parents: 855 diff changeset	59 % The coefficients can either be function handles or grid functions
607c631f175e Small changes to Elastic2dVariable to facilitate adjoing gradient computation. Martin Almquist <malmquist@stanford.edu> parents: 855 diff changeset	60 function obj = Elastic2dVariable(g ,order, lambda, mu, rho, opSet)
687 e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	61 default_arg('opSet',{@sbp.D2Variable, @sbp.D2Variable});
861 607c631f175e Small changes to Elastic2dVariable to facilitate adjoing gradient computation. Martin Almquist <malmquist@stanford.edu> parents: 855 diff changeset	62 default_arg('lambda', @(x,y) 0*x+1);
607c631f175e Small changes to Elastic2dVariable to facilitate adjoing gradient computation. Martin Almquist <malmquist@stanford.edu> parents: 855 diff changeset	63 default_arg('mu', @(x,y) 0*x+1);
607c631f175e Small changes to Elastic2dVariable to facilitate adjoing gradient computation. Martin Almquist <malmquist@stanford.edu> parents: 855 diff changeset	64 default_arg('rho', @(x,y) 0*x+1);
687 e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	65 dim = 2;
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	66
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	67 assert(isa(g, 'grid.Cartesian'))
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	68
861 607c631f175e Small changes to Elastic2dVariable to facilitate adjoing gradient computation. Martin Almquist <malmquist@stanford.edu> parents: 855 diff changeset	69 if isa(lambda, 'function_handle')
607c631f175e Small changes to Elastic2dVariable to facilitate adjoing gradient computation. Martin Almquist <malmquist@stanford.edu> parents: 855 diff changeset	70 lambda = grid.evalOn(g, lambda);
607c631f175e Small changes to Elastic2dVariable to facilitate adjoing gradient computation. Martin Almquist <malmquist@stanford.edu> parents: 855 diff changeset	71 end
607c631f175e Small changes to Elastic2dVariable to facilitate adjoing gradient computation. Martin Almquist <malmquist@stanford.edu> parents: 855 diff changeset	72 if isa(mu, 'function_handle')
607c631f175e Small changes to Elastic2dVariable to facilitate adjoing gradient computation. Martin Almquist <malmquist@stanford.edu> parents: 855 diff changeset	73 mu = grid.evalOn(g, mu);
607c631f175e Small changes to Elastic2dVariable to facilitate adjoing gradient computation. Martin Almquist <malmquist@stanford.edu> parents: 855 diff changeset	74 end
607c631f175e Small changes to Elastic2dVariable to facilitate adjoing gradient computation. Martin Almquist <malmquist@stanford.edu> parents: 855 diff changeset	75 if isa(rho, 'function_handle')
607c631f175e Small changes to Elastic2dVariable to facilitate adjoing gradient computation. Martin Almquist <malmquist@stanford.edu> parents: 855 diff changeset	76 rho = grid.evalOn(g, rho);
607c631f175e Small changes to Elastic2dVariable to facilitate adjoing gradient computation. Martin Almquist <malmquist@stanford.edu> parents: 855 diff changeset	77 end
607c631f175e Small changes to Elastic2dVariable to facilitate adjoing gradient computation. Martin Almquist <malmquist@stanford.edu> parents: 855 diff changeset	78
687 e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	79 m = g.size();
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	80 m_tot = g.N();
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	81
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	82 h = g.scaling();
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	83 lim = g.lim;
729 aa8cf3851de8 Update multiblock.DiffOp to work for systems. Martin Almquist <malmquist@stanford.edu> parents: 727 diff changeset	84 if isempty(lim)
aa8cf3851de8 Update multiblock.DiffOp to work for systems. Martin Almquist <malmquist@stanford.edu> parents: 727 diff changeset	85 x = g.x;
aa8cf3851de8 Update multiblock.DiffOp to work for systems. Martin Almquist <malmquist@stanford.edu> parents: 727 diff changeset	86 lim = cell(length(x),1);
aa8cf3851de8 Update multiblock.DiffOp to work for systems. Martin Almquist <malmquist@stanford.edu> parents: 727 diff changeset	87 for i = 1:length(x)
aa8cf3851de8 Update multiblock.DiffOp to work for systems. Martin Almquist <malmquist@stanford.edu> parents: 727 diff changeset	88 lim{i} = {min(x{i}), max(x{i})};
aa8cf3851de8 Update multiblock.DiffOp to work for systems. Martin Almquist <malmquist@stanford.edu> parents: 727 diff changeset	89 end
aa8cf3851de8 Update multiblock.DiffOp to work for systems. Martin Almquist <malmquist@stanford.edu> parents: 727 diff changeset	90 end
687 e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	91
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	92 % 1D operators
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	93 ops = cell(dim,1);
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	94 for i = 1:dim
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	95 ops{i} = opSet{i}(m(i), lim{i}, order);
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	96 end
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	97
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	98 % Borrowing constants
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	99 for i = 1:dim
882 14fee299ada2 In Elastic2dVariable: Improve notation for borrowing constants. Update interface constant to match corrected derivation. Tests ok. Martin Almquist <malmquist@stanford.edu> parents: 861 diff changeset	100 obj.theta_R{i} = h(i)*ops{i}.borrowing.R.delta_D;
14fee299ada2 In Elastic2dVariable: Improve notation for borrowing constants. Update interface constant to match corrected derivation. Tests ok. Martin Almquist <malmquist@stanford.edu> parents: 861 diff changeset	101 obj.theta_H{i} = h(i)*ops{i}.borrowing.H11;
14fee299ada2 In Elastic2dVariable: Improve notation for borrowing constants. Update interface constant to match corrected derivation. Tests ok. Martin Almquist <malmquist@stanford.edu> parents: 861 diff changeset	102 obj.theta_M{i} = h(i)*ops{i}.borrowing.M.d1;
687 e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	103 end
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	104
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	105 I = cell(dim,1);
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	106 D1 = cell(dim,1);
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	107 D2 = cell(dim,1);
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	108 H = cell(dim,1);
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	109 Hi = cell(dim,1);
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	110 e_l = cell(dim,1);
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	111 e_r = cell(dim,1);
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	112 d1_l = cell(dim,1);
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	113 d1_r = cell(dim,1);
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	114
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	115 for i = 1:dim
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	116 I{i} = speye(m(i));
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	117 D1{i} = ops{i}.D1;
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	118 D2{i} = ops{i}.D2;
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	119 H{i} = ops{i}.H;
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	120 Hi{i} = ops{i}.HI;
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	121 e_l{i} = ops{i}.e_l;
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	122 e_r{i} = ops{i}.e_r;
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	123 d1_l{i} = ops{i}.d1_l;
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	124 d1_r{i} = ops{i}.d1_r;
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	125 end
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	126
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	127 %====== Assemble full operators ========
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	128 LAMBDA = spdiag(lambda);
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	129 obj.LAMBDA = LAMBDA;
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	130 MU = spdiag(mu);
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	131 obj.MU = MU;
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	132 RHO = spdiag(rho);
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	133 obj.RHO = RHO;
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	134 obj.RHOi = inv(RHO);
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	135
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	136 obj.D1 = cell(dim,1);
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	137 obj.D2_lambda = cell(dim,1);
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	138 obj.D2_mu = cell(dim,1);
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	139 obj.e_l = cell(dim,1);
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	140 obj.e_r = cell(dim,1);
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	141 obj.d1_l = cell(dim,1);
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	142 obj.d1_r = cell(dim,1);
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	143
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	144 % D1
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	145 obj.D1{1} = kron(D1{1},I{2});
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	146 obj.D1{2} = kron(I{1},D1{2});
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	147
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	148 % Boundary operators
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	149 obj.e_l{1} = kron(e_l{1},I{2});
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	150 obj.e_l{2} = kron(I{1},e_l{2});
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	151 obj.e_r{1} = kron(e_r{1},I{2});
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	152 obj.e_r{2} = kron(I{1},e_r{2});
963 c75ddd568fcc Turn alpha into a boundary operator. Add properties H_w etc for getBoundaryQuadrature to work. Martin Almquist <malmquist@stanford.edu> parents: 961 diff changeset	153 obj.e_w = obj.e_l{1};
c75ddd568fcc Turn alpha into a boundary operator. Add properties H_w etc for getBoundaryQuadrature to work. Martin Almquist <malmquist@stanford.edu> parents: 961 diff changeset	154 obj.e_e = obj.e_r{1};
c75ddd568fcc Turn alpha into a boundary operator. Add properties H_w etc for getBoundaryQuadrature to work. Martin Almquist <malmquist@stanford.edu> parents: 961 diff changeset	155 obj.e_s = obj.e_l{2};
c75ddd568fcc Turn alpha into a boundary operator. Add properties H_w etc for getBoundaryQuadrature to work. Martin Almquist <malmquist@stanford.edu> parents: 961 diff changeset	156 obj.e_n = obj.e_r{2};
687 e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	157
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	158 obj.d1_l{1} = kron(d1_l{1},I{2});
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	159 obj.d1_l{2} = kron(I{1},d1_l{2});
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	160 obj.d1_r{1} = kron(d1_r{1},I{2});
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	161 obj.d1_r{2} = kron(I{1},d1_r{2});
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	162
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	163 % D2
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	164 for i = 1:dim
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	165 obj.D2_lambda{i} = sparse(m_tot);
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	166 obj.D2_mu{i} = sparse(m_tot);
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	167 end
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	168 ind = grid.funcToMatrix(g, 1:m_tot);
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	169
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	170 for i = 1:m(2)
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	171 D_lambda = D2{1}(lambda(ind(:,i)));
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	172 D_mu = D2{1}(mu(ind(:,i)));
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	173
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	174 p = ind(:,i);
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	175 obj.D2_lambda{1}(p,p) = D_lambda;
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	176 obj.D2_mu{1}(p,p) = D_mu;
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	177 end
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	178
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	179 for i = 1:m(1)
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	180 D_lambda = D2{2}(lambda(ind(i,:)));
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	181 D_mu = D2{2}(mu(ind(i,:)));
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	182
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	183 p = ind(i,:);
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	184 obj.D2_lambda{2}(p,p) = D_lambda;
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	185 obj.D2_mu{2}(p,p) = D_mu;
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	186 end
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	187
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	188 % Quadratures
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	189 obj.H = kron(H{1},H{2});
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	190 obj.Hi = inv(obj.H);
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	191 obj.H_boundary = cell(dim,1);
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	192 obj.H_boundary{1} = H{2};
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	193 obj.H_boundary{2} = H{1};
855 5751262b323b Add 1D quadrature matrices as property in Elastic2dVariable. Martin Almquist <malmquist@stanford.edu> parents: 795 diff changeset	194 obj.H_1D = {H{1}, H{2}};
963 c75ddd568fcc Turn alpha into a boundary operator. Add properties H_w etc for getBoundaryQuadrature to work. Martin Almquist <malmquist@stanford.edu> parents: 961 diff changeset	195 obj.H_w = H{2};
c75ddd568fcc Turn alpha into a boundary operator. Add properties H_w etc for getBoundaryQuadrature to work. Martin Almquist <malmquist@stanford.edu> parents: 961 diff changeset	196 obj.H_e = H{2};
c75ddd568fcc Turn alpha into a boundary operator. Add properties H_w etc for getBoundaryQuadrature to work. Martin Almquist <malmquist@stanford.edu> parents: 961 diff changeset	197 obj.H_s = H{1};
c75ddd568fcc Turn alpha into a boundary operator. Add properties H_w etc for getBoundaryQuadrature to work. Martin Almquist <malmquist@stanford.edu> parents: 961 diff changeset	198 obj.H_n = H{1};
687 e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	199
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	200 % E{i}^T picks out component i.
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	201 E = cell(dim,1);
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	202 I = speye(m_tot,m_tot);
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	203 for i = 1:dim
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	204 e = sparse(dim,1);
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	205 e(i) = 1;
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	206 E{i} = kron(I,e);
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	207 end
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	208 obj.E = E;
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	209
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	210 % Differentiation matrix D (without SAT)
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	211 D2_lambda = obj.D2_lambda;
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	212 D2_mu = obj.D2_mu;
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	213 D1 = obj.D1;
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	214 D = sparse(dimm_tot,dimm_tot);
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	215 d = @kroneckerDelta; % Kronecker delta
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	216 db = @(i,j) 1-d(i,j); % Logical not of Kronecker delta
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	217 for i = 1:dim
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	218 for j = 1:dim
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	219 D = D + E{i}inv(RHO)( d(i,j)D2_lambda{i}E{j}' +...
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	220 db(i,j)D1{i}LAMBDAD1{j}E{j}' ...
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	221 );
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	222 D = D + E{i}inv(RHO)( d(i,j)D2_mu{i}E{j}' +...
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	223 db(i,j)D1{j}MUD1{i}E{j}' + ...
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	224 D2_mu{j}*E{i}' ...
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	225 );
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	226 end
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	227 end
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	228 obj.D = D;
861 607c631f175e Small changes to Elastic2dVariable to facilitate adjoing gradient computation. Martin Almquist <malmquist@stanford.edu> parents: 855 diff changeset	229 %=========================================%'
687 e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	230
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	231 % Numerical traction operators for BC.
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	232 % Because d1 =/= e0^T*D1, the numerical tractions are different
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	233 % at every boundary.
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	234 T_l = cell(dim,1);
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	235 T_r = cell(dim,1);
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	236 tau_l = cell(dim,1);
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	237 tau_r = cell(dim,1);
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	238 % tau^{j}_i = sum_k T^{j}_{ik} u_k
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	239
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	240 d1_l = obj.d1_l;
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	241 d1_r = obj.d1_r;
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	242 e_l = obj.e_l;
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	243 e_r = obj.e_r;
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	244 D1 = obj.D1;
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	245
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	246 % Loop over boundaries
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	247 for j = 1:dim
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	248 T_l{j} = cell(dim,dim);
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	249 T_r{j} = cell(dim,dim);
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	250 tau_l{j} = cell(dim,1);
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	251 tau_r{j} = cell(dim,1);
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	252
958 72cd29107a9a Temporary changes in multiblock.DiffOp. Change traction operators in Elastic2dvariable to be true boundary operators. But adjoint FD conv test fails for dirichlet BC so need to debug! Martin Almquist <malmquist@stanford.edu> parents: 919 diff changeset	253 LAMBDA_l = e_l{j}'LAMBDAe_l{j};
72cd29107a9a Temporary changes in multiblock.DiffOp. Change traction operators in Elastic2dvariable to be true boundary operators. But adjoint FD conv test fails for dirichlet BC so need to debug! Martin Almquist <malmquist@stanford.edu> parents: 919 diff changeset	254 LAMBDA_r = e_r{j}'LAMBDAe_r{j};
72cd29107a9a Temporary changes in multiblock.DiffOp. Change traction operators in Elastic2dvariable to be true boundary operators. But adjoint FD conv test fails for dirichlet BC so need to debug! Martin Almquist <malmquist@stanford.edu> parents: 919 diff changeset	255 MU_l = e_l{j}'MUe_l{j};
72cd29107a9a Temporary changes in multiblock.DiffOp. Change traction operators in Elastic2dvariable to be true boundary operators. But adjoint FD conv test fails for dirichlet BC so need to debug! Martin Almquist <malmquist@stanford.edu> parents: 919 diff changeset	256 MU_r = e_r{j}'MUe_r{j};
72cd29107a9a Temporary changes in multiblock.DiffOp. Change traction operators in Elastic2dvariable to be true boundary operators. But adjoint FD conv test fails for dirichlet BC so need to debug! Martin Almquist <malmquist@stanford.edu> parents: 919 diff changeset	257
72cd29107a9a Temporary changes in multiblock.DiffOp. Change traction operators in Elastic2dvariable to be true boundary operators. But adjoint FD conv test fails for dirichlet BC so need to debug! Martin Almquist <malmquist@stanford.edu> parents: 919 diff changeset	258 [~, n_l] = size(e_l{j});
72cd29107a9a Temporary changes in multiblock.DiffOp. Change traction operators in Elastic2dvariable to be true boundary operators. But adjoint FD conv test fails for dirichlet BC so need to debug! Martin Almquist <malmquist@stanford.edu> parents: 919 diff changeset	259 [~, n_r] = size(e_r{j});
72cd29107a9a Temporary changes in multiblock.DiffOp. Change traction operators in Elastic2dvariable to be true boundary operators. But adjoint FD conv test fails for dirichlet BC so need to debug! Martin Almquist <malmquist@stanford.edu> parents: 919 diff changeset	260
687 e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	261 % Loop over components
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	262 for i = 1:dim
958 72cd29107a9a Temporary changes in multiblock.DiffOp. Change traction operators in Elastic2dvariable to be true boundary operators. But adjoint FD conv test fails for dirichlet BC so need to debug! Martin Almquist <malmquist@stanford.edu> parents: 919 diff changeset	263 tau_l{j}{i} = sparse(n_l, dim*m_tot);
72cd29107a9a Temporary changes in multiblock.DiffOp. Change traction operators in Elastic2dvariable to be true boundary operators. But adjoint FD conv test fails for dirichlet BC so need to debug! Martin Almquist <malmquist@stanford.edu> parents: 919 diff changeset	264 tau_r{j}{i} = sparse(n_r, dim*m_tot);
687 e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	265 for k = 1:dim
813 b374a8aa9246 Correct interface penalty strength in Elastic2dVariable Martin Almquist <malmquist@stanford.edu> parents: 795 diff changeset	266 T_l{j}{i,k} = ...
966 4a85cff53107 Fix bug in Elastic2dVariable that caused crash for m_x neq m_y Martin Almquist <malmquist@stanford.edu> parents: 965 diff changeset	267 -d(i,j)LAMBDA_l(d(i,k)d1_l{j}' + db(i,k)e_l{j}'*D1{k})...
4a85cff53107 Fix bug in Elastic2dVariable that caused crash for m_x neq m_y Martin Almquist <malmquist@stanford.edu> parents: 965 diff changeset	268 -d(j,k)MU_l(d(i,j)d1_l{j}' + db(i,j)e_l{j}'*D1{i})...
958 72cd29107a9a Temporary changes in multiblock.DiffOp. Change traction operators in Elastic2dvariable to be true boundary operators. But adjoint FD conv test fails for dirichlet BC so need to debug! Martin Almquist <malmquist@stanford.edu> parents: 919 diff changeset	269 -d(i,k)MU_ld1_l{j}';
687 e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	270
813 b374a8aa9246 Correct interface penalty strength in Elastic2dVariable Martin Almquist <malmquist@stanford.edu> parents: 795 diff changeset	271 T_r{j}{i,k} = ...
966 4a85cff53107 Fix bug in Elastic2dVariable that caused crash for m_x neq m_y Martin Almquist <malmquist@stanford.edu> parents: 965 diff changeset	272 d(i,j)LAMBDA_r(d(i,k)d1_r{j}' + db(i,k)e_r{j}'*D1{k})...
4a85cff53107 Fix bug in Elastic2dVariable that caused crash for m_x neq m_y Martin Almquist <malmquist@stanford.edu> parents: 965 diff changeset	273 +d(j,k)MU_r(d(i,j)d1_r{j}' + db(i,j)e_r{j}'*D1{i})...
958 72cd29107a9a Temporary changes in multiblock.DiffOp. Change traction operators in Elastic2dvariable to be true boundary operators. But adjoint FD conv test fails for dirichlet BC so need to debug! Martin Almquist <malmquist@stanford.edu> parents: 919 diff changeset	274 +d(i,k)MU_rd1_r{j}';
687 e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	275
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	276 tau_l{j}{i} = tau_l{j}{i} + T_l{j}{i,k}*E{k}';
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	277 tau_r{j}{i} = tau_r{j}{i} + T_r{j}{i,k}*E{k}';
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	278 end
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	279
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	280 end
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	281 end
958 72cd29107a9a Temporary changes in multiblock.DiffOp. Change traction operators in Elastic2dvariable to be true boundary operators. But adjoint FD conv test fails for dirichlet BC so need to debug! Martin Almquist <malmquist@stanford.edu> parents: 919 diff changeset	282
72cd29107a9a Temporary changes in multiblock.DiffOp. Change traction operators in Elastic2dvariable to be true boundary operators. But adjoint FD conv test fails for dirichlet BC so need to debug! Martin Almquist <malmquist@stanford.edu> parents: 919 diff changeset	283 % Transpose T and tau to match boundary operator convention
72cd29107a9a Temporary changes in multiblock.DiffOp. Change traction operators in Elastic2dvariable to be true boundary operators. But adjoint FD conv test fails for dirichlet BC so need to debug! Martin Almquist <malmquist@stanford.edu> parents: 919 diff changeset	284 for i = 1:dim
72cd29107a9a Temporary changes in multiblock.DiffOp. Change traction operators in Elastic2dvariable to be true boundary operators. But adjoint FD conv test fails for dirichlet BC so need to debug! Martin Almquist <malmquist@stanford.edu> parents: 919 diff changeset	285 for j = 1:dim
72cd29107a9a Temporary changes in multiblock.DiffOp. Change traction operators in Elastic2dvariable to be true boundary operators. But adjoint FD conv test fails for dirichlet BC so need to debug! Martin Almquist <malmquist@stanford.edu> parents: 919 diff changeset	286 tau_l{i}{j} = transpose(tau_l{i}{j});
72cd29107a9a Temporary changes in multiblock.DiffOp. Change traction operators in Elastic2dvariable to be true boundary operators. But adjoint FD conv test fails for dirichlet BC so need to debug! Martin Almquist <malmquist@stanford.edu> parents: 919 diff changeset	287 tau_r{i}{j} = transpose(tau_r{i}{j});
72cd29107a9a Temporary changes in multiblock.DiffOp. Change traction operators in Elastic2dvariable to be true boundary operators. But adjoint FD conv test fails for dirichlet BC so need to debug! Martin Almquist <malmquist@stanford.edu> parents: 919 diff changeset	288 for k = 1:dim
72cd29107a9a Temporary changes in multiblock.DiffOp. Change traction operators in Elastic2dvariable to be true boundary operators. But adjoint FD conv test fails for dirichlet BC so need to debug! Martin Almquist <malmquist@stanford.edu> parents: 919 diff changeset	289 T_l{i}{j,k} = transpose(T_l{i}{j,k});
72cd29107a9a Temporary changes in multiblock.DiffOp. Change traction operators in Elastic2dvariable to be true boundary operators. But adjoint FD conv test fails for dirichlet BC so need to debug! Martin Almquist <malmquist@stanford.edu> parents: 919 diff changeset	290 T_r{i}{j,k} = transpose(T_r{i}{j,k});
72cd29107a9a Temporary changes in multiblock.DiffOp. Change traction operators in Elastic2dvariable to be true boundary operators. But adjoint FD conv test fails for dirichlet BC so need to debug! Martin Almquist <malmquist@stanford.edu> parents: 919 diff changeset	291 end
72cd29107a9a Temporary changes in multiblock.DiffOp. Change traction operators in Elastic2dvariable to be true boundary operators. But adjoint FD conv test fails for dirichlet BC so need to debug! Martin Almquist <malmquist@stanford.edu> parents: 919 diff changeset	292 end
72cd29107a9a Temporary changes in multiblock.DiffOp. Change traction operators in Elastic2dvariable to be true boundary operators. But adjoint FD conv test fails for dirichlet BC so need to debug! Martin Almquist <malmquist@stanford.edu> parents: 919 diff changeset	293 end
72cd29107a9a Temporary changes in multiblock.DiffOp. Change traction operators in Elastic2dvariable to be true boundary operators. But adjoint FD conv test fails for dirichlet BC so need to debug! Martin Almquist <malmquist@stanford.edu> parents: 919 diff changeset	294
687 e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	295 obj.T_l = T_l;
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	296 obj.T_r = T_r;
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	297 obj.tau_l = tau_l;
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	298 obj.tau_r = tau_r;
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	299
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	300 % Kroneckered norms and coefficients
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	301 I_dim = speye(dim);
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	302 obj.RHOi_kron = kron(obj.RHOi, I_dim);
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	303 obj.Hi_kron = kron(obj.Hi, I_dim);
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	304
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	305 % Misc.
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	306 obj.m = m;
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	307 obj.h = h;
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	308 obj.order = order;
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	309 obj.grid = g;
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	310 obj.dim = dim;
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	311
965 db3411264b96 Remove use of tensor library in B assembly Martin Almquist <malmquist@stanford.edu> parents: 963 diff changeset	312 % B, used for adjoint optimization
db3411264b96 Remove use of tensor library in B assembly Martin Almquist <malmquist@stanford.edu> parents: 963 diff changeset	313 B = cell(dim, 1);
861 607c631f175e Small changes to Elastic2dVariable to facilitate adjoing gradient computation. Martin Almquist <malmquist@stanford.edu> parents: 855 diff changeset	314 for i = 1:dim
965 db3411264b96 Remove use of tensor library in B assembly Martin Almquist <malmquist@stanford.edu> parents: 963 diff changeset	315 B{i} = cell(m_tot, 1);
db3411264b96 Remove use of tensor library in B assembly Martin Almquist <malmquist@stanford.edu> parents: 963 diff changeset	316 end
db3411264b96 Remove use of tensor library in B assembly Martin Almquist <malmquist@stanford.edu> parents: 963 diff changeset	317
db3411264b96 Remove use of tensor library in B assembly Martin Almquist <malmquist@stanford.edu> parents: 963 diff changeset	318 for i = 1:dim
db3411264b96 Remove use of tensor library in B assembly Martin Almquist <malmquist@stanford.edu> parents: 963 diff changeset	319 for j = 1:m_tot
db3411264b96 Remove use of tensor library in B assembly Martin Almquist <malmquist@stanford.edu> parents: 963 diff changeset	320 B{i}{j} = sparse(m_tot, m_tot);
861 607c631f175e Small changes to Elastic2dVariable to facilitate adjoing gradient computation. Martin Almquist <malmquist@stanford.edu> parents: 855 diff changeset	321 end
607c631f175e Small changes to Elastic2dVariable to facilitate adjoing gradient computation. Martin Almquist <malmquist@stanford.edu> parents: 855 diff changeset	322 end
607c631f175e Small changes to Elastic2dVariable to facilitate adjoing gradient computation. Martin Almquist <malmquist@stanford.edu> parents: 855 diff changeset	323
965 db3411264b96 Remove use of tensor library in B assembly Martin Almquist <malmquist@stanford.edu> parents: 963 diff changeset	324 ind = grid.funcToMatrix(g, 1:m_tot);
db3411264b96 Remove use of tensor library in B assembly Martin Almquist <malmquist@stanford.edu> parents: 963 diff changeset	325
db3411264b96 Remove use of tensor library in B assembly Martin Almquist <malmquist@stanford.edu> parents: 963 diff changeset	326 % Direction 1
db3411264b96 Remove use of tensor library in B assembly Martin Almquist <malmquist@stanford.edu> parents: 963 diff changeset	327 for k = 1:m(1)
db3411264b96 Remove use of tensor library in B assembly Martin Almquist <malmquist@stanford.edu> parents: 963 diff changeset	328 c = sparse(m(1),1);
db3411264b96 Remove use of tensor library in B assembly Martin Almquist <malmquist@stanford.edu> parents: 963 diff changeset	329 c(k) = 1;
db3411264b96 Remove use of tensor library in B assembly Martin Almquist <malmquist@stanford.edu> parents: 963 diff changeset	330 [~, B_1D] = ops{1}.D2(c);
db3411264b96 Remove use of tensor library in B assembly Martin Almquist <malmquist@stanford.edu> parents: 963 diff changeset	331 for l = 1:m(2)
db3411264b96 Remove use of tensor library in B assembly Martin Almquist <malmquist@stanford.edu> parents: 963 diff changeset	332 p = ind(:,l);
db3411264b96 Remove use of tensor library in B assembly Martin Almquist <malmquist@stanford.edu> parents: 963 diff changeset	333 B{1}{(k-1)*m(2) + l}(p, p) = B_1D;
db3411264b96 Remove use of tensor library in B assembly Martin Almquist <malmquist@stanford.edu> parents: 963 diff changeset	334 end
db3411264b96 Remove use of tensor library in B assembly Martin Almquist <malmquist@stanford.edu> parents: 963 diff changeset	335 end
db3411264b96 Remove use of tensor library in B assembly Martin Almquist <malmquist@stanford.edu> parents: 963 diff changeset	336
db3411264b96 Remove use of tensor library in B assembly Martin Almquist <malmquist@stanford.edu> parents: 963 diff changeset	337 % Direction 2
db3411264b96 Remove use of tensor library in B assembly Martin Almquist <malmquist@stanford.edu> parents: 963 diff changeset	338 for k = 1:m(2)
db3411264b96 Remove use of tensor library in B assembly Martin Almquist <malmquist@stanford.edu> parents: 963 diff changeset	339 c = sparse(m(2),1);
db3411264b96 Remove use of tensor library in B assembly Martin Almquist <malmquist@stanford.edu> parents: 963 diff changeset	340 c(k) = 1;
db3411264b96 Remove use of tensor library in B assembly Martin Almquist <malmquist@stanford.edu> parents: 963 diff changeset	341 [~, B_1D] = ops{2}.D2(c);
db3411264b96 Remove use of tensor library in B assembly Martin Almquist <malmquist@stanford.edu> parents: 963 diff changeset	342 for l = 1:m(1)
db3411264b96 Remove use of tensor library in B assembly Martin Almquist <malmquist@stanford.edu> parents: 963 diff changeset	343 p = ind(l,:);
db3411264b96 Remove use of tensor library in B assembly Martin Almquist <malmquist@stanford.edu> parents: 963 diff changeset	344 B{2}{(l-1)*m(2) + k}(p, p) = B_1D;
db3411264b96 Remove use of tensor library in B assembly Martin Almquist <malmquist@stanford.edu> parents: 963 diff changeset	345 end
db3411264b96 Remove use of tensor library in B assembly Martin Almquist <malmquist@stanford.edu> parents: 963 diff changeset	346 end
db3411264b96 Remove use of tensor library in B assembly Martin Almquist <malmquist@stanford.edu> parents: 963 diff changeset	347
db3411264b96 Remove use of tensor library in B assembly Martin Almquist <malmquist@stanford.edu> parents: 963 diff changeset	348 obj.B = B;
db3411264b96 Remove use of tensor library in B assembly Martin Almquist <malmquist@stanford.edu> parents: 963 diff changeset	349
687 e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	350 end
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	351
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	352
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	353 % Closure functions return the operators applied to the own domain to close the boundary
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	354 % Penalty functions return the operators to force the solution. In the case of an interface it returns the operator applied to the other doamin.
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	355 % boundary is a string specifying the boundary e.g. 'l','r' or 'e','w','n','s'.
795 1f6b2fb69225 Revert bcSetup and update bc functions in elastic schemes to be compatible. Martin Almquist <malmquist@stanford.edu> parents: 738 diff changeset	356 % bc is a cell array of component and bc type, e.g. {1, 'd'} for Dirichlet condition
813 b374a8aa9246 Correct interface penalty strength in Elastic2dVariable Martin Almquist <malmquist@stanford.edu> parents: 795 diff changeset	357 % on the first component.
687 e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	358 % data is a function returning the data that should be applied at the boundary.
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	359 % neighbour_scheme is an instance of Scheme that should be interfaced to.
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	360 % neighbour_boundary is a string specifying which boundary to interface to.
795 1f6b2fb69225 Revert bcSetup and update bc functions in elastic schemes to be compatible. Martin Almquist <malmquist@stanford.edu> parents: 738 diff changeset	361 function [closure, penalty] = boundary_condition(obj, boundary, bc, tuning)
738 aa4ef495f1fd Set Dirichlet tuning as parameter in BC function in Elastic2dVariable. Martin Almquist <malmquist@stanford.edu> parents: 734 diff changeset	362 default_arg('tuning', 1.2);
687 e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	363
795 1f6b2fb69225 Revert bcSetup and update bc functions in elastic schemes to be compatible. Martin Almquist <malmquist@stanford.edu> parents: 738 diff changeset	364 assert( iscell(bc), 'The BC type must be a 2x1 cell array' );
1f6b2fb69225 Revert bcSetup and update bc functions in elastic schemes to be compatible. Martin Almquist <malmquist@stanford.edu> parents: 738 diff changeset	365 comp = bc{1};
1f6b2fb69225 Revert bcSetup and update bc functions in elastic schemes to be compatible. Martin Almquist <malmquist@stanford.edu> parents: 738 diff changeset	366 type = bc{2};
729 aa8cf3851de8 Update multiblock.DiffOp to work for systems. Martin Almquist <malmquist@stanford.edu> parents: 727 diff changeset	367
687 e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	368 % j is the coordinate direction of the boundary
726 37d5d69b1a0d Refactor BC code in Elastic2dVariable Martin Almquist <malmquist@stanford.edu> parents: 689 diff changeset	369 j = obj.get_boundary_number(boundary);
970 23d9ca6755be Add getBoundaryQuadrature in Elastic2dVariable. Rename get_boundary_operator -> getBoundaryOperator. Add operators in getBoundaryOperator, with full size so that they work with multiblock.DiffOp.getBoundaryOperator. Martin Almquist <malmquist@stanford.edu> parents: 968 diff changeset	370 [e, T, tau, H_gamma] = obj.getBoundaryOperator({'e','T','tau','H'}, boundary);
687 e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	371
963 c75ddd568fcc Turn alpha into a boundary operator. Add properties H_w etc for getBoundaryQuadrature to work. Martin Almquist <malmquist@stanford.edu> parents: 961 diff changeset	372
687 e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	373 E = obj.E;
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	374 Hi = obj.Hi;
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	375 LAMBDA = obj.LAMBDA;
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	376 MU = obj.MU;
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	377 RHOi = obj.RHOi;
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	378
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	379 dim = obj.dim;
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	380 m_tot = obj.grid.N();
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	381
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	382 % Preallocate
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	383 closure = sparse(dimm_tot, dimm_tot);
795 1f6b2fb69225 Revert bcSetup and update bc functions in elastic schemes to be compatible. Martin Almquist <malmquist@stanford.edu> parents: 738 diff changeset	384 penalty = sparse(dim*m_tot, m_tot/obj.m(j));
687 e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	385
795 1f6b2fb69225 Revert bcSetup and update bc functions in elastic schemes to be compatible. Martin Almquist <malmquist@stanford.edu> parents: 738 diff changeset	386 k = comp;
1f6b2fb69225 Revert bcSetup and update bc functions in elastic schemes to be compatible. Martin Almquist <malmquist@stanford.edu> parents: 738 diff changeset	387 switch type
1f6b2fb69225 Revert bcSetup and update bc functions in elastic schemes to be compatible. Martin Almquist <malmquist@stanford.edu> parents: 738 diff changeset	388
1f6b2fb69225 Revert bcSetup and update bc functions in elastic schemes to be compatible. Martin Almquist <malmquist@stanford.edu> parents: 738 diff changeset	389 % Dirichlet boundary condition
1f6b2fb69225 Revert bcSetup and update bc functions in elastic schemes to be compatible. Martin Almquist <malmquist@stanford.edu> parents: 738 diff changeset	390 case {'D','d','dirichlet','Dirichlet'}
687 e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	391
970 23d9ca6755be Add getBoundaryQuadrature in Elastic2dVariable. Rename get_boundary_operator -> getBoundaryOperator. Add operators in getBoundaryOperator, with full size so that they work with multiblock.DiffOp.getBoundaryOperator. Martin Almquist <malmquist@stanford.edu> parents: 968 diff changeset	392 alpha = obj.getBoundaryOperator('alpha', boundary);
687 e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	393
795 1f6b2fb69225 Revert bcSetup and update bc functions in elastic schemes to be compatible. Martin Almquist <malmquist@stanford.edu> parents: 738 diff changeset	394 % Loop over components that Dirichlet penalties end up on
1f6b2fb69225 Revert bcSetup and update bc functions in elastic schemes to be compatible. Martin Almquist <malmquist@stanford.edu> parents: 738 diff changeset	395 for i = 1:dim
958 72cd29107a9a Temporary changes in multiblock.DiffOp. Change traction operators in Elastic2dvariable to be true boundary operators. But adjoint FD conv test fails for dirichlet BC so need to debug! Martin Almquist <malmquist@stanford.edu> parents: 919 diff changeset	396 C = transpose(T{k,i});
963 c75ddd568fcc Turn alpha into a boundary operator. Add properties H_w etc for getBoundaryQuadrature to work. Martin Almquist <malmquist@stanford.edu> parents: 961 diff changeset	397 A = -e*transpose(alpha{i,k});
958 72cd29107a9a Temporary changes in multiblock.DiffOp. Change traction operators in Elastic2dvariable to be true boundary operators. But adjoint FD conv test fails for dirichlet BC so need to debug! Martin Almquist <malmquist@stanford.edu> parents: 919 diff changeset	398 B = A + e*C;
813 b374a8aa9246 Correct interface penalty strength in Elastic2dVariable Martin Almquist <malmquist@stanford.edu> parents: 795 diff changeset	399 closure = closure + E{i}RHOiHiB'eH_gamma(e'*E{k}' );
795 1f6b2fb69225 Revert bcSetup and update bc functions in elastic schemes to be compatible. Martin Almquist <malmquist@stanford.edu> parents: 738 diff changeset	400 penalty = penalty - E{i}RHOiHiB'e*H_gamma;
813 b374a8aa9246 Correct interface penalty strength in Elastic2dVariable Martin Almquist <malmquist@stanford.edu> parents: 795 diff changeset	401 end
687 e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	402
795 1f6b2fb69225 Revert bcSetup and update bc functions in elastic schemes to be compatible. Martin Almquist <malmquist@stanford.edu> parents: 738 diff changeset	403 % Free boundary condition
1f6b2fb69225 Revert bcSetup and update bc functions in elastic schemes to be compatible. Martin Almquist <malmquist@stanford.edu> parents: 738 diff changeset	404 case {'F','f','Free','free','traction','Traction','t','T'}
958 72cd29107a9a Temporary changes in multiblock.DiffOp. Change traction operators in Elastic2dvariable to be true boundary operators. But adjoint FD conv test fails for dirichlet BC so need to debug! Martin Almquist <malmquist@stanford.edu> parents: 919 diff changeset	405 closure = closure - E{k}RHOiHieH_gamma*tau{k}';
795 1f6b2fb69225 Revert bcSetup and update bc functions in elastic schemes to be compatible. Martin Almquist <malmquist@stanford.edu> parents: 738 diff changeset	406 penalty = penalty + E{k}RHOiHieH_gamma;
687 e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	407
795 1f6b2fb69225 Revert bcSetup and update bc functions in elastic schemes to be compatible. Martin Almquist <malmquist@stanford.edu> parents: 738 diff changeset	408 % Unknown boundary condition
1f6b2fb69225 Revert bcSetup and update bc functions in elastic schemes to be compatible. Martin Almquist <malmquist@stanford.edu> parents: 738 diff changeset	409 otherwise
1f6b2fb69225 Revert bcSetup and update bc functions in elastic schemes to be compatible. Martin Almquist <malmquist@stanford.edu> parents: 738 diff changeset	410 error('No such boundary condition: type = %s',type);
687 e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	411 end
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	412 end
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	413
960 ac566f3dc9b3 Add type to Elastic2dVariable.interface Martin Almquist <malmquist@stanford.edu> parents: 959 diff changeset	414 % type Struct that specifies the interface coupling.
ac566f3dc9b3 Add type to Elastic2dVariable.interface Martin Almquist <malmquist@stanford.edu> parents: 959 diff changeset	415 % Fields:
ac566f3dc9b3 Add type to Elastic2dVariable.interface Martin Almquist <malmquist@stanford.edu> parents: 959 diff changeset	416 % -- tuning: penalty strength, defaults to 1.2
ac566f3dc9b3 Add type to Elastic2dVariable.interface Martin Almquist <malmquist@stanford.edu> parents: 959 diff changeset	417 % -- interpolation: type of interpolation, default 'none'
ac566f3dc9b3 Add type to Elastic2dVariable.interface Martin Almquist <malmquist@stanford.edu> parents: 959 diff changeset	418 function [closure, penalty] = interface(obj,boundary,neighbour_scheme,neighbour_boundary,type)
ac566f3dc9b3 Add type to Elastic2dVariable.interface Martin Almquist <malmquist@stanford.edu> parents: 959 diff changeset	419
ac566f3dc9b3 Add type to Elastic2dVariable.interface Martin Almquist <malmquist@stanford.edu> parents: 959 diff changeset	420 defaultType.tuning = 1.2;
ac566f3dc9b3 Add type to Elastic2dVariable.interface Martin Almquist <malmquist@stanford.edu> parents: 959 diff changeset	421 defaultType.interpolation = 'none';
ac566f3dc9b3 Add type to Elastic2dVariable.interface Martin Almquist <malmquist@stanford.edu> parents: 959 diff changeset	422 default_struct('type', defaultType);
ac566f3dc9b3 Add type to Elastic2dVariable.interface Martin Almquist <malmquist@stanford.edu> parents: 959 diff changeset	423
ac566f3dc9b3 Add type to Elastic2dVariable.interface Martin Almquist <malmquist@stanford.edu> parents: 959 diff changeset	424 switch type.interpolation
ac566f3dc9b3 Add type to Elastic2dVariable.interface Martin Almquist <malmquist@stanford.edu> parents: 959 diff changeset	425 case {'none', ''}
ac566f3dc9b3 Add type to Elastic2dVariable.interface Martin Almquist <malmquist@stanford.edu> parents: 959 diff changeset	426 [closure, penalty] = interfaceStandard(obj,boundary,neighbour_scheme,neighbour_boundary,type);
ac566f3dc9b3 Add type to Elastic2dVariable.interface Martin Almquist <malmquist@stanford.edu> parents: 959 diff changeset	427 case {'op','OP'}
ac566f3dc9b3 Add type to Elastic2dVariable.interface Martin Almquist <malmquist@stanford.edu> parents: 959 diff changeset	428 [closure, penalty] = interfaceNonConforming(obj,boundary,neighbour_scheme,neighbour_boundary,type);
ac566f3dc9b3 Add type to Elastic2dVariable.interface Martin Almquist <malmquist@stanford.edu> parents: 959 diff changeset	429 otherwise
ac566f3dc9b3 Add type to Elastic2dVariable.interface Martin Almquist <malmquist@stanford.edu> parents: 959 diff changeset	430 error('Unknown type of interpolation: %s ', type.interpolation);
ac566f3dc9b3 Add type to Elastic2dVariable.interface Martin Almquist <malmquist@stanford.edu> parents: 959 diff changeset	431 end
ac566f3dc9b3 Add type to Elastic2dVariable.interface Martin Almquist <malmquist@stanford.edu> parents: 959 diff changeset	432 end
ac566f3dc9b3 Add type to Elastic2dVariable.interface Martin Almquist <malmquist@stanford.edu> parents: 959 diff changeset	433
ac566f3dc9b3 Add type to Elastic2dVariable.interface Martin Almquist <malmquist@stanford.edu> parents: 959 diff changeset	434 function [closure, penalty] = interfaceStandard(obj,boundary,neighbour_scheme,neighbour_boundary,type)
963 c75ddd568fcc Turn alpha into a boundary operator. Add properties H_w etc for getBoundaryQuadrature to work. Martin Almquist <malmquist@stanford.edu> parents: 961 diff changeset	435 tuning = type.tuning;
c75ddd568fcc Turn alpha into a boundary operator. Add properties H_w etc for getBoundaryQuadrature to work. Martin Almquist <malmquist@stanford.edu> parents: 961 diff changeset	436
687 e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	437 % u denotes the solution in the own domain
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	438 % v denotes the solution in the neighbour domain
727 6d5953fc090e First implementation of elastic interface Martin Almquist <malmquist@stanford.edu> parents: 726 diff changeset	439 % Operators without subscripts are from the own domain.
6d5953fc090e First implementation of elastic interface Martin Almquist <malmquist@stanford.edu> parents: 726 diff changeset	440
6d5953fc090e First implementation of elastic interface Martin Almquist <malmquist@stanford.edu> parents: 726 diff changeset	441 % Get boundary operators
972 104f0af001e0 Clean up Elastic2dVariable.interfaceStandard Martin Almquist <malmquist@stanford.edu> parents: 970 diff changeset	442 [e, tau] = obj.getBoundaryOperator({'e_tot','tau_tot'}, boundary);
104f0af001e0 Clean up Elastic2dVariable.interfaceStandard Martin Almquist <malmquist@stanford.edu> parents: 970 diff changeset	443 [e_v, tau_v] = neighbour_scheme.getBoundaryOperator({'e_tot','tau_tot'}, neighbour_boundary);
104f0af001e0 Clean up Elastic2dVariable.interfaceStandard Martin Almquist <malmquist@stanford.edu> parents: 970 diff changeset	444 H_gamma = obj.getBoundaryQuadrature(boundary);
727 6d5953fc090e First implementation of elastic interface Martin Almquist <malmquist@stanford.edu> parents: 726 diff changeset	445
6d5953fc090e First implementation of elastic interface Martin Almquist <malmquist@stanford.edu> parents: 726 diff changeset	446 % Operators and quantities that correspond to the own domain only
972 104f0af001e0 Clean up Elastic2dVariable.interfaceStandard Martin Almquist <malmquist@stanford.edu> parents: 970 diff changeset	447 Hi = obj.Hi_kron;
104f0af001e0 Clean up Elastic2dVariable.interfaceStandard Martin Almquist <malmquist@stanford.edu> parents: 970 diff changeset	448 RHOi = obj.RHOi_kron;
727 6d5953fc090e First implementation of elastic interface Martin Almquist <malmquist@stanford.edu> parents: 726 diff changeset	449
972 104f0af001e0 Clean up Elastic2dVariable.interfaceStandard Martin Almquist <malmquist@stanford.edu> parents: 970 diff changeset	450 % Penalty strength operators
104f0af001e0 Clean up Elastic2dVariable.interfaceStandard Martin Almquist <malmquist@stanford.edu> parents: 970 diff changeset	451 alpha_u = 1/4*obj.getBoundaryOperator('alpha_tot', boundary);
104f0af001e0 Clean up Elastic2dVariable.interfaceStandard Martin Almquist <malmquist@stanford.edu> parents: 970 diff changeset	452 alpha_v = 1/4*neighbour_scheme.getBoundaryOperator('alpha_tot', neighbour_boundary);
727 6d5953fc090e First implementation of elastic interface Martin Almquist <malmquist@stanford.edu> parents: 726 diff changeset	453
972 104f0af001e0 Clean up Elastic2dVariable.interfaceStandard Martin Almquist <malmquist@stanford.edu> parents: 970 diff changeset	454 closure = -RHOiHieH_gamma(alpha_u' + alpha_v'e_ve');
104f0af001e0 Clean up Elastic2dVariable.interfaceStandard Martin Almquist <malmquist@stanford.edu> parents: 970 diff changeset	455 penalty = RHOiHieH_gamma(alpha_u'ee_v' + alpha_v');
727 6d5953fc090e First implementation of elastic interface Martin Almquist <malmquist@stanford.edu> parents: 726 diff changeset	456
972 104f0af001e0 Clean up Elastic2dVariable.interfaceStandard Martin Almquist <malmquist@stanford.edu> parents: 970 diff changeset	457 closure = closure - 1/2RHOiHieH_gamma*tau';
104f0af001e0 Clean up Elastic2dVariable.interfaceStandard Martin Almquist <malmquist@stanford.edu> parents: 970 diff changeset	458 penalty = penalty - 1/2RHOiHieH_gamma*tau_v';
727 6d5953fc090e First implementation of elastic interface Martin Almquist <malmquist@stanford.edu> parents: 726 diff changeset	459
972 104f0af001e0 Clean up Elastic2dVariable.interfaceStandard Martin Almquist <malmquist@stanford.edu> parents: 970 diff changeset	460 closure = closure + 1/2RHOiHitauH_gamma*e';
104f0af001e0 Clean up Elastic2dVariable.interfaceStandard Martin Almquist <malmquist@stanford.edu> parents: 970 diff changeset	461 penalty = penalty - 1/2RHOiHitauH_gamma*e_v';
727 6d5953fc090e First implementation of elastic interface Martin Almquist <malmquist@stanford.edu> parents: 726 diff changeset	462
687 e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	463 end
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	464
960 ac566f3dc9b3 Add type to Elastic2dVariable.interface Martin Almquist <malmquist@stanford.edu> parents: 959 diff changeset	465 function [closure, penalty] = interfaceNonConforming(obj,boundary,neighbour_scheme,neighbour_boundary,type)
ac566f3dc9b3 Add type to Elastic2dVariable.interface Martin Almquist <malmquist@stanford.edu> parents: 959 diff changeset	466 error('Non-conforming interfaces not implemented yet.');
ac566f3dc9b3 Add type to Elastic2dVariable.interface Martin Almquist <malmquist@stanford.edu> parents: 959 diff changeset	467 end
ac566f3dc9b3 Add type to Elastic2dVariable.interface Martin Almquist <malmquist@stanford.edu> parents: 959 diff changeset	468
687 e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	469 % Returns the coordinate number and outward normal component for the boundary specified by the string boundary.
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	470 function [j, nj] = get_boundary_number(obj, boundary)
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	471
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	472 switch boundary
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	473 case {'w','W','west','West', 'e', 'E', 'east', 'East'}
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	474 j = 1;
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	475 case {'s','S','south','South', 'n', 'N', 'north', 'North'}
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	476 j = 2;
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	477 otherwise
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	478 error('No such boundary: boundary = %s',boundary);
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	479 end
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	480
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	481 switch boundary
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	482 case {'w','W','west','West','s','S','south','South'}
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	483 nj = -1;
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	484 case {'e', 'E', 'east', 'East','n', 'N', 'north', 'North'}
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	485 nj = 1;
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	486 end
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	487 end
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	488
726 37d5d69b1a0d Refactor BC code in Elastic2dVariable Martin Almquist <malmquist@stanford.edu> parents: 689 diff changeset	489 % Returns the boundary operator op for the boundary specified by the string boundary.
37d5d69b1a0d Refactor BC code in Elastic2dVariable Martin Almquist <malmquist@stanford.edu> parents: 689 diff changeset	490 % op: may be a cell array of strings
970 23d9ca6755be Add getBoundaryQuadrature in Elastic2dVariable. Rename get_boundary_operator -> getBoundaryOperator. Add operators in getBoundaryOperator, with full size so that they work with multiblock.DiffOp.getBoundaryOperator. Martin Almquist <malmquist@stanford.edu> parents: 968 diff changeset	491 % Only operators with name *_tot can be used with multiblock.DiffOp.getBoundaryOperator()
23d9ca6755be Add getBoundaryQuadrature in Elastic2dVariable. Rename get_boundary_operator -> getBoundaryOperator. Add operators in getBoundaryOperator, with full size so that they work with multiblock.DiffOp.getBoundaryOperator. Martin Almquist <malmquist@stanford.edu> parents: 968 diff changeset	492 function [varargout] = getBoundaryOperator(obj, op, boundary)
687 e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	493
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	494 switch boundary
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	495 case {'w','W','west','West', 'e', 'E', 'east', 'East'}
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	496 j = 1;
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	497 case {'s','S','south','South', 'n', 'N', 'north', 'North'}
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	498 j = 2;
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	499 otherwise
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	500 error('No such boundary: boundary = %s',boundary);
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	501 end
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	502
726 37d5d69b1a0d Refactor BC code in Elastic2dVariable Martin Almquist <malmquist@stanford.edu> parents: 689 diff changeset	503 if ~iscell(op)
37d5d69b1a0d Refactor BC code in Elastic2dVariable Martin Almquist <malmquist@stanford.edu> parents: 689 diff changeset	504 op = {op};
37d5d69b1a0d Refactor BC code in Elastic2dVariable Martin Almquist <malmquist@stanford.edu> parents: 689 diff changeset	505 end
37d5d69b1a0d Refactor BC code in Elastic2dVariable Martin Almquist <malmquist@stanford.edu> parents: 689 diff changeset	506
958 72cd29107a9a Temporary changes in multiblock.DiffOp. Change traction operators in Elastic2dvariable to be true boundary operators. But adjoint FD conv test fails for dirichlet BC so need to debug! Martin Almquist <malmquist@stanford.edu> parents: 919 diff changeset	507 for k = 1:length(op)
72cd29107a9a Temporary changes in multiblock.DiffOp. Change traction operators in Elastic2dvariable to be true boundary operators. But adjoint FD conv test fails for dirichlet BC so need to debug! Martin Almquist <malmquist@stanford.edu> parents: 919 diff changeset	508 switch op{k}
726 37d5d69b1a0d Refactor BC code in Elastic2dVariable Martin Almquist <malmquist@stanford.edu> parents: 689 diff changeset	509 case 'e'
37d5d69b1a0d Refactor BC code in Elastic2dVariable Martin Almquist <malmquist@stanford.edu> parents: 689 diff changeset	510 switch boundary
37d5d69b1a0d Refactor BC code in Elastic2dVariable Martin Almquist <malmquist@stanford.edu> parents: 689 diff changeset	511 case {'w','W','west','West','s','S','south','South'}
958 72cd29107a9a Temporary changes in multiblock.DiffOp. Change traction operators in Elastic2dvariable to be true boundary operators. But adjoint FD conv test fails for dirichlet BC so need to debug! Martin Almquist <malmquist@stanford.edu> parents: 919 diff changeset	512 varargout{k} = obj.e_l{j};
726 37d5d69b1a0d Refactor BC code in Elastic2dVariable Martin Almquist <malmquist@stanford.edu> parents: 689 diff changeset	513 case {'e', 'E', 'east', 'East','n', 'N', 'north', 'North'}
958 72cd29107a9a Temporary changes in multiblock.DiffOp. Change traction operators in Elastic2dvariable to be true boundary operators. But adjoint FD conv test fails for dirichlet BC so need to debug! Martin Almquist <malmquist@stanford.edu> parents: 919 diff changeset	514 varargout{k} = obj.e_r{j};
726 37d5d69b1a0d Refactor BC code in Elastic2dVariable Martin Almquist <malmquist@stanford.edu> parents: 689 diff changeset	515 end
970 23d9ca6755be Add getBoundaryQuadrature in Elastic2dVariable. Rename get_boundary_operator -> getBoundaryOperator. Add operators in getBoundaryOperator, with full size so that they work with multiblock.DiffOp.getBoundaryOperator. Martin Almquist <malmquist@stanford.edu> parents: 968 diff changeset	516
23d9ca6755be Add getBoundaryQuadrature in Elastic2dVariable. Rename get_boundary_operator -> getBoundaryOperator. Add operators in getBoundaryOperator, with full size so that they work with multiblock.DiffOp.getBoundaryOperator. Martin Almquist <malmquist@stanford.edu> parents: 968 diff changeset	517 case 'e_tot'
23d9ca6755be Add getBoundaryQuadrature in Elastic2dVariable. Rename get_boundary_operator -> getBoundaryOperator. Add operators in getBoundaryOperator, with full size so that they work with multiblock.DiffOp.getBoundaryOperator. Martin Almquist <malmquist@stanford.edu> parents: 968 diff changeset	518 e = obj.getBoundaryOperator('e', boundary);
23d9ca6755be Add getBoundaryQuadrature in Elastic2dVariable. Rename get_boundary_operator -> getBoundaryOperator. Add operators in getBoundaryOperator, with full size so that they work with multiblock.DiffOp.getBoundaryOperator. Martin Almquist <malmquist@stanford.edu> parents: 968 diff changeset	519 I_dim = speye(obj.dim, obj.dim);
23d9ca6755be Add getBoundaryQuadrature in Elastic2dVariable. Rename get_boundary_operator -> getBoundaryOperator. Add operators in getBoundaryOperator, with full size so that they work with multiblock.DiffOp.getBoundaryOperator. Martin Almquist <malmquist@stanford.edu> parents: 968 diff changeset	520 varargout{k} = kron(e, I_dim);
23d9ca6755be Add getBoundaryQuadrature in Elastic2dVariable. Rename get_boundary_operator -> getBoundaryOperator. Add operators in getBoundaryOperator, with full size so that they work with multiblock.DiffOp.getBoundaryOperator. Martin Almquist <malmquist@stanford.edu> parents: 968 diff changeset	521
726 37d5d69b1a0d Refactor BC code in Elastic2dVariable Martin Almquist <malmquist@stanford.edu> parents: 689 diff changeset	522 case 'd'
37d5d69b1a0d Refactor BC code in Elastic2dVariable Martin Almquist <malmquist@stanford.edu> parents: 689 diff changeset	523 switch boundary
37d5d69b1a0d Refactor BC code in Elastic2dVariable Martin Almquist <malmquist@stanford.edu> parents: 689 diff changeset	524 case {'w','W','west','West','s','S','south','South'}
958 72cd29107a9a Temporary changes in multiblock.DiffOp. Change traction operators in Elastic2dvariable to be true boundary operators. But adjoint FD conv test fails for dirichlet BC so need to debug! Martin Almquist <malmquist@stanford.edu> parents: 919 diff changeset	525 varargout{k} = obj.d1_l{j};
726 37d5d69b1a0d Refactor BC code in Elastic2dVariable Martin Almquist <malmquist@stanford.edu> parents: 689 diff changeset	526 case {'e', 'E', 'east', 'East','n', 'N', 'north', 'North'}
958 72cd29107a9a Temporary changes in multiblock.DiffOp. Change traction operators in Elastic2dvariable to be true boundary operators. But adjoint FD conv test fails for dirichlet BC so need to debug! Martin Almquist <malmquist@stanford.edu> parents: 919 diff changeset	527 varargout{k} = obj.d1_r{j};
726 37d5d69b1a0d Refactor BC code in Elastic2dVariable Martin Almquist <malmquist@stanford.edu> parents: 689 diff changeset	528 end
970 23d9ca6755be Add getBoundaryQuadrature in Elastic2dVariable. Rename get_boundary_operator -> getBoundaryOperator. Add operators in getBoundaryOperator, with full size so that they work with multiblock.DiffOp.getBoundaryOperator. Martin Almquist <malmquist@stanford.edu> parents: 968 diff changeset	529
726 37d5d69b1a0d Refactor BC code in Elastic2dVariable Martin Almquist <malmquist@stanford.edu> parents: 689 diff changeset	530 case 'T'
37d5d69b1a0d Refactor BC code in Elastic2dVariable Martin Almquist <malmquist@stanford.edu> parents: 689 diff changeset	531 switch boundary
37d5d69b1a0d Refactor BC code in Elastic2dVariable Martin Almquist <malmquist@stanford.edu> parents: 689 diff changeset	532 case {'w','W','west','West','s','S','south','South'}
958 72cd29107a9a Temporary changes in multiblock.DiffOp. Change traction operators in Elastic2dvariable to be true boundary operators. But adjoint FD conv test fails for dirichlet BC so need to debug! Martin Almquist <malmquist@stanford.edu> parents: 919 diff changeset	533 varargout{k} = obj.T_l{j};
726 37d5d69b1a0d Refactor BC code in Elastic2dVariable Martin Almquist <malmquist@stanford.edu> parents: 689 diff changeset	534 case {'e', 'E', 'east', 'East','n', 'N', 'north', 'North'}
958 72cd29107a9a Temporary changes in multiblock.DiffOp. Change traction operators in Elastic2dvariable to be true boundary operators. But adjoint FD conv test fails for dirichlet BC so need to debug! Martin Almquist <malmquist@stanford.edu> parents: 919 diff changeset	535 varargout{k} = obj.T_r{j};
726 37d5d69b1a0d Refactor BC code in Elastic2dVariable Martin Almquist <malmquist@stanford.edu> parents: 689 diff changeset	536 end
970 23d9ca6755be Add getBoundaryQuadrature in Elastic2dVariable. Rename get_boundary_operator -> getBoundaryOperator. Add operators in getBoundaryOperator, with full size so that they work with multiblock.DiffOp.getBoundaryOperator. Martin Almquist <malmquist@stanford.edu> parents: 968 diff changeset	537
726 37d5d69b1a0d Refactor BC code in Elastic2dVariable Martin Almquist <malmquist@stanford.edu> parents: 689 diff changeset	538 case 'tau'
37d5d69b1a0d Refactor BC code in Elastic2dVariable Martin Almquist <malmquist@stanford.edu> parents: 689 diff changeset	539 switch boundary
37d5d69b1a0d Refactor BC code in Elastic2dVariable Martin Almquist <malmquist@stanford.edu> parents: 689 diff changeset	540 case {'w','W','west','West','s','S','south','South'}
958 72cd29107a9a Temporary changes in multiblock.DiffOp. Change traction operators in Elastic2dvariable to be true boundary operators. But adjoint FD conv test fails for dirichlet BC so need to debug! Martin Almquist <malmquist@stanford.edu> parents: 919 diff changeset	541 varargout{k} = obj.tau_l{j};
726 37d5d69b1a0d Refactor BC code in Elastic2dVariable Martin Almquist <malmquist@stanford.edu> parents: 689 diff changeset	542 case {'e', 'E', 'east', 'East','n', 'N', 'north', 'North'}
958 72cd29107a9a Temporary changes in multiblock.DiffOp. Change traction operators in Elastic2dvariable to be true boundary operators. But adjoint FD conv test fails for dirichlet BC so need to debug! Martin Almquist <malmquist@stanford.edu> parents: 919 diff changeset	543 varargout{k} = obj.tau_r{j};
855 5751262b323b Add 1D quadrature matrices as property in Elastic2dVariable. Martin Almquist <malmquist@stanford.edu> parents: 795 diff changeset	544 end
970 23d9ca6755be Add getBoundaryQuadrature in Elastic2dVariable. Rename get_boundary_operator -> getBoundaryOperator. Add operators in getBoundaryOperator, with full size so that they work with multiblock.DiffOp.getBoundaryOperator. Martin Almquist <malmquist@stanford.edu> parents: 968 diff changeset	545
23d9ca6755be Add getBoundaryQuadrature in Elastic2dVariable. Rename get_boundary_operator -> getBoundaryOperator. Add operators in getBoundaryOperator, with full size so that they work with multiblock.DiffOp.getBoundaryOperator. Martin Almquist <malmquist@stanford.edu> parents: 968 diff changeset	546 case 'tau_tot'
23d9ca6755be Add getBoundaryQuadrature in Elastic2dVariable. Rename get_boundary_operator -> getBoundaryOperator. Add operators in getBoundaryOperator, with full size so that they work with multiblock.DiffOp.getBoundaryOperator. Martin Almquist <malmquist@stanford.edu> parents: 968 diff changeset	547 [e, tau] = obj.getBoundaryOperator({'e', 'tau'}, boundary);
23d9ca6755be Add getBoundaryQuadrature in Elastic2dVariable. Rename get_boundary_operator -> getBoundaryOperator. Add operators in getBoundaryOperator, with full size so that they work with multiblock.DiffOp.getBoundaryOperator. Martin Almquist <malmquist@stanford.edu> parents: 968 diff changeset	548
23d9ca6755be Add getBoundaryQuadrature in Elastic2dVariable. Rename get_boundary_operator -> getBoundaryOperator. Add operators in getBoundaryOperator, with full size so that they work with multiblock.DiffOp.getBoundaryOperator. Martin Almquist <malmquist@stanford.edu> parents: 968 diff changeset	549 I_dim = speye(obj.dim, obj.dim);
23d9ca6755be Add getBoundaryQuadrature in Elastic2dVariable. Rename get_boundary_operator -> getBoundaryOperator. Add operators in getBoundaryOperator, with full size so that they work with multiblock.DiffOp.getBoundaryOperator. Martin Almquist <malmquist@stanford.edu> parents: 968 diff changeset	550 e_tot = kron(e, I_dim);
23d9ca6755be Add getBoundaryQuadrature in Elastic2dVariable. Rename get_boundary_operator -> getBoundaryOperator. Add operators in getBoundaryOperator, with full size so that they work with multiblock.DiffOp.getBoundaryOperator. Martin Almquist <malmquist@stanford.edu> parents: 968 diff changeset	551 E = obj.E;
23d9ca6755be Add getBoundaryQuadrature in Elastic2dVariable. Rename get_boundary_operator -> getBoundaryOperator. Add operators in getBoundaryOperator, with full size so that they work with multiblock.DiffOp.getBoundaryOperator. Martin Almquist <malmquist@stanford.edu> parents: 968 diff changeset	552 tau_tot = (e_tot'E{1}e*tau{1}')';
23d9ca6755be Add getBoundaryQuadrature in Elastic2dVariable. Rename get_boundary_operator -> getBoundaryOperator. Add operators in getBoundaryOperator, with full size so that they work with multiblock.DiffOp.getBoundaryOperator. Martin Almquist <malmquist@stanford.edu> parents: 968 diff changeset	553 for i = 2:obj.dim
23d9ca6755be Add getBoundaryQuadrature in Elastic2dVariable. Rename get_boundary_operator -> getBoundaryOperator. Add operators in getBoundaryOperator, with full size so that they work with multiblock.DiffOp.getBoundaryOperator. Martin Almquist <malmquist@stanford.edu> parents: 968 diff changeset	554 tau_tot = tau_tot + (e_tot'E{i}e*tau{i}')';
23d9ca6755be Add getBoundaryQuadrature in Elastic2dVariable. Rename get_boundary_operator -> getBoundaryOperator. Add operators in getBoundaryOperator, with full size so that they work with multiblock.DiffOp.getBoundaryOperator. Martin Almquist <malmquist@stanford.edu> parents: 968 diff changeset	555 end
23d9ca6755be Add getBoundaryQuadrature in Elastic2dVariable. Rename get_boundary_operator -> getBoundaryOperator. Add operators in getBoundaryOperator, with full size so that they work with multiblock.DiffOp.getBoundaryOperator. Martin Almquist <malmquist@stanford.edu> parents: 968 diff changeset	556 varargout{k} = tau_tot;
23d9ca6755be Add getBoundaryQuadrature in Elastic2dVariable. Rename get_boundary_operator -> getBoundaryOperator. Add operators in getBoundaryOperator, with full size so that they work with multiblock.DiffOp.getBoundaryOperator. Martin Almquist <malmquist@stanford.edu> parents: 968 diff changeset	557
963 c75ddd568fcc Turn alpha into a boundary operator. Add properties H_w etc for getBoundaryQuadrature to work. Martin Almquist <malmquist@stanford.edu> parents: 961 diff changeset	558 case 'H'
c75ddd568fcc Turn alpha into a boundary operator. Add properties H_w etc for getBoundaryQuadrature to work. Martin Almquist <malmquist@stanford.edu> parents: 961 diff changeset	559 varargout{k} = obj.H_boundary{j};
970 23d9ca6755be Add getBoundaryQuadrature in Elastic2dVariable. Rename get_boundary_operator -> getBoundaryOperator. Add operators in getBoundaryOperator, with full size so that they work with multiblock.DiffOp.getBoundaryOperator. Martin Almquist <malmquist@stanford.edu> parents: 968 diff changeset	560
861 607c631f175e Small changes to Elastic2dVariable to facilitate adjoing gradient computation. Martin Almquist <malmquist@stanford.edu> parents: 855 diff changeset	561 case 'alpha'
882 14fee299ada2 In Elastic2dVariable: Improve notation for borrowing constants. Update interface constant to match corrected derivation. Tests ok. Martin Almquist <malmquist@stanford.edu> parents: 861 diff changeset	562 % alpha = alpha(i,j) is the penalty strength for displacement BC.
970 23d9ca6755be Add getBoundaryQuadrature in Elastic2dVariable. Rename get_boundary_operator -> getBoundaryOperator. Add operators in getBoundaryOperator, with full size so that they work with multiblock.DiffOp.getBoundaryOperator. Martin Almquist <malmquist@stanford.edu> parents: 968 diff changeset	563 e = obj.getBoundaryOperator('e', boundary);
963 c75ddd568fcc Turn alpha into a boundary operator. Add properties H_w etc for getBoundaryQuadrature to work. Martin Almquist <malmquist@stanford.edu> parents: 961 diff changeset	564
861 607c631f175e Small changes to Elastic2dVariable to facilitate adjoing gradient computation. Martin Almquist <malmquist@stanford.edu> parents: 855 diff changeset	565 tuning = 1.2;
607c631f175e Small changes to Elastic2dVariable to facilitate adjoing gradient computation. Martin Almquist <malmquist@stanford.edu> parents: 855 diff changeset	566 LAMBDA = obj.LAMBDA;
607c631f175e Small changes to Elastic2dVariable to facilitate adjoing gradient computation. Martin Almquist <malmquist@stanford.edu> parents: 855 diff changeset	567 MU = obj.MU;
607c631f175e Small changes to Elastic2dVariable to facilitate adjoing gradient computation. Martin Almquist <malmquist@stanford.edu> parents: 855 diff changeset	568
607c631f175e Small changes to Elastic2dVariable to facilitate adjoing gradient computation. Martin Almquist <malmquist@stanford.edu> parents: 855 diff changeset	569 dim = obj.dim;
919 e30aaa4a3e09 Bugfix in Elastic2dVariable.get_boundary_ops Martin Almquist <malmquist@stanford.edu> parents: 882 diff changeset	570 theta_R = obj.theta_R{j};
e30aaa4a3e09 Bugfix in Elastic2dVariable.get_boundary_ops Martin Almquist <malmquist@stanford.edu> parents: 882 diff changeset	571 theta_H = obj.theta_H{j};
e30aaa4a3e09 Bugfix in Elastic2dVariable.get_boundary_ops Martin Almquist <malmquist@stanford.edu> parents: 882 diff changeset	572 theta_M = obj.theta_M{j};
861 607c631f175e Small changes to Elastic2dVariable to facilitate adjoing gradient computation. Martin Almquist <malmquist@stanford.edu> parents: 855 diff changeset	573
919 e30aaa4a3e09 Bugfix in Elastic2dVariable.get_boundary_ops Martin Almquist <malmquist@stanford.edu> parents: 882 diff changeset	574 a_lambda = dim/theta_H + 1/theta_R;
e30aaa4a3e09 Bugfix in Elastic2dVariable.get_boundary_ops Martin Almquist <malmquist@stanford.edu> parents: 882 diff changeset	575 a_mu_i = 2/theta_M;
e30aaa4a3e09 Bugfix in Elastic2dVariable.get_boundary_ops Martin Almquist <malmquist@stanford.edu> parents: 882 diff changeset	576 a_mu_ij = 2/theta_H + 1/theta_R;
861 607c631f175e Small changes to Elastic2dVariable to facilitate adjoing gradient computation. Martin Almquist <malmquist@stanford.edu> parents: 855 diff changeset	577
607c631f175e Small changes to Elastic2dVariable to facilitate adjoing gradient computation. Martin Almquist <malmquist@stanford.edu> parents: 855 diff changeset	578 d = @kroneckerDelta; % Kronecker delta
607c631f175e Small changes to Elastic2dVariable to facilitate adjoing gradient computation. Martin Almquist <malmquist@stanford.edu> parents: 855 diff changeset	579 db = @(i,j) 1-d(i,j); % Logical not of Kronecker delta
958 72cd29107a9a Temporary changes in multiblock.DiffOp. Change traction operators in Elastic2dvariable to be true boundary operators. But adjoint FD conv test fails for dirichlet BC so need to debug! Martin Almquist <malmquist@stanford.edu> parents: 919 diff changeset	580 alpha = cell(obj.dim, obj.dim);
72cd29107a9a Temporary changes in multiblock.DiffOp. Change traction operators in Elastic2dvariable to be true boundary operators. But adjoint FD conv test fails for dirichlet BC so need to debug! Martin Almquist <malmquist@stanford.edu> parents: 919 diff changeset	581
959 c226fb8c2b8a Bugfix in Elastic2dVariable.getBoundaryOps, alpha. Adjoint FD conv works now! Martin Almquist <malmquist@stanford.edu> parents: 958 diff changeset	582 alpha_func = @(i,j) tuning( d(i,j) a_lambda*LAMBDA ...
919 e30aaa4a3e09 Bugfix in Elastic2dVariable.get_boundary_ops Martin Almquist <malmquist@stanford.edu> parents: 882 diff changeset	583 + d(i,j)* a_mu_i*MU ...
e30aaa4a3e09 Bugfix in Elastic2dVariable.get_boundary_ops Martin Almquist <malmquist@stanford.edu> parents: 882 diff changeset	584 + db(i,j)a_mu_ijMU );
959 c226fb8c2b8a Bugfix in Elastic2dVariable.getBoundaryOps, alpha. Adjoint FD conv works now! Martin Almquist <malmquist@stanford.edu> parents: 958 diff changeset	585 for i = 1:obj.dim
c226fb8c2b8a Bugfix in Elastic2dVariable.getBoundaryOps, alpha. Adjoint FD conv works now! Martin Almquist <malmquist@stanford.edu> parents: 958 diff changeset	586 for l = 1:obj.dim
963 c75ddd568fcc Turn alpha into a boundary operator. Add properties H_w etc for getBoundaryQuadrature to work. Martin Almquist <malmquist@stanford.edu> parents: 961 diff changeset	587 alpha{i,l} = d(i,l)alpha_func(i,j)e;
958 72cd29107a9a Temporary changes in multiblock.DiffOp. Change traction operators in Elastic2dvariable to be true boundary operators. But adjoint FD conv test fails for dirichlet BC so need to debug! Martin Almquist <malmquist@stanford.edu> parents: 919 diff changeset	588 end
72cd29107a9a Temporary changes in multiblock.DiffOp. Change traction operators in Elastic2dvariable to be true boundary operators. But adjoint FD conv test fails for dirichlet BC so need to debug! Martin Almquist <malmquist@stanford.edu> parents: 919 diff changeset	589 end
919 e30aaa4a3e09 Bugfix in Elastic2dVariable.get_boundary_ops Martin Almquist <malmquist@stanford.edu> parents: 882 diff changeset	590
958 72cd29107a9a Temporary changes in multiblock.DiffOp. Change traction operators in Elastic2dvariable to be true boundary operators. But adjoint FD conv test fails for dirichlet BC so need to debug! Martin Almquist <malmquist@stanford.edu> parents: 919 diff changeset	591 varargout{k} = alpha;
970 23d9ca6755be Add getBoundaryQuadrature in Elastic2dVariable. Rename get_boundary_operator -> getBoundaryOperator. Add operators in getBoundaryOperator, with full size so that they work with multiblock.DiffOp.getBoundaryOperator. Martin Almquist <malmquist@stanford.edu> parents: 968 diff changeset	592
23d9ca6755be Add getBoundaryQuadrature in Elastic2dVariable. Rename get_boundary_operator -> getBoundaryOperator. Add operators in getBoundaryOperator, with full size so that they work with multiblock.DiffOp.getBoundaryOperator. Martin Almquist <malmquist@stanford.edu> parents: 968 diff changeset	593 case 'alpha_tot'
23d9ca6755be Add getBoundaryQuadrature in Elastic2dVariable. Rename get_boundary_operator -> getBoundaryOperator. Add operators in getBoundaryOperator, with full size so that they work with multiblock.DiffOp.getBoundaryOperator. Martin Almquist <malmquist@stanford.edu> parents: 968 diff changeset	594 % alpha = alpha(i,j) is the penalty strength for displacement BC.
23d9ca6755be Add getBoundaryQuadrature in Elastic2dVariable. Rename get_boundary_operator -> getBoundaryOperator. Add operators in getBoundaryOperator, with full size so that they work with multiblock.DiffOp.getBoundaryOperator. Martin Almquist <malmquist@stanford.edu> parents: 968 diff changeset	595 [e, e_tot, alpha] = obj.getBoundaryOperator({'e', 'e_tot', 'alpha'}, boundary);
23d9ca6755be Add getBoundaryQuadrature in Elastic2dVariable. Rename get_boundary_operator -> getBoundaryOperator. Add operators in getBoundaryOperator, with full size so that they work with multiblock.DiffOp.getBoundaryOperator. Martin Almquist <malmquist@stanford.edu> parents: 968 diff changeset	596 E = obj.E;
23d9ca6755be Add getBoundaryQuadrature in Elastic2dVariable. Rename get_boundary_operator -> getBoundaryOperator. Add operators in getBoundaryOperator, with full size so that they work with multiblock.DiffOp.getBoundaryOperator. Martin Almquist <malmquist@stanford.edu> parents: 968 diff changeset	597 [m, n] = size(alpha{1,1});
23d9ca6755be Add getBoundaryQuadrature in Elastic2dVariable. Rename get_boundary_operator -> getBoundaryOperator. Add operators in getBoundaryOperator, with full size so that they work with multiblock.DiffOp.getBoundaryOperator. Martin Almquist <malmquist@stanford.edu> parents: 968 diff changeset	598 alpha_tot = sparse(mobj.dim, nobj.dim);
23d9ca6755be Add getBoundaryQuadrature in Elastic2dVariable. Rename get_boundary_operator -> getBoundaryOperator. Add operators in getBoundaryOperator, with full size so that they work with multiblock.DiffOp.getBoundaryOperator. Martin Almquist <malmquist@stanford.edu> parents: 968 diff changeset	599 for i = 1:obj.dim
23d9ca6755be Add getBoundaryQuadrature in Elastic2dVariable. Rename get_boundary_operator -> getBoundaryOperator. Add operators in getBoundaryOperator, with full size so that they work with multiblock.DiffOp.getBoundaryOperator. Martin Almquist <malmquist@stanford.edu> parents: 968 diff changeset	600 for l = 1:obj.dim
23d9ca6755be Add getBoundaryQuadrature in Elastic2dVariable. Rename get_boundary_operator -> getBoundaryOperator. Add operators in getBoundaryOperator, with full size so that they work with multiblock.DiffOp.getBoundaryOperator. Martin Almquist <malmquist@stanford.edu> parents: 968 diff changeset	601 alpha_tot = alpha_tot + (e_tot'E{i}ealpha{i,l}'E{l}')';
23d9ca6755be Add getBoundaryQuadrature in Elastic2dVariable. Rename get_boundary_operator -> getBoundaryOperator. Add operators in getBoundaryOperator, with full size so that they work with multiblock.DiffOp.getBoundaryOperator. Martin Almquist <malmquist@stanford.edu> parents: 968 diff changeset	602 end
23d9ca6755be Add getBoundaryQuadrature in Elastic2dVariable. Rename get_boundary_operator -> getBoundaryOperator. Add operators in getBoundaryOperator, with full size so that they work with multiblock.DiffOp.getBoundaryOperator. Martin Almquist <malmquist@stanford.edu> parents: 968 diff changeset	603 end
23d9ca6755be Add getBoundaryQuadrature in Elastic2dVariable. Rename get_boundary_operator -> getBoundaryOperator. Add operators in getBoundaryOperator, with full size so that they work with multiblock.DiffOp.getBoundaryOperator. Martin Almquist <malmquist@stanford.edu> parents: 968 diff changeset	604 varargout{k} = alpha_tot;
23d9ca6755be Add getBoundaryQuadrature in Elastic2dVariable. Rename get_boundary_operator -> getBoundaryOperator. Add operators in getBoundaryOperator, with full size so that they work with multiblock.DiffOp.getBoundaryOperator. Martin Almquist <malmquist@stanford.edu> parents: 968 diff changeset	605
726 37d5d69b1a0d Refactor BC code in Elastic2dVariable Martin Almquist <malmquist@stanford.edu> parents: 689 diff changeset	606 otherwise
958 72cd29107a9a Temporary changes in multiblock.DiffOp. Change traction operators in Elastic2dvariable to be true boundary operators. But adjoint FD conv test fails for dirichlet BC so need to debug! Martin Almquist <malmquist@stanford.edu> parents: 919 diff changeset	607 error(['No such operator: operator = ' op{k}]);
726 37d5d69b1a0d Refactor BC code in Elastic2dVariable Martin Almquist <malmquist@stanford.edu> parents: 689 diff changeset	608 end
687 e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	609 end
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	610
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	611 end
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	612
970 23d9ca6755be Add getBoundaryQuadrature in Elastic2dVariable. Rename get_boundary_operator -> getBoundaryOperator. Add operators in getBoundaryOperator, with full size so that they work with multiblock.DiffOp.getBoundaryOperator. Martin Almquist <malmquist@stanford.edu> parents: 968 diff changeset	613 function H = getBoundaryQuadrature(obj, boundary)
23d9ca6755be Add getBoundaryQuadrature in Elastic2dVariable. Rename get_boundary_operator -> getBoundaryOperator. Add operators in getBoundaryOperator, with full size so that they work with multiblock.DiffOp.getBoundaryOperator. Martin Almquist <malmquist@stanford.edu> parents: 968 diff changeset	614 switch boundary
23d9ca6755be Add getBoundaryQuadrature in Elastic2dVariable. Rename get_boundary_operator -> getBoundaryOperator. Add operators in getBoundaryOperator, with full size so that they work with multiblock.DiffOp.getBoundaryOperator. Martin Almquist <malmquist@stanford.edu> parents: 968 diff changeset	615 case {'w','W','west','West', 'e', 'E', 'east', 'East'}
23d9ca6755be Add getBoundaryQuadrature in Elastic2dVariable. Rename get_boundary_operator -> getBoundaryOperator. Add operators in getBoundaryOperator, with full size so that they work with multiblock.DiffOp.getBoundaryOperator. Martin Almquist <malmquist@stanford.edu> parents: 968 diff changeset	616 j = 1;
23d9ca6755be Add getBoundaryQuadrature in Elastic2dVariable. Rename get_boundary_operator -> getBoundaryOperator. Add operators in getBoundaryOperator, with full size so that they work with multiblock.DiffOp.getBoundaryOperator. Martin Almquist <malmquist@stanford.edu> parents: 968 diff changeset	617 case {'s','S','south','South', 'n', 'N', 'north', 'North'}
23d9ca6755be Add getBoundaryQuadrature in Elastic2dVariable. Rename get_boundary_operator -> getBoundaryOperator. Add operators in getBoundaryOperator, with full size so that they work with multiblock.DiffOp.getBoundaryOperator. Martin Almquist <malmquist@stanford.edu> parents: 968 diff changeset	618 j = 2;
23d9ca6755be Add getBoundaryQuadrature in Elastic2dVariable. Rename get_boundary_operator -> getBoundaryOperator. Add operators in getBoundaryOperator, with full size so that they work with multiblock.DiffOp.getBoundaryOperator. Martin Almquist <malmquist@stanford.edu> parents: 968 diff changeset	619 otherwise
23d9ca6755be Add getBoundaryQuadrature in Elastic2dVariable. Rename get_boundary_operator -> getBoundaryOperator. Add operators in getBoundaryOperator, with full size so that they work with multiblock.DiffOp.getBoundaryOperator. Martin Almquist <malmquist@stanford.edu> parents: 968 diff changeset	620 error('No such boundary: boundary = %s',boundary);
23d9ca6755be Add getBoundaryQuadrature in Elastic2dVariable. Rename get_boundary_operator -> getBoundaryOperator. Add operators in getBoundaryOperator, with full size so that they work with multiblock.DiffOp.getBoundaryOperator. Martin Almquist <malmquist@stanford.edu> parents: 968 diff changeset	621 end
23d9ca6755be Add getBoundaryQuadrature in Elastic2dVariable. Rename get_boundary_operator -> getBoundaryOperator. Add operators in getBoundaryOperator, with full size so that they work with multiblock.DiffOp.getBoundaryOperator. Martin Almquist <malmquist@stanford.edu> parents: 968 diff changeset	622 H = obj.H_boundary{j};
23d9ca6755be Add getBoundaryQuadrature in Elastic2dVariable. Rename get_boundary_operator -> getBoundaryOperator. Add operators in getBoundaryOperator, with full size so that they work with multiblock.DiffOp.getBoundaryOperator. Martin Almquist <malmquist@stanford.edu> parents: 968 diff changeset	623 I_dim = speye(obj.dim, obj.dim);
23d9ca6755be Add getBoundaryQuadrature in Elastic2dVariable. Rename get_boundary_operator -> getBoundaryOperator. Add operators in getBoundaryOperator, with full size so that they work with multiblock.DiffOp.getBoundaryOperator. Martin Almquist <malmquist@stanford.edu> parents: 968 diff changeset	624 H = kron(H, I_dim);
23d9ca6755be Add getBoundaryQuadrature in Elastic2dVariable. Rename get_boundary_operator -> getBoundaryOperator. Add operators in getBoundaryOperator, with full size so that they work with multiblock.DiffOp.getBoundaryOperator. Martin Almquist <malmquist@stanford.edu> parents: 968 diff changeset	625 end
23d9ca6755be Add getBoundaryQuadrature in Elastic2dVariable. Rename get_boundary_operator -> getBoundaryOperator. Add operators in getBoundaryOperator, with full size so that they work with multiblock.DiffOp.getBoundaryOperator. Martin Almquist <malmquist@stanford.edu> parents: 968 diff changeset	626
687 e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	627 function N = size(obj)
734 eebe24a636c7 Make Elastic2dVariable.size account for components. Martin Almquist <malmquist@stanford.edu> parents: 730 diff changeset	628 N = obj.dim*prod(obj.m);
687 e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	629 end
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	630 end
e8fc3aa1faf6 Rename elastic scheme. Martin Almquist <malmquist@stanford.edu> parents: diff changeset	631 end

Mercurial > repos > public > sbplib

annotate +scheme/Elastic2dVariable.m @ 972:104f0af001e0 feature/poroelastic