Mercurial > repos > public > sbplib

comparison +scheme/Schrodinger2dCurve.m @ 710:023a95f63950 feature/quantumTriangles

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Remove stubid print-out
author Ylva Rydin <ylva.rydin@telia.com>
date Tue, 21 Nov 2017 18:12:00 +0100
parents f004b9e9d17a
children
comparison
equal deleted inserted replaced
709:f004b9e9d17a 710:023a95f63950
211 % neighbour_boundary is a string specifying which boundary to interface to. 211 % neighbour_boundary is a string specifying which boundary to interface to.
212 function [closure, penalty,closureHamiltonian,penaltyHamiltonian] = boundary_condition(obj, boundary,~) 212 function [closure, penalty,closureHamiltonian,penaltyHamiltonian] = boundary_condition(obj, boundary,~)
213 [e, d_n, d_t, coeff_t, coeff_n s, halfnorm_inv_n, halfnorm_inv_t, halfnorm_t,g] = obj.get_boundary_ops(boundary); 213 [e, d_n, d_t, coeff_t, coeff_n s, halfnorm_inv_n, halfnorm_inv_t, halfnorm_t,g] = obj.get_boundary_ops(boundary);
214 a_t = @(t) spdiag(coeff_t(t)); 214 a_t = @(t) spdiag(coeff_t(t));
215 a_n = @(t) spdiag(coeff_n(t)); 215 a_n = @(t) spdiag(coeff_n(t));
216 boundary
217 216
218 F = @(t)(s * a_n(t)*d_n' + s * a_t(t) *d_t')'; 217 F = @(t)(s * a_n(t)*d_n' + s * a_t(t) *d_t')';
219 tau1 = 1; 218 tau1 = 1;
220 a = @(t)spdiag(g(t)); 219 a = @(t)spdiag(g(t));
221 tau2 = @(t) (1*s*a(t))/2; 220 tau2 = @(t) (1*s*a(t))/2;