Bonjour,

Dans le but de résoudre un "grand" nombre d'edp qui ne diffèrent entre elles que par leurs seconds membres, j'envisage de résoudre les N systèmes linéaires engendrés par une méthode directe de type factorisation LU ou autre.
Ceci offrirait la possibilité de n'effectuer qu'une seule factorisation de matrice pour la résolution des N systèmes linéaires.

La possibilité de conserver la factorisation après résolution du système linéaire pourrait aussi être intéressante pour la résolution d'un problème instationnaire à coefficients constants si l'on discrétise le temps avec un pas constant.


Existe-t-il dans Castem des procédures fonctionnant de cette manière ou bien permettant de le faire(factorisation et résolution à partir de la factorisation)?


Merci d'avance pour vos réponses.