Problème de pointeurs avec un tableau

**Tom_Tom_RoF** · 05/10/2007, 16h11

Bonjour,

J'essaye de mettre les coordonees d'un poyedre compose de triangles dans deux tableaux.
Le premier contient les coordonees de chaque sommet et le second indique les sommets respectifs pour chaque face.
Je pensais utiliser des pointeurs dans le second tableau pour indiquer quel sont les points a utiliser.
Malheuresement

le compilateur (absoft) refuse cette utilisation des pointeurs:
"the pointer-object in a pointer assignment statement must have the pointer attribute."

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
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
 
 
real, target, dimension(3,4) :: vertex_1 ! coordonnee des points
real, target, dimension(3,4) :: vertex_2 
 
real, pointer, dimension (3,4) :: face_1 ! points pour chaque face
real, pointer, dimension (3,4) :: face_2
 
!vertex 1 - ployedron 1
Vertex_1 (1,1) = 0
Vertex_1 (2,1) = 5
Vertex_1 (3,1) = 5
 
!vertex 2 - ployedron 1
vertex_1 (1,2) = 0
vertex_1 (2,2) = 0
vertex_1 (3,2) = 7
 
!vertex 3 - ployedron 1
vertex_1 (1,3) = 5
vertex_1 (2,3) = 0
vertex_1 (3,3) = 0
 
!vertex 4 - ployedron 1
vertex_1 (1,4) = 0
vertex_1 (2,4) = 0
vertex_1 (3,4) = 2
 
!vertex 1 - ployedron 2
Vertex_2 (1,1) = 0
Vertex_2 (2,1) = 7
Vertex_2 (3,1) = 5
 
!vertex 2 - ployedron 2
vertex_2 (1,2) = 0
vertex_2 (2,2) = 3
vertex_2 (3,2) = 7
 
!vertex 3 - ployedron 2
vertex_2 (1,3) = 5
vertex_2 (2,3) = 3
vertex_2 (3,3) = 0
 
!vertex 4 - ployedron 2
vertex_2 (1,4) = 0
vertex_2 (2,4) = 3
vertex_2 (3,4) = 2
 
!face 1 - ployedron 1
face_1 (1,1) => Vertex_1 (1,:) ! ca pose probleme a partir d'ici
face_1 (2,1) => vertex_1 (2,:)
face_1 (3,1) => vertex_1 (3,:)
 
!face 2 - ployedron 1
face_1 (1,2) => Vertex_1 (1,:)
face_1 (2,2) => vertex_1 (3,:)
face_1 (3,2) => vertex_1 (4,:)
 
!face 3 - ployedron 1
face_1 (1,3) => Vertex_1 (1,:)
face_1 (2,3) => vertex_1 (2,:)
face_1 (3,3) => vertex_1 (4,:)
 
!face 4 - ployedron 1
face_1 (1,4) => vertex_1 (2,:)
face_1 (2,4) => vertex_1 (3,:)
face_1 (3,4) => vertex_1 (4,:)
 
!face 1 - ployedron 2
face_2 (1,1) => vertex_2 (1,:)
face_2 (2,1) => vertex_2 (2,:)
face_2 (3,1) => vertex_2 (3,:)
 
!face 2 - ployedron 2
face_2 (1,2) => vertex_2 (1,:)
face_2 (2,2) => vertex_2 (3,:)
face_2 (3,2) => vertex_2 (4,:)
 
!face 3 - ployedron 2
face_2 (1,3) => vertex_2 (1,:)
face_2 (2,3) => vertex_2 (2,:)
face_2 (3,3) => vertex_2 (4,:)
 
!face 4 - ployedron 2
face_2 (1,4) => vertex_2 (2,:)
face_2 (2,4) => vertex_2 (3,:)
face_2 (3,4) => vertex_2 (4,:)
 
!pour apres l'utiliser ici:
!tester chaque arrete du premier polyedre avec chaque face du second polyedre
do m = 1, 4 ! 4 est le nombre d'arretes du poly 1
 
	do i = 1, 4 ! 4 est le nombre de faces du poly 2
 
		call ray_intersecting_polygon (face_2 (1,i), face_2 (2,i), face_2 (3,i), face_1 (1,m), face_1(2,m), collision)
		call ray_intersecting_polygon (face_2 (1,i), face_2 (2,i), face_2 (3,i), face_1 (2,m), face_1(3,m), collision)
		call ray_intersecting_polygon (face_2 (1,i), face_2 (2,i), face_2 (3,i), face_1 (1,m), face_1(3,m), collision)
 
	enddo
 
enddo

Comment resoudre mon probleme?

Merci.

**FR119492** · 05/10/2007, 16h27

Salut!

