Discussion :

[solo12]

Bonjour à tous;
j'ai essayer d'executer un programme en fortran qui me fait des icrementations en 3 directions càd un maillage dans un cylindre creux,les résultats que j'obtient sont toujours des 0 je n'arrive pas à trouver où est ce l'erreur pouvez vous m'aider?

voici mon prog:

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
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
 
	DIMENSION  X(100,100,100),Y(100,100,100),Z(100,100,100)
 
 
		pi=4.0*atan(1.0)
		 ni=80;       nj=80;      nk=80 
	m=ni+1
	n=nj+1	   
	zp=1.5	            
 
	 zl=5.
	r1=8.
	r0=4.
	dth=2*pi/float(nj)
 
	do 100 k=1,nk
 
	   do 200 j=1,n
 
	   x(1,j,k)=(r0+(r1-r0)*(j-1)/float(nj))
	y(1,j,k)=0.0
	x(m,j,k)=x(1,j,k)
	y(m,j,k)=y(1,j,k)
  200 continue
      do 300 i=1,m
	theta=(1-i)*dth+2.0*pi
	x(i,1,k)=r0*cos(theta)
	y(i,1,k)=r0*sin(theta)
	x(i,n,k)=r1*cos(theta)
	y(i,n,k)=r1*sin(theta)
  300 
 
	continue
	z(1,1,1)=0.0
	z(i,j,k)=zl*((k-1)/nk)**zp
 
  100	continue
 
 
	    do 110 i=1,m
  110   write(*,*)x(i,j,k)
       do 120 j=1,n
  120  write(*,*) y(i,j,k)
        do 130 k=1,nk
  130  write(*,*)z(i,j,k) 
    	WRITE(*,*)x(i,j,k),y(i,j,k),z(i,j,k)

[yogitetradim]

bonjour,

une piste (pas vérifiée...désolé) : ta déclaration de tableau semble ne pas indiquer le type de X,Y et Z. Ainsi, si le tableau est enregistré comme un tableau d'entiers, il se pourrait que tes 0 correspondent à la partie entière de ton résultat.

Quelquechose comme celà devrait être plus précis :

Code :

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

1
2
 
real , dimension(1:100,1:100,1:100) :: x,y,z

De plus, un conseil utilise le" Implicit none" , ça permet de bien vérifier que toutes tes variables ont bien été déclarées, ça pourrait permettre d'identifier des causes d'erreurs dans ton programme.
En espérant que ceci ait pu t'aider ;-)

[gregory2526]

il me semble que le type implite de ces tableaux est REAL
(si je ne me trompe pas, c'est les variables dont la première lettre est entre i et n qui sont des INTEGER)

mais effectivement, l'utilisation de "IMPLICIT NONE" permet d'une part de bien contrôler le type de toutes les variables, mais aussi d'éviter les erreurs de frappe dans le nom d'une variable.

Pour ton problème, tu devrai afficher les résultats intermédiaires (en réduisant les ni, nj et nk pour que ce soit plus simple à débugger). Tu verrai ainsi a quel moment le calcul devient faux.

[solo12]

Bonjour à tous;
merci yogitetradim et gregory2526 de votre aide oui j'ai déclaré mes variables comme real ca a marché pour les deux premières composantes mais la troisième n'est pas calculé

[solo12]

Bonjour à tous;
merci yogitetradim et gregory2526 de votre aide oui j'ai déclaré mes variables comme real ca a marché pour les deux premières composantes mais la troisième n'est pas calculé (le prog est toujours executable)

[FR119492]

Salut !

Lorsque tu testes un programme, affiche les résultats intermédiaires en format E pour distinguer les valeurs très petites de zéro.

D'autre part, on obtient souvent de bons résultats avec la méthode de la bissection: Tu affiches des valeurs intermédiaires à peu près au milieu de ton programme. Si ces résultats sont faux, l'erreur se trouve dans la première moitié; s'ils sont justes,c'est dans la seconde. Puis tu recommences avec la moitié qui contient l'erreur, et ainsi de suite; ça ne marche pas toujours, mais assez souvent et c'est assez rapide.

Jean-Marc Blanc