Mercurial > repos > public > sbplib

comparison operator_def/assemble_opTest.m @ 830:e15a667ffde2 feature/operator_files

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Add quadrature stencils and operator assembly. Add tests
author Vidar Stiernström <vidar.stiernstrom@it.uu.se>
date Mon, 10 Sep 2018 19:27:27 +0200
parents f82da6644f42
children
comparison
equal deleted inserted replaced
829:e0913772dc1c 830:e15a667ffde2
40 verifyEqual(testCase,e_l,op.e_l,'AbsTol',1e-10) 40 verifyEqual(testCase,e_l,op.e_l,'AbsTol',1e-10)
41 verifyEqual(testCase,e_r,op.e_r,'AbsTol',1e-10) 41 verifyEqual(testCase,e_r,op.e_r,'AbsTol',1e-10)
42 verifyEqual(testCase,d1_l,op.d1_l,'AbsTol',1e-10) 42 verifyEqual(testCase,d1_l,op.d1_l,'AbsTol',1e-10)
43 verifyEqual(testCase,d1_r,op.d1_r,'AbsTol',1e-10) 43 verifyEqual(testCase,d1_r,op.d1_r,'AbsTol',1e-10)
44 end 44 end
45
46 function TestAssembleH(testCase)
47 m = 10;
48 op = sbp.D2Standard(m,{0 3},2);
49 h = op.h;
50 [H, Hinv] = assemble_H('H_2',h,m);
51 verifyEqual(testCase,H,op.H,'AbsTol',1e-10)
52 verifyEqual(testCase,Hinv,op.HI,'AbsTol',1e-10)
53 end