Mercurial > repos > public > sbplib

comparison +multiblock/Laplace.m @ 776:23ad69a347dd feature/grids

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Remove multiplication with jacobian where it shouldn't be
author Jonatan Werpers <jonatan@werpers.com>
date Mon, 23 Jul 2018 09:32:56 -0700
parents e05465aa2e25
children
comparison
equal deleted inserted replaced
775:c596122374df 776:23ad69a347dd
19 obj.order = order; 19 obj.order = order;
20 obj.mbDiffOp = multiblock.DiffOp(@scheme.LaplaceCurvilinear, obj.grid, order, {a,b,opGen}); 20 obj.mbDiffOp = multiblock.DiffOp(@scheme.LaplaceCurvilinear, obj.grid, order, {a,b,opGen});
21 21
22 obj.D = obj.mbDiffOp.D; 22 obj.D = obj.mbDiffOp.D;
23 obj.J = obj.jacobian(); 23 obj.J = obj.jacobian();
24 obj.H = obj.mbDiffOp.H * obj.jacobian(); 24 obj.H = obj.mbDiffOp.H;
25 end 25 end
26 26
27 function s = size(obj) 27 function s = size(obj)
28 s = size(obj.mbDiffOp); 28 s = size(obj.mbDiffOp);
29 end 29 end