Mercurial > repos > public > sbplib

comparison +scheme/Laplace1d.m @ 1083:a22b23098021 feature/laplace1dUpdate

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Divide interface penalty strength by 4, to make it sharp. Stability verified by eigenvalue computations.
author Martin Almquist <malmquist@stanford.edu>
date Tue, 05 Mar 2019 12:41:50 -0800
parents a20fc67e9ac0
children 33c378e508d2
comparison
equal deleted inserted replaced
1082:a20fc67e9ac0 1083:a22b23098021
102 gamm_u = obj.gamm; 102 gamm_u = obj.gamm;
103 gamm_v = neighbour_scheme.gamm; 103 gamm_v = neighbour_scheme.gamm;
104 104
105 tuning = 1.1; 105 tuning = 1.1;
106 106
107 tau1 = -(a_u/gamm_u + a_v/gamm_v) * tuning; 107 tau1 = -1/4*(a_u/gamm_u + a_v/gamm_v) * tuning;
108 tau2 = 1/2*a_u; 108 tau2 = 1/2*a_u;
109 sig1 = -1/2; 109 sig1 = -1/2;
110 sig2 = 0; 110 sig2 = 0;
111 111
112 tau = tau1*e_u + tau2*d_u; 112 tau = tau1*e_u + tau2*d_u;