Discussion :

[feynman]

Salut,
j'ai un petit problème dans l'instruction if.
J'ai une matrice de taille NxN, avec si i=N/2 et (en même temps) j=n/2+1, alors A(i,j)=0., sinon on a que A(i,j)=1./(i-N/2.+k*(N/2-1-j)) avec k c'est le complexe i .
J'ai eu la place 0 à la place que je veux, mais j'ai eu des 0 dans la dernière ligne et la troisième colonne !
J'ai essayé de trouver l'erreur mais ça me semble vraiment bizarre, pouvez-vous m'aider svp ?

Voilà la matrice que j'ai eu avec N=4

Code :

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

1
2
3
4
5
  A=
  -0.40000001  -0.20000000  -0.50000000  -0.50000000   0.00000000   0.00000000   0.50000000  -0.50000000
  -0.50000000   0.00000000  -1.00000000   0.00000000   0.00000000   0.00000000   1.00000000   0.00000000
  -0.40000001   0.20000000  -0.50000000   0.50000000   0.00000000   0.00000000   0.50000000   0.50000000
   0.00000000   0.00000000   0.00000000   0.00000000   0.00000000   0.00000000   0.00000000   0.00000000

Et voilà mon code

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
 
      program fft
      integer i,j,N
      parameter(N=4)
      complex A(0:N-1,0:N-1),ii
      parameter(ii=(0.,1.))
 
 
 
 
      do i=0,N-1
         do j=0,N-1
            if((i.ne.(int(real(N)/2.))).and.
&                     (j.ne.(int(real(N)/2.)+1))) then
               A(i,j)=1./(real(i)-real(N)/2.+ii*(real(N)/2.-1.-real(j)))
             else
               A(i,j)=(0.,0.)
            endif
         enddo
      enddo
 
      print 10,((A(i,j),i=0,N-1),j=0,N-1)
10    format('A=',8(/,1X,8(F12.8,1X)))
 
      end

[FR119492]

Salut !

kekseksa: int(real(N)/2. ???

Jean-Marc Blanc

[feynman]

Salut !

kekseksa: int(real(N)/2. ???

Jean-Marc Blanc

c est la fonction int c est a dire la partie entière d un nombre reel .

[feynman]

Moi je veux mettre un element de la matrice nul et non pas des colones et des lignes, mon programme rend nul la ligne et la colone qui contient cet element

[VITALTH]

Bonjour,
Tu veux

Code :

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

1
2
3
4
5
si conditon1= vrai et si condition2= vrai  alors
   instructIon1
sinon (donc complément de : condition1 = vrai et condition2=vrai)
   instruction2
fin si

Or le complément de condition1 = vrai et condition2= vrai c'est condition1= faux ou condition2 = faux alors que tu as écrit condition1=faux et condition2=faux.

Dans ton code là où tu as écrit :

Code :

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

1
2
3
4
5
6
if((i.ne.(int(real(N)/2.))).and.
&                     (j.ne.(int(real(N)/2.)+1))) then
A(i,j)=1./(real(i)-real(N)/2.+ii*(real(N)/2.-1.-real(j)))
else
      A(i,j)=(0.,0.)
endif

J'aurais écrit

Code :

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

1
2
3
4
5
6
7
if((i.ne.(int(real(N)/2.))).or.
&                     (j.ne.(int(real(N)/2.)+1))) then
 
A(i,j)=1./(real(i)-real(N)/2.+ii*(real(N)/2.-1.-real(j)))
else
      A(i,j)=(0.,0.)
endif

ou

Code :

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

1
2
3
4
5
6
if((i.eq.(int(real(N)/2.))).and.
&                     (j.eq.(int(real(N)/2.)+1))) then
A(i,j)=(0.,0.)
else
A(i,j)=1./(real(i)-real(N)/2.+ii*(real(N)/2.-1.-real(j)))
      endif

Ou alors je n'ai pas bien compris ta requête

[gregory2526]

int(real(N)/2.

c est la fonction int c est a dire la partie entière d un nombre reel .

tu peux expliquer là ?

N est un entier
tu le transforme en réel pour le diviser
et tu reprend la partie entière ???

ça me parrait un peu étrange tout de même...

[feynman]

Moi je veux creer la matrice suivante : Avec N etant un nombre paire
si i=N/2 et (en même temps) j=N/2+1, alors A(i,j)=0., sinon on a que A(i,j)=1./(i-N/2.+k*(N/2-1-j))

En d autre terme je veux seulement l element A(N/2,N/2-1)=0.

[feynman]

voila la matrice but que je veux

Code :

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

1
2
3
4
5
6
         A  =
 
  - 0.4 - 0.2i    - 0.5 - 0.5i  - i      0.5 - 0.5i  
  - 0.5           - 1.            0      1.          
  - 0.4 + 0.2i    - 0.5 + 0.5i    i      0.5 + 0.5i  
  - 0.25 + 0.25i  - 0.2 + 0.4i    i      0.2 + 0.4i

[feynman]

Bonjour,
Tu veux
si conditon1= vrai et si condition2= vrai alors
instruciton1
sinon (donc complément de : condition1 = vrai et condition2=vrai)
instruction2
fin si

Or le complément de condition1 = vrai et condition2= vrai c'est condition1= faux ou condition2 = faux alors que tu as écrit condition1=faux et condition2=faux.

Dans ton code là où tu as écrit :
if((i.ne.(int(real(N)/2.))).and.
& (j.ne.(int(real(N)/2.)+1))) then
A(i,j)=1./(real(i)-real(N)/2.+ii*(real(N)/2.-1.-real(j)))
else
A(i,j)=(0.,0.)
endif

J'aurais écrit
if((i.ne.(int(real(N)/2.))).or.
& (j.ne.(int(real(N)/2.)+1))) then

A(i,j)=1./(real(i)-real(N)/2.+ii*(real(N)/2.-1.-real(j)))
else
A(i,j)=(0.,0.)
endif


ou

if((i.eq.(int(real(N)/2.))).and.
& (j.eq.(int(real(N)/2.)+1))) then
A(i,j)=(0.,0.)
else
A(i,j)=1./(real(i)-real(N)/2.+ii*(real(N)/2.-1.-real(j)))
endif

Ou alors je n'ai pas bien compris ta requête

Merci VITALTH , mais j ai obtenu avec tes proposition la matrice suivante :

Code :

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

1
2
3
4
5
6
 
A=
  -0.40000001  -0.20000000  -0.50000000  -0.50000000   0.00000000  -1.00000000   0.50000000  -0.50000000
  -0.50000000   0.00000000  -1.00000000   0.00000000          nan          nan   1.00000000   0.00000000
  -0.40000001   0.20000000  -0.50000000   0.50000000   0.00000000   1.00000000   0.50000000   0.50000000
  -0.25000000   0.25000000  -0.20000000   0.40000001   0.00000000   0.00000000   0.20000000   0.40000001

[feynman]

sayez c est résolu , merci beacoups pour votre aide,
vraiment sans votre aide et patience je ne pourrai y arriver