Mercurial > repos > public > sbplib
comparison +scheme/LaplaceCurvilinear.m @ 743:f4595f14d696 feature/utux2D
Change schemes to work for special coefficients.
author | Martin Almquist <malmquist@stanford.edu> |
---|---|
date | Mon, 21 May 2018 23:11:34 -0700 |
parents | 07f8311374c6 |
children | 14b093a344eb |
comparison
equal
deleted
inserted
replaced
720:07f8311374c6 | 743:f4595f14d696 |
---|---|
55 if b ~=1 | 55 if b ~=1 |
56 error('Not implemented yet') | 56 error('Not implemented yet') |
57 end | 57 end |
58 | 58 |
59 % assert(isa(g, 'grid.Curvilinear')) | 59 % assert(isa(g, 'grid.Curvilinear')) |
60 if isa(a, 'function_handle') | |
61 a = grid.evalOn(g, a); | |
62 a = spdiag(a); | |
63 end | |
60 | 64 |
61 m = g.size(); | 65 m = g.size(); |
62 m_u = m(1); | 66 m_u = m(1); |
63 m_v = m(2); | 67 m_v = m(2); |
64 m_tot = g.N(); | 68 m_tot = g.N(); |