Mercurial > repos > public > sbplib_julia

comparison test/SbpOperators/volumeops/derivatives/first_derivative_test.jl @ 2096:5af7534e5b3c feature/sbp_operators/laplace_curvilinear tip

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Merge default
author Jonatan Werpers <jonatan@werpers.com>
date Mon, 02 Mar 2026 15:58:27 +0100
parents e21c295fb2da
children
comparison
equal deleted inserted replaced
2082:87157cfca640 2096:5af7534e5b3c
14 g₁ = equidistant_grid(0., 1., 11) 14 g₁ = equidistant_grid(0., 1., 11)
15 g₂ = equidistant_grid((0.,1.), (1.,3.), 11, 14) 15 g₂ = equidistant_grid((0.,1.), (1.,3.), 11, 14)
16 16
17 @test first_derivative(g₁, stencil_set) isa LazyTensor{Float64,1,1} 17 @test first_derivative(g₁, stencil_set) isa LazyTensor{Float64,1,1}
18 @test first_derivative(g₂, stencil_set, 2) isa LazyTensor{Float64,2,2} 18 @test first_derivative(g₂, stencil_set, 2) isa LazyTensor{Float64,2,2}
19
20 @test first_derivative(g₁, stencil_set) == first_derivative(g₁, stencil_set, 1)
21 @test_throws DomainError(3, "Derivative direction must be 1.") first_derivative(g₁, stencil_set, 3)
22 @test_throws DomainError(3, "Derivative direction must be in 1:2.") first_derivative(g₂, stencil_set, 3)
19 23
20 interior_stencil = CenteredStencil(-1,0,1) 24 interior_stencil = CenteredStencil(-1,0,1)
21 closure_stencils = [Stencil(-1,1, center=1)] 25 closure_stencils = [Stencil(-1,1, center=1)]
22 26
23 @test first_derivative(g₁, interior_stencil, closure_stencils) isa LazyTensor{Float64,1,1} 27 @test first_derivative(g₁, interior_stencil, closure_stencils) isa LazyTensor{Float64,1,1}