Publicité
+ Répondre à la discussion
Affichage des résultats 1 à 2 sur 2
  1. #1
    Nouveau Membre du Club
    Homme Profil pro Thomas Chauve
    Étudiant
    Inscrit en
    décembre 2012
    Messages
    37
    Détails du profil
    Informations personnelles :
    Nom : Homme Thomas Chauve
    Localisation : France, Isère (Rhône Alpes)

    Informations professionnelles :
    Activité : Étudiant

    Informations forums :
    Inscription : décembre 2012
    Messages : 37
    Points : 30
    Points
    30

    Par défaut Cherche Subroutine lapack A*x=B

    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

  2. #2
    Rédacteur/Modérateur

    Homme Profil pro Jean-Marc Blanc
    Comme retraité, des masses
    Inscrit en
    avril 2007
    Messages
    2 961
    Détails du profil
    Informations personnelles :
    Nom : Homme Jean-Marc Blanc
    Âge : 74
    Localisation : Suisse

    Informations professionnelles :
    Activité : Comme retraité, des masses
    Secteur : Industrie

    Informations forums :
    Inscription : avril 2007
    Messages : 2 961
    Points : 5 047
    Points
    5 047

    Par défaut

    Salut!
    Ta matrice est tridiagonale. Il faudrait donc utiliser DGTSV.
    Jean-Marc Blanc
    Calcul numérique de processus industriels
    Formation, conseil, développement

    Point n'est besoin d'espérer pour entreprendre, ni de réussir pour persévérer. (Guillaume le Taiseux)

+ Répondre à la discussion
Cette discussion est résolue.

Liens sociaux

Règles de messages

  • Vous ne pouvez pas créer de nouvelles discussions
  • Vous ne pouvez pas envoyer des réponses
  • Vous ne pouvez pas envoyer des pièces jointes
  • Vous ne pouvez pas modifier vos messages
  •