Three-dimensional fourier transform processing method for shared memory scalar parallel computer

A computer readable recording medium storing a program for controlling a shared memory scalar parallel computer to realize a three-dimensional Fourier transform for use with the computer, which stores, in a continuous area, data in a direction in which i of input data x (i, j, k) is continuously arranged, where an index of a first dimension is i, an index of a second dimension is j, and an index of a third dimension is k in high-speed access to the continuous data area, according to operations comprising:

• loading into a cache each unit of a two-dimensional array data x (i, j, k*) having the indexes i and j and having a value ≧0 as the index k*;
• performing a Fourier transform of the first dimension in a direction in which the index i changes for a multiplicity in the j direction with the i continuously changed based upon a j-radix;
• performing a Fourier transform of the second dimension in a direction in which the index j changes for a multiplicity in the i direction while prioritizing an arithmetic operation in a direction in which the i continuously changes based upon an i-radix;
• performing a Fourier transform of the third dimension in a direction in which the index k changes based upon a k-radix;
• and
• outputting a calculated analysis according to the loading of the data into the cache and the Fourier transforms of the first, second, and/or third data dimensions loaded in the cache.

A shared memory scalar parallel computing apparatus executing a three-dimensional Fourier transform by storing in a continuous area, data in a direction in which i of input data x (i, j, k) is continuously arranged, where an index of a first dimension is i, an index of a second dimension is j, and an index of a third dimension is k in high-speed access to the continuous data area, the apparatus comprising:

• a loading unit loading into each unit of a two-dimensional array data x (i, j, k*) having the indexes i and j and having a value ≧0 as the index k*;
• a first dimension transform unit performing a Fourier transform of the first dimension in a direction in which the index i changes for a multiplicity in the j direction with the i continuously changed based upon the j-radix;
• a second dimension transform unit performing a Fourier transform of the second dimension in a direction in which the index j changes for a multiplicity in the i direction while prioritizing an arithmetic operation in a direction in which the i continuously changes based upon an i-radix;
• a third dimension transform unit performing a Fourier transform of the third dimension in a direction in which the index k changes based upon a k-radix;
• and
• an output unit outputting a calculated analysis according to the loading of the data into the cache and the Fourier transforms of the first, second, and/or third data dimensions loaded in the cache.