TP3 - Aides : 3. 2 - Les termes constants / t

TP3 - Aides : Etape 3 : schéma en temps : 3. 2 - Les termes constants / t

3. 2 - Les termes constants / t

construire les termes indépendants du temps :

Les matrices $M$ , $A=\rho c_ vM +\Delta t \theta K$ et $B=A-\Delta t K$ et le vecteur $Fc$ sont indépendants du temps, on construit donc les matrices et $Fc$ avant la boucle en temps.

A=rc*M+dt*theta*K ; B=A-dt*K ;
si les matrices $M$ et $K$ sont construites Morse, $A$ et $B$ sont stockées automatiquement Morse,
si $M$ et $K$ sont construites pleines, B=sparse(A-dt*K) définie $B$ en stockage Morse (gain de place).
Les noeuds concernés par les CL étant les mêmes $\forall t$ , on modifie $A$ avant la boucle en temps pour prendre en compte les CL.
(attention : si $A$ est en stockage plein, utiliser ApleinCLdi du TP1, si $A$ est en stockage morse, utiliser AmorseCLdi du TP1).
La décomposition de la matrice $A$ en $A=LU$ est faite une seule fois, avant la boucle en temps (gain de temps) : A=lufact(A) ;
Remarque : si $A$ est stockée pleine, passer en stockage Morse par A=lufact(sparse(A)) ; (gain de place)
On sauvegarde la partie constante du second membre avec la CL de Neumann (constante en temps) dans Fc=M*fs.