18#include <grass/config.h>
20#if defined(HAVE_FFTW_H) || defined(HAVE_DFFTW_H) || defined(HAVE_FFTW3_H)
23#include <grass/gmath.h>
25#include <grass/glocale.h>
43int del2g(
double *img[2],
int size,
double w)
48 fft(FORWARD, img, size * size, size, size);
50 g[0] = (
double *)G_malloc(size * size *
sizeof(
double));
51 g[1] = (
double *)G_malloc(size * size *
sizeof(
double));
56 G_message(_(
" taking FFT of del**2 g..."));
57 fft(FORWARD,
g, size * size, size, size);
60 G_message(_(
" multiplying transforms..."));
64 fft(INVERSE, img, size * size, size, size);
int getg(double w, double *g[2], int size)
void G_message(const char *msg,...)
Print a message to stderr.
int G_math_complex_mult(double *v1[2], int size1, double *v2[2], int size2, double *v3[2], int size3)