Mercurial > repos > public > sbplib_julia
comparison test/ext/sparse_arrays_test.jl @ 1699:3e9c3986930d feature/lazy_tensors/sparse_conversions
Add extensions for SparseArrays and SparseArrayKit that allow conversion of a LazyTensor
| author | Jonatan Werpers <jonatan@werpers.com> |
|---|---|
| date | Mon, 02 Sep 2024 15:35:54 +0200 |
| parents | |
| children | 1dd64b46ca2a |
comparison
equal
deleted
inserted
replaced
| 1676:5a1f51b4e3d9 | 1699:3e9c3986930d |
|---|---|
| 1 using Test | |
| 2 | |
| 3 using Sbplib | |
| 4 using Sbplib.Grids | |
| 5 using Sbplib.SbpOperators | |
| 6 using Sbplib.RegionIndices | |
| 7 | |
| 8 using SparseArrays | |
| 9 using Tokens | |
| 10 | |
| 11 | |
| 12 @testset "SparseArray" begin | |
| 13 g = equidistant_grid((0,0),(1,2), 20,30) | |
| 14 stencil_set = read_stencil_set(sbp_operators_path()*"standard_diagonal.toml"; order=4) | |
| 15 | |
| 16 | |
| 17 @testset let Δ = laplace(g, stencil_set), M = sparse(Δ) | |
| 18 @test ndims(M) == 2 | |
| 19 @test size(M) == (20*30,20*30) | |
| 20 | |
| 21 v = rand(size(g)...) | |
| 22 | |
| 23 Mv = M*reshape(v,:) | |
| 24 @test Mv ≈ reshape(Δ*v,:) | |
| 25 end | |
| 26 | |
| 27 @testset let dₙ = normal_derivative(g, stencil_set,CartesianBoundary{1,Lower}()), M = sparse(dₙ) | |
| 28 @test ndims(M) == 2 | |
| 29 @test size(M) == (30,20*30) | |
| 30 | |
| 31 v = rand(size(g)...) | |
| 32 Mv = M*reshape(v,:) | |
| 33 @test Mv ≈ reshape(dₙ*v,:) | |
| 34 end | |
| 35 end |