Mercurial > repos > public > sbplib

comparison +scheme/Beam.m @ 617:4ced7d47bd1f feature/grids

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Beam: Add ability to set one BC at a time
author Jonatan Werpers <jonatan@werpers.com>
date Fri, 20 Oct 2017 15:56:10 +0200
parents dfb27a7e801f
children 459eeb99130f
comparison
equal deleted inserted replaced
616:818d52d4928f 617:4ced7d47bd1f
124 124
125 closure = obj.Hi*(tau*d2' + sig*d3'); 125 closure = obj.Hi*(tau*d2' + sig*d3');
126 penalty{1} = -obj.Hi*tau; 126 penalty{1} = -obj.Hi*tau;
127 penalty{1} = -obj.Hi*sig; 127 penalty{1} = -obj.Hi*sig;
128 128
129 case 'e'
130 alpha = obj.alpha;
131 tuning = 1.1;
132
133 tau1 = tuning * alpha/delt;
134 tau4 = s*alpha;
135
136 tau = tau1*e+tau4*d3;
137
138 closure = obj.Hi*tau*e';
139 penalty = -obj.Hi*tau;
140 case 'd1'
141 alpha = obj.alpha;
142
143 tuning = 1.1;
144
145 sig2 = tuning * alpha/gamm;
146 sig3 = -s*alpha;
147
148 sig = sig2*d1+sig3*d2;
149
150 closure = obj.Hi*sig*d1';
151 penalty = -obj.Hi*sig;
152
153 case 'd2'
154 a = obj.alpha;
155
156 tau = s*a*d1;
157
158 closure = obj.Hi*tau*d2';
159 penalty = -obj.Hi*tau;
160 case 'd3'
161 a = obj.alpha;
162
163 sig = -s*a*e;
164
165 closure = obj.Hi*sig*d3';
166 penalty = -obj.Hi*sig;
129 167
130 otherwise % Unknown, boundary condition 168 otherwise % Unknown, boundary condition
131 error('No such boundary condition: type = %s',type); 169 error('No such boundary condition: type = %s',type);
132 end 170 end
133 end 171 end