Chris@10: Previous: Fortran Examples, Chris@10: Up: Calling FFTW from Legacy Fortran Chris@10:
In this section, we discuss how one can import/export FFTW wisdom Chris@10: (saved plans) to/from a Fortran program; we assume that the reader is Chris@10: already familiar with wisdom, as described in Words of Wisdom-Saving Plans. Chris@10: Chris@10:
The basic problem is that is difficult to (portably) pass files and
Chris@10: strings between Fortran and C, so we cannot provide a direct Fortran
Chris@10: equivalent to the fftw_export_wisdom_to_file, etcetera,
Chris@10: functions. Fortran interfaces are provided for the functions
Chris@10: that do not take file/string arguments, however:
Chris@10: dfftw_import_system_wisdom, dfftw_import_wisdom,
Chris@10: dfftw_export_wisdom, and dfftw_forget_wisdom.
Chris@10:
Chris@10:
Chris@10:
So, for example, to import the system-wide wisdom, you would do: Chris@10: Chris@10:
integer isuccess Chris@10: call dfftw_import_system_wisdom(isuccess) Chris@10:Chris@10:
As usual, the C return value is turned into a first parameter;
Chris@10: isuccess is non-zero on success and zero on failure (e.g. if
Chris@10: there is no system wisdom installed).
Chris@10:
Chris@10:
If you want to import/export wisdom from/to an arbitrary file or
Chris@10: elsewhere, you can employ the generic dfftw_import_wisdom and
Chris@10: dfftw_export_wisdom functions, for which you must supply a
Chris@10: subroutine to read/write one character at a time. The FFTW package
Chris@10: contains an example file doc/f77_wisdom.f demonstrating how to
Chris@10: implement import_wisdom_from_file and
Chris@10: export_wisdom_to_file subroutines in this way. (These routines
Chris@10: cannot be compiled into the FFTW library itself, lest all FFTW-using
Chris@10: programs be required to link with the Fortran I/O library.)
Chris@10:
Chris@10:
Chris@10: