Mercurial > repos > public > sbplib_julia

comparison test/SbpOperators/volumeops/derivatives/second_derivative_test.jl @ 1158:c94a12327737 refactor/sbpoperators/inflation

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Disregard review comments about multi-d tests. After discussion on Discord we concluded that removing multi-d tests for volume operators would make them too implementation dependent.
author Jonatan Werpers <jonatan@werpers.com>
date Tue, 29 Nov 2022 21:46:54 +0100
parents f1bb1b6d85dd
children 7d52c4835d15
comparison
equal deleted inserted replaced
1157:03c217c50d7c 1158:c94a12327737
4 using Sbplib.Grids 4 using Sbplib.Grids
5 using Sbplib.LazyTensors 5 using Sbplib.LazyTensors
6 6
7 import Sbplib.SbpOperators.VolumeOperator 7 import Sbplib.SbpOperators.VolumeOperator
8 8
9 # Review: 9 # TODO: Refactor these test to look more like the tests in first_derivative_test.jl.
10 # If we test LazyTensor.inflation for e.g 1D-3D general tensors then we should only need to test 10
11 # the 1D second derivative.
12 @testset "SecondDerivative" begin 11 @testset "SecondDerivative" begin
13 operator_path = sbp_operators_path()*"standard_diagonal.toml" 12 operator_path = sbp_operators_path()*"standard_diagonal.toml"
14 stencil_set = read_stencil_set(operator_path; order=4) 13 stencil_set = read_stencil_set(operator_path; order=4)
15 inner_stencil = parse_stencil(stencil_set["D2"]["inner_stencil"]) 14 inner_stencil = parse_stencil(stencil_set["D2"]["inner_stencil"])
16 closure_stencils = parse_stencil.(stencil_set["D2"]["closure_stencils"]) 15 closure_stencils = parse_stencil.(stencil_set["D2"]["closure_stencils"])