|
Chris@42
|
1 c Copyright (c) 2003, 2007-14 Matteo Frigo
|
|
Chris@42
|
2 c Copyright (c) 2003, 2007-14 Massachusetts Institute of Technology
|
|
Chris@42
|
3 c
|
|
Chris@42
|
4 c This program is free software; you can redistribute it and/or modify
|
|
Chris@42
|
5 c it under the terms of the GNU General Public License as published by
|
|
Chris@42
|
6 c the Free Software Foundation; either version 2 of the License, or
|
|
Chris@42
|
7 c (at your option) any later version.
|
|
Chris@42
|
8 c
|
|
Chris@42
|
9 c This program is distributed in the hope that it will be useful,
|
|
Chris@42
|
10 c but WITHOUT ANY WARRANTY; without even the implied warranty of
|
|
Chris@42
|
11 c MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
|
|
Chris@42
|
12 c GNU General Public License for more details.
|
|
Chris@42
|
13 c
|
|
Chris@42
|
14 c You should have received a copy of the GNU General Public License
|
|
Chris@42
|
15 c along with this program; if not, write to the Free Software
|
|
Chris@42
|
16 c Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
|
|
Chris@42
|
17 c
|
|
Chris@42
|
18 cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
|
|
Chris@42
|
19 c
|
|
Chris@42
|
20 c This is an example implementation of Fortran wisdom export/import
|
|
Chris@42
|
21 c to/from a Fortran unit (file), exploiting the generic
|
|
Chris@42
|
22 c dfftw_export_wisdom/dfftw_import_wisdom functions.
|
|
Chris@42
|
23 c
|
|
Chris@42
|
24 c We cannot compile this file into the FFTW library itself, lest all
|
|
Chris@42
|
25 c FFTW-calling programs be required to link to the Fortran I/O
|
|
Chris@42
|
26 c libraries.
|
|
Chris@42
|
27 c
|
|
Chris@42
|
28 cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
|
|
Chris@42
|
29
|
|
Chris@42
|
30 c Strictly speaking, the '$' format specifier, which allows us to
|
|
Chris@42
|
31 c write a character without a trailing newline, is not standard F77.
|
|
Chris@42
|
32 c However, it seems to be a nearly universal extension.
|
|
Chris@42
|
33 subroutine write_char(c, iunit)
|
|
Chris@42
|
34 character c
|
|
Chris@42
|
35 integer iunit
|
|
Chris@42
|
36 write(iunit,321) c
|
|
Chris@42
|
37 321 format(a,$)
|
|
Chris@42
|
38 end
|
|
Chris@42
|
39
|
|
Chris@42
|
40 subroutine export_wisdom_to_file(iunit)
|
|
Chris@42
|
41 integer iunit
|
|
Chris@42
|
42 external write_char
|
|
Chris@42
|
43 call dfftw_export_wisdom(write_char, iunit)
|
|
Chris@42
|
44 end
|
|
Chris@42
|
45
|
|
Chris@42
|
46 c Fortran 77 does not have any portable way to read an arbitrary
|
|
Chris@42
|
47 c file one character at a time. The best alternative seems to be to
|
|
Chris@42
|
48 c read a whole line into a buffer, since for fftw-exported wisdom we
|
|
Chris@42
|
49 c can bound the line length. (If the file contains longer lines,
|
|
Chris@42
|
50 c then the lines will be truncated and the wisdom import should
|
|
Chris@42
|
51 c simply fail.) Ugh.
|
|
Chris@42
|
52 subroutine read_char(ic, iunit)
|
|
Chris@42
|
53 integer ic
|
|
Chris@42
|
54 integer iunit
|
|
Chris@42
|
55 character*256 buf
|
|
Chris@42
|
56 save buf
|
|
Chris@42
|
57 integer ibuf
|
|
Chris@42
|
58 data ibuf/257/
|
|
Chris@42
|
59 save ibuf
|
|
Chris@42
|
60 if (ibuf .lt. 257) then
|
|
Chris@42
|
61 ic = ichar(buf(ibuf:ibuf))
|
|
Chris@42
|
62 ibuf = ibuf + 1
|
|
Chris@42
|
63 return
|
|
Chris@42
|
64 endif
|
|
Chris@42
|
65 read(iunit,123,end=666) buf
|
|
Chris@42
|
66 ic = ichar(buf(1:1))
|
|
Chris@42
|
67 ibuf = 2
|
|
Chris@42
|
68 return
|
|
Chris@42
|
69 666 ic = -1
|
|
Chris@42
|
70 ibuf = 257
|
|
Chris@42
|
71 123 format(a256)
|
|
Chris@42
|
72 end
|
|
Chris@42
|
73
|
|
Chris@42
|
74 subroutine import_wisdom_from_file(isuccess, iunit)
|
|
Chris@42
|
75 integer isuccess
|
|
Chris@42
|
76 integer iunit
|
|
Chris@42
|
77 external read_char
|
|
Chris@42
|
78 call dfftw_import_wisdom(isuccess, read_char, iunit)
|
|
Chris@42
|
79 end
|
