中文字幕日韩精品一区二区免费_精品一区二区三区国产精品无卡在_国精品无码专区一区二区三区_国产αv三级中文在线

高斯消元的Python源碼規(guī)模排旋轉(zhuǎn)

如下的資料是關于高斯消元的Python規(guī)模排旋轉(zhuǎn)的內(nèi)容,希望能對大家有些好處。

成都網(wǎng)絡公司-成都網(wǎng)站建設公司創(chuàng)新互聯(lián)建站十年經(jīng)驗成就非凡,專業(yè)從事成都網(wǎng)站建設、成都做網(wǎng)站,成都網(wǎng)頁設計,成都網(wǎng)頁制作,軟文營銷1元廣告等。十年來已成功提供全面的成都網(wǎng)站建設方案,打造行業(yè)特色的成都網(wǎng)站建設案例,建站熱線:13518219792,我們期待您的來電!

''' x = gaussPivot(a,b,tol=1.0e-9).
    Solves [a]{x} =  by Gauss elimination with
    scaled row pivoting
'''    
from numpy import zeros,argmax,dot
import swap 
import error

def gaussPivot(a,b,tol=1.0e-12):
    n = len(b)

  # Set up scale factors
    s = zeros(n)
    for i in range(n):
        s[i] = max(abs(a[i,:]))

    for k in range(0,n-1):

      # Row interchange, if needed
        p = argmax(abs(a[k:n,k])/s[k:n]) + k
        if abs(a[p,k]) < tol: error.err('Matrix is singular')
        if p != k:
            swap.swapRows(b,k,p)
            swap.swapRows(s,k,p)
            swap.swapRows(a,k,p)

      # Elimination
        for i in range(k+1,n):
            if a[i,k] != 0.0:
                lam = a[i,k]/a[k,k]
    if abs(a[n-1,n-1]) < tol: error.err('Matrix is singular')

  # Back substitution
    b[n-1] = b[n-1]/a[n-1,n-1]
    for k in range(n-2,-1,-1):
        b[k] = (b[k] - dot(a[k,k+1:n],b[k+1:n]))/a[k,k]
    return b

分享文章:高斯消元的Python源碼規(guī)模排旋轉(zhuǎn)
瀏覽路徑:http://m.rwnh.cn/article30/igposo.html

成都網(wǎng)站建設公司_創(chuàng)新互聯(lián),為您提供網(wǎng)站改版、建站公司品牌網(wǎng)站建設、網(wǎng)站維護小程序開發(fā)、網(wǎng)站排名

廣告

聲明:本網(wǎng)站發(fā)布的內(nèi)容(圖片、視頻和文字)以用戶投稿、用戶轉(zhuǎn)載內(nèi)容為主,如果涉及侵權請盡快告知,我們將會在第一時間刪除。文章觀點不代表本網(wǎng)站立場,如需處理請聯(lián)系客服。電話:028-86922220;郵箱:631063699@qq.com。內(nèi)容未經(jīng)允許不得轉(zhuǎn)載,或轉(zhuǎn)載時需注明來源: 創(chuàng)新互聯(lián)

h5響應式網(wǎng)站建設
如皋市| 安达市| 西盟| 茂名市| 福泉市| 青神县| 神农架林区| 夏津县| 黄骅市| 荃湾区| 涞源县| 勐海县| 三门县| 喜德县| 阳谷县| 百色市| 九龙城区| 邹城市| 郯城县| 沐川县| 宜兴市| 天水市| 连平县| 保定市| 深圳市| 保亭| 临西县| 丹寨县| 霍林郭勒市| 潮州市| 万荣县| 兴山县| 饶平县| 博白县| 新安县| 东兴市| 义乌市| 千阳县| 六安市| 南靖县| 西乡县|