blocksparsezeros with Float64 default

Author: mtfishman

Project.toml

@@ -1,7 +1,7 @@
 name = "BlockSparseArrays"
 uuid = "2c9a651f-6452-4ace-a6ac-809f4280fbb4"
 authors = ["ITensor developers <[email protected]> and contributors"]
-version = "0.7.11"
+version = "0.7.12"
 
 [deps]
 Adapt = "79e6a3ab-5dfb-504d-930d-738a2a938a0e"

src/blocksparsearray/blocksparsearray.jl

@@ -250,6 +250,9 @@ end
 function blocksparsezeros(elt::Type, axes...)
   return BlockSparseArray{elt}(undef, axes...)
 end
+function blocksparsezeros(axes...)
+  return blocksparsezeros(Float64, axes...)
+end
 function blocksparsezeros(::BlockType{A}, axes...) where {A<:AbstractArray}
   # TODO: Use:
   # ```julia

test/test_basics.jl

@@ -159,6 +159,15 @@ arrayts = (Array, JLArray)
     end
 
     # Convenient constructors.
+    a = blocksparsezeros([2, 3], [2, 3])
+    @test iszero(a)
+    @test iszero(blockstoredlength(a))
+    @test a isa BlockSparseMatrix{Float64,Matrix{Float64}}
+    @test blocktype(a) == Matrix{Float64}
+    @test blockstype(a) <: SparseMatrixDOK{Matrix{Float64}}
+    @test blocksize(a) == (2, 2)
+    @test blocksizes(a) == [(2, 2) (2, 3); (3, 2) (3, 3)]
+
     a = blocksparsezeros(elt, [2, 3], [2, 3])
     @test iszero(a)
     @test iszero(blockstoredlength(a))