Discussion :

[dva2tlse]

Bonjour, bonsoir le forum,
ce post fait suite à ces autres : http://www.developpez.net/forums/d13...n-c-voire-cpp/ et http://forum.ubuntu-fr.org/viewtopic.php?id=1310981, où je demandais quel langage de programmation choisir entre le fortran et le C voire le C++. J'ai finalement choisi de coder en fortran, parce que c'est le langage que je connais le mieux, et que j'envisage par la suite que la plupart du travail soit effectuée en parallèle sur plusieurs processeurs. (il y a des serveurs de calcul qui tournent déjà sur quarante proc's dans ma boite, pour d'autres applications)

Le travail en question est la traduction d'un fichier d'entrée d'environ deux millions de lignes, en une suite de valeurs numériques déterminées d'après les informations de chaque ligne, et avec d'autres infos que l'on récupère dans deux autres fichiers d'entrée de quelques centaines de lignes chacun; on m'a informé sur ce forum ou sur celui de linux-ubuntu où j'ai exposé le problème aussi, que le travail en question s'apparentait à un "jointure" de bases de données.

La première fois où le programme a tourné correctement jusqu'au bout, il a mis environ vingt minutes; suite aux discussions sur les forums, j'en suis arrivé à deux minutes environ (en améliorant sensiblement la logique de l'algo et les I/O), et aujourd'hui vendredi, en passant d'une vieille HP9000 à bout de souffle, à un bon PC sous linux avec un (seul) processeur itanium, j'en suis à 21 secondes.

Petite question, j'utilise Gfortran pour compiler mon programme, et il râle quand des indices de tableaux sont eux-mêmes des variables indicées.
Par exemple, je ne peux pas utiliser Ccvt1(iFScvt, M(iFScvt), T(iFScvt)), mais je dois passer par des variables intemédiaires,Mtmp=M(iFScvt), et Ttmp=T(iFScvt), pour pouvoir finalement accéder à Ccvt1(iFScvt, Mtmp , Ttmp). Y a-t'il un moyen de contourner cela ?

Merci de me donner toutes les idées, même farfelues, que vous inspire ce post', car n'étant pas informaticien mais mécanicien (en calcul de structures), mon imagination et mes connaissances sont assez pauvres dans ces domaines,
merci encore,

[Matthieu Brucher]

que j'envisage par la suite que la plupart du travail soit effectuée en parallèle sur plusieurs processeurs. (il y a des serveurs de calcul qui tournent déjà sur quarante proc's dans ma boite, pour d'autres applications)

Rien à voir avec le reste de ton message, mais c'est plus simple de faire ça en C ou en C++ qu'en Fortran (et je sais de quoi je parle).

[dva2tlse]

Bonjour Matthieu,
merci de ton conseil, mais je connais pas mal le fortran alors que je ne suis que débutant en C, donc je vais faire ce que je connais, et cela devrait m'être beaucoup plus facile pour les nombreux tableaux que j'ai à utiliser parce que je galère avec les pointeurs...
à +

[Sylvain Bergeron]

Petite question, j'utilise Gfortran pour compiler mon programme, et il râle quand des indices de tableaux sont eux-mêmes des variables indicées.
Par exemple, je ne peux pas utiliser Ccvt1(iFScvt, M(iFScvt), T(iFScvt)), mais je dois passer par des variables intemédiaires,Mtmp=M(iFScvt), et Ttmp=T(iFScvt), pour pouvoir finalement accéder à Ccvt1(iFScvt, Mtmp , Ttmp). Y a-t'il un moyen de contourner cela ?

Donne nous le message d'erreur exacte parce qu'une expression de la forme Tableau(Tableau(indice)) est tout à fait conforme.

[Matthieu Brucher]

Raler = erreur ou warning ? C'est peut-etre un warning disant qu'il ne peut pas vectoriser ce genre de cas ?

[Sylvain Bergeron]

Vectoriser au sens pipeliner ? Possible. Si c'est le cas, c'est un warning très soft. Des situations / constructions limitant la vectorisation sont très fréquentes et relèvent plus de l'analyse / suggestion que du warning...

[Matthieu Brucher]

Vectoriser au sens pipeliner ? Possible. Si c'est le cas, c'est un warning très soft. Des situations / constructions limitant la vectorisation sont très fréquentes et relèvent plus de l'analyse / suggestion que du warning...

Tout a fait d'accord.

[dva2tlse]

Bonjour, je viens de m'apercevoir ce matin en essayant de sélectionner les lignes qui faisaient planter la compilation, que c'était impossible.
Donc ça ne plante plus quand l'indice d'un tableau est lui même une valeur indicée, ce qui est une des choses qui rend intéressant de travailler en fortran.
Donc "Résolu", merci

[dva2tlse]

@Matthieu et Sylvain, je ne suis plus au boulot, donc demain dès que j'arrive je mets en commentaire les deux ou trois petites modif's de rien du tout que j'avais faites pour que la compil' aille jusqu'au bout. (Il le semble me souvenir qu'elle ne se terminait pas)
Merci et à demain,