TP1 - Partie 1 - Aides : 4. 2

TP1 - Partie 1 - Aides : Etape 4 : résolution : 4. 2 - Stockage morse

4. 2 - Stockage morse

Vous disposez maintenant d’un (petit) code résolvant l’équation de la chaleur stationnaire 1D avec des conditions aux limites de Dirichlet.

Par contre le stockage de la matrice du système n’est pas optimal. Un stockage morse serait mieux ...

Pour cela, écrire les sous-programmes de

construction la matrice $K$ en complétant les zones grises de Kmorse1DP1 (mettre avec Kplein1DP1 dans un même fichier, par exemple K1DP1.sce)

$\includegraphics[width=9.5cm]{./prog-images/Kmorse1DP1.png}$
modification la matrice $K$ pour les CL en complétant les zones grises de AmorseCLdi (mettre cette function dans AFCL.sce)
$\includegraphics[width=9.5cm]{./prog-images/AmorseCLdi.png}$
résolution de $KX=F$ en complétant les zones grises du fichier chaleur1DP1_morse.sci
$\includegraphics[width=7cm]{./prog-images/chaleur1DP1dimo.png}$

A tester de la même manière que précédemment :

exec(’K1DP1.sce’,-1) ; exec(’AFCL.sce’,-1) ;

exec(’chaleur1D1P_morse.sci’,-1) ;

...

T=chaleur1D1P_di_morse(nx,x,kc,cldi) ; T

fin de l’aide sur l’étape 4