Bibliothèque Python permettant de Gauss-Seidel Solveur Itératif?

Est-il une algèbre linéaire bibliothèque qui implémente itérative de Gauss-Seidel pour résoudre des systèmes linéaires? Ou peut-être un préconditionnés gradient solveur?

Grâce

EDIT: En fin de compte j'ai utilisé une sorte de brut mais bonne façon de le résoudre. Comme j'ai eu à créer la matrice A (pour Ax=b) de toute façon, j'ai partitionné de la matrice comme 

A = M - N

avec 

 M = (D + L) and N = -U

où D est la diagonale, L est triangulaire inférieure de la section, et U triangulaire supérieure de la section. Puis

Pinv = scipy.linalg.inv(M)
x_k_1 = np.dot(Pinv,np.dot(N,x_k)) + np.dot(Pinv,b)

Également fait quelques tests de convergence. Elle fonctionne.

OriginalL'auteur Ivan | 2011-04-11

  1. 3

    Il est un exemple de mise en œuvre où ils vont à travers différentes implémentations ici:

    http://www.scipy.org/PerformancePython

    un autre exemple ici:

    http://www.dur.ac.uk/physics.astrolab/py_source/kiusalaas/v1_with_numpy/gaussSeidel.py

    Vous pourriez également vouloir vérifier les liaisons python pour petsc:

    http://code.google.com/p/petsc4py/

    Merci pour l'entrée. Dans la fin ne ma propre solveur.
    Malheureusement, les deux premiers liens sont rompus, et le dernier de la date.

    OriginalL'auteur JoshAdel

  2. 0
    from sage.all import *

a=matrix(QQ,[[12,3,-5],[1,5,3],[3,7,13]])

b=matrix(QQ,[[1],[28],[76]])


x=[]
r=len(a[0])
i=0
while(i<r):
    li=raw_input('give approximate solution :')
    h=eval(li)
    x.append(h)
    i=i+1


def gausseidel():
    tol=0.00001;l=0;itern=10
    fnd=0

    while(l<itern and fnd==0):
        j=0
        while(j<r):
            temp=b[j,0]
            k=0
            while(k<r):
                if (j!=k):
                    temp-=a[j,k]*x[k]
                k=k+1
            c=temp/a[j,j]

            if (abs(x[j]-c)<=tol):
                fnd=1
                break
            x[j]=c
            j=j+1
        l=l+1
        if l==itern:
            print "Iterations are over"

    for p in range(r):
        print 'x',p,'=',float(x[p])


gausseidel()
    J'ai formaté votre code pour vous (à l'aide de la {} bouton). Merci de modifier votre réponse à expliquer d'où sage.all vient et pourquoi vous l'avez utilisé.
    sage est le général scientifique package: sagemath.org. Merci pour la réponse.

    OriginalL'auteur bokks