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
Partager