Using Float128 in RandomMatrices

New to Julia
1

Hi guys:

I am using the package RandomMatrices to generate Haar random unitary matrices. But recently I want to generate unitary matrices with higher precision. So I found the package Quadmath which provides Float128 object. But how to combines these two packages together? How to generate a Haar random unitary matrix with Float128? By the way, is there a way to use things like randn(Float128,(N,N))? It seems that this doesn’t work for Float128 in Quadmath

2

I’m not sure if Quadmath has the necessary support, but you can do this with DoubleFloats (which give you roughly 100 bits of precision and are usually about 2x faster).

julia> using DoubleFloats
julia> randn(Double64, 10, 10)
10×10 Matrix{Double64}:
...
1 Like
3

Thanks this sounds good.