LU decomposition of symetric pentadiagonal matrix in Python,,''' d,e,f =
文章由Byrx.net分享于2019-03-23 09:03:32
LU decomposition of symetric pentadiagonal matrix in Python,,''' d,e,f =
''' d,e,f = LUdecomp5(d,e,f). LU decomposition of symetric pentadiagonal matrix [f\e\d\e\f]. On output {d},{e} and {f} are the diagonals of the decomposed matrix. x = LUsolve5(d,e,f,b). Solves [f\e\d\e\f]{x} = {b}, where {d}, {e} and {f} are the vectors returned from LUdecomp5. '''def LUdecomp5(d,e,f): n = len(d) for k in range(n-2): lam = e[k]/d[k] d[k+1] = d[k+1] - lam*e[k] e[k+1] = e[k+1] - lam*f[k] e[k] = lam lam = f[k]/d[k] d[k+2] = d[k+2] - lam*f[k] f[k] = lam lam = e[n-2]/d[n-2] d[n-1] = d[n-1] - lam*e[n-2] e[n-2] = lam return d,e,fdef LUsolve5(d,e,f,b): n = len(d) b[1] = b[1] - e[0]*b[0] for k in range(2,n): b[k] = b[k] - e[k-1]*b[k-1] - f[k-2]*b[k-2] b[n-1] = b[n-1]/d[n-1] b[n-2] = b[n-2]/d[n-2] - e[n-2]*b[n-1] for k in range(n-3,-1,-1): b[k] = b[k]/d[k] - e[k]*b[k+1] - f[k]*b[k+2] return b
相关内容
- python中使用尾递归代码范例,python尾递归范例,# This p
- python提取url中的域名和端口号,pythonurl,import urlli
- python执行get提交的操作,pythonget提交,import sys,
- django获得用户ip地址,django获得ip,def get_clie
- python 装饰器记录日志,python日志,from time im
- Firefly 开源的游戏服务器端框架,firefly开源,Firefly是免
- 在VIM中使用GOOGLE进行搜索或者翻译,vimgoogle,一个在VI
- easy_install 错误Setup script exited with error: Unable to find vcv
- python简单的函数定义和使用范例,python函数范例,def c
- 找出数组中出现频率超过一半的数,数组超过一半,有一
评论关闭