brightway2 中的技术圈矩阵有多大

Question

在 Brightway2 中有多大的技术圈矩阵仍然是可逆的？据我了解 here Brightway 使用 Pardiso 库的包装器来加速矩阵求逆。与标准 scipy 库相比，Pardiso 确实快得多。但是当我测试它时，我运行遇到了大于 10000 行的矩阵的内存问题。要反转矩阵，我需要定义一个 numpy 数组，它在我的笔记本电脑上不能超过 1000000 行。

from pypardiso import spsolve
import numpy as np
ar=csc_matrix(np.eye(10000))
%time spsolve(ar,np.eye(10000))

这里是限制，还是我可以绕过 numpy 数组的使用？

(PS：我用标准 scipy 求解器进行了测试，我可以反转 10^5 行的稀疏矩阵，但速度很慢。

Answer 1

我知道 Adrian Haas 使用的稀疏矩阵超过 100 000 rows/columns。唯一的限制应该是你机器的内存，而不是软件本身。

您可以在合理的时间内逐列创建逆向，并获得更好的数值稳定性和 pardiso 速度的好处：

In [1]: from brightway2 import *

In [2]: import pyprind

In [3]: from time import time

In [4]: db = Database("ecoinvent 3.5 cutoff")

In [5]: def invert(database):
   ...:    lca = LCA({database.random(): 1})
   ...:    lca.lci()
   ...:    for act in pyprind.prog_bar(database):
   ...:        lca.redo_lci({act: 1})
   ...:

In [6]: start = time(); invert(db); print(time() - start)
0%                          100%
[##############################] | ETA: 00:00:00
Total time elapsed: 00:03:21
202.16850423812866

然而，根据我的经验，很少有人真正需要技术圈的逆向。

brightway2 中的技术圈矩阵有多大

how big technosphere matrices can be in brightway2

python-3.x

brightway