Mercurial > repos > public > sbplib_julia

comparison test/ext/sparse_arrays_test.jl @ 1725:8317252e4535

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Merge feature/lazy_tensors/sparse_conversions
author Jonatan Werpers <jonatan@werpers.com>
date Thu, 05 Sep 2024 23:16:39 +0200
parents 7274ec5050a0
children 471a948cd2b2
comparison
equal deleted inserted replaced
1719:d81d0660b80d 1725:8317252e4535
1 using Test
2
3 using Sbplib
4 using Sbplib.Grids
5 using Sbplib.SbpOperators
6
7 using SparseArrays
8 using Tokens
9
10
11 @testset "SparseArray" begin
12 g = equidistant_grid((0,0),(1,2), 20,30)
13 stencil_set = read_stencil_set(sbp_operators_path()*"standard_diagonal.toml"; order=4)
14
15
16 @testset let Δ = laplace(g, stencil_set), M = sparse(Δ)
17 @test ndims(M) == 2
18 @test size(M) == (20*30,20*30)
19
20 v = rand(size(g)...)
21
22 Mv = M*reshape(v,:)
23 @test Mv ≈ reshape(Δ*v,:)
24 end
25
26 @testset let dₙ = normal_derivative(g, stencil_set,CartesianBoundary{1,LowerBoundary}()), M = sparse(dₙ)
27 @test ndims(M) == 2
28 @test size(M) == (30,20*30)
29
30 v = rand(size(g)...)
31 Mv = M*reshape(v,:)
32 @test Mv ≈ reshape(dₙ*v,:)
33 end
34 end