Utilisation de subroutine Fortran 77 en Fortran 90

Bonsoir,

Je travaille actuellement sur un code Fortran 90 (pour une résolution de système linéaire). Ce code fonctionne bien, mais reste lent.
Comme le système à résoudre est très creux, je cherche à faire évoluer le code pour qu'il puisse travailler avec des matrices creuses.

Après quelques recherches, j'ai récupéré des sources dans ce but. Elles sont écrites en Fortran 77, et on y trouve différentes subroutines à appeler.

Comment est ce que je peux m'y prendre pour pouvoir les utiliser depuis mon code Fortran 90 ?

Merci d'avance

Tu n'as rien d'autre à faire que de les compiler. À part quelques exceptions près, le Fortran 77 est un sous-ensemble du Fortran 90/95. Les compilateurs Fortran 90 supportent toujours le sous-ensemble Fortran 77.

Et la différence de format (fixe ou libre) n'en est pas une. Le Fortran 77 ne permettait que le format fixe, alors que depuis Fortran 90, on peut utiliser au choix le format fixe ou le format libre. La seule contrainte est qu'un fichier ne peut pas contenir des sections en format fixe et des sections en format variable. Tu devras donc avoir au moins 2 fichiers sources.

Merci pour la réponse rapide.

Les deux versions de Fortran sont bien dans des fichiers séparés, donc pas de soucis de ce côté là.

Dans le fichier principal (celui en Fortran 90), quand j'appelais des routines d'autres modules Fortran 90, je devais mettre un use en début de fichier. N'y a-t-il pas une syntaxe similaire à respecter ?

Use ? Tu ne mets pas un use pour pouvoir utiliser des routines F90. Tu mets un use si tu veux utiliser des routines membres d'un module. Un module est une organisation du code possible en F90, mais qui n'est pas obligatoire (bien que fortement suggéré).

Dans ton cas, comme tes routines F77 ne sont (probablement) pas dans un module, tu n'as pas à faire un use module pour les utiliser.

Cela dit, tu pourrais peut-être "emballer" tes routines F77 dans un module pour faciliter leur utilisation. Il y a cependant un risque : les routines de résolution numérique ont souvent en argument des vecteurs (1D) auxquels des tableaux (2D) sont associés et ce genre de pirouette résiste assez mal aux interfaces explicites générées par les modules...