Merge branch 'main' of https://github.com/quantumstatistics/Lehmann.jl

Author: LittleBug

.DS_Store

@@ -0,0 +1,11 @@
+using JLD2
+# using CodecZlib
+
+mutable struct A
+    m::Matrix{Float64}
+    n::Matrix{Float64}
+end
+
+N = 100
+a = A(rand(N, N), rand(N, N))
+save("test2.jld2", Dict("a" => a), compress = true)

experiment/test.jl

@@ -0,0 +1,7 @@
+[deps]
+CodecZlib = "944b1d66-785c-5afd-91f1-9de20f533193"
+FastGaussQuadrature = "442a2c76-b920-505d-bb47-c5924d526838"
+JLD2 = "033835bb-8acc-5ee8-8aae-3f567f8a3819"
+Lehmann = "95bf888a-8996-4655-9f35-1c0506bdfefe"
+Printf = "de0858da-6303-5e67-8744-51eddeeeb8d7"
+Test = "8dfed614-e22c-5e08-85e1-65c5234f0b40"

test/Project.toml

@@ -247,3 +247,18 @@ end
     @test dlr_load.Λ ≈ dlr.Λ
     rm(file)
 end
+
+@testset "JLD2 IO" begin
+    dlr = DLRGrid(isFermi = true, beta = 10.0)
+    save("test.jld2", Dict("dlr" => dlr), compress = true)
+    dlr_new = load("test.jld2")["dlr"]
+    println(dlr_new)
+    @test dlr.isFermi == dlr_new.isFermi
+    @test dlr.Euv ≈ dlr_new.Euv
+    @test dlr.β ≈ dlr_new.β
+    @test dlr.τ ≈ dlr_new.τ
+    @test dlr.n ≈ dlr_new.n
+    @test dlr.ωn ≈ dlr_new.ωn
+    @test dlr.symmetry == dlr_new.symmetry
+    rm("test.jld2")
+end

test/dlr.jl

@@ -2,6 +2,8 @@ using Base: BufferStream
 using Lehmann
 using Test
 using Printf
+using JLD2
+using CodecZlib #it is important to import CodecZlib explicitly. Otherwise, the code may try to dynamically load this package, and sometimes leads to error
 
 include("spectral.jl")
 include("dlr.jl")