Le calcul numérique de haute performance - Université de Laval

Une racine carrée est 10 à 20 fois plus
longue à calculer qu’une addition ou une
multiplication.
Solution: considérer la fonction � comme
étant une fonction de r 2 plutôt que r.
On crée de nouvelles tables d’interpolation
dans lesquelles la fonction est évaluée à
des intervalles constants en r 2 .
Nouvel incrément: �r 2
chi(0:npt)
0
1
2
3
4
5
6
7
8
9
npt
�(0)
�(�r 2 )
�(2�r 2 )
�(3�r 2 )
�(4�r 2 )
�(5�r 2 )
�(6�r 2 )
�(7�r 2 )
�(8�r 2 )
�(9�r 2 )
...
�(npt*�r 2 )
0
1
2
3
4
5
6
7
8
9
npt-1
dchi(0:npt-1)
�(�r 2 )��(0)
�(2�r 2 )��(�r 2 )
�(3�r 2 )��(2�r 2 )
�(4�r 2 )��(3�r 2 )
�(5�r 2 )��(4�r 2 )
�(6�r 2 )��(5�r 2 )
�(7�r 2 )-�(6�r 2 )
�(8�r 2 )��(7�r 2 )
�(9�r 2 )��(8�r 2 )
�(10�r 2 )��(9�r 2 )
...
�(npt*�r 2 )-�[(npt-1)*�r 2 ]