4 * This source code is part of
8 * GROningen MAchine for Chemical Simulations
11 * Copyright (c) 1991-2001, University of Groningen, The Netherlands
12 * This program is free software; you can redistribute it and/or
13 * modify it under the terms of the GNU General Public License
14 * as published by the Free Software Foundation; either version 2
15 * of the License, or (at your option) any later version.
17 * If you want to redistribute modifications, please consider that
18 * scientific software is very special. Version control is crucial -
19 * bugs must be traceable. We will be happy to consider code for
20 * inclusion in the official distribution, but derived work must not
21 * be called official GROMACS. Details are found in the README & COPYING
22 * files - if they are missing, get the official version at www.gromacs.org.
24 * To help us fund GROMACS development, we humbly ask that you cite
25 * the papers on the package - you can find them in the top README file.
27 * For more info, check our website at http://www.gromacs.org
30 * GROup of MAchos and Cynical Suckers
47 int main(int argc
,char *argv
[])
49 int mmm
[] = { 8, 10, 12, 15, 16, 18, 20, 24, 25, 27, 30, 32, 36, 40,
50 45, 48, 50, 54, 60, 64, 72, 75, 80, 81, 90, 100 };
51 int nnn
[] = { 24, 32, 48, 60, 72, 84, 96 };
52 #define NNN asize(nnn)
55 int i
,j
,n
,nit
,ntot
,n3
,rsize
;
60 static bool bOptFFT
= FALSE
;
62 static int nitfac
= 1;
64 { "-opt", FALSE
, etBOOL
, {&bOptFFT
},
66 { "-np", FALSE
, etINT
, {&nnode
},
68 { "-itfac", FALSE
, etINT
, {&nitfac
},
69 "Multiply number of iterations by this" }
71 static t_filenm fnm
[] = {
72 { efLOG
, "-g", "fft", ffWRITE
},
73 { efXVG
, "-o", "fft", ffWRITE
}
75 #define NFILE asize(fnm)
77 cr
= init_par(&argc
,&argv
);
79 CopyRight(stdout
,argv
[0]);
80 parse_common_args(&argc
,argv
,
81 PCA_CAN_SET_DEFFNM
| (MASTER(cr
) ? 0 : PCA_QUIET
),
82 TRUE
,NFILE
,fnm
,asize(pa
),pa
,0,NULL
,0,NULL
);
83 open_log(ftp2fn(efLOG
,NFILE
,fnm
),cr
);
89 for(i
=0; (i
<NNN
); i
++) {
103 fprintf(stderr
,"\r3D FFT (%s precision) %3d^3, niter %3d ",
104 (rsize
== 8) ? "Double" : "Single",n
,nit
);
106 g
= mk_fftgrid(stdlog
,(nnode
> 1),n
,n
,n
,bOptFFT
);
112 for(j
=0; (j
<nit
); j
++) {
113 gmxfft3D(g
,FFTW_FORWARD
,cr
);
114 gmxfft3D(g
,FFTW_BACKWARD
,cr
);
117 rt
[i
] = time(NULL
)-start
;
126 fprintf(stderr
,"\n");
127 fp
=xvgropen(ftp2fn(efXVG
,NFILE
,fnm
),
128 "FFT timings","n^3","t (s)");
129 for(i
=0; (i
<NNN
); i
++) {
130 n3
= 2*niter
[i
]*nnn
[i
]*nnn
[i
]*nnn
[i
];
131 fprintf(fp
,"%10d %10g\n",nnn
[i
],rt
[i
]/(2*niter
[i
]));