Mercurial > repos > public > sbplib

annotate +scheme/Hypsyst3dCurve.m @ 351:7cc3d5bd3692 feature/hypsyst

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Curve linear code in 3D converges for manufactures solution
author Ylva Rydin <ylva.rydin@telia.com>
date Thu, 03 Nov 2016 22:03:09 +0100
parents 5d5652fe826a
children dbac99d2c318
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
350
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
1 classdef Hypsyst3dCurve < scheme.Scheme
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
2 properties
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
3 m % Number of points in each direction, possibly a vector
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
4 n %size of system
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
5 h % Grid spacing
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
6 X, Y, Z% Values of x and y for each grid point
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
7 Yx, Zx, Xy, Zy, Xz, Yz %Grid values for boundary surfaces
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
8
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
9 xi,eta,zeta
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
10 Xi, Eta, Zeta
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
11
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
12 Eta_xi, Zeta_xi, Xi_eta, Zeta_eta, Xi_zeta, Eta_zeta
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
13
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
14 X_xi, X_eta, X_zeta,Y_xi,Y_eta,Y_zeta,Z_xi,Z_eta,Z_zeta
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
15 Aev
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
16
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
17 metric_terms
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
18
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
19 order % Order accuracy for the approximation
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
20
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
21 D % non-stabalized scheme operator
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
22 Aevaluated, Bevaluated, Cevaluated, Eevaluated
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
23 Ahat, Bhat, Chat, E
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
24 A,B,C
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
25
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
26 J, Ji
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
27
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
28 H % Discrete norm
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
29 % Norms in the x, y and z directions
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
30 Hxii,Hetai,Hzetai, Hzi % Kroneckerd norms. 1'*Hx*v corresponds to integration in the x dir.
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
31 I_xi,I_eta,I_zeta, I_N,onesN
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
32 e_w, e_e, e_s, e_n, e_b, e_t
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
33 index_w, index_e,index_s,index_n, index_b, index_t
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
34 params %parameters for the coeficient matrice
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
35 end
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
36
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
37
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
38 methods
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
39 function obj = Hypsyst3dCurve(m, order, A, B,C, E, params,ti)
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
40 xilim ={0 1};
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
41 etalim = {0 1};
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
42 zetalim = {0 1};
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
43
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
44 if length(m) == 1
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
45 m = [m m m];
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
46 end
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
47 m_xi = m(1);
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
48 m_eta = m(2);
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
49 m_zeta=m(3);
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
50 m_tot=m_xi*m_eta*m_zeta;
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
51 obj.params = params;
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
52 obj.n = length(A(obj,0,0,0));
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
53
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
54 obj.m=m;
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
55
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
56 obj.order=order;
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
57 obj.onesN=ones(obj.n);
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
58 ops_xi = sbp.D2Standard(m_xi,xilim,order);
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
59 ops_eta = sbp.D2Standard(m_eta,etalim,order);
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
60 ops_zeta = sbp.D2Standard(m_zeta,zetalim,order);
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
61
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
62 obj.xi = ops_xi.x;
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
63 obj.eta = ops_eta.x;
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
64 obj.zeta = ops_zeta.x;
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
65
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
66 obj.Xi = kr(obj.xi,ones(m_eta,1),ones(m_zeta,1));%% Que pasa?
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
67 obj.Eta = kr(ones(m_xi,1),obj.eta,ones(m_zeta,1));
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
68 obj.Zeta = kr(ones(m_xi,1),ones(m_eta,1),obj.zeta);
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
69
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
70 obj.Eta_xi=kr(obj.eta,ones(m_xi,1));
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
71 obj.Zeta_xi=kr(ones(m_eta,1),obj.zeta);
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
72
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
73 obj.Xi_eta=kr(obj.xi,ones(m_zeta,1));
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
74 obj.Zeta_eta=kr(ones(m_xi,1),obj.zeta);
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
75
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
76 obj.Xi_zeta=kr(obj.xi,ones(m_eta,1));
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
77 obj.Eta_zeta=kr(ones(m_zeta,1),obj.eta);
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
78
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
79 [X,Y,Z] = ti.map(obj.Xi,obj.Eta,obj.Zeta);
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
80 obj.X=X;
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
81 obj.Y=Y;
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
82 obj.Z=Z;
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
83
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
84 I_n = eye(obj.n);
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
85 I_xi = speye(m_xi);
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
86 obj.I_xi = I_xi;
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
87 I_eta = speye(m_eta);
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
88 obj.I_eta = I_eta;
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
89 I_zeta = speye(m_zeta);
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
90 obj.I_zeta = I_zeta;
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
91
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
92
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
93 O_xi=ones(m_xi,1);
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
94 O_eta=ones(m_eta,1);
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
95 O_zeta=ones(m_zeta,1);
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
96
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
97 D1_xi = kr(ops_xi.D1, I_eta,I_zeta);
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
98 obj.Hxii = kr(I_n, ops_xi.HI, I_eta,I_zeta);
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
99 D1_eta = kr(I_xi, ops_eta.D1,I_zeta);
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
100 obj.Hetai = kr(I_n, I_xi, ops_eta.HI,I_zeta);
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
101 D1_zeta = kr(I_xi, I_eta,ops_zeta.D1);
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
102 obj.Hzetai = kr(I_n, I_xi,I_eta, ops_zeta.HI);
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
103 obj.h=[ops_xi.h ops_eta.h ops_zeta.h];
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
104
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
105 obj.e_w = kr(I_n, ops_xi.e_l, I_eta,I_zeta);
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
106 obj.e_e = kr(I_n, ops_xi.e_r, I_eta,I_zeta);
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
107 obj.e_s = kr(I_n, I_xi, ops_eta.e_l,I_zeta);
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
108 obj.e_n = kr(I_n, I_xi, ops_eta.e_r,I_zeta);
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
109 obj.e_b = kr(I_n, I_xi, I_eta, ops_zeta.e_l);
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
110 obj.e_t = kr(I_n, I_xi, I_eta, ops_zeta.e_r);
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
111
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
112 obj.A=A;
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
113 obj.B=B;
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
114 obj.C=C;
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
115
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
116 obj.X_xi=D1_xi*X;
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
117 obj.X_eta=D1_eta*X;
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
118 obj.X_zeta=D1_zeta*X;
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
119 obj.Y_xi=D1_xi*Y;
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
120 obj.Y_eta=D1_eta*Y;
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
121 obj.Y_zeta=D1_zeta*Y;
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
122 obj.Z_xi=D1_xi*Z;
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
123 obj.Z_eta=D1_eta*Z;
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
124 obj.Z_zeta=D1_zeta*Z;
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
125
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
126 D1_xi=kr(I_n,D1_xi);
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
127 D1_eta=kr(I_n,D1_eta);
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
128 D1_zeta=kr(I_n,D1_zeta);
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
129
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
130 obj.index_w=(kr(ops_xi.e_l, O_eta,O_zeta)==1);
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
131 obj.index_e=(kr(ops_xi.e_r, O_eta,O_zeta)==1);
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
132 obj.index_s=(kr(O_xi, ops_eta.e_l,O_zeta)==1);
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
133 obj.index_n=(kr(O_xi, ops_eta.e_r,O_zeta)==1);
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
134 obj.index_b=(kr(O_xi, O_eta, ops_zeta.e_l)==1);
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
135 obj.index_t=(kr(O_xi, O_eta, ops_zeta.e_r)==1);
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
136
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
137
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
138 obj.Ahat=@transform_coefficient_matrix;
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
139 obj.Bhat=@transform_coefficient_matrix;
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
140 obj.Chat=@transform_coefficient_matrix;
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
141 obj.E=@(obj,x,y,z,~,~,~,~,~,~)E(obj,x,y,z);
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
142
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
143 obj.Aevaluated = obj.evaluateCoefficientMatrix(obj.Ahat,obj.X, obj.Y,obj.Z, obj.X_eta,obj.X_zeta,obj.Y_eta,obj.Y_zeta,obj.Z_eta,obj.Z_zeta);
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
144 obj.Bevaluated = obj.evaluateCoefficientMatrix(obj.Bhat,obj.X, obj.Y,obj.Z, obj.X_zeta,obj.X_xi,obj.Y_zeta,obj.Y_xi,obj.Z_zeta,obj.Z_xi);
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
145 obj.Cevaluated = obj.evaluateCoefficientMatrix(obj.Chat,obj.X,obj.Y,obj.Z, obj.X_xi,obj.X_eta,obj.Y_xi,obj.Y_eta,obj.Z_xi,obj.Z_eta);
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
146 obj.Eevaluated = obj.evaluateCoefficientMatrix(obj.E, obj.X, obj.Y,obj.Z,[],[],[],[],[],[]);
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
147
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
148 obj.J=obj.X_xi.*obj.Y_eta.*obj.Z_zeta...
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
149 +obj.X_zeta.*obj.Y_xi.*obj.Z_eta...
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
150 +obj.X_eta.*obj.Y_zeta.*obj.Z_xi...
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
151 -obj.X_xi.*obj.Y_zeta.*obj.Z_eta...
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
152 -obj.X_eta.*obj.Y_xi.*obj.Z_zeta...
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
153 -obj.X_zeta.*obj.Y_eta.*obj.Z_xi;
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
154
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
155 obj.Ji =kr(I_n,spdiags(1./obj.J,0,m_tot,m_tot));
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
156
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
157 obj.D=obj.Ji*(-obj.Aevaluated*D1_xi-obj.Bevaluated*D1_eta -obj.Cevaluated*D1_zeta)-obj.Eevaluated;
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
158 end
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
159
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
160 function [ret]=transform_coefficient_matrix(obj,x,y,z,x_1,x_2,y_1,y_2,z_1,z_2)
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
161 ret=obj.A(obj,x,y,z).*(y_1.*z_2-z_1.*y_2);
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
162 ret=ret+obj.B(obj,x,y,z).*(x_2.*z_1-x_1.*z_2);
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
163 ret=ret+obj.C(obj,x,y,z).*(x_1.*y_2-x_2.*y_1);
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
164 end
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
165
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
166
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
167 % Closure functions return the opertors applied to the own doamin to close the boundary
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
168 % Penalty functions return the opertors to force the solution. In the case of an interface it returns the operator applied to the other doamin.
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
169 % boundary is a string specifying the boundary e.g. 'l','r' or 'e','w','n','s'.
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
170 % type is a string specifying the type of boundary condition if there are several.
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
171 % data is a function returning the data that should be applied at the boundary.
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
172 function [closure, penalty] = boundary_condition(obj,boundary,type,L)
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
173 default_arg('type','char');
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
174 BM=boundary_matrices(obj,boundary);
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
175
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
176 switch type
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
177 case{'c','char'}
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
178 [closure,penalty]=boundary_condition_char(obj,BM);
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
179 case{'general'}
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
180 [closure,penalty]=boundary_condition_general(obj,BM,boundary,L);
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
181 otherwise
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
182 error('No such boundary condition')
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
183 end
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
184 end
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
185
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
186 function [closure, penalty] = interface(obj,boundary,neighbour_scheme,neighbour_boundary)
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
187 error('An interface function does not exist yet');
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
188 end
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
189
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
190 function N = size(obj)
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
191 N = obj.m;
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
192 end
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
193
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
194 function [ret] = evaluateCoefficientMatrix(obj,mat, X, Y, Z , x_1 , x_2 , y_1 , y_2 , z_1 , z_2)
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
195 params=obj.params;
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
196 side=max(length(X),length(Y));
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
197 if isa(mat,'function_handle')
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
198 [rows,cols]=size(mat(obj,0,0,0,0,0,0,0,0,0));
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
199 x_1=kr(obj.onesN,x_1);
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
200 x_2=kr(obj.onesN,x_2);
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
201 y_1=kr(obj.onesN,y_1);
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
202 y_2=kr(obj.onesN,y_2);
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
203 z_1=kr(obj.onesN,z_1);
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
204 z_2=kr(obj.onesN,z_2);
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
205 matVec=mat(obj,X',Y',Z',x_1',x_2',y_1',y_2',z_1',z_2');
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
206 matVec=sparse(matVec);
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
207 else
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
208 matVec=mat;
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
209 [rows,cols]=size(matVec);
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
210 side=max(length(X),length(Y));
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
211 cols=cols/side;
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
212 end
351
7cc3d5bd3692 Curve linear code in 3D converges for manufactures solution
Ylva Rydin <ylva.rydin@telia.com>
parents: 350
diff changeset
213 ret=cell(rows,cols);
7cc3d5bd3692 Curve linear code in 3D converges for manufactures solution
Ylva Rydin <ylva.rydin@telia.com>
parents: 350
diff changeset
214
350
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
215
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
216 for ii=1:rows
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
217 for jj=1:cols
351
7cc3d5bd3692 Curve linear code in 3D converges for manufactures solution
Ylva Rydin <ylva.rydin@telia.com>
parents: 350
diff changeset
218 ret{ii,jj}=diag(matVec(ii,(jj-1)*side+1:jj*side));
350
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
219 end
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
220 end
351
7cc3d5bd3692 Curve linear code in 3D converges for manufactures solution
Ylva Rydin <ylva.rydin@telia.com>
parents: 350
diff changeset
221
7cc3d5bd3692 Curve linear code in 3D converges for manufactures solution
Ylva Rydin <ylva.rydin@telia.com>
parents: 350
diff changeset
222 ret=cell2mat(ret);
350
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
223 end
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
224
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
225
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
226 function [BM]=boundary_matrices(obj,boundary)
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
227 params=obj.params;
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
228 BM.boundary=boundary;
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
229 switch boundary
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
230 case {'w','W','west'}
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
231 BM.e_=obj.e_w;
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
232 mat=obj.Ahat;
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
233 BM.boundpos='l';
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
234 BM.Hi=obj.Hxii;
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
235 BM.index=obj.index_w;
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
236 BM.x_1=obj.X_eta(BM.index);
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
237 BM.x_2=obj.X_zeta(BM.index);
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
238 BM.y_1=obj.Y_eta(BM.index);
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
239 BM.y_2=obj.Y_zeta(BM.index);
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
240 BM.z_1=obj.Z_eta(BM.index);
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
241 BM.z_2=obj.Z_zeta(BM.index);
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
242 case {'e','E','east'}
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
243 BM.e_=obj.e_e;
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
244 mat=obj.Ahat;
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
245 BM.boundpos='r';
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
246 BM.Hi=obj.Hxii;
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
247 BM.index=obj.index_e;
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
248 BM.x_1=obj.X_eta(BM.index);
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
249 BM.x_2=obj.X_zeta(BM.index);
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
250 BM.y_1=obj.Y_eta(BM.index);
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
251 BM.y_2=obj.Y_zeta(BM.index);
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
252 BM.z_1=obj.Z_eta(BM.index);
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
253 BM.z_2=obj.Z_zeta(BM.index);
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
254 case {'s','S','south'}
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
255 BM.e_=obj.e_s;
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
256 mat=obj.Bhat;
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
257 BM.boundpos='l';
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
258 BM.Hi=obj.Hetai;
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
259 BM.index=obj.index_s;
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
260 BM.x_1=obj.X_zeta(BM.index);
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
261 BM.x_2=obj.X_xi(BM.index);
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
262 BM.y_1=obj.Y_zeta(BM.index);
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
263 BM.y_2=obj.Y_xi(BM.index);
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
264 BM.z_1=obj.Z_zeta(BM.index);
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
265 BM.z_2=obj.Z_xi(BM.index);
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
266 case {'n','N','north'}
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
267 BM.e_=obj.e_n;
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
268 mat=obj.Bhat;
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
269 BM.boundpos='r';
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
270 BM.Hi=obj.Hetai;
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
271 BM.index=obj.index_n;
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
272 BM.x_1=obj.X_zeta(BM.index);
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
273 BM.x_2=obj.X_xi(BM.index);
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
274 BM.y_1=obj.Y_zeta(BM.index);
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
275 BM.y_2=obj.Y_xi(BM.index);
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
276 BM.z_1=obj.Z_zeta(BM.index);
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
277 BM.z_2=obj.Z_xi(BM.index);
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
278 case{'b','B','Bottom'}
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
279 BM.e_=obj.e_b;
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
280 mat=obj.Chat;
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
281 BM.boundpos='l';
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
282 BM.Hi=obj.Hzetai;
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
283 BM.index=obj.index_b;
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
284 BM.x_1=obj.X_xi(BM.index);
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
285 BM.x_2=obj.X_eta(BM.index);
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
286 BM.y_1=obj.Y_xi(BM.index);
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
287 BM.y_2=obj.Y_eta(BM.index);
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
288 BM.z_1=obj.Z_xi(BM.index);
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
289 BM.z_2=obj.Z_eta(BM.index);
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
290 case{'t','T','Top'}
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
291 BM.e_=obj.e_t;
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
292 mat=obj.Chat;
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
293 BM.boundpos='r';
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
294 BM.Hi=obj.Hzetai;
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
295 BM.index=obj.index_t;
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
296 BM.x_1=obj.X_xi(BM.index);
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
297 BM.x_2=obj.X_eta(BM.index);
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
298 BM.y_1=obj.Y_xi(BM.index);
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
299 BM.y_2=obj.Y_eta(BM.index);
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
300 BM.z_1=obj.Z_xi(BM.index);
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
301 BM.z_2=obj.Z_eta(BM.index);
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
302 end
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
303 [BM.V,BM.Vi,BM.D,signVec]=obj.matrixDiag(mat,obj.X(BM.index),obj.Y(BM.index),obj.Z(BM.index),...
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
304 BM.x_1,BM.x_2,BM.y_1,BM.y_2,BM.z_1,BM.z_2);
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
305 BM.side=sum(BM.index);
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
306 BM.pos=signVec(1); BM.zeroval=signVec(2); BM.neg=signVec(3);
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
307 end
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
308
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
309
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
310 function [closure, penalty]=boundary_condition_char(obj,BM)
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
311 side = BM.side;
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
312 pos = BM.pos;
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
313 neg = BM.neg;
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
314 zeroval=BM.zeroval;
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
315 V = BM.V;
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
316 Vi = BM.Vi;
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
317 Hi=BM.Hi;
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
318 D=BM.D;
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
319 e_=BM.e_;
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
320
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
321 switch BM.boundpos
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
322 case {'l'}
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
323 tau=sparse(obj.n*side,pos);
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
324 Vi_plus=Vi(1:pos,:);
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
325 tau(1:pos,:)=-abs(D(1:pos,1:pos));
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
326 closure=Hi*e_*V*tau*Vi_plus*e_';
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
327 penalty=-Hi*e_*V*tau*Vi_plus;
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
328 case {'r'}
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
329 tau=sparse(obj.n*side,neg);
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
330 tau((pos+zeroval)+1:obj.n*side,:)=-abs(D((pos+zeroval)+1:obj.n*side,(pos+zeroval)+1:obj.n*side));
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
331 Vi_minus=Vi((pos+zeroval)+1:obj.n*side,:);
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
332 closure=Hi*e_*V*tau*Vi_minus*e_';
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
333 penalty=-Hi*e_*V*tau*Vi_minus;
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
334 end
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
335 end
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
336
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
337
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
338 function [closure,penalty]=boundary_condition_general(obj,BM,boundary,L)
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
339 side = BM.side;
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
340 pos = BM.pos;
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
341 neg = BM.neg;
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
342 zeroval=BM.zeroval;
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
343 V = BM.V;
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
344 Vi = BM.Vi;
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
345 Hi=BM.Hi;
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
346 D=BM.D;
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
347 e_=BM.e_;
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
348 index=BM.index;
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
349
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
350 switch BM.boundary
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
351 case{'b','B','bottom'}
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
352 Ji_vec=diag(obj.Ji);
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
353 Ji=diag(Ji_vec(index));
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
354 Zeta_x=Ji*(obj.Y_xi(index).*obj.Z_eta(index)-obj.Z_xi(index).*obj.Y_eta(index));
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
355 Zeta_y=Ji*(obj.X_eta(index).*obj.Z_xi(index)-obj.X_xi(index).*obj.Z_eta(index));
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
356 Zeta_z=Ji*(obj.X_xi(index).*obj.Y_eta(index)-obj.Y_xi(index).*obj.X_eta(index));
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
357
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
358 L=obj.evaluateCoefficientMatrix(L,Zeta_x,Zeta_y,Zeta_z,[],[],[],[],[],[]);
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
359 end
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
360
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
361 switch BM.boundpos
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
362 case {'l'}
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
363 tau=sparse(obj.n*side,pos);
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
364 Vi_plus=Vi(1:pos,:);
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
365 Vi_minus=Vi(pos+zeroval+1:obj.n*side,:);
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
366 V_plus=V(:,1:pos);
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
367 V_minus=V(:,(pos+zeroval)+1:obj.n*side);
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
368
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
369 tau(1:pos,:)=-abs(D(1:pos,1:pos));
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
370 R=-inv(L*V_plus)*(L*V_minus);
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
371 closure=Hi*e_*V*tau*(Vi_plus-R*Vi_minus)*e_';
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
372 penalty=-Hi*e_*V*tau*inv(L*V_plus)*L;
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
373 case {'r'}
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
374 tau=sparse(obj.n*side,neg);
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
375 tau((pos+zeroval)+1:obj.n*side,:)=-abs(D((pos+zeroval)+1:obj.n*side,(pos+zeroval)+1:obj.n*side));
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
376 Vi_plus=Vi(1:pos,:);
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
377 Vi_minus=Vi((pos+zeroval)+1:obj.n*side,:);
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
378
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
379 V_plus=V(:,1:pos);
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
380 V_minus=V(:,(pos+zeroval)+1:obj.n*side);
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
381 R=-inv(L*V_minus)*(L*V_plus);
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
382 closure=Hi*e_*V*tau*(Vi_minus-R*Vi_plus)*e_';
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
383 penalty=-Hi*e_*V*tau*inv(L*V_minus)*L;
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
384 end
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
385 end
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
386
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
387
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
388 function [V,Vi, D,signVec]=matrixDiag(obj,mat,x,y,z,x_1,x_2,y_1,y_2,z_1,z_2)
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
389 params=obj.params;
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
390 eps=10^(-10);
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
391 if(sum(abs(x_1))>eps)
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
392 syms x_1s
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
393 else
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
394 x_1s=0;
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
395 end
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
396
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
397 if(sum(abs(x_2))>eps)
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
398 syms x_2s;
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
399 else
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
400 x_2s=0;
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
401 end
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
402
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
403
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
404 if(sum(abs(y_1))>eps)
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
405 syms y_1s
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
406 else
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
407 y_1s=0;
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
408 end
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
409
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
410 if(sum(abs(y_2))>eps)
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
411 syms y_2s;
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
412 else
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
413 y_2s=0;
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
414 end
351
7cc3d5bd3692 Curve linear code in 3D converges for manufactures solution
Ylva Rydin <ylva.rydin@telia.com>
parents: 350
diff changeset
415
350
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
416 if(sum(abs(z_1))>eps)
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
417 syms z_1s
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
418 else
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
419 z_1s=0;
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
420 end
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
421
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
422 if(sum(abs(z_2))>eps)
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
423 syms z_2s;
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
424 else
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
425 z_2s=0;
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
426 end
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
427
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
428 syms xs ys zs
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
429 [V, D]=eig(mat(obj,xs,ys,zs,x_1s,x_2s,y_1s,y_2s,z_1s,z_2s));
351
7cc3d5bd3692 Curve linear code in 3D converges for manufactures solution
Ylva Rydin <ylva.rydin@telia.com>
parents: 350
diff changeset
430 Vi=inv(V);
7cc3d5bd3692 Curve linear code in 3D converges for manufactures solution
Ylva Rydin <ylva.rydin@telia.com>
parents: 350
diff changeset
431 % syms x_1s x_2s y_1s y_2s z_1s z_2s
7cc3d5bd3692 Curve linear code in 3D converges for manufactures solution
Ylva Rydin <ylva.rydin@telia.com>
parents: 350
diff changeset
432 xs=x;
7cc3d5bd3692 Curve linear code in 3D converges for manufactures solution
Ylva Rydin <ylva.rydin@telia.com>
parents: 350
diff changeset
433 ys=y;
7cc3d5bd3692 Curve linear code in 3D converges for manufactures solution
Ylva Rydin <ylva.rydin@telia.com>
parents: 350
diff changeset
434 zs=z;
7cc3d5bd3692 Curve linear code in 3D converges for manufactures solution
Ylva Rydin <ylva.rydin@telia.com>
parents: 350
diff changeset
435 x_1s=x_1;
7cc3d5bd3692 Curve linear code in 3D converges for manufactures solution
Ylva Rydin <ylva.rydin@telia.com>
parents: 350
diff changeset
436 x_2s=x_2;
7cc3d5bd3692 Curve linear code in 3D converges for manufactures solution
Ylva Rydin <ylva.rydin@telia.com>
parents: 350
diff changeset
437 y_1s=y_1;
7cc3d5bd3692 Curve linear code in 3D converges for manufactures solution
Ylva Rydin <ylva.rydin@telia.com>
parents: 350
diff changeset
438 y_2s=y_2;
7cc3d5bd3692 Curve linear code in 3D converges for manufactures solution
Ylva Rydin <ylva.rydin@telia.com>
parents: 350
diff changeset
439 z_1s=z_1;
7cc3d5bd3692 Curve linear code in 3D converges for manufactures solution
Ylva Rydin <ylva.rydin@telia.com>
parents: 350
diff changeset
440 z_2s=z_2;
350
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
441
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
442 side=max(length(x),length(y));
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
443 Dret=zeros(obj.n,side*obj.n);
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
444 Vret=zeros(obj.n,side*obj.n);
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
445 Viret=zeros(obj.n,side*obj.n);
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
446
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
447 for ii=1:obj.n
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
448 for jj=1:obj.n
351
7cc3d5bd3692 Curve linear code in 3D converges for manufactures solution
Ylva Rydin <ylva.rydin@telia.com>
parents: 350
diff changeset
449 Dret(jj,(ii-1)*side+1:side*ii)=eval(D(jj,ii));
7cc3d5bd3692 Curve linear code in 3D converges for manufactures solution
Ylva Rydin <ylva.rydin@telia.com>
parents: 350
diff changeset
450 Vret(jj,(ii-1)*side+1:side*ii)=eval(V(jj,ii));
7cc3d5bd3692 Curve linear code in 3D converges for manufactures solution
Ylva Rydin <ylva.rydin@telia.com>
parents: 350
diff changeset
451 Viret(jj,(ii-1)*side+1:side*ii)=eval(Vi(jj,ii));
350
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
452 end
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
453 end
351
7cc3d5bd3692 Curve linear code in 3D converges for manufactures solution
Ylva Rydin <ylva.rydin@telia.com>
parents: 350
diff changeset
454
350
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
455 D=sparse(Dret);
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
456 V=sparse(Vret);
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
457 Vi=sparse(Viret);
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
458 V=obj.evaluateCoefficientMatrix(V,x,y,z,x_1,x_2,y_1,y_2,z_1,z_2);
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
459 D=obj.evaluateCoefficientMatrix(D,x,y,z,x_1,x_2,y_1,y_2,z_1,z_2);
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
460 Vi=obj.evaluateCoefficientMatrix(Vi,x,y,z,x_1,x_2,y_1,y_2,z_1,z_2);
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
461 DD=diag(D);
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
462
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
463 poseig=(DD>0);
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
464 zeroeig=(DD==0);
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
465 negeig=(DD<0);
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
466
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
467 D=diag([DD(poseig); DD(zeroeig); DD(negeig)]);
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
468 V=[V(:,poseig) V(:,zeroeig) V(:,negeig)];
351
7cc3d5bd3692 Curve linear code in 3D converges for manufactures solution
Ylva Rydin <ylva.rydin@telia.com>
parents: 350
diff changeset
469 Vi=[Vi(poseig,:); Vi(zeroeig,:); Vi(negeig,:)];
350
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
470 signVec=[sum(poseig),sum(zeroeig),sum(negeig)];
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
471 end
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
472 end
5d5652fe826a A commit before I try resolving the performance issues
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
473 end