Mercurial > repos > public > sbplib

comparison +scheme/Schrodinger1dCurve.m @ 692:dd3e9a0f5032 feature/quantumTriangles

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
added jacobian to stichshemes
author Ylva Rydin <ylva.rydin@telia.com>
date Fri, 15 Sep 2017 09:21:29 +0200
parents 32a24485f3e8
children 40299ae5fabe
comparison
equal deleted inserted replaced
691:527a7b954f26 692:dd3e9a0f5032
138 138
139 function [closure, penalty] = interface(obj,boundary,neighbour_scheme,neighbour_boundary) 139 function [closure, penalty] = interface(obj,boundary,neighbour_scheme,neighbour_boundary)
140 % u denotes the solution in the own domain 140 % u denotes the solution in the own domain
141 % v denotes the solution in the neighbour domain 141 % v denotes the solution in the neighbour domain
142 [e_u,d_u,s_u,p_u] = obj.get_boundary_ops(boundary); 142 [e_u,d_u,s_u,p_u] = obj.get_boundary_ops(boundary);
143 [e_v,d_v,s_v,p_v] = neighbour_scheme.get_boundary_ops(neighbour_boundary); 143 [e_v,d_v] = neighbour_scheme.get_boundary_ops(neighbour_boundary);
144 144
145 a1 = s_u* 1/2 * 1i ; 145 a1 = s_u* 1/2 * 1i ;
146 b1 = -s_u* 1/2 * 1i; 146 b1 = -s_u* 1/2 * 1i;
147 gamma = @(a) s_u*a(p_u,p_u)/2*e_u; 147 gamma = @(a) s_u*a(p_u,p_u)/2*e_u;
148 148