Mercurial > repos > public > sbplib_julia

comparison test/SbpOperators/volumeops/derivatives/second_derivative_test.jl @ 1529:43aaf710463e refactor/equidistant_grid/signature

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Change to signature of equidistant_grid to same style as many array methods. See for example Array{T}(undef, dims...), zeros(T, dims...), fill(a, dims...) and more.
author Jonatan Werpers <jonatan@werpers.com>
date Thu, 11 Apr 2024 22:31:04 +0200
parents 356ec6a72974
children 471a948cd2b2
comparison
equal deleted inserted replaced
1528:d641798539c2 1529:43aaf710463e
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)