Clojure中的分数矩阵幂

Question

我正在尝试使用clojure中的过渡矩阵。在将年度债券评级转换矩阵转换为季度时，我需要0.25次方矩阵。

在python中，我们将fractional_matrix_power from scipy作为：

>>> from scipy.linalg import fractional_matrix_power
>>> a = np.array([[1.0, 3.0], [1.0, 4.0]])
>>> b = fractional_matrix_power(a, 0.5)
>>> b
array([[ 0.75592895,  1.13389342],
       [ 0.37796447,  1.88982237]])

在搜索Incanter和Parallel Colt时，我还没有找到任何内容。浏览javadocs和谷歌搜索没有帮助，但可能“矩阵”和“权力”过于通用，无法深入到我想要的内容。

我真的需要transcode a python或R函数，还是有一些很酷的colt doc网站我不知道？

Answer 1

您可能已经看过the core.matrix library for Clojure。它专注于基本线性代数，并且没有（现在！）正如您所描述的矩阵指数函数。

但是，如果您将上述算法从Python转码为Clojure，我相信它会成为core.matrix的一个受欢迎的补充。

你还应该研究the Neanderthal project，它专注于使用GPU进行矩阵运算。