Mercurial > repos > public > sbplib_julia

comparison test/SbpOperators/volumeops/derivatives/second_derivative_test.jl @ 1594:d68d02dd882f feature/boundary_conditions

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Merge with default
author Vidar Stiernström <vidar.stiernstrom@gmail.com>
date Sat, 25 May 2024 16:07:10 -0700
parents 43aaf710463e
children 471a948cd2b2
comparison
equal deleted inserted replaced
1591:615eeb6e662e 1594:d68d02dd882f
13 stencil_set = read_stencil_set(operator_path; order=4) 13 stencil_set = read_stencil_set(operator_path; order=4)
14 inner_stencil = parse_stencil(stencil_set["D2"]["inner_stencil"]) 14 inner_stencil = parse_stencil(stencil_set["D2"]["inner_stencil"])
15 closure_stencils = parse_stencil.(stencil_set["D2"]["closure_stencils"]) 15 closure_stencils = parse_stencil.(stencil_set["D2"]["closure_stencils"])
16 Lx = 3.5 16 Lx = 3.5
17 Ly = 3. 17 Ly = 3.
18 g_1D = equidistant_grid(121, 0.0, Lx) 18 g_1D = equidistant_grid(0.0, Lx, 121)
19 g_2D = equidistant_grid((121,123), (0.0, 0.0), (Lx, Ly)) 19 g_2D = equidistant_grid((0.0, 0.0), (Lx, Ly), 121, 123)
20 20
21 @testset "Constructors" begin 21 @testset "Constructors" begin
22 @testset "1D" begin 22 @testset "1D" begin
23 Dₓₓ = second_derivative(g_1D, stencil_set) 23 Dₓₓ = second_derivative(g_1D, stencil_set)
24 @test Dₓₓ == second_derivative(g_1D, inner_stencil, closure_stencils) 24 @test Dₓₓ == second_derivative(g_1D, inner_stencil, closure_stencils)