Pendant longtemps, les pointeurs ont été inconnus en Fortran; il semblent qu'ils n'ont été introduits que pour satisfaire les besoins de ceux qui les avaient utilisés dans d'autres langages, mais on peut très bien s'en passer. Je me souviens d'un cours du prof. Ch. Rapin, de l'EPF-Lausanne, qui traitait toutes les structures d'information (piles, queues, listes diverses) à l'aide de tableaux et d'indices qui étaient des entiers.

Pour ton problème, j'essaierais par exemple
Real*8 Vertex(1000,3)
Integer Face(2000,3)
en admettant que tu limites ton programme à 1000 sommets et 2000 faces.

Bonne chance
Jean-Marc Blanc

**feynman** · 05/10/2007, 16h31

Oui c est vrai Jean Marc,
Dans tout les cours que j ai vu jusqu a present en Fortran, j ai pas vu la notion de Pointeurs la dedant,
tu connais une reference qui traite les pointeur en fortran?
merci

**FR119492** · 05/10/2007, 16h42

Salut!

Pendant longtemps, les pointeurs ont été inconnus en Fortran; il semblent qu'ils n'ont été introduits que pour satisfaire les besoins de ceux qui les avaient utilisés dans d'autres langages, mais on peut très bien s'en passer. Je me souviens d'un cours du prof. Ch. Rapin, de l'EPF-Lausanne, qui traitait toutes les structures d'information (piles, queues, listes diverses) à l'aide de tableaux et d'indices qui étaient des entiers.

Pour ton problème, j'essaierais par exemple
Real*8 Vertex(1000,3)
Integer Face(2000,3)
en admettant que tu limites ton programme à 1000 sommets et 2000 facec.

Bonne chance
Jean-Marc Blanc

**Tom_Tom_RoF** · 05/10/2007, 16h45

Envoyé par feynman

Oui c est vrai Jean Marc,
Dans tout les cours que j ai vu jusqu a present en Fortran, j ai pas vu la notion de Pointeurs la dedant,
tu connais une reference qui traite les pointeur en fortran?
merci

Tu pourra trouver ca dans les cours de IDRIS cours 95

http://gershwin.ens.fr/vdaniel/Doc-L...rs.html#Chap-7

Mais j'ai pas compris grand chose. Il doit me manquer quelques connaissances.

**Sylvain Bergeron** · 05/10/2007, 17h00

Sans discuter de la pertinence ou non d'utiliser des pointeurs, ton problème semble surtout être un problème de rang. Dans ton code :

Code :

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

1
2
3
4
5
6
 
real, target, dimension(3,4) :: vertex_1 ! coordonnee des points
...
real, pointer, dimension (3,4) :: face_1 ! points pour chaque face
...
face_1 (1,1) => Vertex_1 (1,:)

Tu essaie d'utiliser un pointeur réel scalaire pour pointer sur un vecteur...

Correction

Mauvais diagnostic (voir la réponse suivante de SB). Comme on ne peut pas avoir de vecteur de pointeur, on ne peut pas avoir Pointeur(indice) => ...

**Tom_Tom_RoF** · 05/10/2007, 17h10

Envoyé par Sylvain Bergeron

Tu essaie d'utiliser un pointeur réel scalaire pour pointer sur un vecteur...

Oui, c'est exactement ca mon probleme!!!

Donc si je pointe sur l'adresse et non pas sur la valeur ca devrai fonctionner.
Je vais essayer de faire ca.

Merci de toutes vos reponses.

**Sylvain Bergeron** · 05/10/2007, 17h22

Envoyé par Tom_Tom_RoF

Donc si je pointe sur l'adresse et non pas sur la valeur ca devrai fonctionner.

Ça sent le programmeur C... Un pointeur Fortran, c'est plus qu'une adresse.

Mais, Ah, j'ai trouvé où ça coince :

En Fortran, tu ne peux pas avoir un vecteur de pointeur. Ce n'est pas prévu par le langage. La solution est de passer par des types :

Code :

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

1
2
3
4
5
6
7
8
9
10
 
type TPtrPoint
   real, pointer :: Ptr(:)
end type
type (TPtrPoint) :: Face(n) ! J'ai maintenant un vecteur de n pointeur (pointant sur des points 3D)
 
...
 
Face(k)%Ptr => Vertex(j,:)
...

**Tom_Tom_RoF** · 05/10/2007, 17h34

Envoyé par Sylvain Bergeron

Ça sent le programmeur C

Ce serai plutot un programmeur qui a juste approcher le C et qui doit faire un gros projet en fortran!!

Sinon ta solution est a mon avis ce qu'il me faut.
Reste a bien tout comprendre et a l'interger a mon code.

Merci du coup de main!

PS: je metterai en resolu quand je serai sur que c'est le cas.

PPS: a a l'air bon!!! Merci encore

Problème de pointeurs avec un tableau

Fortran

Discussions similaires

Partager

Partager