Discussion :

[pascale_92]

J'ai encapsulé un sous programme fortran dans un programme C. Le problème subsistant est le passage de tableaux multidimensionnels en C à des tableaux fortran.

[DaZumba]

C'est penible a faire. Les tableaux multi-dimensionnels ne sont pas classes de la meme facon en memoire en Fortran et en C.
Les tableaux Fortran sont en column-major, i.e. la premiere dimension varie la plus rapidement. Les tableaux C sont en row-major, i.e. la dimension la plus a droite varie le plus rapidement.

Le plus simple est certainement de considerer les tableaux comme ayant une seule dimension (linearisation) et d'utiliser les formules suivantes pour acceder aux elements (exemple en deux dimensions):
En Fortran:

Code :

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

1
2
 
index  = (( col  - 1) * Nrow ) + row

ou les indices vont de 1 a Nrow*Ncol.

En C:

Code :

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

1
2
 
index  = ( row  * Ncol ) + col

ou les indices vont de 0 a Nrow*Ncol - 1.

[pascale_92]

A plusieurs dimensions on fait comment ?

[DaZumba]

Ben, il est assez simple de generaliser les formules ci-dessus.
Par exemple, pour un tableau a trois dimensions N1, N2, N3, pour acceder a l'element i,j,k, on fait en Fortran:

Code :

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

1
2
 
index =  ((k-1) * N2 + j - 1) * N1 + i

avec 1 <= i <= N1, 1 <= j <= N2, 1 <= k <= N3.

en C:

Code :

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

1
2
 
index = (i * N2 + j)*N3 + k;

avec 0 <= i < N1, 0<= j < N2, 0 <= k < N3.

Mais le plus simple est tout de meme de ne jamais utiliser de tableaux de plus de 1 dimension entre C et Fortran, et d'utiliser une formule commune pour 'simuler' un tableau multi-dimensionnel dans les deux langages (en faisant attention au decalage de 1 dans les indices).

[pascale_92]

Il me semble plus simple de passer d'un tableau fortran à un tableau C simplement en inversant les indices non?

J'ai testé en dimension 3...

[DaZumba]

Faire des tests est effectivement necessaire. Apres, cela depend de qui appelle qui. En effet, tu peux imposer au C de suivre une organisation 'a la Fortran'.
Mais personnellement, je ne m'embete pas et je travaille uniquement avec des tableaux a une dimension. Les tableaux multi-dimensionnels sont linearises en jouant avec les indices.