Discussion :

[Bob_Tom]

Bonjour,

je cherche a résouldre A*x=B, pour le moment j'utilise DGESV mais qui le fait pour des matrice dense (quelconque). le truc c'est que ma matrice A a plein de 0

# # 0 ..................................
# # # 0 ...............................
0 # # # 0 ............................
0 0 # # # 0 .........................
0 0 0 # # # 0 ......................

............0 # # # 0 ...............

...........................0 # # # 0
.............................0 # # #
...............................0 # #

En d'autre mot seul aij aij-1 et aij+1 ne sont pas nul

je sais qu'il existe des routines pour résoudre les systèmes avec des matrices symétriques, antisymétriques, définie positive

je souhaite diminuer la mémoire utiliser de mon programme, augmenter le vitesse d'exécution et diminué la propagation d'erreur surtout

merci de votre aide

[FR119492]

Salut!
Ta matrice est tridiagonale. Il faudrait donc utiliser DGTSV.
Jean-Marc Blanc