MSE12 版 (精华区)

发信人: chumsdock (微笑服务), 信区: MSE12       
标  题: 代码……
发信站: BBS 听涛站 (Sat Oct 29 22:36:22 2005), 转信

from scipy import *
from RandomArray import *
from Numeric import *
import pickle
from jacobi import *

if __name__ == '__main__':
    n=150
    seed()

    V=random((n,n,n))
    r=randint(0,n,(12,))

    B=jacobi(V,r)

    filev=open('temp11.dat','wb')
    pickle.dump(B,filev)
    filev.close()

三维雅克比的那段程序
Fortran编的jacobi子程比较难看，不贴了
--
    OK 5.4


※ 来源:·BBS 听涛站 tingtao.net·[FROM: 59.66.203.181]