sbplib: +scheme/Hypsyst3d.m annotate

author	Ylva Rydin <ylva.rydin@telia.com>
date	Thu, 13 Oct 2016 09:34:30 +0200
parents
children	5d5652fe826a

rev	line source
349 cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	1 classdef Hypsyst3d < scheme.Scheme
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	2 properties
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	3 m % Number of points in each direction, possibly a vector
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	4 n %size of system
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	5 h % Grid spacing
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	6 x, y, z % Grid
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	7 X, Y, Z% Values of x and y for each grid point
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	8 Yx, Zx, Xy, Zy, Xz, Yz %Grid values for boundary surfaces
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	9 order % Order accuracy for the approximation
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	10
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	11 D % non-stabalized scheme operator
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	12 A, B, C, E
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	13
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	14 H % Discrete norm
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	15 % Norms in the x, y and z directions
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	16 Hxi,Hyi, Hzi % Kroneckerd norms. 1'Hxv corresponds to integration in the x dir.
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	17 I_x,I_y, I_z, I_N
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	18 e_w, e_e, e_s, e_n, e_b, e_t
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	19 params %parameters for the coeficient matrice
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	20 end
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	21
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	22
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	23 methods
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	24 function obj = Hypsyst3d(m, lim, order, A, B,C, E, params)
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	25 default_arg('E', [])
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	26 xlim = lim{1};
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	27 ylim = lim{2};
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	28 zlim = lim{3};
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	29
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	30 if length(m) == 1
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	31 m = [m m m];
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	32 end
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	33
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	34 obj.A=A;
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	35 obj.B=B;
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	36 obj.C=C;
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	37 obj.E=E;
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	38 m_x = m(1);
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	39 m_y = m(2);
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	40 m_z=m(3);
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	41 obj.params = params;
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	42
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	43 ops_x = sbp.D2Standard(m_x,xlim,order);
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	44 ops_y = sbp.D2Standard(m_y,ylim,order);
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	45 ops_z = sbp.D2Standard(m_z,zlim,order);
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	46
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	47 obj.x = ops_x.x;
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	48 obj.y = ops_y.x;
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	49 obj.z = ops_z.x;
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	50
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	51 obj.X = kr(obj.x,ones(m_y,1),ones(m_z,1));%% Que pasa?
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	52 obj.Y = kr(ones(m_x,1),obj.y,ones(m_z,1));
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	53 obj.Z = kr(ones(m_x,1),ones(m_y,1),obj.z);
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	54
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	55 obj.Yx=kr(obj.y,ones(m_z,1));
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	56 obj.Zx=kr(ones(m_y,1),obj.z);
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	57
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	58 obj.Xy=kr(obj.x,ones(m_z,1));
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	59 obj.Zy=kr(ones(m_x,1),obj.z);
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	60
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	61 obj.Xz=kr(obj.x,ones(m_y,1));
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	62 obj.Yz=kr(ones(m_z,1),obj.y);
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	63
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	64 Aevaluated = obj.evaluateCoefficientMatrix(A, obj.X, obj.Y,obj.Z);
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	65 Bevaluated = obj.evaluateCoefficientMatrix(B, obj.X, obj.Y,obj.Z);
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	66 Cevaluated = obj.evaluateCoefficientMatrix(C, obj.X, obj.Y,obj.Z);
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	67 Eevaluated = obj.evaluateCoefficientMatrix(E, obj.X, obj.Y,obj.Z);
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	68
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	69 obj.n = length(A(obj.params,0,0,0));
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	70
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	71 I_n = eye(obj.n);
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	72 I_x = speye(m_x);
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	73 obj.I_x = I_x;
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	74 I_y = speye(m_y);
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	75 obj.I_y = I_y;
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	76 I_z = speye(m_z);
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	77 obj.I_z = I_z;
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	78
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	79
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	80 D1_x = kr(I_n, ops_x.D1, I_y,I_z);
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	81 obj.Hxi = kr(I_n, ops_x.HI, I_y,I_z);
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	82 D1_y = kr(I_n, I_x, ops_y.D1,I_z);
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	83 obj.Hyi = kr(I_n, I_x, ops_y.HI,I_z);
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	84 D1_z = kr(I_n, I_x, I_y,ops_z.D1);
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	85 obj.Hzi = kr(I_n, I_x,I_y, ops_y.HI);
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	86
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	87 obj.e_w = kr(I_n, ops_x.e_l, I_y,I_z);
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	88 obj.e_e = kr(I_n, ops_x.e_r, I_y,I_z);
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	89 obj.e_s = kr(I_n, I_x, ops_y.e_l,I_z);
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	90 obj.e_n = kr(I_n, I_x, ops_y.e_r,I_z);
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	91 obj.e_b = kr(I_n, I_x, I_y, ops_z.e_l);
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	92 obj.e_t = kr(I_n, I_x, I_y, ops_z.e_r);
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	93
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	94 obj.m=m;
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	95 obj.h=[ops_x.h ops_y.h ops_x.h];
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	96 obj.order=order;
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	97
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	98 obj.D=-AevaluatedD1_x-BevaluatedD1_y-Cevaluated*D1_z-Eevaluated;
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	99 end
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	100
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	101 % Closure functions return the opertors applied to the own doamin to close the boundary
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	102 % Penalty functions return the opertors to force the solution. In the case of an interface it returns the operator applied to the other doamin.
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	103 % boundary is a string specifying the boundary e.g. 'l','r' or 'e','w','n','s'.
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	104 % type is a string specifying the type of boundary condition if there are several.
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	105 % data is a function returning the data that should be applied at the boundary.
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	106 function [closure, penalty] = boundary_condition(obj,boundary,type,L)
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	107 default_arg('type','char');
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	108 BM=boundary_matrices(obj,boundary);
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	109
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	110 switch type
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	111 case{'c','char'}
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	112 [closure,penalty]=boundary_condition_char(obj,BM);
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	113 case{'general'}
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	114 [closure,penalty]=boundary_condition_general(obj,BM,boundary,L);
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	115 otherwise
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	116 error('No such boundary condition')
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	117 end
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	118 end
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	119
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	120 function [closure, penalty] = interface(obj,boundary,neighbour_scheme,neighbour_boundary)
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	121 error('An interface function does not exist yet');
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	122 end
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	123
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	124 function N = size(obj)
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	125 N = obj.m;
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	126 end
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	127
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	128 function [ret] = evaluateCoefficientMatrix(obj, mat, X, Y, Z)
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	129 params=obj.params;
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	130 side=max(length(X),length(Y));
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	131 if isa(mat,'function_handle')
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	132 [rows,cols]=size(mat(params,0,0,0));
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	133 matVec=mat(params,X',Y',Z');
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	134 matVec=sparse(matVec);
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	135 else
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	136 matVec=mat;
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	137 [rows,cols]=size(matVec);
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	138 side=max(length(X),length(Y));
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	139 cols=cols/side;
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	140 end
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	141 ret=kron(ones(rows,cols),speye(side));
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	142
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	143 for ii=1:rows
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	144 for jj=1:cols
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	145 ret((ii-1)side+1:iiside,(jj-1)side+1:jjside)=diag(matVec(ii,(jj-1)side+1:jjside));
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	146 end
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	147 end
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	148 end
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	149
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	150
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	151 function [BM]=boundary_matrices(obj,boundary)
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	152 params=obj.params;
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	153
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	154 switch boundary
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	155 case {'w','W','west'}
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	156 BM.e_=obj.e_w;
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	157 mat=obj.A;
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	158 BM.boundpos='l';
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	159 BM.Hi=obj.Hxi;
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	160 [BM.V,BM.Vi,BM.D,signVec]=obj.matrixDiag(mat,obj.X(1),obj.Yx,obj.Zx);
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	161 BM.side=length(obj.Yx);
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	162 case {'e','E','east'}
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	163 BM.e_=obj.e_e;
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	164 mat=obj.A;
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	165 BM.boundpos='r';
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	166 BM.Hi=obj.Hxi;
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	167 [BM.V,BM.Vi,BM.D,signVec]=obj.matrixDiag(mat,obj.X(end),obj.Yx,obj.Zx);
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	168 BM.side=length(obj.Yx);
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	169 case {'s','S','south'}
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	170 BM.e_=obj.e_s;
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	171 mat=obj.B;
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	172 BM.boundpos='l';
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	173 BM.Hi=obj.Hyi;
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	174 [BM.V,BM.Vi,BM.D,signVec]=obj.matrixDiag(mat,obj.Xy,obj.Y(1),obj.Zy);
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	175 BM.side=length(obj.Xy);
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	176 case {'n','N','north'}
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	177 BM.e_=obj.e_n;
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	178 mat=obj.B;
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	179 BM.boundpos='r';
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	180 BM.Hi=obj.Hyi;
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	181 [BM.V,BM.Vi,BM.D,signVec]=obj.matrixDiag(mat,obj.Xy,obj.Y(end),obj.Zy);
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	182 BM.side=length(obj.Xy);
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	183 case{'b','B','Bottom'}
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	184 BM.e_=obj.e_b;
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	185 mat=obj.C;
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	186 BM.boundpos='l';
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	187 BM.Hi=obj.Hzi;
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	188 [BM.V,BM.Vi,BM.D,signVec]=obj.matrixDiag(mat,obj.Xz,obj.Yz,obj.Z(1));
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	189 BM.side=length(obj.Xz);
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	190 case{'t','T','Top'}
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	191 BM.e_=obj.e_t;
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	192 mat=obj.C;
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	193 BM.boundpos='r';
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	194 BM.Hi=obj.Hzi;
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	195 [BM.V,BM.Vi,BM.D,signVec]=obj.matrixDiag(mat,obj.Xz,obj.Yz,obj.Z(end));
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	196 BM.side=length(obj.Xz);
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	197 end
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	198
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	199 BM.pos=signVec(1); BM.zeroval=signVec(2); BM.neg=signVec(3);
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	200 end
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	201
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	202
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	203 function [closure, penalty]=boundary_condition_char(obj,BM)
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	204 side = BM.side;
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	205 pos = BM.pos;
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	206 neg = BM.neg;
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	207 zeroval=BM.zeroval;
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	208 V = BM.V;
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	209 Vi = BM.Vi;
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	210 Hi=BM.Hi;
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	211 D=BM.D;
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	212 e_=BM.e_;
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	213
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	214 switch BM.boundpos
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	215 case {'l'}
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	216 tau=sparse(obj.n*side,pos);
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	217 Vi_plus=Vi(1:pos,:);
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	218 tau(1:pos,:)=-abs(D(1:pos,1:pos));
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	219 closure=Hie_VtauVi_plus*e_';
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	220 penalty=-Hie_VtauVi_plus;
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	221 case {'r'}
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	222 tau=sparse(obj.n*side,neg);
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	223 tau((pos+zeroval)+1:obj.nside,:)=-abs(D((pos+zeroval)+1:obj.nside,(pos+zeroval)+1:obj.n*side));
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	224 Vi_minus=Vi((pos+zeroval)+1:obj.n*side,:);
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	225 closure=Hie_VtauVi_minus*e_';
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	226 penalty=-Hie_VtauVi_minus;
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	227 end
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	228 end
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	229
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	230
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	231 function [closure,penalty]=boundary_condition_general(obj,BM,boundary,L)
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	232 side = BM.side;
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	233 pos = BM.pos;
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	234 neg = BM.neg;
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	235 zeroval=BM.zeroval;
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	236 V = BM.V;
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	237 Vi = BM.Vi;
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	238 Hi=BM.Hi;
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	239 D=BM.D;
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	240 e_=BM.e_;
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	241 switch boundary
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	242 case {'w','W','west'}
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	243 L=obj.evaluateCoefficientMatrix(L,obj.x(1),obj.Yx,obj.Zx);
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	244 case {'e','E','east'}
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	245 L=obj.evaluateCoefficientMatrix(L,obj.x(end),obj.Yx,obj.Zx);
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	246 case {'s','S','south'}
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	247 L=obj.evaluateCoefficientMatrix(L,obj.Xy,obj.y(1),obj.Zy);
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	248 case {'n','N','north'}
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	249 L=obj.evaluateCoefficientMatrix(L,obj.Xy,obj.y(end),obj.Zy);
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	250 case {'b','B','bottom'}
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	251 L=obj.evaluateCoefficientMatrix(L,obj.Xz,obj.Yz,obj.z(1));
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	252 case {'t','T','top'}
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	253 L=obj.evaluateCoefficientMatrix(L,obj.Xz,obj.Yz,obj.z(end));
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	254 end
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	255
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	256 switch BM.boundpos
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	257 case {'l'}
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	258 tau=sparse(obj.n*side,pos);
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	259 Vi_plus=Vi(1:pos,:);
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	260 Vi_minus=Vi(pos+zeroval+1:obj.n*side,:);
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	261 V_plus=V(:,1:pos);
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	262 V_minus=V(:,(pos+zeroval)+1:obj.n*side);
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	263
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	264 tau(1:pos,:)=-abs(D(1:pos,1:pos));
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	265 R=-inv(LV_plus)(L*V_minus);
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	266 closure=Hie_Vtau(Vi_plus-RVi_minus)e_';
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	267 penalty=-Hie_Vtauinv(LV_plus)L;
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	268 case {'r'}
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	269 tau=sparse(obj.n*side,neg);
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	270 tau((pos+zeroval)+1:obj.nside,:)=-abs(D((pos+zeroval)+1:obj.nside,(pos+zeroval)+1:obj.n*side));
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	271 Vi_plus=Vi(1:pos,:);
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	272 Vi_minus=Vi((pos+zeroval)+1:obj.n*side,:);
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	273
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	274 V_plus=V(:,1:pos);
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	275 V_minus=V(:,(pos+zeroval)+1:obj.n*side);
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	276 R=-inv(LV_minus)(L*V_plus);
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	277 closure=Hie_Vtau(Vi_minus-RVi_plus)e_';
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	278 penalty=-Hie_Vtauinv(LV_minus)L;
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	279 end
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	280 end
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	281
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	282
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	283 function [V,Vi, D,signVec]=matrixDiag(obj,mat,x,y,z)
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	284 params=obj.params;
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	285 syms xs ys zs
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	286 [V, D]=eig(mat(params,xs,ys,zs));
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	287 xs=x;
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	288 ys=y;
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	289 zs=z;
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	290
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	291
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	292 side=max(length(x),length(y));
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	293 Dret=zeros(obj.n,side*obj.n);
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	294 Vret=zeros(obj.n,side*obj.n);
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	295 for ii=1:obj.n
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	296 for jj=1:obj.n
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	297 Dret(jj,(ii-1)side+1:sideii)=eval(D(jj,ii));
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	298 Vret(jj,(ii-1)side+1:sideii)=eval(V(jj,ii));
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	299 end
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	300 end
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	301
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	302 D=sparse(Dret);
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	303 V=sparse(Vret);
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	304 V=obj.evaluateCoefficientMatrix(V,x,y,z);
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	305 D=obj.evaluateCoefficientMatrix(D,x,y,z);
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	306 DD=diag(D);
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	307
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	308 poseig=(DD>0);
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	309 zeroeig=(DD==0);
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	310 negeig=(DD<0);
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	311
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	312 D=diag([DD(poseig); DD(zeroeig); DD(negeig)]);
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	313 V=[V(:,poseig) V(:,zeroeig) V(:,negeig)];
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	314 Vi=inv(V);
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	315 signVec=[sum(poseig),sum(zeroeig),sum(negeig)];
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	316 end
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	317 end
cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	318 end

349

cd6a29ab3746 A 3D is added and an attempt to imlement 3D transfinit interpolation has been initialized