function - raiz - problemas resueltos de programación en fortran 95 pdf

Cómo alias un nombre de función en Fortran

Question

No estoy seguro si el título está bien puesto. Sugerencias bienvenidas.

Esto es lo que quiero hacer. Verifique una condición y luego decida qué función usar en un bucle. Por ejemplo:

if (a < 0) then
    loop_func = func1
else
    loop_func = func2
endif

Entonces puedo usar loop_func como un puntero al escribir mi ciclo. Ambas funciones toman exactamente las mismas entradas, y son diferentes enfoques para abordar el problema en función del valor de a . Esto me permitirá tener solo un bloque de código, en lugar de dos bloques casi idénticos. Esto también podría aplicarse a las subrutinas.

¿Alguna idea de cómo esto podría implementarse?

Gracias.

Answer 1

La mayoría de las implementaciones de Fortran no tienen una forma estándar de manipular punteros de función o punteros de procedimiento. Sin embargo, Fortran 2003 y posteriores tienen algo. (Consulte la página 6 de this ).

Para la situación dada, esto funcionará bastante bien en su lugar:

 function func1 (p1, p2, etc)
 ...  as you have it already
 end

 function func2 (p1, p2, etc)
 ...  as you have it already
 end

 function funcselect (a, p1, p2, etc)
     if (a < 0) then
          x = func1 (p1, p2, etc)
     else
          x = func2 (p1, p2, etc)
     endif
 end

Luego simplemente llame a funcselect con el parámetro adicional en lugar de lo que habría hecho con loop_func .

Answer 2

Sí, Fortran tiene punteros de procedimiento, por lo que puede, en efecto, alias un nombre de función. Aquí hay un ejemplo de código que asigna al puntero de función "f_ptr" una función u otra. A partir de entonces, el programa puede usar "f_ptr" y la función seleccionada será invocada.

module ExampleFuncs

   implicit none

contains

function f1 (x)
  real :: f1
  real, intent (in) :: x

  f1 = 2.0 * x

  return
end function f1


function f2 (x)
   real :: f2
   real, intent (in) :: x

   f2 = 3.0 * x**2

   return
end function f2

end module ExampleFuncs


program test_func_ptrs

    use ExampleFuncs
    implicit none

   abstract interface
      function func (z)
         real :: func
         real, intent (in) :: z
      end function func
   end interface

   procedure (func), pointer :: f_ptr => null ()

   real :: input

   write (*, ''( / "Input test value: ")'', advance="no" )
   read (*, *) input

   if ( input < 0 ) then
      f_ptr => f1
   else
      f_ptr => f2
   end if

   write (*, ''(/ "evaluate function: ", ES14.4 )'' )  f_ptr (input)

   stop

end program test_func_ptrs