https://github.com/JuliaStats/PDMats.jl is a highly optimized Julia linear algebra library for computations regarding positive-definite matrices. Since covariance matrices of Gaussian quantum states are necessarily positive-definite, then having support from this package might be worthwhile. I would be interested to see how large the performance benefits are, specifically for the use case of Gabs.