L'instruction if

Version imprimable

L'instruction if

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:

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:

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

26/03/2008, 23h17
FR119492

Salut !

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

Jean-Marc Blanc
27/03/2008, 07h57
feynman

Citation:

Envoyé par FR119492

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 .
27/03/2008, 09h11
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

Bonjour,
Tu veux

Code:

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:

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:

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

Code:

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

27/03/2008, 09h34
gregory2526

Citation:

int(real(N)/2.

Citation:

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...
27/03/2008, 09h47
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.

voila la matrice but que je veux

Code:

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

Citation:

Envoyé par VITALTH

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:

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

27/03/2008, 10h00
feynman

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