scipy中的float128 SVD

Question

我有一个float128类型的ndarray。我想计算此矩阵的SVD，但是当我使用sp.linalg.svd时，返回的数组的类型为float64。是否有一种简单的方法来使计算在float128中完成？

如果没有，是否可以直接调用LAPACK函数？

我还看到了mpmath，因为没有简单的方法可以使用它，所以我可能不得不使用它。有没有一种简单的方法可以在我的numpy数组上使用此函数？我无法将float128转换为mpmath矩阵。