484 static MagickBooleanType ForwardFourier(const FourierInfo *fourier_info,  in ForwardFourier()  argument
527   magnitude_info=AcquireVirtualMemory((size_t) fourier_info->width,  in ForwardFourier()
528     fourier_info->height*sizeof(*magnitude_pixels));  in ForwardFourier()
529   phase_info=AcquireVirtualMemory((size_t) fourier_info->width,  in ForwardFourier()
530     fourier_info->height*sizeof(*phase_pixels));  in ForwardFourier()
543   (void) memset(magnitude_pixels,0,fourier_info->width*  in ForwardFourier()
544     fourier_info->height*sizeof(*magnitude_pixels));  in ForwardFourier()
546   (void) memset(phase_pixels,0,fourier_info->width*  in ForwardFourier()
547     fourier_info->height*sizeof(*phase_pixels));  in ForwardFourier()
548   status=ForwardQuadrantSwap(fourier_info->width,fourier_info->height,  in ForwardFourier()
551     status=ForwardQuadrantSwap(fourier_info->width,fourier_info->height,phase,  in ForwardFourier()
553   CorrectPhaseLHS(fourier_info->width,fourier_info->height,phase_pixels);  in ForwardFourier()
554   if (fourier_info->modulus != MagickFalse)  in ForwardFourier()
557       for (y=0L; y < (ssize_t) fourier_info->height; y++)  in ForwardFourier()
558         for (x=0L; x < (ssize_t) fourier_info->width; x++)  in ForwardFourier()
567   for (y=0L; y < (ssize_t) fourier_info->height; y++)  in ForwardFourier()
569     q=GetCacheViewAuthenticPixels(magnitude_view,0L,y,fourier_info->width,1UL,  in ForwardFourier()
573     for (x=0L; x < (ssize_t) fourier_info->width; x++)  in ForwardFourier()
575       switch (fourier_info->channel)  in ForwardFourier()
619   for (y=0L; y < (ssize_t) fourier_info->height; y++)  in ForwardFourier()
621     q=GetCacheViewAuthenticPixels(phase_view,0L,y,fourier_info->width,1UL,  in ForwardFourier()
625     for (x=0L; x < (ssize_t) fourier_info->width; x++)  in ForwardFourier()
627       switch (fourier_info->channel)  in ForwardFourier()
674 static MagickBooleanType ForwardFourierTransform(FourierInfo *fourier_info,  in ForwardFourierTransform()  argument
710   source_info=AcquireVirtualMemory((size_t) fourier_info->width,  in ForwardFourierTransform()
711     fourier_info->height*sizeof(*source_pixels));  in ForwardFourierTransform()
719   memset(source_pixels,0,fourier_info->width*fourier_info->height*  in ForwardFourierTransform()
723   for (y=0L; y < (ssize_t) fourier_info->height; y++)  in ForwardFourierTransform()
725     p=GetCacheViewVirtualPixels(image_view,0L,y,fourier_info->width,1UL,  in ForwardFourierTransform()
729     for (x=0L; x < (ssize_t) fourier_info->width; x++)  in ForwardFourierTransform()
731       switch (fourier_info->channel)  in ForwardFourierTransform()
765   forward_info=AcquireVirtualMemory((size_t) fourier_info->width,  in ForwardFourierTransform()
766     (fourier_info->height/2+1)*sizeof(*forward_pixels));  in ForwardFourierTransform()
778   fftw_r2c_plan=fftw_plan_dft_r2c_2d(fourier_info->width,fourier_info->height,  in ForwardFourierTransform()
793       gamma=PerceptibleReciprocal((double) fourier_info->width*  in ForwardFourierTransform()
794         fourier_info->height);  in ForwardFourierTransform()
795       for (y=0L; y < (ssize_t) fourier_info->height; y++)  in ForwardFourierTransform()
796         for (x=0L; x < (ssize_t) fourier_info->center; x++)  in ForwardFourierTransform()
811   if (fourier_info->modulus != MagickFalse)  in ForwardFourierTransform()
812     for (y=0L; y < (ssize_t) fourier_info->height; y++)  in ForwardFourierTransform()
813       for (x=0L; x < (ssize_t) fourier_info->center; x++)  in ForwardFourierTransform()
820     for (y=0L; y < (ssize_t) fourier_info->height; y++)  in ForwardFourierTransform()
821       for (x=0L; x < (ssize_t) fourier_info->center; x++)  in ForwardFourierTransform()
840     fourier_info;  in ForwardFourierTransformChannel()  local
849   fourier_info.width=image->columns;  in ForwardFourierTransformChannel()
850   fourier_info.height=image->rows;  in ForwardFourierTransformChannel()
855       fourier_info.width=(extent & 0x01) == 1 ? extent+1UL : extent;  in ForwardFourierTransformChannel()
857   fourier_info.height=fourier_info.width;  in ForwardFourierTransformChannel()
858   fourier_info.center=(ssize_t) (fourier_info.width/2L)+1L;  in ForwardFourierTransformChannel()
859   fourier_info.channel=channel;  in ForwardFourierTransformChannel()
860   fourier_info.modulus=modulus;  in ForwardFourierTransformChannel()
861   magnitude_info=AcquireVirtualMemory((size_t) fourier_info.width,  in ForwardFourierTransformChannel()
862     (fourier_info.height/2+1)*sizeof(*magnitude_pixels));  in ForwardFourierTransformChannel()
863   phase_info=AcquireVirtualMemory((size_t) fourier_info.width,  in ForwardFourierTransformChannel()
864     (fourier_info.height/2+1)*sizeof(*phase_pixels));  in ForwardFourierTransformChannel()
878   status=ForwardFourierTransform(&fourier_info,image,magnitude_pixels,  in ForwardFourierTransformChannel()
881     status=ForwardFourier(&fourier_info,fourier_image,magnitude_pixels,  in ForwardFourierTransformChannel()
1089 static MagickBooleanType InverseFourier(FourierInfo *fourier_info,  in InverseFourier()  argument
1123   magnitude_info=AcquireVirtualMemory((size_t) fourier_info->width,  in InverseFourier()
1124     fourier_info->height*sizeof(*magnitude_pixels));  in InverseFourier()
1125   phase_info=AcquireVirtualMemory((size_t) fourier_info->width,  in InverseFourier()
1126     fourier_info->height*sizeof(*phase_pixels));  in InverseFourier()
1127   inverse_info=AcquireVirtualMemory((size_t) fourier_info->width,  in InverseFourier()
1128     (fourier_info->height/2+1)*sizeof(*inverse_pixels));  in InverseFourier()
1149   for (y=0L; y < (ssize_t) fourier_info->height; y++)  in InverseFourier()
1151     p=GetCacheViewVirtualPixels(magnitude_view,0L,y,fourier_info->width,1UL,  in InverseFourier()
1155     for (x=0L; x < (ssize_t) fourier_info->width; x++)  in InverseFourier()
1157       switch (fourier_info->channel)  in InverseFourier()
1191   status=InverseQuadrantSwap(fourier_info->width,fourier_info->height,  in InverseFourier()
1193   (void) memcpy(magnitude_pixels,inverse_pixels,fourier_info->height*  in InverseFourier()
1194     fourier_info->center*sizeof(*magnitude_pixels));  in InverseFourier()
1197   for (y=0L; y < (ssize_t) fourier_info->height; y++)  in InverseFourier()
1199     p=GetCacheViewVirtualPixels(phase_view,0,y,fourier_info->width,1,  in InverseFourier()
1203     for (x=0L; x < (ssize_t) fourier_info->width; x++)  in InverseFourier()
1205       switch (fourier_info->channel)  in InverseFourier()
1238   if (fourier_info->modulus != MagickFalse)  in InverseFourier()
1241       for (y=0L; y < (ssize_t) fourier_info->height; y++)  in InverseFourier()
1242         for (x=0L; x < (ssize_t) fourier_info->width; x++)  in InverseFourier()
1250   CorrectPhaseLHS(fourier_info->width,fourier_info->height,phase_pixels);  in InverseFourier()
1252     status=InverseQuadrantSwap(fourier_info->width,fourier_info->height,  in InverseFourier()
1254   (void) memcpy(phase_pixels,inverse_pixels,fourier_info->height*  in InverseFourier()
1255     fourier_info->center*sizeof(*phase_pixels));  in InverseFourier()
1261   if (fourier_info->modulus != MagickFalse)  in InverseFourier()
1262     for (y=0L; y < (ssize_t) fourier_info->height; y++)  in InverseFourier()
1263        for (x=0L; x < (ssize_t) fourier_info->center; x++)  in InverseFourier()
1275     for (y=0L; y < (ssize_t) fourier_info->height; y++)  in InverseFourier()
1276       for (x=0L; x < (ssize_t) fourier_info->center; x++)  in InverseFourier()
1291 static MagickBooleanType InverseFourierTransform(FourierInfo *fourier_info,  in InverseFourierTransform()  argument
1319   source_info=AcquireVirtualMemory((size_t) fourier_info->width,  in InverseFourierTransform()
1320     fourier_info->height*sizeof(*source_pixels));  in InverseFourierTransform()
1338       gamma=PerceptibleReciprocal((double) fourier_info->width*  in InverseFourierTransform()
1339         fourier_info->height);  in InverseFourierTransform()
1340       for (y=0L; y < (ssize_t) fourier_info->height; y++)  in InverseFourierTransform()
1341         for (x=0L; x < (ssize_t) fourier_info->center; x++)  in InverseFourierTransform()
1355   fftw_c2r_plan=fftw_plan_dft_c2r_2d(fourier_info->width,fourier_info->height,  in InverseFourierTransform()
1361   for (y=0L; y < (ssize_t) fourier_info->height; y++)  in InverseFourierTransform()
1365     q=GetCacheViewAuthenticPixels(image_view,0L,y,fourier_info->width >  in InverseFourierTransform()
1366       image->columns ? image->columns : fourier_info->width,1UL,exception);  in InverseFourierTransform()
1369     for (x=0L; x < (ssize_t) fourier_info->width; x++)  in InverseFourierTransform()
1372         switch (fourier_info->channel)  in InverseFourierTransform()
1425     fourier_info;  in InverseFourierTransformChannel()  local
1433   fourier_info.width=magnitude_image->columns;  in InverseFourierTransformChannel()
1434   fourier_info.height=magnitude_image->rows;  in InverseFourierTransformChannel()
1441       fourier_info.width=(extent & 0x01) == 1 ? extent+1UL : extent;  in InverseFourierTransformChannel()
1443   fourier_info.height=fourier_info.width;  in InverseFourierTransformChannel()
1444   fourier_info.center=(ssize_t) (fourier_info.width/2L)+1L;  in InverseFourierTransformChannel()
1445   fourier_info.channel=channel;  in InverseFourierTransformChannel()
1446   fourier_info.modulus=modulus;  in InverseFourierTransformChannel()
1447   inverse_info=AcquireVirtualMemory((size_t) fourier_info.width,  in InverseFourierTransformChannel()
1448     (fourier_info.height/2+1)*sizeof(*inverse_pixels));  in InverseFourierTransformChannel()
1457   status=InverseFourier(&fourier_info,magnitude_image,phase_image,  in InverseFourierTransformChannel()
1460     status=InverseFourierTransform(&fourier_info,inverse_pixels,fourier_image,  in InverseFourierTransformChannel()