# lusolve

solveur de système linéaire creux

## Séquence d'appel

    x=lusolve(hand,b)
                x=lusolve(A,b)

## Paramètres

b

  matrice réelle pleine

A

  matrice réelle carrée creuse inversible

hand

  pointeur vers des facteurs L,U creux déjà calculés (sortie de lufact)

x

  matrice réelle pleine

## Description

`x=lusolve(hand,b)` résout le système linéaire `A*x = b`, où `hand` est obtenu par un appel préliminaire à `lufact` :`[hand,rk]=lufact(A)`.

`x=lusolve(A,b)` résout le système linéaire creux `A*x = b`

## Exemples

    non_zeros=[1,2,3,4];rows_cols=[1,1;2,2;3,3;4,4];
    sp=sparse(rows_cols,non_zeros);
    [h,rk]=lufact(sp);x=lusolve(h,[1;1;1;1]);ludel(h)
    rk,sp*x

    non_zeros=[1,2,3,4];rows_cols=[1,1;2,2;3,3;4,4];
    sp=sparse(rows_cols,non_zeros);
    x=lusolve(sp,-ones(4,1));
    sp*x

## Voir aussi

- sparse
- lufact
- slash
- backslash