Prepare version 0.3.0

SymSparseMatrixCSR only stores upper triangle, instead of lower,
in order to interact with GridapPardiso package

Author: victorsndvg

Project.toml

@@ -1,7 +1,7 @@
 name = "SparseMatricesCSR"
 uuid = "a0a7dd2c-ebf4-11e9-1f05-cf50bc540ca1"
 authors = ["Víctor Sande <[email protected]>"]
-version = "0.2.1"
+version = "0.3.0"
 
 [deps]
 LinearAlgebra = "37e2e46d-f89d-539d-b4ee-838fcccc9c8e"

src/SymSparseMatrixCSR.jl

@@ -14,7 +14,7 @@ end
 
 show(io::IO, A::SymSparseMatrixCSR) = show(io, A.lowertrian)
 size(A::SymSparseMatrixCSR) = size(A.lowertrian)
-getindex(A::SymSparseMatrixCSR, x::Integer, y::Integer) = getindex(A.lowertrian,min(x,y),max(x,y))
+getindex(A::SymSparseMatrixCSR, x::Integer, y::Integer) = getindex(A.lowertrian,max(x,y),min(x,y))
 
 """
     nnz(S::SymSparseMatrixCSR)
@@ -88,9 +88,10 @@ end
     function push_coo!(::Type{SymSparseMatrixCSR},I,J,V,ik,jk,vk) 
 
 Inserts entries in COO vectors for further building a SymSparseMatrixCSR.
+It stores only the upper triangle, ignoring entries with (ik>jk) coordinates.
 """
 function push_coo!(::Type{SymSparseMatrixCSR},I::Vector,J::Vector,V::Vector,ik::Integer,jk::Integer,vk::Number) 
-    (ik<jk) && return
+    (ik>jk) && return
     (push!(I, ik), push!(J, jk), push!(V, vk))
 end
 

test/SymSparseMatrixCSR.jl

@@ -11,7 +11,7 @@
             push_coo!(SymSparseMatrixCSR,I,J,V,ik,jk,vk)
         end
         finalize_coo!(SymSparseMatrixCSR,I,J,V,maxrows, maxcols)
-        SYMCSC = Symmetric(sparse(I, J, V, maxrows, maxcols),:L)
+        SYMCSC = Symmetric(sparse(I, J, V, maxrows, maxcols),:U)
         SYMCSR = symsparsecsr(I, J, V, maxrows, maxcols)
 
         @test size(SYMCSC)==size(SYMCSR)