五对角矩阵工具箱。
pentap的Python项目详细描述
欢迎使用Pentapy
目的
pentapy是python中处理五对角矩阵的工具箱。
安装
可以通过pip安装包。 在windows上,您可以安装WinPython来获取 python和pip正在运行。
pip install pentapy
参考文献
求解器基于ptrans-i和ptrans-ii算法 由Askar et al. 2015提出。
示例
求解五对角线性方程组
这是一个用五对角矩阵求解les的例子。
importnumpyasnpfrompentapyimportsolvesize=1000# create a flattened pentadiagonal matrixM_flat=(np.random.random((5,size))-0.5)*1e-5V=np.random.random(size)*1e5# solve the LESX=solve(M_flat,V,is_flat=True)# create the corresponding matrix for checkingM=(np.diag(M_flat[0,:-2],2)+np.diag(M_flat[1,:-1],1)+np.diag(M_flat[2,:],0)+np.diag(M_flat[3,1:],-1)+np.diag(M_flat[4,2:],-2))# calculate the errorprint(np.max(np.abs(np.dot(M,X)-V)))
这应该提供如下信息:
4.257890395820141e-08
要求:
联系人
你可以通过info@geostat-framework.org联系我们。
许可证
GPL©2019年