Passer un vrai tableau si un ensemble complexe devrait

Question

J'ai un sous-programme Fortran qui attend un ensemble complexe comme

subroutine foo(cnumbers, n) 
    integer :: n 
    complex :: cnumbers(n) 
    ... 
end subroutine foo 

et plus tard, je veux l'appeler comme

real :: rnumbers(40) 
... 
call foo(rnumbers, 20) 

Cependant, je reçois l'erreur du compilateur:

error #6633: The type of the actual argument differs from the type of the dummy argument.

Bien sûr, cela est compréhensible car un tableau réel n'est pas un ensemble complexe. Mais il doit y avoir un moyen de le faire fonctionner.

Parce que si le sous-programme foo et l'appel de foo sont différents modules et sont écrits dans différents fichiers Fortran, le compilateur ne pas se plaindre, et tout fonctionne très bien.

Est-ce que quelqu'un sait comment le faire fonctionner? Comment passer un tableau réel si un tableau complexe est attendu?

Answer 1

Vous pouvez utiliser transfer(rnumbers, ...) pour convertir le type (un tableau temporaire est susceptible d'être créé) ou utiliser l'équivalence pour éviter

real :: rnumbers(40) 
    complex :: cnumbers(20) 
    equivalence (rnumbers, cnumbers) 

    set the value of rnumbers 

    call foo(cnumbers, 20) 

Si vous avez besoin des tableaux affectables l'équivalence ne fonctionnera pas.

Vous pouvez également utiliser un sous-programme externe et placer le compilateur sur l'interface et simplement passer le tableau réel au lieu du complexe. Ce n'est pas conforme aux normes, mais il est parfois utilisé. Voir aussi Gfortran complex actual to real dummy argument