问题:A是正方形,满级,稀疏和带状.在Matlab中有太多的元素存储为单个矩阵(至少~4.6*10 18,理想情况下~10 40,两者都超过最大数组大小.编辑:A存储为稀疏,问题不是内存有限,但元素数量有限).因此,我必须将其存储为较小数组(行/对角线/列/块)的集合.
寻找:一种解决Ax = b的方法,A给定为较小数组的集合.理想情况下,在Matlab中,但不是必须的.
或者,如果不是在Matlab中:也许有一个程序可以存储和解决这么大的A？
到目前为止找到:方法如果A是三/五对角线,但我的A有N个对角线.还发现了一些关于将A分区为块的问题,但是找不到用这些块解决线性系统的方法.
ps系统是64位.
感谢大家!