Section d'un tableau en paramètre d'une fonction

Version imprimable

06/04/2007, 23h28
ZeLL

Section d'un tableau en paramètre d'une fonction
Bonjour à tous,

Voila mon problème,
J'ai pour global un tableau "tres" grand, permetant de gerer toute les données de mon application:

Code:

REAL WS(COM_SIZE)

Dans mon main, j'apelle une fonction foo:

Code:

CALL FOO(WS(OFFSET))

Dans ma fonction FOO, ce paramètre corespond a un tableau a 5 dimensions:
Code:

1 2 3 4 SUBROUTINE FOO(FIELDS) /*input/output*/ REAL FIELDS(MXLON,NLAT,NVER,NF,NLAYERS)
La shape de FIELDS corespond a une signification logique au niveau des données.
...La tout va bien dans le meilleur des mondes

Maintenant je veut appeller une autre fonction a partir de mon main,

Code:

CALL BAR(WS(OFFSET))

le meme tableau en fait, mais pas toutes les données,
du point de vu logique je veut passer:

Code:

FIELDS(1:NLON,NLAT,NVER,NF,NLAYERS)

avec NLON<MXLON

Ce qui corespont pour mon tableaux WS, (de rang 1) a une sous section non contigus.
Je peut, par l'intermediaire d'un vecteur d'indice, faire une indexation directe
(i.e. je fait un vecteur avec tous les indices que je veux voir dans ma sous section)
Le probleme, c'est qu'un tel sous-tableau ne peut pas ﬁgurer à gauche d’un signe d’affectation, donc le passage en parametre ne marche pas non plus.

De plus j'ai lu dans un tutorial:

Citation:

Sections non régulières en argument de procédures
Si on passe une section non régulière en argument d’appel d’une
procédure, il faut savoir que :
• c’est une copie contiguë en mémoire qui est passée par le
compilateur,
• l’argument muet correspondant de la procédure ne doit pas avoir la
vocation INTENT(inout) ou INTENT(out) ; autrement dit, en
retour de la procédure, il n’y a pas mise-à-jour du tableau “père” de la
section irrégulière.

Toutes les remarques seront les bienvenues
12/04/2007, 14h38
DaZumba
Je te conseille de travailler sur une copie locale dans la fonction FOO, avec recopie des elements necessaires avant l'appel a BAR. Ainsi, dans FOO et BAR tu auras
Code:

1 2 REAL L_FIELDS(NLON,NLAT,NVER,NF,NLAYERS)
ce qui simplifie grandement les choses.
12/04/2007, 15h34
ZeLL

Merci pour ta reponse,

C'est actuellement ce que je fais, et ca marche tres bien.
Mais le problème c'est que c'est du code scientific, et j'ai besoin des meilleurs
performances possibles.
L'utilisation d'un tableau "swap", bien que tres pratique, me coute bien trop cher en temps et en memoire (jusqu'a 1E6 element de type double precisions).
16/04/2007, 08h00
ZeLL

Pas d'idées?
J commence à croire que c'est pas possible :-(
j trouve rien sur le net, meme avec du f95

J'ai trouvé !

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
26
 
program Test
...
REAL WS(COM_SIZE)
...
call InterfaceBar(WS(Offset))
...
end
subroutine InterfaceBar(FIELDS)
   /*input/output*/
    REAL FIELDS(MXLON,NLAT,NVER,NF,NLAYERS)
 
   interface
      subroutine Bar(a)
         real, intent(inout) :: a(:,:,:,:,:)
      end subroutine
   end interface
 
   call Bar(Fields(1:NLON,NLAT,NVER,NF,NLAYERS))
   end subroutine
 
subroutine Bar(Fields)
  implicit none
  real,intent(inout) :: Fields(:,:,:,:,:)
 
...

18/04/2007, 14h49
ZeLL

Je n'ai qu'une chose à dire, bien joué:king:

Fallait y penser, ca parait tellement trivial quand la réponse apparait comme
par magie.

Le seul hic c'est que dans mon projet ca marche pas, mais la fonction que j'appelle est un petit peu plus compliquée.
En fait Bar(), vas juste enregistrer la variable pour chaque process (programation parallele) selon une distribution donnée. Puis a l appelle d'une autre fonction, les variable enregistrée (comme FIELDS par example) sont sauvegarder a travers le grid via un logicielle appeller Internet Backplane Protocol, et c'est la que ca marche pas.
Bref je suis presque sur que ta solution fonctionne a merveille, et si ca marche pas chez moi, c'est que mon programme est assez complex ( de plus Bar est ecrit en c ).
Je suis maintenant le seul a pouvoir m'aider :aie:

En tout cas merci Sylvain pour ta soluce,