朱莉娅的大型矩阵幂

Question

我正在处理非常大的图及其相应的加权邻接矩阵，并且需要将这些大矩阵采用类似的大幂（即，将矩阵提高到成千上万的幂）。

我遇到的问题是矩阵的元素很快变得太大，计算机无法处理，我想知道如何解决这个问题。

以前有人解决过这样的问题吗（将矩阵提高为大数），您如何解决它们？

我知道Python的numpy可以处理这些计算。 Julia中有一个类似的库也可以做到这一点吗？

Answer 1

您可以将类型转换为BigFloat：

julia> A = [1.5 2 -4; 3 -1 -6; -10 2.3 4]
3×3 Array{Float64,2}:
   1.5   2.0  -4.0
   3.0  -1.0  -6.0
 -10.0   2.3   4.0

julia> (BigFloat.(A))^32000
3×3 Array{BigFloat,2}:
  4.16164e+31019   8.71351e+31017  -3.22788e+31019
  4.60207e+31019   9.63565e+31017  -3.56949e+31019
 -5.83403e+31019  -1.22151e+31018   4.52503e+31019