Discussion :

[Merel]

salut à tous,

Dans le cadre de mes travaux concernant le problèmes inverses, je dois effectuer un calcul simple qui me pose beaucoup de problèmes...

J'ai 20412 données à inverser qui bloquent complètement Matlab.
Je dois en effet créer une matrice de 20412x20412 éléments ce qui représente en gros 3Go de mémoire. Je peux réduire la place à 1.5Go en passant par des classes single. Mais le problème reste le même.

J'obtiens le message d'erruer "maximum variable size..." et quand ca n'est pas lui (sur une matrice 20412x750) c'est le "out of memory..." que j'obtiens.

J'ai essayé sur un matlab 7.1 R14 SP3 avec vista 32bits 4Go de ram
et sur un matlab 7.1 R14 SP3 avec un linux 2.6 4Go de ram également

Je dois avoir ces variables en mémoires :

1 matrice 20412x20412
1 matrice 20412x750
3 vecteurs 750
1 matrice 750x750

Est-ce que quelqu'un à éventuellement une idée ? Passer en Fortran, C ?

[Jerome Briot]

Comme souvent dans ces problèmes, la limite n'est pas donnée par la quantité de mémoire disponible mais par le nombre d'adresse que le système d'exploitation peut gérer.

Dans ton cas, si tu veux disposer de toutes les matrices entières en même temps, il te faut une machine 64bits.

Voir les liens donnés dans la : Out of memory. Type HELP MEMORY for your options.

[Merel]

Merci Dut pour cette réponse,

Ca confirme mes craintes... J'ai lu la faq concernant ce type de problème. C'est pour ca que je cherchais une autre méthode... Opter pour du C ou du fortran etc... plutôt que de dénicher une machine 64 bits... La ca sera mission impossible...

[Jerome Briot]

Je dois avoir ces variables en mémoires :

1 matrice 20412x20412
1 matrice 20412x750
3 vecteurs 750
1 matrice 750x750

N'y a-t-il pas moyen d'utiliser une méthode itérative ne nécessitant que des morceaux de ces matrices plutôt que de vouloir utiliser une méthode frontale ?

[Merel]

Hmm, ca me parait compliqué puisque je dois résoudre une équation qui utilise tous ces éléments.

En gros, et avec des lettres, la résolution d'un système inverse par moindre carré conduit à une équation de ce type

résultat = A+inv(B*C*B+D)*B*C*(E*F-B*A)

Et ca, je ne vois pas comment le faire par morceaux. Je peux juste jouer sur des paramètres de grilles et sur les données pour réduire la taille des matrices. J'ai testé en divisant par deux le nombre de données et ca plante toujours...

Mais bon dans l'idéal, j'aurais besoin de ces résultats pour valider ou invalider une méthode.

[Jerome Briot]

résultat = A+inv(B*C*B+D)*B*C*(E*F-B*A)

Et ca, je ne vois pas comment le faire par morceaux.

Tu peux toujours poser ta question ici : Algorithmes > Mathématiques (en mentionnant cette discussion aussi dans ton message)

[Merel]

Ouep, c'est une bonne idée mais je ne pense pas que l'équation puisse être modifiée. Je posterai quand même après quelques tests. J'ai pu trouver deux machines 64 bits. La première avec un linux et l'autre avec windows.

Sur la machine linux la création de toutes les variables passe nikel. Par contre l'inversion de la grosse matrice 20412x20412 plante matlab. Aucun message d'erreur, juste une fermeture pure et simple du truc.

Je testerais ce soir le 64 bits windows. Il a plus de ram (16Go). J'espère que ca passera. Merci ;-)

[FR119492]

Salut!

l'inversion de la grosse matrice 20412x20412

Pourquoi calculer l'inverse d'une telle matrice?
Jean-Marc Blanc