Discussion :

[feynman]

Bonjour,
J'ai besoin de votre avis sur la façon de programmer l'équation en attachement.

J'ai écrit le code mais je doute quand même de sa validité vu que je n'obtiens pas la bonne réponse.

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
 
         do f=0,TT
           call matbar(2*MM,mu)
           call matwise(2*MM,TR,mu,D)
           call dft2d(2*MM,D,Dhat)
           call dft2d(2*MM,Gh,Ghat)
           call matwise(2*MM,Ghat,Dhat,DD)
           call idft2d(2*MM,DD,DDinvhat)
           do i=0,(2*MM-1)
              do j=0,(2*MM-1)
                 mu(i,j)=III(i,j)-h**2*DDinvhat(i,j)
              enddo
           enddo
        enddo

Avec matbar qui est une subroutine qui transforme une matrice en son conjugué complexe, matwise pour le produit de Frobinus, dft2d pour la transformée de Fourier discrète, idft2d pour la transformée de Fourier discrète inverse.

Alors merci de me donner votre avis de la programmation de cette formule et si c'est possible vos idées ou vos suggestions.

Merci de nouveau

[feynman]

et alors aucune suggestion?