如何在Julia中调用LAPACK代码（cpbtrf）

Question

我正在尝试将我现有的Python代码翻译成Julia，我需要计算一个带状复杂矩阵的Cholesky分解。正确的LAPACK例程是cpbtrf（当前由SciPy调用的例程），我正努力让它在Julia中工作。

我不确定要提供什么额外的细节，我对Julia很新，我确信我做的事情很愚蠢。 LAPACK调用在info变量中返回1，表示某些东西不是正定的，但我知道它是（SciPy很高兴地分解相同的矩阵）。

=IF(A1=B1,CONCATENATE(A1,B1),"")

编辑：只是为了澄清，这是一个骨架的例子。我正在编写的代码处理10 ^ 5或更大的矩阵，并且可能需要五角形，六角形，七对角矩阵等等。我需要一个特定于带状的算法。

Answer 1

除LAPACK子程序外，它都是正确的。您使用的是128位复数，因此应使用:zpbtrf_代替:cpbtrf_。