From e6f96a475be5be5838f0a8e335c7f0d96e1ff54f Mon Sep 17 00:00:00 2001 From: amartin Date: Fri, 26 Nov 2021 21:24:59 +1100 Subject: [PATCH] Using CSR format in PoissonTests.jl --- Project.toml | 2 +- test/PoissonTests.jl | 4 +++- 2 files changed, 4 insertions(+), 2 deletions(-) diff --git a/Project.toml b/Project.toml index 836d6e0..c8d0d9b 100644 --- a/Project.toml +++ b/Project.toml @@ -21,7 +21,7 @@ GridapDistributed = "0.2" MPI = "0.14, 0.15, 0.16, 0.17, 0.18, 0.19" PETSc_jll = "3.13" PartitionedArrays = "0.2.4" -SparseMatricesCSR = "0.6.1" +SparseMatricesCSR = "0.6.6" julia = "1.3" [extras] diff --git a/test/PoissonTests.jl b/test/PoissonTests.jl index aeff2c4..9487c62 100644 --- a/test/PoissonTests.jl +++ b/test/PoissonTests.jl @@ -1,3 +1,4 @@ +using SparseMatricesCSR using Gridap using Gridap.Algebra using Gridap.FESpaces @@ -35,7 +36,8 @@ function main(parts) a(u,v) = ∫( ∇(v)⋅∇(u) )dΩ l(v) = ∫( v*f )dΩ + ∫( v*g )dΓn - op = AffineFEOperator(a,l,U,V) + assem=SparseMatrixAssembler(SparseMatrixCSR{0,PetscScalar,PetscInt},Vector{Float64},U,V) + op = AffineFEOperator(a,l,U,V,assem) ls = PETScLinearSolver() fels = LinearFESolver(ls)