build.jl

# Use this script to generate grid files and save them into a given folder
using Lehmann
using Printf

# rtol = [-6, -8, -10, -12, -14]
# rtol = [-7, -9, -11, -13]
# Λ = [100, 1000, 10000, 100000, 1e6, 1e7, 1e8, 1e9, 1e10, 1e11, 1e12, 1e13, 1e14]
# Λ = [1e11, 1e12, 1e13, 1e14]
rtol = [-12]
Λ = [1000000,]
algorithm = :functional
# folder = "./basis/"
folder = "./"

for lambda in Λ
    for err in rtol
        if lambda <= 1e8 #the universal grid beyond 1e8 suffers from the accuracy loss, so that we turn it off
            DLRGrid(1.0, lambda, 10.0^err, true, :none, rebuild = true, algorithm = algorithm, folder = folder, verbose = true)
        end
        DLRGrid(1.0, lambda, 10.0^err, true, :ph, rebuild = true, algorithm = algorithm, folder = folder, verbose = true)
        DLRGrid(1.0, lambda, 10.0^err, true, :pha, rebuild = true, algorithm = algorithm, folder = folder, verbose = true)
    end
end