inline void StopWatch::Stop()
{
if (!Running) return;
#if (MFEM_TIMER_TYPE == 0)
user_time += ( std::clock() - start_utime );
#elif (MFEM_TIMER_TYPE == 1)
clock_t curr_rtime, curr_utime, curr_stime;
Current(&curr_rtime, &curr_utime, &curr_stime);
real_time += ( curr_rtime - start_rtime );
user_time += ( curr_utime - start_utime );
syst_time += ( curr_stime - start_stime );
#elif (MFEM_TIMER_TYPE == 2)
struct timespec curr_rtime, curr_utime;
GetRealTime(curr_rtime);
GetUserTime(curr_utime);
real_time.tv_sec += ( curr_rtime.tv_sec - start_rtime.tv_sec );
real_time.tv_nsec += ( curr_rtime.tv_nsec - start_rtime.tv_nsec );
user_time.tv_sec += ( curr_utime.tv_sec - start_utime.tv_sec );
user_time.tv_nsec += ( curr_utime.tv_nsec - start_utime.tv_nsec );
#elif (MFEM_TIMER_TYPE == 3)
LARGE_INTEGER curr_rtime;
QueryPerformanceCounter(&curr_rtime);
real_time.QuadPart += (curr_rtime.QuadPart - start_rtime.QuadPart);
#endif
Running = 0;
}
bool
Search::
out_of_time()
{
int const& timelimit = m_options.search.timeout;
if (timelimit > 0) {
double elapsed_time = GetUserTime();
if (elapsed_time > timelimit) {
VERBOSE(1,"Decoding is out of time (" << elapsed_time << ","
<< timelimit << ")" << std::endl);
interrupted_flag = 1;
return true;
}
}
int const& segment_timelimit = m_options.search.segment_timeout;
if (segment_timelimit > 0) {
double elapsed_time = m_timer.get_elapsed_time();
if (elapsed_time > segment_timelimit) {
VERBOSE(1,"Decoding for segment is out of time (" << elapsed_time << ","
<< segment_timelimit << ")" << std::endl);
interrupted_flag = 1;
return true;
}
}
return false;
}
inline double StopWatch::UserTime()
{
#if (MFEM_TIMER_TYPE == 0)
std::clock_t utime = user_time;
if (Running)
utime += (std::clock() - start_utime);
return (double)(utime) / CLOCKS_PER_SEC;
#elif (MFEM_TIMER_TYPE == 1)
clock_t curr_rtime, curr_utime, curr_stime;
clock_t utime = user_time;
if (Running)
{
Current(&curr_rtime, &curr_utime, &curr_stime);
utime += (curr_utime - start_utime);
}
return (double)(utime) / my_CLK_TCK;
#elif (MFEM_TIMER_TYPE == 2)
if (Running)
{
struct timespec curr_utime;
GetUserTime(curr_utime);
return ((user_time.tv_sec + (curr_utime.tv_sec - start_utime.tv_sec)) +
1e-9*(user_time.tv_nsec +
(curr_utime.tv_nsec - start_utime.tv_nsec)));
}
else
{
return user_time.tv_sec + 1e-9*user_time.tv_nsec;
}
#elif (MFEM_TIMER_TYPE == 3)
return RealTime();
#endif
}
bool
Search::
out_of_time()
{
int const& timelimit = m_options.search.timeout;
if (!timelimit) return false;
double elapsed_time = GetUserTime();
if (elapsed_time <= timelimit) return false;
VERBOSE(1,"Decoding is out of time (" << elapsed_time << ","
<< timelimit << ")" << std::endl);
interrupted_flag = 1;
return true;
}
inline void StopWatch::Start()
{
if (Running) return;
#if (MFEM_TIMER_TYPE == 0)
start_utime = std::clock();
#elif (MFEM_TIMER_TYPE == 1)
Current(&start_rtime, &start_utime, &start_stime);
#elif (MFEM_TIMER_TYPE == 2)
GetRealTime(start_rtime);
GetUserTime(start_utime);
#elif (MFEM_TIMER_TYPE == 3)
QueryPerformanceCounter(&start_rtime);
#endif
Running = 1;
}
ULONGLONG CThreadInfo::GetTotalTime()
{
ULONGLONG nTime = 0;
if ( IsValid() )
{
ULARGE_INTEGER ulKernel;
ULARGE_INTEGER ulUser;
memcpy( &ulKernel, &(GetKernelTime()), sizeof(FILETIME) );
memcpy( &ulUser, &(GetUserTime() ), sizeof(FILETIME) );
nTime = ulKernel.QuadPart + ulUser.QuadPart;
}
return nTime;
}
inline void StopWatch::Clear()
{
#if (MFEM_TIMER_TYPE == 0)
user_time = 0;
if (Running)
start_utime = std::clock();
#elif (MFEM_TIMER_TYPE == 1)
real_time = user_time = syst_time = 0;
if (Running)
Current(&start_rtime, &start_utime, &start_stime);
#elif (MFEM_TIMER_TYPE == 2)
real_time.tv_sec = user_time.tv_sec = 0;
real_time.tv_nsec = user_time.tv_nsec = 0;
if (Running)
{
GetRealTime(start_rtime);
GetUserTime(start_utime);
}
#elif (MFEM_TIMER_TYPE == 3)
real_time.QuadPart = 0;
if (Running)
QueryPerformanceCounter(&start_rtime);
#endif
}
double MIDASTimer::GetElapsedTime() {
return (GetUserTime() - start_time);
}
void MIDASTimer::StartTiming() {start_time = GetUserTime();}
MagickExport MagickBooleanType IdentifyImage(Image *image,FILE *file,
const MagickBooleanType verbose,ExceptionInfo *exception)
{
char
color[MaxTextExtent],
format[MaxTextExtent],
key[MaxTextExtent];
ChannelFeatures
*channel_features;
ChannelStatistics
*channel_statistics;
ColorspaceType
colorspace;
const char
*artifact,
*name,
*property,
*registry,
*value;
const MagickInfo
*magick_info;
double
elapsed_time,
user_time;
ImageType
type;
MagickBooleanType
ping;
register const Quantum
*p;
register ssize_t
i,
x;
size_t
distance,
scale;
ssize_t
y;
assert(image != (Image *) NULL);
assert(image->signature == MagickSignature);
if (image->debug != MagickFalse)
(void) LogMagickEvent(TraceEvent,GetMagickModule(),"%s",image->filename);
if (file == (FILE *) NULL)
file=stdout;
*format='\0';
elapsed_time=GetElapsedTime(&image->timer);
user_time=GetUserTime(&image->timer);
GetTimerInfo(&image->timer);
if (verbose == MagickFalse)
{
/*
Display summary info about the image.
*/
if (*image->magick_filename != '\0')
if (LocaleCompare(image->magick_filename,image->filename) != 0)
(void) FormatLocaleFile(file,"%s=>",image->magick_filename);
if ((GetPreviousImageInList(image) == (Image *) NULL) &&
(GetNextImageInList(image) == (Image *) NULL) &&
(image->scene == 0))
(void) FormatLocaleFile(file,"%s ",image->filename);
else
(void) FormatLocaleFile(file,"%s[%.20g] ",image->filename,(double)
image->scene);
(void) FormatLocaleFile(file,"%s ",image->magick);
if ((image->magick_columns != 0) || (image->magick_rows != 0))
if ((image->magick_columns != image->columns) ||
(image->magick_rows != image->rows))
(void) FormatLocaleFile(file,"%.20gx%.20g=>",(double)
image->magick_columns,(double) image->magick_rows);
(void) FormatLocaleFile(file,"%.20gx%.20g ",(double) image->columns,
(double) image->rows);
if ((image->page.width != 0) || (image->page.height != 0) ||
(image->page.x != 0) || (image->page.y != 0))
(void) FormatLocaleFile(file,"%.20gx%.20g%+.20g%+.20g ",(double)
image->page.width,(double) image->page.height,(double) image->page.x,
(double) image->page.y);
(void) FormatLocaleFile(file,"%.20g-bit ",(double) image->depth);
if (image->type != UndefinedType)
(void) FormatLocaleFile(file,"%s ",CommandOptionToMnemonic(
MagickTypeOptions,(ssize_t) image->type));
if (image->storage_class == DirectClass)
{
(void) FormatLocaleFile(file,"DirectClass ");
if (image->total_colors != 0)
{
(void) FormatMagickSize(image->total_colors,MagickFalse,format);
(void) FormatLocaleFile(file,"%s ",format);
}
}
else
if (image->total_colors <= image->colors)
(void) FormatLocaleFile(file,"PseudoClass %.20gc ",(double)
image->colors);
else
(void) FormatLocaleFile(file,"PseudoClass %.20g=>%.20gc ",(double)
image->total_colors,(double) image->colors);
if (image->error.mean_error_per_pixel != 0.0)
(void) FormatLocaleFile(file,"%.20g/%f/%fdb ",(double)
(image->error.mean_error_per_pixel+0.5),
image->error.normalized_mean_error,
image->error.normalized_maximum_error);
if (GetBlobSize(image) != 0)
{
(void) FormatMagickSize(GetBlobSize(image),MagickFalse,format);
(void) FormatLocaleFile(file,"%s ",format);
}
(void) FormatLocaleFile(file,"%0.3fu %lu:%02lu.%03lu",user_time,
(unsigned long) (elapsed_time/60.0),(unsigned long) floor(fmod(
elapsed_time,60.0)),(unsigned long) (1000.0*(elapsed_time-
floor(elapsed_time))));
(void) FormatLocaleFile(file,"\n");
(void) fflush(file);
return(ferror(file) != 0 ? MagickFalse : MagickTrue);
}
/*
Display verbose info about the image.
*/
p=GetVirtualPixels(image,0,0,1,1,exception);
ping=p == (const Quantum *) NULL ? MagickTrue : MagickFalse;
type=GetImageType(image,exception);
(void) SignatureImage(image,exception);
(void) FormatLocaleFile(file,"Image: %s\n",image->filename);
if (*image->magick_filename != '\0')
if (LocaleCompare(image->magick_filename,image->filename) != 0)
{
char
filename[MaxTextExtent];
GetPathComponent(image->magick_filename,TailPath,filename);
(void) FormatLocaleFile(file," Base filename: %s\n",filename);
}
magick_info=GetMagickInfo(image->magick,exception);
if ((magick_info == (const MagickInfo *) NULL) ||
(*GetMagickDescription(magick_info) == '\0'))
(void) FormatLocaleFile(file," Format: %s\n",image->magick);
else
(void) FormatLocaleFile(file," Format: %s (%s)\n",image->magick,
GetMagickDescription(magick_info));
(void) FormatLocaleFile(file," Class: %s\n",CommandOptionToMnemonic(
MagickClassOptions,(ssize_t) image->storage_class));
(void) FormatLocaleFile(file," Geometry: %.20gx%.20g%+.20g%+.20g\n",(double)
image->columns,(double) image->rows,(double) image->tile_offset.x,(double)
image->tile_offset.y);
if ((image->magick_columns != 0) || (image->magick_rows != 0))
if ((image->magick_columns != image->columns) ||
(image->magick_rows != image->rows))
(void) FormatLocaleFile(file," Base geometry: %.20gx%.20g\n",(double)
image->magick_columns,(double) image->magick_rows);
if ((image->resolution.x != 0.0) && (image->resolution.y != 0.0))
{
(void) FormatLocaleFile(file," Resolution: %gx%g\n",image->resolution.x,
image->resolution.y);
(void) FormatLocaleFile(file," Print size: %gx%g\n",(double)
image->columns/image->resolution.x,(double) image->rows/
image->resolution.y);
}
(void) FormatLocaleFile(file," Units: %s\n",CommandOptionToMnemonic(
MagickResolutionOptions,(ssize_t) image->units));
(void) FormatLocaleFile(file," Type: %s\n",CommandOptionToMnemonic(
MagickTypeOptions,(ssize_t) type));
if (image->type != UndefinedType)
(void) FormatLocaleFile(file," Base type: %s\n",CommandOptionToMnemonic(
MagickTypeOptions,(ssize_t) image->type));
(void) FormatLocaleFile(file," Endianess: %s\n",CommandOptionToMnemonic(
MagickEndianOptions,(ssize_t) image->endian));
/*
Detail channel depth and extrema.
*/
(void) FormatLocaleFile(file," Colorspace: %s\n",CommandOptionToMnemonic(
MagickColorspaceOptions,(ssize_t) image->colorspace));
channel_statistics=(ChannelStatistics *) NULL;
channel_features=(ChannelFeatures *) NULL;
colorspace=image->colorspace;
scale=1;
if (ping == MagickFalse)
{
size_t
depth;
channel_statistics=GetImageStatistics(image,exception);
artifact=GetImageArtifact(image,"identify:features");
if (artifact != (const char *) NULL)
{
distance=StringToUnsignedLong(artifact);
channel_features=GetImageFeatures(image,distance,exception);
}
depth=GetImageDepth(image,exception);
if (image->depth == depth)
(void) FormatLocaleFile(file," Depth: %.20g-bit\n",(double)
image->depth);
else
(void) FormatLocaleFile(file," Depth: %.20g/%.20g-bit\n",(double)
image->depth,(double) depth);
(void) FormatLocaleFile(file," Channel depth:\n");
if (IsImageGray(image,exception) != MagickFalse)
colorspace=GRAYColorspace;
switch (colorspace)
{
case RGBColorspace:
default:
{
(void) FormatLocaleFile(file," red: %.20g-bit\n",(double)
channel_statistics[RedPixelChannel].depth);
(void) FormatLocaleFile(file," green: %.20g-bit\n",(double)
channel_statistics[GreenPixelChannel].depth);
(void) FormatLocaleFile(file," blue: %.20g-bit\n",(double)
channel_statistics[BluePixelChannel].depth);
break;
}
case CMYKColorspace:
{
(void) FormatLocaleFile(file," cyan: %.20g-bit\n",(double)
channel_statistics[CyanPixelChannel].depth);
(void) FormatLocaleFile(file," magenta: %.20g-bit\n",(double)
channel_statistics[MagentaPixelChannel].depth);
(void) FormatLocaleFile(file," yellow: %.20g-bit\n",(double)
channel_statistics[YellowPixelChannel].depth);
(void) FormatLocaleFile(file," black: %.20g-bit\n",(double)
channel_statistics[BlackPixelChannel].depth);
break;
}
case GRAYColorspace:
{
(void) FormatLocaleFile(file," gray: %.20g-bit\n",(double)
channel_statistics[GrayPixelChannel].depth);
break;
}
}
if (image->matte != MagickFalse)
(void) FormatLocaleFile(file," alpha: %.20g-bit\n",(double)
channel_statistics[AlphaPixelChannel].depth);
scale=1;
if (image->depth <= MAGICKCORE_QUANTUM_DEPTH)
scale=QuantumRange/((size_t) QuantumRange >> ((size_t)
MAGICKCORE_QUANTUM_DEPTH-image->depth));
}
/*
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% %
% %
% %
% M a i n %
% %
% %
% %
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
%
*/
int main(int argc,char **argv)
{
char
*option;
double
elapsed_time,
user_time;
ExceptionInfo
*exception;
ImageInfo
*image_info;
MagickBooleanType
regard_warnings,
status;
register long
i;
TimerInfo
timer;
unsigned long
iterations;
MagickCoreGenesis(*argv,MagickTrue);
exception=AcquireExceptionInfo();
iterations=1;
status=MagickFalse;
regard_warnings=MagickFalse;
for (i=1; i < (long) (argc-1); i++)
{
option=argv[i];
if ((strlen(option) == 1) || ((*option != '-') && (*option != '+')))
continue;
if (LocaleCompare("bench",option+1) == 0)
iterations=(unsigned long) atol(argv[++i]);
if (LocaleCompare("debug",option+1) == 0)
(void) SetLogEventMask(argv[++i]);
if (LocaleCompare("regard-warnings",option+1) == 0)
regard_warnings=MagickTrue;
}
GetTimerInfo(&timer);
for (i=0; i < (long) iterations; i++)
{
image_info=AcquireImageInfo();
status=ConvertImageCommand(image_info,argc,argv,(char **) NULL,exception);
if (exception->severity != UndefinedException)
{
if ((exception->severity > ErrorException) ||
(regard_warnings != MagickFalse))
status=MagickTrue;
CatchException(exception);
}
image_info=DestroyImageInfo(image_info);
}
if (iterations > 1)
{
elapsed_time=GetElapsedTime(&timer);
user_time=GetUserTime(&timer);
(void) fprintf(stderr,"Performance: %lui %gips %0.3fu %ld:%02ld\n",
iterations,1.0*iterations/elapsed_time,user_time,(long)
(elapsed_time/60.0+0.5),(long) ceil(fmod(elapsed_time,60.0)));
}
exception=DestroyExceptionInfo(exception);
MagickCoreTerminus();
return(status == MagickFalse ? 0 : 1);
}
/*
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% %
% %
% %
% I d e n t i f y I m a g e %
% %
% %
% %
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
% IdentifyImage() identifies an image by printing its attributes to the file.
% Attributes include the image width, height, size, and others.
%
% The format of the IdentifyImage method is:
%
% MagickBooleanType IdentifyImage(Image *image,FILE *file,
% const MagickBooleanType verbose)
%
% A description of each parameter follows:
%
% o image: The image.
%
% o file: The file, typically stdout.
%
% o verbose: A value other than zero prints more detailed information
% about the image.
%
*/
MagickExport MagickBooleanType IdentifyImage(Image *image,FILE *file,
const MagickBooleanType verbose)
{
#define IdentifyFormat " %s:\n Min: " QuantumFormat \
" (%g)\n Max: " QuantumFormat " (%g)\n" \
" Mean: %g (%g)\n Standard deviation: %g (%g)\n"
char
color[MaxTextExtent],
format[MaxTextExtent],
key[MaxTextExtent];
ColorspaceType
colorspace;
const char
*property,
*value;
const MagickInfo
*magick_info;
const PixelPacket
*pixels;
double
elapsed_time,
user_time;
ExceptionInfo
*exception;
Image
*p;
ImageType
type;
long
y;
MagickBooleanType
ping;
register long
i,
x;
unsigned long
scale;
assert(image != (Image *) NULL);
assert(image->signature == MagickSignature);
if (image->debug != MagickFalse)
(void) LogMagickEvent(TraceEvent,GetMagickModule(),"%s",image->filename);
if (file == (FILE *) NULL)
file=stdout;
*format='\0';
elapsed_time=GetElapsedTime(&image->timer);
user_time=GetUserTime(&image->timer);
GetTimerInfo(&image->timer);
if (verbose == MagickFalse)
{
/*
Display summary info about the image.
*/
if (*image->magick_filename != '\0')
if (LocaleCompare(image->magick_filename,image->filename) != 0)
(void) fprintf(file,"%s=>",image->magick_filename);
if ((GetPreviousImageInList(image) == (Image *) NULL) &&
(GetNextImageInList(image) == (Image *) NULL) && (image->scene == 0))
(void) fprintf(file,"%s ",image->filename);
else
(void) fprintf(file,"%s[%lu] ",image->filename,image->scene);
(void) fprintf(file,"%s ",image->magick);
if ((image->magick_columns != 0) || (image->magick_rows != 0))
if ((image->magick_columns != image->columns) ||
(image->magick_rows != image->rows))
(void) fprintf(file,"%lux%lu=>",image->magick_columns,
image->magick_rows);
(void) fprintf(file,"%lux%lu ",image->columns,image->rows);
if ((image->page.width != 0) || (image->page.height != 0) ||
(image->page.x != 0) || (image->page.y != 0))
(void) fprintf(file,"%lux%lu%+ld%+ld ",image->page.width,
image->page.height,image->page.x,image->page.y);
if (image->storage_class == DirectClass)
{
(void) fprintf(file,"DirectClass ");
if (image->total_colors != 0)
{
(void) FormatMagickSize(image->total_colors,format);
(void) fprintf(file,"%s ",format);
}
}
else
if (image->total_colors <= image->colors)
(void) fprintf(file,"PseudoClass %luc ",image->colors);
else
(void) fprintf(file,"PseudoClass %lu=>%luc ",image->total_colors,
image->colors);
(void) fprintf(file,"%lu-bit ",image->depth);
if (image->error.mean_error_per_pixel != 0.0)
(void) fprintf(file,"%ld/%f/%fdb ",
(long) (image->error.mean_error_per_pixel+0.5),
image->error.normalized_mean_error,
image->error.normalized_maximum_error);
if (GetBlobSize(image) != 0)
{
(void) FormatMagickSize(GetBlobSize(image),format);
(void) fprintf(file,"%s ",format);
}
if (elapsed_time > 0.06)
(void) fprintf(file,"%0.3fu %ld:%02ld",user_time,
(long) (elapsed_time/60.0+0.5),(long) ceil(fmod(elapsed_time,60.0)));
(void) fprintf(file,"\n");
(void) fflush(file);
return(ferror(file) != 0 ? MagickFalse : MagickTrue);
}
/*
Display verbose info about the image.
*/
exception=AcquireExceptionInfo();
pixels=AcquireImagePixels(image,0,0,1,1,exception);
exception=DestroyExceptionInfo(exception);
ping=pixels == (const PixelPacket *) NULL ? MagickTrue : MagickFalse;
type=GetImageType(image,&image->exception);
(void) SignatureImage(image);
(void) fprintf(file,"Image: %s\n",image->filename);
if (*image->magick_filename != '\0')
if (LocaleCompare(image->magick_filename,image->filename) != 0)
{
char
filename[MaxTextExtent];
GetPathComponent(image->magick_filename,TailPath,filename);
(void) fprintf(file," Base filename: %s\n",filename);
}
magick_info=GetMagickInfo(image->magick,&image->exception);
if ((magick_info == (const MagickInfo *) NULL) ||
(*GetMagickDescription(magick_info) == '\0'))
(void) fprintf(file," Format: %s\n",image->magick);
else
(void) fprintf(file," Format: %s (%s)\n",image->magick,
GetMagickDescription(magick_info));
(void) fprintf(file," Class: %s\n",
MagickOptionToMnemonic(MagickClassOptions,(long) image->storage_class));
(void) fprintf(file," Geometry: %lux%lu%+ld%+ld\n",image->columns,
image->rows,image->tile_offset.x,image->tile_offset.y);
if ((image->magick_columns != 0) || (image->magick_rows != 0))
if ((image->magick_columns != image->columns) ||
(image->magick_rows != image->rows))
(void) fprintf(file," Base geometry: %lux%lu\n",image->magick_columns,
image->magick_rows);
(void) fprintf(file," Type: %s\n",MagickOptionToMnemonic(MagickTypeOptions,
(long) type));
(void) fprintf(file," Endianess: %s\n",MagickOptionToMnemonic(
MagickEndianOptions,(long) image->endian));
/*
Detail channel depth and extrema.
*/
colorspace=image->colorspace;
if (IsGrayImage(image,&image->exception) != MagickFalse)
colorspace=GRAYColorspace;
(void) fprintf(file," Colorspace: %s\n",
MagickOptionToMnemonic(MagickColorspaceOptions,(long) colorspace));
if (ping == MagickFalse)
{
ChannelStatistics
*channel_statistics;
(void) fprintf(file," Depth: %lu-bit\n",GetImageDepth(image,
&image->exception));
channel_statistics=GetImageChannelStatistics(image,&image->exception);
(void) fprintf(file," Channel depth:\n");
switch (colorspace)
{
case RGBColorspace:
default:
{
(void) fprintf(file," Red: %lu-bit\n",
channel_statistics[RedChannel].depth);
(void) fprintf(file," Green: %lu-bit\n",
channel_statistics[GreenChannel].depth);
(void) fprintf(file," Blue: %lu-bit\n",
channel_statistics[BlueChannel].depth);
if (image->matte != MagickFalse)
(void) fprintf(file," Alpha: %lu-bit\n",
channel_statistics[OpacityChannel].depth);
break;
}
case CMYKColorspace:
{
(void) fprintf(file," Cyan: %lu-bit\n",
channel_statistics[CyanChannel].depth);
(void) fprintf(file," Magenta: %lu-bit\n",
channel_statistics[MagentaChannel].depth);
(void) fprintf(file," Yellow: %lu-bit\n",
channel_statistics[YellowChannel].depth);
(void) fprintf(file," Black: %lu-bit\n",
channel_statistics[BlackChannel].depth);
if (image->matte != MagickFalse)
(void) fprintf(file," Alpha: %lu-bit\n",
channel_statistics[OpacityChannel].depth);
break;
}
case GRAYColorspace:
{
(void) fprintf(file," Gray: %lu-bit\n",
channel_statistics[GrayChannel].depth);
if (image->matte != MagickFalse)
(void) fprintf(file," Alpha: %lu-bit\n",
channel_statistics[OpacityChannel].depth);
break;
}
}
scale=QuantumRange/((unsigned long) QuantumRange >> ((unsigned long)
MAGICKCORE_QUANTUM_DEPTH-channel_statistics[AllChannels].depth));
(void) fprintf(file," Channel statistics:\n");
switch (colorspace)
{
case RGBColorspace:
default:
{
(void) fprintf(file,IdentifyFormat,"Red",(Quantum)
(channel_statistics[RedChannel].minima/scale),(double)
channel_statistics[RedChannel].minima/(double) QuantumRange,
(Quantum) (channel_statistics[RedChannel].maxima/scale),(double)
channel_statistics[RedChannel].maxima/(double) QuantumRange,
channel_statistics[RedChannel].mean/(double) scale,
channel_statistics[RedChannel].mean/(double) QuantumRange,
channel_statistics[RedChannel].standard_deviation/(double) scale,
channel_statistics[RedChannel].standard_deviation/(double)
QuantumRange);
(void) fprintf(file,IdentifyFormat,"Green",(Quantum)
(channel_statistics[GreenChannel].minima/scale),(double)
channel_statistics[GreenChannel].minima/(double) QuantumRange,
(Quantum) (channel_statistics[GreenChannel].maxima/scale),(double)
channel_statistics[GreenChannel].maxima/(double) QuantumRange,
channel_statistics[GreenChannel].mean/(double) scale,
channel_statistics[GreenChannel].mean/(double) QuantumRange,
channel_statistics[GreenChannel].standard_deviation/(double) scale,
channel_statistics[GreenChannel].standard_deviation/(double)
QuantumRange);
(void) fprintf(file,IdentifyFormat,"Blue",(Quantum)
(channel_statistics[BlueChannel].minima/scale),(double)
channel_statistics[BlueChannel].minima/(double) QuantumRange,
(Quantum) (channel_statistics[BlueChannel].maxima/scale),(double)
channel_statistics[BlueChannel].maxima/(double) QuantumRange,
channel_statistics[BlueChannel].mean/(double) scale,
channel_statistics[BlueChannel].mean/(double) QuantumRange,
channel_statistics[BlueChannel].standard_deviation/(double) scale,
channel_statistics[BlueChannel].standard_deviation/(double)
QuantumRange);
break;
}
case CMYKColorspace:
{
(void) fprintf(file,IdentifyFormat,"Cyan",(Quantum)
(channel_statistics[CyanChannel].minima/scale),(double)
channel_statistics[CyanChannel].minima/(double) QuantumRange,
(Quantum) (channel_statistics[CyanChannel].maxima/scale),(double)
channel_statistics[CyanChannel].maxima/(double) QuantumRange,
channel_statistics[CyanChannel].mean/(double) scale,
channel_statistics[CyanChannel].mean/(double) QuantumRange,
channel_statistics[CyanChannel].standard_deviation/(double) scale,
channel_statistics[CyanChannel].standard_deviation/(double)
QuantumRange);
(void) fprintf(file,IdentifyFormat,"Magenta",(Quantum)
(channel_statistics[MagentaChannel].minima/scale),(double)
channel_statistics[MagentaChannel].minima/(double) QuantumRange,
(Quantum) (channel_statistics[MagentaChannel].maxima/scale),(double)
channel_statistics[MagentaChannel].maxima/(double) QuantumRange,
channel_statistics[MagentaChannel].mean/(double) scale,
channel_statistics[MagentaChannel].mean/(double) QuantumRange,
channel_statistics[MagentaChannel].standard_deviation/(double)
scale,channel_statistics[MagentaChannel].standard_deviation/(double)
QuantumRange);
(void) fprintf(file,IdentifyFormat,"Yellow",(Quantum)
(channel_statistics[YellowChannel].minima/scale),(double)
channel_statistics[YellowChannel].minima/(double) QuantumRange,
(Quantum) (channel_statistics[YellowChannel].maxima/scale),(double)
channel_statistics[YellowChannel].maxima/(double) QuantumRange,
channel_statistics[YellowChannel].mean/(double) scale,
channel_statistics[YellowChannel].mean/(double) QuantumRange,
channel_statistics[YellowChannel].standard_deviation/(double) scale,
channel_statistics[YellowChannel].standard_deviation/(double)
QuantumRange);
(void) fprintf(file,IdentifyFormat,"Black",(Quantum)
(channel_statistics[BlackChannel].minima/scale),(double)
channel_statistics[BlackChannel].minima/(double) QuantumRange,
(Quantum) (channel_statistics[BlackChannel].maxima/scale),(double)
channel_statistics[BlackChannel].maxima/(double) QuantumRange,
channel_statistics[BlackChannel].mean/(double) scale,
channel_statistics[BlackChannel].mean/(double) QuantumRange,
channel_statistics[BlackChannel].standard_deviation/(double) scale,
channel_statistics[BlackChannel].standard_deviation/(double)
QuantumRange);
break;
}
case GRAYColorspace:
{
(void) fprintf(file,IdentifyFormat,"Gray",(Quantum)
(channel_statistics[GrayChannel].minima/scale),(double)
channel_statistics[GrayChannel].minima/(double) QuantumRange,
(Quantum) (channel_statistics[GrayChannel].maxima/scale),(double)
channel_statistics[GrayChannel].maxima/(double) QuantumRange,
channel_statistics[GrayChannel].mean/(double) scale,
channel_statistics[GrayChannel].mean/(double) QuantumRange,
channel_statistics[GrayChannel].standard_deviation/(double) scale,
channel_statistics[GrayChannel].standard_deviation/(double)
QuantumRange);
break;
}
}
if (image->matte != MagickFalse)
(void) fprintf(file,IdentifyFormat,"Opacity",(Quantum)
(channel_statistics[OpacityChannel].minima/scale),(double)
channel_statistics[OpacityChannel].minima/(double) QuantumRange,
(Quantum) (channel_statistics[OpacityChannel].maxima/scale),(double)
channel_statistics[OpacityChannel].maxima/(double) QuantumRange,
channel_statistics[OpacityChannel].mean/(double) scale,
channel_statistics[OpacityChannel].mean/(double) QuantumRange,
channel_statistics[OpacityChannel].standard_deviation/(double) scale,
channel_statistics[OpacityChannel].standard_deviation/(double)
QuantumRange);
channel_statistics=(ChannelStatistics *) RelinquishMagickMemory(
channel_statistics);
if (colorspace == CMYKColorspace)
(void) fprintf(file," Total ink density: %.0f%%\n",100.0*
GetImageTotalInkDensity(image)/(double) QuantumRange);
x=0;
p=NewImageList();
if (image->matte != MagickFalse)
{
register const IndexPacket
*indexes;
register const PixelPacket
*p;
p=(PixelPacket *) NULL;
indexes=(IndexPacket *) NULL;
for (y=0; y < (long) image->rows; y++)
{
p=AcquireImagePixels(image,0,y,image->columns,1,&image->exception);
if (p == (const PixelPacket *) NULL)
break;
indexes=AcquireIndexes(image);
for (x=0; x < (long) image->columns; x++)
{
if (p->opacity == (Quantum) TransparentOpacity)
break;
p++;
}
if (x < (long) image->columns)
break;
}
if ((x < (long) image->columns) || (y < (long) image->rows))
{
char
tuple[MaxTextExtent];
MagickPixelPacket
pixel;
GetMagickPixelPacket(image,&pixel);
SetMagickPixelPacket(image,p,indexes+x,&pixel);
(void) QueryMagickColorname(image,&pixel,SVGCompliance,
MagickFalse,tuple,&image->exception);
(void) fprintf(file," Alpha: %s ",tuple);
(void) QueryMagickColorname(image,&pixel,SVGCompliance,MagickTrue,
tuple,&image->exception);
(void) fprintf(file," %s\n",tuple);
}
}
if ((ping == MagickFalse) &&
(IsHistogramImage(image,&image->exception) != MagickFalse))
{
(void) fprintf(file," Histogram:\n");
(void) GetNumberColors(image,file,&image->exception);
}
}
int main(int argc,char **argv)
{
#define DestroyValidate() \
{ \
timer=DestroyTimerInfo(timer); \
image_info=DestroyImageInfo(image_info); \
exception=DestroyExceptionInfo(exception); \
}
#define ThrowValidateException(asperity,tag,option) \
{ \
(void) ThrowMagickException(exception,GetMagickModule(),asperity,tag,"`%s'", \
option); \
CatchException(exception); \
DestroyValidate(); \
return(MagickFalse); \
}
char
output_filename[MaxTextExtent],
reference_filename[MaxTextExtent],
*option;
double
elapsed_time,
user_time;
ExceptionInfo
*exception;
Image
*reference_image;
ImageInfo
*image_info;
MagickBooleanType
regard_warnings,
status;
register long
i;
TimerInfo
*timer;
unsigned long
fail,
iterations,
tests;
ValidateType
type;
/*
Validate the ImageMagick image processing suite.
*/
MagickCoreGenesis(*argv,MagickFalse);
iterations=1;
status=MagickFalse;
type=AllValidate;
regard_warnings=MagickFalse;
exception=AcquireExceptionInfo();
image_info=AcquireImageInfo();
(void) CopyMagickString(image_info->filename,ReferenceFilename,MaxTextExtent);
timer=AcquireTimerInfo();
GetTimerInfo(timer);
for (i=1; i < (long) argc; i++)
{
option=argv[i];
if (IsMagickOption(option) == MagickFalse)
{
(void) CopyMagickString(image_info->filename,option,MaxTextExtent);
continue;
}
switch (*(option+1))
{
case 'b':
{
if (LocaleCompare("bench",option+1) == 0)
{
iterations=(unsigned long) atol(argv[++i]);
break;
}
ThrowValidateException(OptionError,"UnrecognizedOption",option)
}
case 'd':
{
if (LocaleCompare("debug",option+1) == 0)
{
(void) SetLogEventMask(argv[++i]);
break;
}
ThrowValidateException(OptionError,"UnrecognizedOption",option)
}
case 'h':
{
if (LocaleCompare("help",option+1) == 0)
{
(void) ValidateUsage();
return(0);
}
ThrowValidateException(OptionError,"UnrecognizedOption",option)
}
case 'l':
{
if (LocaleCompare("log",option+1) == 0)
{
if (*option != '+')
(void) SetLogFormat(argv[i+1]);
break;
}
ThrowValidateException(OptionError,"UnrecognizedOption",option)
}
case 'r':
{
if (LocaleCompare("regard-warnings",option+1) == 0)
{
regard_warnings=MagickTrue;
break;
}
ThrowValidateException(OptionError,"UnrecognizedOption",option)
}
case 'v':
{
if (LocaleCompare("validate",option+1) == 0)
{
long
validate;
if (*option == '+')
break;
i++;
if (i == (long) argc)
ThrowValidateException(OptionError,"MissingArgument",option);
validate=ParseMagickOption(MagickValidateOptions,MagickFalse,
argv[i]);
if (validate < 0)
ThrowValidateException(OptionError,"UnrecognizedValidateType",
argv[i]);
type=(ValidateType) validate;
break;
}
if ((LocaleCompare("version",option+1) == 0) ||
(LocaleCompare("-version",option+1) == 0))
{
(void) fprintf(stdout,"Version: %s\n",
GetMagickVersion((unsigned long *) NULL));
(void) fprintf(stdout,"Copyright: %s\n\n",GetMagickCopyright());
return(0);
}
ThrowValidateException(OptionError,"UnrecognizedOption",option)
}
default:
ThrowValidateException(OptionError,"UnrecognizedOption",option)
}
}
reference_image=ReadImage(image_info,exception);
tests=0;
fail=0;
if (reference_image == (Image *) NULL)
fail++;
else
{
if (LocaleCompare(image_info->filename,ReferenceFilename) == 0)
(void) CopyMagickString(reference_image->magick,ReferenceImageFormat,
MaxTextExtent);
(void) AcquireUniqueFilename(reference_filename);
(void) AcquireUniqueFilename(output_filename);
(void) CopyMagickString(reference_image->filename,reference_filename,
MaxTextExtent);
status=WriteImage(image_info,reference_image);
InheritException(exception,&reference_image->exception);
reference_image=DestroyImage(reference_image);
if (status == MagickFalse)
fail++;
else
{
(void) fprintf(stdout,"Version: %s\n",
GetMagickVersion((unsigned long *) NULL));
(void) fprintf(stdout,"Copyright: %s\n\n",
GetMagickCopyright());
(void) fprintf(stdout,"ImageMagick Validation Suite (%s)\n\n",
MagickOptionToMnemonic(MagickValidateOptions,(long) type));
if ((type & CompareValidate) != 0)
tests+=ValidateCompareCommand(image_info,reference_filename,
output_filename,&fail,exception);
if ((type & CompositeValidate) != 0)
tests+=ValidateCompositeCommand(image_info,reference_filename,
output_filename,&fail,exception);
if ((type & ConvertValidate) != 0)
tests+=ValidateConvertCommand(image_info,reference_filename,
output_filename,&fail,exception);
if ((type & FormatsInMemoryValidate) != 0)
tests+=ValidateImageFormatsInMemory(image_info,reference_filename,
output_filename,&fail,exception);
if ((type & FormatsOnDiskValidate) != 0)
tests+=ValidateImageFormatsOnDisk(image_info,reference_filename,
output_filename,&fail,exception);
if ((type & IdentifyValidate) != 0)
tests+=ValidateIdentifyCommand(image_info,reference_filename,
output_filename,&fail,exception);
if ((type & ImportExportValidate) != 0)
tests+=ValidateImportExportPixels(image_info,reference_filename,
output_filename,&fail,exception);
if ((type & MontageValidate) != 0)
tests+=ValidateMontageCommand(image_info,reference_filename,
output_filename,&fail,exception);
if ((type & StreamValidate) != 0)
tests+=ValidateStreamCommand(image_info,reference_filename,
output_filename,&fail,exception);
(void) fprintf(stdout,"validation suite: %lu tests; %lu passed; "
"%lu failed.\n",tests,tests-fail,fail);
}
(void) RelinquishUniqueFileResource(output_filename);
(void) RelinquishUniqueFileResource(reference_filename);
}
if (exception->severity != UndefinedException)
CatchException(exception);
if (iterations > 1)
{
elapsed_time=GetElapsedTime(timer);
user_time=GetUserTime(timer);
(void) fprintf(stderr,"Performance: %lui %gips %0.3fu %ld:%02ld\n",
iterations,1.0*iterations/elapsed_time,user_time,(long)
(elapsed_time/60.0+0.5),(long) ceil(fmod(elapsed_time,60.0)));
}
DestroyValidate();
MagickCoreTerminus();
return(fail == 0 ? 0 : 1);
}
static int WizardMain(int argc,char **argv)
{
char
*option;
double
elapsed_time,
user_time;
ExceptionInfo
*exception;
register ssize_t
i;
TimerInfo
*timer;
unsigned int
iterations;
WizardBooleanType
regard_warnings,
status;
WizardsToolkitGenesis(*argv);
exception=AcquireExceptionInfo();
iterations=1;
status=WizardTrue;
regard_warnings=WizardFalse;
for (i=1; i < (ssize_t) (argc-1); i++)
{
option=argv[i];
if ((strlen(option) == 1) || ((*option != '-') && (*option != '+')))
continue;
if (LocaleCompare("bench",option+1) == 0)
iterations=(unsigned int) StringToUnsignedLong(argv[++i]);
if (LocaleCompare("debug",option+1) == 0)
(void) SetLogEventMask(argv[++i]);
if (LocaleCompare("regard-warnings",option+1) == 0)
regard_warnings=WizardTrue;
}
timer=(TimerInfo *) NULL;
if (iterations > 1)
timer=AcquireTimerInfo();
for (i=0; i < (ssize_t) iterations; i++)
{
status=KeyringCommand(argc,argv,exception);
if ((status == WizardFalse) ||
(GetExceptionSeverity(exception) != UndefinedException))
{
if ((GetExceptionSeverity(exception) < ErrorException) &&
(regard_warnings == WizardFalse))
status=WizardTrue;
CatchException(exception);
}
}
if (iterations > 1)
{
elapsed_time=GetElapsedTime(timer);
user_time=GetUserTime(timer);
(void) fprintf(stderr,"Performance: %ui %gips %0.3fu %.20g:%02g.%03g\n",
iterations,1.0*iterations/elapsed_time,user_time,(double)
(elapsed_time/60.0+0.5),floor(fmod(elapsed_time,60.0)),(double)
(1000.0*(elapsed_time-floor(elapsed_time))+0.5));
timer=DestroyTimerInfo(timer);
}
exception=DestroyExceptionInfo(exception);
WizardsToolkitTerminus();
return(status == WizardFalse ? 1 : 0);
}
void TimeOnePfRFFT(int count, int fft_log_size, float signal_value,
int signal_type) {
struct AlignedPtr* x_aligned;
struct AlignedPtr* y_aligned;
struct AlignedPtr* z_aligned;
struct AlignedPtr* y_tmp_aligned;
float* x;
struct ComplexFloat* y;
OMX_F32* z;
float* y_true;
float* y_tmp;
int n;
int fft_size;
struct timeval start_time;
struct timeval end_time;
double elapsed_time;
PFFFT_Setup *s;
struct SnrResult snr_forward;
struct SnrResult snr_inverse;
fft_size = 1 << fft_log_size;
x_aligned = AllocAlignedPointer(32, sizeof(*x) * fft_size);
y_aligned = AllocAlignedPointer(32, sizeof(*y) * (fft_size + 2));
z_aligned = AllocAlignedPointer(32, sizeof(*z) * fft_size);
y_tmp_aligned = AllocAlignedPointer(32, sizeof(*y_tmp) * (fft_size + 2));
y_true = (float*) malloc(sizeof(*y_true) * 2 * fft_size);
x = x_aligned->aligned_pointer_;
y = y_aligned->aligned_pointer_;
z = z_aligned->aligned_pointer_;
y_tmp = y_tmp_aligned->aligned_pointer_;
s = pffft_new_setup(fft_size, PFFFT_REAL);
if (!s) {
fprintf(stderr, "TimeOnePfRFFT: Could not initialize structure for order %d\n",
fft_log_size);
}
GenerateRealFloatSignal(x, (struct ComplexFloat*) y_true, fft_size, signal_type, signal_value);
if (do_forward_test) {
GetUserTime(&start_time);
for (n = 0; n < count; ++n) {
pffft_transform_ordered(s, (float*)x, (float*)y, NULL, PFFFT_FORWARD);
}
GetUserTime(&end_time);
elapsed_time = TimeDifference(&start_time, &end_time);
/*
* Arrange the output of the FFT to match the expected output.
*/
y[fft_size / 2].Re = y[0].Im;
y[fft_size / 2].Im = 0;
y[0].Im = 0;
CompareComplexFloat(&snr_forward, (OMX_FC32*) y, (OMX_FC32*) y_true, fft_size / 2 + 1);
PrintResult("Forward PFFFT RFFT", fft_log_size, elapsed_time, count, snr_forward.complex_snr_);
if (verbose >= 255) {
printf("FFT Actual:\n");
DumpArrayComplexFloat("y", fft_size / 2 + 1, (OMX_FC32*) y);
printf("FFT Expected:\n");
DumpArrayComplexFloat("true", fft_size / 2 + 1, (OMX_FC32*) y_true);
}
}
if (do_inverse_test) {
float scale = 1.0 / fft_size;
/* Copy y_true to true, but arrange the values according to what rdft wants. */
memcpy(y_tmp, y_true, sizeof(y_tmp[0]) * fft_size);
y_tmp[1] = y_true[fft_size / 2];
GetUserTime(&start_time);
for (n = 0; n < count; ++n) {
int m;
pffft_transform_ordered(s, (float*)y_tmp, (float*)z, NULL, PFFFT_BACKWARD);
/*
* Need to include cost of scaling the inverse
*/
ScaleVector(z, fft_size, fft_size);
}
GetUserTime(&end_time);
elapsed_time = TimeDifference(&start_time, &end_time);
CompareFloat(&snr_inverse, (OMX_F32*) z, (OMX_F32*) x, fft_size);
PrintResult("Inverse PFFFT RFFT", fft_log_size, elapsed_time, count, snr_inverse.complex_snr_);
if (verbose >= 255) {
printf("IFFT Actual:\n");
DumpArrayFloat("z", fft_size, z);
printf("IFFT Expected:\n");
DumpArrayFloat("x", fft_size, x);
}
}
FreeAlignedPointer(x_aligned);
FreeAlignedPointer(y_aligned);
FreeAlignedPointer(z_aligned);
FreeAlignedPointer(y_tmp_aligned);
pffft_destroy_setup(s);
free(y_true);
}
/* Argument s16s32:
* S32: Calculate RFFT16 with 32 bit complex FFT;
* otherwise: Calculate RFFT16 with 16 bit complex FFT.
*/
void TimeOneRFFT16(int count, int fft_log_size, float signal_value,
int signal_type, s16_s32 s16s32) {
OMX_S16* x;
OMX_S32* y;
OMX_S16* z;
OMX_S32* y_true;
OMX_F32* xr;
OMX_F32* yrTrue;
struct AlignedPtr* x_aligned;
struct AlignedPtr* y_aligned;
struct AlignedPtr* z_aligned;
struct AlignedPtr* y_trueAligned;
struct AlignedPtr* xr_aligned;
struct AlignedPtr* yr_true_aligned;
OMX_S16* temp16;
OMX_S32* temp32;
OMX_INT n, fft_spec_buffer_size;
OMXResult status;
OMXFFTSpec_R_S16 * fft_fwd_spec = NULL;
OMXFFTSpec_R_S16 * fft_inv_spec = NULL;
int fft_size;
struct timeval start_time;
struct timeval end_time;
double elapsed_time;
int scaleFactor;
fft_size = 1 << fft_log_size;
scaleFactor = fft_log_size;
x_aligned = AllocAlignedPointer(32, sizeof(*x) * fft_size);
y_aligned = AllocAlignedPointer(32, sizeof(*y) * (fft_size + 2));
z_aligned = AllocAlignedPointer(32, sizeof(*z) * fft_size);
y_trueAligned = AllocAlignedPointer(32, sizeof(*y_true) * (fft_size + 2));
xr_aligned = AllocAlignedPointer(32, sizeof(*xr) * fft_size);
yr_true_aligned = AllocAlignedPointer(32, sizeof(*yrTrue) * (fft_size + 2));
x = x_aligned->aligned_pointer_;
y = y_aligned->aligned_pointer_;
z = z_aligned->aligned_pointer_;
y_true = y_trueAligned->aligned_pointer_;
xr = xr_aligned->aligned_pointer_;
yrTrue = yr_true_aligned->aligned_pointer_;
temp16 = (OMX_S16*) malloc(sizeof(*temp16) * fft_size);
temp32 = (OMX_S32*) malloc(sizeof(*temp32) * fft_size);
GenerateRFFT16Signal(x, (OMX_SC32*) y_true, fft_size, signal_type,
signal_value);
/*
* Generate a real version so we can measure scaling costs
*/
GenerateRealFloatSignal(xr, (OMX_FC32*) yrTrue, fft_size, signal_type,
signal_value);
if(s16s32 == S32) {
status = omxSP_FFTGetBufSize_R_S16S32(fft_log_size, &fft_spec_buffer_size);
fft_fwd_spec = malloc(fft_spec_buffer_size);
fft_inv_spec = malloc(fft_spec_buffer_size);
status = omxSP_FFTInit_R_S16S32(fft_fwd_spec, fft_log_size);
status = omxSP_FFTInit_R_S16S32(fft_inv_spec, fft_log_size);
}
else {
status = omxSP_FFTGetBufSize_R_S16(fft_log_size, &fft_spec_buffer_size);
fft_fwd_spec = malloc(fft_spec_buffer_size);
fft_inv_spec = malloc(fft_spec_buffer_size);
status = omxSP_FFTInit_R_S16(fft_fwd_spec, fft_log_size);
status = omxSP_FFTInit_R_S16(fft_inv_spec, fft_log_size);
}
if (do_forward_test) {
if (include_conversion) {
int k;
float factor = -1;
GetUserTime(&start_time);
for (k = 0; k < count; ++k) {
/*
* Spend some time computing the max of the signal, and then scaling it.
*/
for (n = 0; n < fft_size; ++n) {
if (abs(xr[n]) > factor) {
factor = abs(xr[n]);
}
}
factor = 32767 / factor;
for (n = 0; n < fft_size; ++n) {
temp16[n] = factor * xr[n];
}
if(s16s32 == S32) {
status = omxSP_FFTFwd_RToCCS_S16S32_Sfs(x, y,
(OMXFFTSpec_R_S16S32*)fft_fwd_spec, (OMX_INT) scaleFactor);
}
else {
status = omxSP_FFTFwd_RToCCS_S16_Sfs(x, (OMX_S16*)y,
(OMXFFTSpec_R_S16*)fft_fwd_spec, (OMX_INT) scaleFactor);
}
/*
* Now spend some time converting the fixed-point FFT back to float.
*/
factor = 1 / factor;
for (n = 0; n < fft_size + 2; ++n) {
xr[n] = y[n] * factor;
}
}
GetUserTime(&end_time);
} else {
float factor = -1;
GetUserTime(&start_time);
for (n = 0; n < count; ++n) {
if(s16s32 == S32) {
status = omxSP_FFTFwd_RToCCS_S16S32_Sfs(x, y,
(OMXFFTSpec_R_S16S32*)fft_fwd_spec, (OMX_INT) scaleFactor);
}
else {
status = omxSP_FFTFwd_RToCCS_S16_Sfs(x, (OMX_S16*)y,
(OMXFFTSpec_R_S16*)fft_fwd_spec, (OMX_INT) scaleFactor);
}
}
GetUserTime(&end_time);
}
elapsed_time = TimeDifference(&start_time, &end_time);
if(s16s32 == S32) {
PrintResultNoSNR("Forward RFFT16 (with S32)",
fft_log_size, elapsed_time, count);
}
else {
PrintResultNoSNR("Forward RFFT16 (with S16)",
fft_log_size, elapsed_time, count);
}
}
if (do_inverse_test) {
if (include_conversion) {
int k;
float factor = -1;
GetUserTime(&start_time);
for (k = 0; k < count; ++k) {
/*
* Spend some time scaling the FFT signal to fixed point.
*/
for (n = 0; n < fft_size; ++n) {
if (abs(yrTrue[n]) > factor) {
factor = abs(yrTrue[n]);
}
}
for (n = 0; n < fft_size; ++n) {
temp32[n] = factor * yrTrue[n];
}
if(s16s32 == S32) {
status = omxSP_FFTInv_CCSToR_S32S16_Sfs(y, z,
(OMXFFTSpec_R_S16S32*)fft_inv_spec, 0);
}
else {
status = omxSP_FFTInv_CCSToR_S16_Sfs((OMX_S16*)y, z,
(OMXFFTSpec_R_S16*)fft_inv_spec, 0);
}
/*
* Spend some time converting the result back to float
*/
factor = 1 / factor;
for (n = 0; n < fft_size; ++n) {
xr[n] = factor * z[n];
}
}
GetUserTime(&end_time);
} else {
GetUserTime(&start_time);
for (n = 0; n < count; ++n) {
if(s16s32 == S32) {
status = omxSP_FFTInv_CCSToR_S32S16_Sfs(y, z,
(OMXFFTSpec_R_S16S32*)fft_inv_spec, 0);
}
else {
status = omxSP_FFTInv_CCSToR_S16_Sfs((OMX_S16*)y, z,
(OMXFFTSpec_R_S16*)fft_inv_spec, 0);
}
}
GetUserTime(&end_time);
}
elapsed_time = TimeDifference(&start_time, &end_time);
if(s16s32 == S32) {
PrintResultNoSNR("Inverse RFFT16 (with S32)",
fft_log_size, elapsed_time, count);
}
else {
PrintResultNoSNR("Inverse RFFT16 (with S16)",
fft_log_size, elapsed_time, count);
}
}
FreeAlignedPointer(x_aligned);
FreeAlignedPointer(y_aligned);
FreeAlignedPointer(z_aligned);
FreeAlignedPointer(y_trueAligned);
FreeAlignedPointer(xr_aligned);
FreeAlignedPointer(yr_true_aligned);
free(fft_fwd_spec);
free(fft_inv_spec);
}
void TimeOneSC16FFT(int count, int fft_log_size, float signal_value,
int signal_type) {
OMX_SC16* x;
OMX_SC16* y;
OMX_SC16* z;
struct AlignedPtr* x_aligned;
struct AlignedPtr* y_aligned;
struct AlignedPtr* z_aligned;
OMX_SC16* y_true;
OMX_SC16* temp16a;
OMX_SC16* temp16b;
OMX_INT n, fft_spec_buffer_size;
OMXResult status;
OMXFFTSpec_C_SC16 * fft_fwd_spec = NULL;
OMXFFTSpec_C_SC16 * fft_inv_spec = NULL;
int fft_size;
struct timeval start_time;
struct timeval end_time;
double elapsed_time;
struct SnrResult snr_forward;
struct SnrResult snr_inverse;
fft_size = 1 << fft_log_size;
x_aligned = AllocAlignedPointer(32, sizeof(*x) * fft_size);
y_aligned = AllocAlignedPointer(32, sizeof(*y) * fft_size);
z_aligned = AllocAlignedPointer(32, sizeof(*z) * fft_size);
y_true = (OMX_SC16*) malloc(sizeof(*y_true) * fft_size);
temp16a = (OMX_SC16*) malloc(sizeof(*temp16a) * fft_size);
temp16b = (OMX_SC16*) malloc(sizeof(*temp16b) * fft_size);
x = x_aligned->aligned_pointer_;
y = y_aligned->aligned_pointer_;
z = z_aligned->aligned_pointer_;
generateSC16Signal(x, y_true, fft_size, signal_type, signal_value);
status = omxSP_FFTGetBufSize_C_SC16(fft_log_size, &fft_spec_buffer_size);
fft_fwd_spec = (OMXFFTSpec_C_SC16*) malloc(fft_spec_buffer_size);
fft_inv_spec = (OMXFFTSpec_C_SC16*) malloc(fft_spec_buffer_size);
status = omxSP_FFTInit_C_SC16(fft_fwd_spec, fft_log_size);
status = omxSP_FFTInit_C_SC16(fft_inv_spec, fft_log_size);
if (do_forward_test) {
if (include_conversion) {
int k;
float factor = -1;
GetUserTime(&start_time);
for (k = 0; k < count; ++k) {
for (n = 0; n < fft_size; ++n) {
if (abs(x[n].Re) > factor) {
factor = abs(x[n].Re);
}
if (abs(x[n].Im) > factor) {
factor = abs(x[n].Im);
}
}
factor = ((1 << 15) - 1) / factor;
for (n = 0; n < fft_size; ++n) {
temp16a[n].Re = factor * x[n].Re;
temp16a[n].Im = factor * x[n].Im;
}
omxSP_FFTFwd_CToC_SC16_Sfs(x, y, fft_fwd_spec, 0);
factor = 1 / factor;
for (n = 0; n < fft_size; ++n) {
temp16b[n].Re = y[n].Re * factor;
temp16b[n].Im = y[n].Im * factor;
}
}
GetUserTime(&end_time);
} else {
GetUserTime(&start_time);
for (n = 0; n < count; ++n) {
omxSP_FFTFwd_CToC_SC16_Sfs(x, y, fft_fwd_spec, 0);
}
GetUserTime(&end_time);
}
elapsed_time = TimeDifference(&start_time, &end_time);
CompareComplex16(&snr_forward, y, y_true, fft_size);
PrintResult("Forward SC16 FFT", fft_log_size, elapsed_time, count, snr_forward.complex_snr_);
}
if (do_inverse_test) {
if (include_conversion) {
int k;
float factor = -1;
GetUserTime(&start_time);
for (k = 0; k < count; ++k) {
for (n = 0; n < fft_size; ++n) {
if (abs(x[n].Re) > factor) {
factor = abs(x[n].Re);
}
if (abs(x[n].Im) > factor) {
factor = abs(x[n].Im);
}
}
factor = ((1 << 15) - 1) / factor;
for (n = 0; n < fft_size; ++n) {
temp16a[n].Re = factor * x[n].Re;
temp16a[n].Im = factor * x[n].Im;
}
status = omxSP_FFTInv_CToC_SC16_Sfs(y, z, fft_inv_spec, 0);
factor = 1 / factor;
for (n = 0; n < fft_size; ++n) {
temp16b[n].Re = y[n].Re * factor;
temp16b[n].Im = y[n].Im * factor;
}
}
GetUserTime(&end_time);
} else {
GetUserTime(&start_time);
for (n = 0; n < count; ++n) {
status = omxSP_FFTInv_CToC_SC16_Sfs(y, z, fft_inv_spec, 0);
}
GetUserTime(&end_time);
}
elapsed_time = TimeDifference(&start_time, &end_time);
CompareComplex16(&snr_inverse, z, x, fft_size);
PrintResult("Inverse SC16 FFT", fft_log_size, elapsed_time, count, snr_inverse.complex_snr_);
}
FreeAlignedPointer(x_aligned);
FreeAlignedPointer(y_aligned);
FreeAlignedPointer(z_aligned);
free(temp16a);
free(temp16b);
free(fft_fwd_spec);
free(fft_inv_spec);
}
void TimeOneFloatFFT(int count, int fft_log_size, float signal_value,
int signal_type) {
struct AlignedPtr* x_aligned;
struct AlignedPtr* y_aligned;
struct AlignedPtr* z_aligned;
struct AlignedPtr* y_true_aligned;
struct ComplexFloat* x;
struct ComplexFloat* y;
OMX_FC32* z;
struct ComplexFloat* y_true;
OMX_INT n, fft_spec_buffer_size;
OMXFFTSpec_C_FC32 * fft_fwd_spec = NULL;
OMXFFTSpec_C_FC32 * fft_inv_spec = NULL;
int fft_size;
struct timeval start_time;
struct timeval end_time;
double elapsed_time;
struct SnrResult snr_forward;
struct SnrResult snr_inverse;
fft_size = 1 << fft_log_size;
x_aligned = AllocAlignedPointer(32, sizeof(*x) * fft_size);
y_aligned = AllocAlignedPointer(32, sizeof(*y) * (fft_size + 2));
z_aligned = AllocAlignedPointer(32, sizeof(*z) * fft_size);
y_true_aligned = AllocAlignedPointer(32, sizeof(*z) * fft_size);
x = x_aligned->aligned_pointer_;
y = y_aligned->aligned_pointer_;
z = z_aligned->aligned_pointer_;
y_true = y_true_aligned->aligned_pointer_;
GenerateTestSignalAndFFT(x, y_true, fft_size, signal_type, signal_value, 0);
omxSP_FFTGetBufSize_C_FC32(fft_log_size, &fft_spec_buffer_size);
fft_fwd_spec = (OMXFFTSpec_C_FC32*) malloc(fft_spec_buffer_size);
fft_inv_spec = (OMXFFTSpec_C_FC32*) malloc(fft_spec_buffer_size);
omxSP_FFTInit_C_FC32(fft_fwd_spec, fft_log_size);
omxSP_FFTInit_C_FC32(fft_inv_spec, fft_log_size);
if (do_forward_test) {
GetUserTime(&start_time);
for (n = 0; n < count; ++n) {
FORWARD_FLOAT_FFT((OMX_FC32*) x, (OMX_FC32*) y, fft_fwd_spec);
}
GetUserTime(&end_time);
elapsed_time = TimeDifference(&start_time, &end_time);
CompareComplexFloat(&snr_forward, (OMX_FC32*) y, (OMX_FC32*) y_true, fft_size);
PrintResult("Forward Float FFT", fft_log_size, elapsed_time, count, snr_forward.complex_snr_);
}
if (do_inverse_test) {
GetUserTime(&start_time);
for (n = 0; n < count; ++n) {
INVERSE_FLOAT_FFT((OMX_FC32*) y_true, z, fft_inv_spec);
}
GetUserTime(&end_time);
elapsed_time = TimeDifference(&start_time, &end_time);
CompareComplexFloat(&snr_inverse, (OMX_FC32*) z, (OMX_FC32*) x, fft_size);
PrintResult("Inverse Float FFT", fft_log_size, elapsed_time, count, snr_inverse.complex_snr_);
}
FreeAlignedPointer(x_aligned);
FreeAlignedPointer(y_aligned);
FreeAlignedPointer(z_aligned);
FreeAlignedPointer(y_true_aligned);
free(fft_fwd_spec);
free(fft_inv_spec);
}
void TimeOneFloatRFFT(int count, int fft_log_size, float signal_value,
int signal_type) {
OMX_F32* x; /* Source */
OMX_F32* y; /* Transform */
OMX_F32* z; /* Inverse transform */
OMX_F32* y_true; /* True FFT */
struct AlignedPtr* x_aligned;
struct AlignedPtr* y_aligned;
struct AlignedPtr* z_aligned;
struct AlignedPtr* y_true_aligned;
OMX_INT n, fft_spec_buffer_size;
OMXResult status;
OMXFFTSpec_R_F32 * fft_fwd_spec = NULL;
OMXFFTSpec_R_F32 * fft_inv_spec = NULL;
int fft_size;
struct timeval start_time;
struct timeval end_time;
double elapsed_time;
struct SnrResult snr_forward;
struct SnrResult snr_inverse;
fft_size = 1 << fft_log_size;
x_aligned = AllocAlignedPointer(32, sizeof(*x) * fft_size);
/* The transformed value is in CCS format and is has fft_size + 2 values */
y_aligned = AllocAlignedPointer(32, sizeof(*y) * (fft_size + 2));
z_aligned = AllocAlignedPointer(32, sizeof(*z) * fft_size);
y_true_aligned = AllocAlignedPointer(32, sizeof(*z) * (fft_size + 2));
x = x_aligned->aligned_pointer_;
y = y_aligned->aligned_pointer_;
z = z_aligned->aligned_pointer_;
y_true = y_true_aligned->aligned_pointer_;
GenerateRealFloatSignal(x, (OMX_FC32*) y_true, fft_size, signal_type,
signal_value);
status = omxSP_FFTGetBufSize_R_F32(fft_log_size, &fft_spec_buffer_size);
fft_fwd_spec = (OMXFFTSpec_R_F32*) malloc(fft_spec_buffer_size);
fft_inv_spec = (OMXFFTSpec_R_F32*) malloc(fft_spec_buffer_size);
status = omxSP_FFTInit_R_F32(fft_fwd_spec, fft_log_size);
status = omxSP_FFTInit_R_F32(fft_inv_spec, fft_log_size);
if (do_forward_test) {
GetUserTime(&start_time);
for (n = 0; n < count; ++n) {
FORWARD_FLOAT_RFFT(x, y, fft_fwd_spec);
}
GetUserTime(&end_time);
elapsed_time = TimeDifference(&start_time, &end_time);
CompareComplexFloat(&snr_forward, (OMX_FC32*) y, (OMX_FC32*) y_true, fft_size / 2 + 1);
PrintResult("Forward Float RFFT", fft_log_size, elapsed_time, count, snr_forward.complex_snr_);
}
if (do_inverse_test) {
GetUserTime(&start_time);
for (n = 0; n < count; ++n) {
INVERSE_FLOAT_RFFT(y_true, z, fft_inv_spec);
}
GetUserTime(&end_time);
elapsed_time = TimeDifference(&start_time, &end_time);
CompareFloat(&snr_inverse, (OMX_F32*) z, (OMX_F32*) x, fft_size);
PrintResult("Inverse Float RFFT", fft_log_size, elapsed_time, count, snr_inverse.complex_snr_);
}
FreeAlignedPointer(x_aligned);
FreeAlignedPointer(y_aligned);
FreeAlignedPointer(z_aligned);
free(fft_fwd_spec);
free(fft_inv_spec);
}
void TimeOneNE10FFT(int count, int fft_log_size, float signal_value,
int signal_type) {
struct AlignedPtr* x_aligned;
struct AlignedPtr* y_aligned;
struct AlignedPtr* z_aligned;
struct ComplexFloat* x;
struct ComplexFloat* y;
OMX_FC32* z;
struct ComplexFloat* y_true;
int n;
ne10_result_t status;
ne10_fft_cfg_float32_t fft_fwd_spec;
int fft_size;
struct timeval start_time;
struct timeval end_time;
double elapsed_time;
struct SnrResult snr_forward;
struct SnrResult snr_inverse;
fft_size = 1 << fft_log_size;
x_aligned = AllocAlignedPointer(32, sizeof(*x) * fft_size);
y_aligned = AllocAlignedPointer(32, sizeof(*y) * 2 * fft_size);
z_aligned = AllocAlignedPointer(32, sizeof(*z) * 2 * fft_size);
y_true = (struct ComplexFloat*) malloc(sizeof(*y_true) * fft_size);
x = x_aligned->aligned_pointer_;
y = y_aligned->aligned_pointer_;
z = z_aligned->aligned_pointer_;
GenerateTestSignalAndFFT(x, y_true, fft_size, signal_type, signal_value, 0);
fft_fwd_spec = ne10_fft_alloc_c2c_float32(fft_size);
if (!fft_fwd_spec) {
fprintf(stderr, "NE10 FFT: Cannot initialize FFT structure for order %d\n", fft_log_size);
return;
}
if (do_forward_test) {
GetUserTime(&start_time);
for (n = 0; n < count; ++n) {
ne10_fft_c2c_1d_float32_neon((ne10_fft_cpx_float32_t *) y,
(ne10_fft_cpx_float32_t *) x,
fft_fwd_spec,
0);
}
GetUserTime(&end_time);
elapsed_time = TimeDifference(&start_time, &end_time);
CompareComplexFloat(&snr_forward, (OMX_FC32*) y, (OMX_FC32*) y_true, fft_size);
PrintResult("Forward NE10 FFT", fft_log_size, elapsed_time, count, snr_forward.complex_snr_);
if (verbose >= 255) {
printf("Input data:\n");
DumpArrayComplexFloat("x", fft_size, (OMX_FC32*) x);
printf("FFT Actual:\n");
DumpArrayComplexFloat("y", fft_size, (OMX_FC32*) y);
printf("FFT Expected:\n");
DumpArrayComplexFloat("true", fft_size, (OMX_FC32*) y_true);
}
}
if (do_inverse_test) {
GetUserTime(&start_time);
for (n = 0; n < count; ++n) {
ne10_fft_c2c_1d_float32_neon((ne10_fft_cpx_float32_t *) z,
(ne10_fft_cpx_float32_t *) y_true,
fft_fwd_spec,
1);
}
GetUserTime(&end_time);
elapsed_time = TimeDifference(&start_time, &end_time);
CompareComplexFloat(&snr_inverse, (OMX_FC32*) z, (OMX_FC32*) x, fft_size);
PrintResult("Inverse NE10 FFT", fft_log_size, elapsed_time, count, snr_inverse.complex_snr_);
if (verbose >= 255) {
printf("Input data:\n");
DumpArrayComplexFloat("y", fft_size, (OMX_FC32*) y_true);
printf("IFFT Actual:\n");
DumpArrayComplexFloat("z", fft_size, z);
printf("IFFT Expected:\n");
DumpArrayComplexFloat("x", fft_size, (OMX_FC32*) x);
}
}
ne10_fft_destroy_c2c_float32(fft_fwd_spec);
FreeAlignedPointer(x_aligned);
FreeAlignedPointer(y_aligned);
FreeAlignedPointer(z_aligned);
free(y_true);
}
void TimeOneRFFT32(int count, int fft_log_size, float signal_value,
int signal_type) {
OMX_S32* x;
OMX_S32* y;
OMX_S32* z;
OMX_S32* y_true;
OMX_F32* xr;
OMX_F32* yrTrue;
struct AlignedPtr* x_aligned;
struct AlignedPtr* y_aligned;
struct AlignedPtr* z_aligned;
struct AlignedPtr* y_true_aligned;
OMX_S32* temp1;
OMX_S32* temp2;
OMX_INT n, fft_spec_buffer_size;
OMXResult status;
OMXFFTSpec_R_S16S32 * fft_fwd_spec = NULL;
OMXFFTSpec_R_S16S32 * fft_inv_spec = NULL;
int fft_size;
struct timeval start_time;
struct timeval end_time;
double elapsed_time;
int scaleFactor = 0;
struct SnrResult snr_forward;
struct SnrResult snr_inverse;
fft_size = 1 << fft_log_size;
x_aligned = AllocAlignedPointer(32, sizeof(*x) * fft_size);
y_aligned = AllocAlignedPointer(32, sizeof(*y) * (fft_size + 2));
z_aligned = AllocAlignedPointer(32, sizeof(*z) * fft_size);
y_true_aligned = AllocAlignedPointer(32, sizeof(*y_true) * (fft_size + 2));
x = x_aligned->aligned_pointer_;
y = y_aligned->aligned_pointer_;
z = z_aligned->aligned_pointer_;
y_true = y_true_aligned->aligned_pointer_;
if (verbose > 3) {
printf("x = %p\n", (void*)x);
printf("y = %p\n", (void*)y);
printf("z = %p\n", (void*)z);
}
xr = (OMX_F32*) malloc(sizeof(*x) * fft_size);
yrTrue = (OMX_F32*) malloc(sizeof(*y) * (fft_size + 2));
temp1 = (OMX_S32*) malloc(sizeof(*temp1) * fft_size);
temp2 = (OMX_S32*) malloc(sizeof(*temp2) * (fft_size + 2));
GenerateRFFT32Signal(x, (OMX_SC32*) y_true, fft_size, signal_type,
signal_value);
if (verbose > 63) {
printf("Signal\n");
printf("n\tx[n]\n");
for (n = 0; n < fft_size; ++n) {
printf("%4d\t%d\n", n, x[n]);
}
}
status = omxSP_FFTGetBufSize_R_S32(fft_log_size, &fft_spec_buffer_size);
if (verbose > 3) {
printf("fft_spec_buffer_size = %d\n", fft_spec_buffer_size);
}
fft_fwd_spec = (OMXFFTSpec_R_S32*) malloc(fft_spec_buffer_size);
fft_inv_spec = (OMXFFTSpec_R_S32*) malloc(fft_spec_buffer_size);
status = omxSP_FFTInit_R_S32(fft_fwd_spec, fft_log_size);
if (status) {
printf("Failed to init forward FFT: status = %d\n", status);
}
status = omxSP_FFTInit_R_S32(fft_inv_spec, fft_log_size);
if (status) {
printf("Failed to init backward FFT: status = %d\n", status);
}
if (do_forward_test) {
if (include_conversion) {
int k;
float factor = -1;
GetUserTime(&start_time);
for (k = 0; k < count; ++k) {
/*
* Spend some time computing the max of the signal, and then scaling it.
*/
for (n = 0; n < fft_size; ++n) {
if (abs(xr[n]) > factor) {
factor = abs(xr[n]);
}
}
factor = (1 << 20) / factor;
for (n = 0; n < fft_size; ++n) {
temp1[n] = factor * xr[n];
}
status = omxSP_FFTFwd_RToCCS_S32_Sfs(x, y, fft_fwd_spec,
(OMX_INT) scaleFactor);
/*
* Now spend some time converting the fixed-point FFT back to float.
*/
factor = 1 / factor;
for (n = 0; n < fft_size + 2; ++n) {
xr[n] = y[n] * factor;
}
}
GetUserTime(&end_time);
} else {
float factor = -1;
GetUserTime(&start_time);
for (n = 0; n < count; ++n) {
status = omxSP_FFTFwd_RToCCS_S32_Sfs(x, y, fft_fwd_spec,
(OMX_INT) scaleFactor);
}
GetUserTime(&end_time);
}
elapsed_time = TimeDifference(&start_time, &end_time);
CompareComplex32(&snr_forward, (OMX_SC32*) y, (OMX_SC32*) y_true, fft_size / 2 + 1);
PrintResult("Forward RFFT32", fft_log_size, elapsed_time, count, snr_forward.complex_snr_);
}
if (do_inverse_test) {
if (include_conversion) {
int k;
float factor = -1;
GetUserTime(&start_time);
for (k = 0; k < count; ++k) {
/*
* Spend some time scaling the FFT signal to fixed point.
*/
for (n = 0; n < fft_size + 2; ++n) {
if (abs(yrTrue[n]) > factor) {
factor = abs(yrTrue[n]);
}
}
for (n = 0; n < fft_size + 2; ++n) {
temp2[n] = factor * yrTrue[n];
}
status = omxSP_FFTInv_CCSToR_S32_Sfs(y_true, z, fft_inv_spec, 0);
/*
* Spend some time converting the result back to float
*/
factor = 1 / factor;
for (n = 0; n < fft_size; ++n) {
xr[n] = factor * z[n];
}
}
GetUserTime(&end_time);
} else {
GetUserTime(&start_time);
for (n = 0; n < count; ++n) {
status = omxSP_FFTInv_CCSToR_S32_Sfs(y_true, z, fft_inv_spec, 0);
}
GetUserTime(&end_time);
}
elapsed_time = TimeDifference(&start_time, &end_time);
CompareReal32(&snr_inverse, z, x, fft_size);
PrintResult("Inverse RFFT32", fft_log_size, elapsed_time, count, snr_inverse.real_snr_);
}
FreeAlignedPointer(x_aligned);
FreeAlignedPointer(y_aligned);
FreeAlignedPointer(z_aligned);
FreeAlignedPointer(y_true_aligned);
free(fft_fwd_spec);
free(fft_inv_spec);
}
void TimeOnePfFFT(int count, int fft_log_size, float signal_value,
int signal_type) {
struct AlignedPtr* x_aligned;
struct AlignedPtr* y_aligned;
struct AlignedPtr* z_aligned;
struct AlignedPtr* y_true_aligned;
struct ComplexFloat* x;
struct ComplexFloat* y;
OMX_FC32* z;
struct ComplexFloat* y_true;
int n;
int fft_size;
struct timeval start_time;
struct timeval end_time;
double elapsed_time;
PFFFT_Setup *s;
struct SnrResult snr_forward;
struct SnrResult snr_inverse;
fft_size = 1 << fft_log_size;
x_aligned = AllocAlignedPointer(32, sizeof(*x) * fft_size);
y_aligned = AllocAlignedPointer(32, sizeof(*y) * (fft_size + 2));
z_aligned = AllocAlignedPointer(32, sizeof(*z) * fft_size);
y_true_aligned = AllocAlignedPointer(32, sizeof(*y_true) * fft_size);
x = x_aligned->aligned_pointer_;
y = y_aligned->aligned_pointer_;
z = z_aligned->aligned_pointer_;
y_true = y_true_aligned->aligned_pointer_;
s = pffft_new_setup(fft_size, PFFFT_COMPLEX);
if (!s) {
fprintf(stderr, "TimeOnePfFFT: Could not initialize structure for order %d\n",
fft_log_size);
}
GenerateTestSignalAndFFT(x, y_true, fft_size, signal_type, signal_value, 0);
if (do_forward_test) {
GetUserTime(&start_time);
for (n = 0; n < count; ++n) {
pffft_transform_ordered(s, (float*)x, (float*)y, NULL, PFFFT_FORWARD);
}
GetUserTime(&end_time);
elapsed_time = TimeDifference(&start_time, &end_time);
CompareComplexFloat(&snr_forward, (OMX_FC32*) y, (OMX_FC32*) y_true, fft_size);
PrintResult("Forward PFFFT FFT", fft_log_size, elapsed_time, count, snr_forward.complex_snr_);
if (verbose >= 255) {
printf("FFT Actual:\n");
DumpArrayComplexFloat("y", fft_size, (OMX_FC32*) y);
printf("FFT Expected:\n");
DumpArrayComplexFloat("true", fft_size, (OMX_FC32*) y_true);
}
}
if (do_inverse_test) {
float scale = 1.0 / fft_size;
memcpy(y, y_true, sizeof(*y) * (fft_size + 2));
GetUserTime(&start_time);
for (n = 0; n < count; ++n) {
int m;
pffft_transform_ordered(s, (float*)y_true, (float*)z, NULL, PFFFT_BACKWARD);
/*
* Need to include cost of scaling the inverse
*/
ScaleVector((OMX_F32*) z, 2 * fft_size, fft_size);
}
GetUserTime(&end_time);
elapsed_time = TimeDifference(&start_time, &end_time);
CompareComplexFloat(&snr_inverse, (OMX_FC32*) z, (OMX_FC32*) x, fft_size);
PrintResult("Inverse PFFFT FFT", fft_log_size, elapsed_time, count, snr_inverse.complex_snr_);
if (verbose >= 255) {
printf("IFFT Actual:\n");
DumpArrayComplexFloat("z", fft_size, z);
printf("IFFT Expected:\n");
DumpArrayComplexFloat("x", fft_size, (OMX_FC32*) x);
}
}
FreeAlignedPointer(x_aligned);
FreeAlignedPointer(y_aligned);
FreeAlignedPointer(z_aligned);
FreeAlignedPointer(y_true_aligned);
pffft_destroy_setup(s);
}
void TimeOneNE10RFFT(int count, int fft_log_size, float signal_value,
int signal_type) {
OMX_F32* x; /* Source */
OMX_FC32* y; /* Transform */
OMX_F32* z; /* Inverse transform */
OMX_F32* temp;
OMX_F32* y_true; /* True FFT */
struct AlignedPtr* x_aligned;
struct AlignedPtr* y_aligned;
struct AlignedPtr* z_aligned;
int n;
ne10_result_t status;
ne10_fft_r2c_cfg_float32_t fft_fwd_spec;
int fft_size;
struct timeval start_time;
struct timeval end_time;
double elapsed_time;
struct SnrResult snr_forward;
struct SnrResult snr_inverse;
fft_size = 1 << fft_log_size;
x_aligned = AllocAlignedPointer(32, sizeof(*x) * 4 * fft_size);
/* The transformed value is in CCS format and is has fft_size + 2 values */
y_aligned = AllocAlignedPointer(32, sizeof(*y) * (4 * fft_size + 2));
z_aligned = AllocAlignedPointer(32, sizeof(*z) * 4 * fft_size);
x = x_aligned->aligned_pointer_;
y = y_aligned->aligned_pointer_;
z = z_aligned->aligned_pointer_;
y_true = (OMX_F32*) malloc(sizeof(*y_true) * (fft_size + 2));
GenerateRealFloatSignal(x, (struct ComplexFloat*) y_true, fft_size, signal_type,
signal_value);
fft_fwd_spec = ne10_fft_alloc_r2c_float32(fft_size);
if (!fft_fwd_spec) {
fprintf(stderr, "NE10 RFFT: Cannot initialize FFT structure for order %d\n", fft_log_size);
return;
}
if (do_forward_test) {
GetUserTime(&start_time);
for (n = 0; n < count; ++n) {
ne10_fft_r2c_1d_float32_neon((ne10_fft_cpx_float32_t *) y,
x,
fft_fwd_spec);
}
GetUserTime(&end_time);
elapsed_time = TimeDifference(&start_time, &end_time);
CompareComplexFloat(&snr_forward, (OMX_FC32*) y, (OMX_FC32*) y_true, fft_size / 2 + 1);
PrintResult("Forward NE10 RFFT", fft_log_size, elapsed_time, count, snr_forward.complex_snr_);
if (verbose >= 255) {
printf("FFT Actual:\n");
DumpArrayComplexFloat("y", fft_size / 2 + 1, y);
printf("FFT Expected:\n");
DumpArrayComplexFloat("true", fft_size / 2 + 1, (OMX_FC32*) y_true);
}
}
if (do_inverse_test) {
// Ne10 FFTs destroy the input.
GetUserTime(&start_time);
for (n = 0; n < count; ++n) {
//memcpy(y, y_true, (fft_size >> 1) * sizeof(*y));
// The inverse appears not to be working.
ne10_fft_c2r_1d_float32_neon(z,
(ne10_fft_cpx_float32_t *) y_true,
fft_fwd_spec);
}
GetUserTime(&end_time);
elapsed_time = TimeDifference(&start_time, &end_time);
CompareFloat(&snr_inverse, (OMX_F32*) z, (OMX_F32*) x, fft_size);
PrintResult("Inverse NE10 RFFT", fft_log_size, elapsed_time, count, snr_inverse.complex_snr_);
if (verbose >= 255) {
printf("IFFT Actual:\n");
DumpArrayFloat("z", fft_size, z);
printf("IFFT Expected:\n");
DumpArrayFloat("x", fft_size, x);
}
}
ne10_fft_destroy_r2c_float32(fft_fwd_spec);
FreeAlignedPointer(x_aligned);
FreeAlignedPointer(y_aligned);
FreeAlignedPointer(z_aligned);
}
/*
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% %
% %
% %
% D e s c r i b e I m a g e %
% %
% %
% %
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
% DescribeImage() describes an image by printing its attributes to the file.
% Attributes include the image width, height, size, and others.
%
% The format of the DescribeImage method is:
%
% void DescribeImage(Image *image,FILE *file,const MagickBool verbose)
%
% A description of each parameter follows:
%
% o image: The image.
%
% o file: The file, typically stdout.
%
% o verbose: A value other than zero prints more detailed information
% about the image. Values greater than one enable counting the number of
% colors in the image.
%
%
*/
MagickExport MagickPassFail DescribeImage(Image *image,FILE *file,
const MagickBool verbose)
{
char
color[MaxTextExtent],
format[MaxTextExtent];
const unsigned char
*profile;
size_t
profile_length;
const ImageAttribute
*attribute;
const MagickInfo
*magick_info;
double
elapsed_time,
user_time;
unsigned long
columns,
rows;
magick_int64_t
pixels_per_second;
Image
*p;
long
y;
register long
i,
x;
unsigned long
count;
assert(image != (Image *) NULL);
assert(image->signature == MagickSignature);
assert(file != (FILE *) NULL);
elapsed_time=GetElapsedTime(&image->timer);
user_time=GetUserTime(&image->timer);
GetTimerInfo(&image->timer);
if (!verbose)
{
/*
Display summary info about the image.
*/
if (*image->magick_filename != '\0')
if (LocaleCompare(image->magick_filename,image->filename) != 0)
(void) fprintf(file,"%.1024s=>",image->magick_filename);
if ((image->previous == (Image *) NULL) &&
(image->next == (Image *) NULL) && (image->scene == 0))
(void) fprintf(file,"%.1024s ",image->filename);
else
(void) fprintf(file,"%.1024s[%lu] ",image->filename,image->scene);
(void) fprintf(file,"%.1024s ",image->magick);
columns=image->columns;
rows=image->rows;
if ((image->magick_columns != 0) || (image->magick_rows != 0))
if ((image->magick_columns != image->columns) ||
(image->magick_rows != image->rows))
{
columns=image->magick_columns;
rows=image->magick_rows;
(void) fprintf(file,"%lux%lu=>",image->magick_columns,
image->magick_rows);
}
(void) fprintf(file,"%lux%lu%+ld%+ld ",image->columns,image->rows,
image->page.x,image->page.y);
if (image->storage_class == DirectClass)
{
(void) fprintf(file,"DirectClass ");
if (image->total_colors != 0)
{
FormatSize(image->total_colors,format);
(void) fprintf(file,"%.1024s ",format);
}
}
else if (image->total_colors <= image->colors)
(void) fprintf(file,"PseudoClass %uc ",image->colors);
else
{
(void) fprintf(file,"PseudoClass %lu=>%uc ",image->total_colors,
image->colors);
(void) fprintf(file,"%ld/%.6f/%.6fe ",
(long) image->error.mean_error_per_pixel,
image->error.normalized_mean_error,
image->error.normalized_maximum_error);
}
(void) fprintf(file,"%u-bit ",image->depth);
if (GetBlobSize(image) != 0)
{
FormatSize(GetBlobSize(image),format);
(void) fprintf(file,"%.1024s ",format);
}
(void) fprintf(file,"%0.3fu %ld:%02ld",user_time,
(long) (elapsed_time/60.0),
(long) ceil(fmod(elapsed_time,60.0)));
/*
Only display pixel read rate if the time accumulated is at
least six times the timer's resolution (typically 0.01 on
Unix).
*/
if (elapsed_time >= GetTimerResolution()*6)
{
pixels_per_second=(magick_int64_t) ((double) rows*columns/ elapsed_time);
FormatSize(pixels_per_second,format);
(void) fprintf(file," (%s pixels/s)",format);
}
(void) fprintf(file,"\n");
return (ferror(file) ? MagickFail : MagickPass);
}
/*
Display verbose info about the image.
*/
(void) SignatureImage(image);
if (verbose > 1)
image->total_colors=GetNumberColors(image,(FILE *) NULL,&image->exception);
(void) fprintf(file,"Image: %.1024s\n",image->filename);
magick_info=GetMagickInfo(image->magick,&image->exception);
if ((magick_info == (const MagickInfo *) NULL) ||
(*magick_info->description == '\0'))
(void) fprintf(file," Format: %.1024s\n",image->magick);
else
(void) fprintf(file," Format: %.1024s (%.1024s)\n",image->magick,
magick_info->description);
(void) fprintf(file," Geometry: %lux%lu\n",image->columns,image->rows);
if (image->storage_class == DirectClass)
(void) fprintf(file," Class: DirectClass\n");
else
(void) fprintf(file," Class: PseudoClass\n");
if ((image->magick_columns != 0) || (image->magick_rows != 0))
if ((image->magick_columns != image->columns) ||
(image->magick_rows != image->rows))
(void) fprintf(file," Base geometry: %lux%lu\n",image->magick_columns,
image->magick_rows);
(void) fprintf(file," Type: ");
switch (GetImageType(image,&image->exception))
{
case BilevelType:
(void) fprintf(file,"bilevel");
break;
case GrayscaleType:
(void) fprintf(file,"grayscale");
break;
case GrayscaleMatteType:
(void) fprintf(file,"grayscale with transparency");
break;
case PaletteType:
(void) fprintf(file,"palette");
break;
case PaletteMatteType:
(void) fprintf(file,"palette with transparency");
break;
case TrueColorType:
(void) fprintf(file,"true color");
break;
case TrueColorMatteType:
(void) fprintf(file,"true color with transparency");
break;
case ColorSeparationType:
(void) fprintf(file,"color separated");
break;
case ColorSeparationMatteType:
(void) fprintf(file,"color separated with transparency");
break;
default:
(void) fprintf(file,"undefined");
break;
}
(void) fprintf(file,"\n");
(void) fprintf(file," Depth: %lu bits-per-pixel component\n",
GetImageDepth(image,&image->exception));
(void) fprintf(file," Channel Depths:\n");
if (image->colorspace == CMYKColorspace)
{
(void) fprintf(file," Cyan: %u bits\n",
GetImageChannelDepth(image,CyanChannel,&image->exception));
(void) fprintf(file," Magenta: %u bits\n",
GetImageChannelDepth(image,MagentaChannel,&image->exception));
(void) fprintf(file," Yellow: %u bits\n",
GetImageChannelDepth(image,YellowChannel,&image->exception));
(void) fprintf(file," Black: %u bits\n",
GetImageChannelDepth(image,BlackChannel,&image->exception));
}
else if ((IsGrayColorspace(image->colorspace)) ||
(image->is_grayscale == True))
{
(void) fprintf(file," Gray: %u bits\n",
GetImageChannelDepth(image,RedChannel,&image->exception));
}
else
{
(void) fprintf(file," Red: %u bits\n",
GetImageChannelDepth(image,RedChannel,&image->exception));
(void) fprintf(file," Green: %u bits\n",
GetImageChannelDepth(image,GreenChannel,&image->exception));
(void) fprintf(file," Blue: %u bits\n",
GetImageChannelDepth(image,BlueChannel,&image->exception));
}
if (image->matte)
(void) fprintf(file," Opacity: %u bits\n",
GetImageChannelDepth(image,OpacityChannel,&image->exception));
(void) fprintf(file," Channel Statistics:\n");
{
ImageStatistics
statistics;
(void) GetImageStatistics(image,&statistics,&image->exception);
if (image->colorspace == CMYKColorspace)
{
(void) fprintf(file," Cyan:\n");
(void) fprintf(file," Minimum: %13.02lf (%1.4f)\n",
MaxRGB*statistics.red.minimum,
statistics.red.minimum);
(void) fprintf(file," Maximum: %13.02lf (%1.4f)\n",
MaxRGB*statistics.red.maximum,
statistics.red.maximum);
(void) fprintf(file," Mean: %13.02lf (%1.4f)\n",
MaxRGB*statistics.red.mean,
statistics.red.mean);
(void) fprintf(file," Standard Deviation: %13.02lf (%1.4f)\n",
MaxRGB*statistics.red.standard_deviation,
statistics.red.standard_deviation);
(void) fprintf(file," Magenta:\n");
(void) fprintf(file," Minimum: %13.02lf (%1.4f)\n",
MaxRGB*statistics.green.minimum,
statistics.green.minimum);
(void) fprintf(file," Maximum: %13.02lf (%1.4f)\n",
MaxRGB*statistics.green.maximum,
statistics.green.maximum);
(void) fprintf(file," Mean: %13.02lf (%1.4f)\n",
MaxRGB*statistics.green.mean,
statistics.green.mean);
(void) fprintf(file," Standard Deviation: %13.02lf (%1.4f)\n",
MaxRGB*statistics.green.standard_deviation,
statistics.green.standard_deviation);
(void) fprintf(file," Yellow:\n");
(void) fprintf(file," Minimum: %13.02lf (%1.4f)\n",
MaxRGB*statistics.blue.minimum,
statistics.blue.minimum);
(void) fprintf(file," Maximum: %13.02lf (%1.4f)\n",
MaxRGB*statistics.blue.maximum,
statistics.blue.maximum);
(void) fprintf(file," Mean: %13.02lf (%1.4f)\n",
MaxRGB*statistics.blue.mean,
statistics.blue.mean);
(void) fprintf(file," Standard Deviation: %13.02lf (%1.4f)\n",
MaxRGB*statistics.blue.standard_deviation,
statistics.blue.standard_deviation);
(void) fprintf(file," Black:\n");
(void) fprintf(file," Minimum: %13.02lf (%1.4f)\n",
MaxRGB*statistics.opacity.minimum,
statistics.opacity.minimum);
(void) fprintf(file," Maximum: %13.02lf (%1.4f)\n",
MaxRGB*statistics.opacity.maximum,
statistics.opacity.maximum);
(void) fprintf(file," Mean: %13.02lf (%1.4f)\n",
MaxRGB*statistics.opacity.mean,
statistics.opacity.mean);
(void) fprintf(file," Standard Deviation: %13.02lf (%1.4f)\n",
MaxRGB*statistics.opacity.standard_deviation,
statistics.opacity.standard_deviation);
/*
if (image->matte)
(void) fprintf(file," Opacity:\n");
*/
}
else if ((IsGrayColorspace(image->colorspace)) ||
(image->is_grayscale == True))
{
(void) fprintf(file," Gray:\n");
(void) fprintf(file," Minimum: %13.02lf (%1.4f)\n",
MaxRGB*statistics.red.minimum,
statistics.red.minimum);
(void) fprintf(file," Maximum: %13.02lf (%1.4f)\n",
MaxRGB*statistics.red.maximum,
statistics.red.maximum);
(void) fprintf(file," Mean: %13.02lf (%1.4f)\n",
MaxRGB*statistics.red.mean,
statistics.red.mean);
(void) fprintf(file," Standard Deviation: %13.02lf (%1.4f)\n",
MaxRGB*statistics.red.standard_deviation,
statistics.red.standard_deviation);
if (image->matte)
{
(void) fprintf(file," Opacity:\n");
(void) fprintf(file," Minimum: %13.02lf (%1.4f)\n",
MaxRGB*statistics.opacity.minimum,
statistics.opacity.minimum);
(void) fprintf(file," Maximum: %13.02lf (%1.4f)\n",
MaxRGB*statistics.opacity.maximum,
statistics.opacity.maximum);
(void) fprintf(file," Mean: %13.02lf (%1.4f)\n",
MaxRGB*statistics.opacity.mean,
statistics.opacity.mean);
(void) fprintf(file," Standard Deviation: %13.02lf (%1.4f)\n",
MaxRGB*statistics.opacity.standard_deviation,
statistics.opacity.standard_deviation);
}
}
else
{
(void) fprintf(file," Red:\n");
(void) fprintf(file," Minimum: %13.02lf (%1.4f)\n",
MaxRGB*statistics.red.minimum,
statistics.red.minimum);
(void) fprintf(file," Maximum: %13.02lf (%1.4f)\n",
MaxRGB*statistics.red.maximum,
statistics.red.maximum);
(void) fprintf(file," Mean: %13.02lf (%1.4f)\n",
MaxRGB*statistics.red.mean,
statistics.red.mean);
(void) fprintf(file," Standard Deviation: %13.02lf (%1.4f)\n",
MaxRGB*statistics.red.standard_deviation,
statistics.red.standard_deviation);
(void) fprintf(file," Green:\n");
(void) fprintf(file," Minimum: %13.02lf (%1.4f)\n",
MaxRGB*statistics.green.minimum,
statistics.green.minimum);
(void) fprintf(file," Maximum: %13.02lf (%1.4f)\n",
MaxRGB*statistics.green.maximum,
statistics.green.maximum);
(void) fprintf(file," Mean: %13.02lf (%1.4f)\n",
MaxRGB*statistics.green.mean,
statistics.green.mean);
(void) fprintf(file," Standard Deviation: %13.02lf (%1.4f)\n",
MaxRGB*statistics.green.standard_deviation,
statistics.green.standard_deviation);
(void) fprintf(file," Blue:\n");
(void) fprintf(file," Minimum: %13.02lf (%1.4f)\n",
MaxRGB*statistics.blue.minimum,
statistics.blue.minimum);
(void) fprintf(file," Maximum: %13.02lf (%1.4f)\n",
MaxRGB*statistics.blue.maximum,
statistics.blue.maximum);
(void) fprintf(file," Mean: %13.02lf (%1.4f)\n",
MaxRGB*statistics.blue.mean,
statistics.blue.mean);
(void) fprintf(file," Standard Deviation: %13.02lf (%1.4f)\n",
MaxRGB*statistics.blue.standard_deviation,
statistics.blue.standard_deviation);
if (image->matte)
{
(void) fprintf(file," Opacity:\n");
(void) fprintf(file," Minimum: %13.02lf (%1.4f)\n",
MaxRGB*statistics.opacity.minimum,
statistics.opacity.minimum);
(void) fprintf(file," Maximum: %13.02lf (%1.4f)\n",
MaxRGB*statistics.opacity.maximum,
statistics.opacity.maximum);
(void) fprintf(file," Mean: %13.02lf (%1.4f)\n",
MaxRGB*statistics.opacity.mean,
statistics.opacity.mean);
(void) fprintf(file," Standard Deviation: %13.02lf (%1.4f)\n",
MaxRGB*statistics.opacity.standard_deviation,
statistics.opacity.standard_deviation);
}
}
}
x=0;
p=(Image *) NULL;
if ((image->matte && (strcmp(image->magick,"GIF") != 0)) || image->taint)
{
char
tuple[MaxTextExtent];
register const PixelPacket
*p;
p=(PixelPacket *) NULL;
for (y=0; y < (long) image->rows; y++)
{
p=AcquireImagePixels(image,0,y,image->columns,1,&image->exception);
if (p == (const PixelPacket *) NULL)
break;
for (x=0; x < (long) image->columns; x++)
{
if (p->opacity == TransparentOpacity)
break;
p++;
}
if (x < (long) image->columns)
break;
}
if ((x < (long) image->columns) || (y < (long) image->rows))
{
GetColorTuple(p,image->depth,image->matte,False,tuple);
(void) fprintf(file," Opacity: %.1024s\t",tuple);
GetColorTuple(p,image->depth,image->matte,True,tuple);
(void) fprintf(file," %.1024s\n",tuple);
}
}
if (image->storage_class == DirectClass)
{
if (image->total_colors != 0)
(void) fprintf(file," Colors: %lu\n",image->total_colors);
}
else
{
if (image->total_colors <= image->colors)
(void) fprintf(file," Colors: %u\n",image->colors);
else
(void) fprintf(file," Colors: %lu=>%u\n",image->total_colors,
image->colors);
}
if (image->storage_class == DirectClass)
{
if (image->total_colors < 1024)
if (verbose > 1)
(void) GetNumberColors(image,file,&image->exception);
}
else
{
char
name[MaxTextExtent];
register PixelPacket
*p;
/*
Display image colormap.
*/
p=image->colormap;
for (i=0; i < (long) image->colors; i++)
{
char
tuple[MaxTextExtent];
GetColorTuple(p,image->depth,image->matte,False,tuple);
(void) fprintf(file," %lu: %.1024s",i,tuple);
(void) fprintf(file,"\t");
(void) QueryColorname(image,p,SVGCompliance,name,&image->exception);
(void) fprintf(file," %.1024s",name);
(void) fprintf(file,"\n");
p++;
}
}
if (image->error.mean_error_per_pixel != 0.0)
(void) fprintf(file," Mean Exception Per Pixel: %ld\n",
(long) image->error.mean_error_per_pixel);
if (image->error.normalized_mean_error != 0.0)
(void) fprintf(file," Normalized Mean Exception: %g\n",
image->error.normalized_mean_error);
if (image->error.normalized_maximum_error != 0.0)
(void) fprintf(file," Normalized Maximum Exception: %gn",
image->error.normalized_maximum_error);
if (image->rendering_intent == SaturationIntent)
(void) fprintf(file," Rendering-Intent: saturation\n");
else if (image->rendering_intent == PerceptualIntent)
(void) fprintf(file," Rendering-Intent: perceptual\n");
else if (image->rendering_intent == AbsoluteIntent)
(void) fprintf(file," Rendering-Intent: absolute\n");
else if (image->rendering_intent == RelativeIntent)
(void) fprintf(file," Rendering-Intent: relative\n");
if (image->gamma != 0.0)
(void) fprintf(file," Gamma: %g\n",image->gamma);
if ((image->chromaticity.red_primary.x != 0.0) ||
(image->chromaticity.green_primary.x != 0.0) ||
(image->chromaticity.blue_primary.x != 0.0) ||
(image->chromaticity.white_point.x != 0.0))
{
/*
Display image chromaticity.
*/
(void) fprintf(file," Chromaticity:\n");
(void) fprintf(file," red primary: (%g,%g)\n",
image->chromaticity.red_primary.x,image->chromaticity.red_primary.y);
(void) fprintf(file," green primary: (%g,%g)\n",
image->chromaticity.green_primary.x,
image->chromaticity.green_primary.y);
(void) fprintf(file," blue primary: (%g,%g)\n",
image->chromaticity.blue_primary.x,image->chromaticity.blue_primary.y);
(void) fprintf(file," white point: (%g,%g)\n",
image->chromaticity.white_point.x,image->chromaticity.white_point.y);
}
if ((image->tile_info.width*image->tile_info.height) != 0)
(void) fprintf(file," Tile geometry: %lux%lu%+ld%+ld\n",
image->tile_info.width,image->tile_info.height,image->tile_info.x,
image->tile_info.y);
if ((image->x_resolution != 0.0) && (image->y_resolution != 0.0))
{
/*
Display image resolution.
*/
(void) fprintf(file," Resolution: %gx%g",image->x_resolution,
image->y_resolution);
if (image->units == UndefinedResolution)
(void) fprintf(file," pixels\n");
else if (image->units == PixelsPerInchResolution)
(void) fprintf(file," pixels/inch\n");
else if (image->units == PixelsPerCentimeterResolution)
(void) fprintf(file," pixels/centimeter\n");
else
(void) fprintf(file,"\n");
}
FormatSize(GetBlobSize(image),format);
(void) fprintf(file," Filesize: %.1024s\n",format);
fprintf(file," Interlace: %s\n",
InterlaceTypeToString(image->interlace == UndefinedInterlace ?
NoInterlace : image->interlace));
(void) fprintf(file," Orientation: %s\n", OrientationTypeToString(image->orientation));
(void) QueryColorname(image,&image->background_color,SVGCompliance,color,
&image->exception);
(void) fprintf(file," Background Color: %.1024s\n",color);
(void) QueryColorname(image,&image->border_color,SVGCompliance,color,
&image->exception);
(void) fprintf(file," Border Color: %.1024s\n",color);
(void) QueryColorname(image,&image->matte_color,SVGCompliance,color,
&image->exception);
(void) fprintf(file," Matte Color: %.1024s\n",color);
if ((image->page.width != 0) && (image->page.height != 0))
(void) fprintf(file," Page geometry: %lux%lu%+ld%+ld\n",image->page.width,
image->page.height,image->page.x,image->page.y);
(void) fprintf(file," Compose: %s\n",
CompositeOperatorToString(image->compose));
(void) fprintf(file," Dispose: ");
switch (image->dispose)
{
case UndefinedDispose:
(void) fprintf(file,"Undefined\n");
break;
case NoneDispose:
(void) fprintf(file,"None\n");
break;
case BackgroundDispose:
(void) fprintf(file,"Background\n");
break;
case PreviousDispose:
(void) fprintf(file,"Previous\n");
break;
default:
(void) fprintf(file,"\n");
break;
}
if (image->delay != 0)
(void) fprintf(file," Delay: %lu\n",image->delay);
if (image->iterations != 1)
(void) fprintf(file," Iterations: %lu\n",image->iterations);
p=image;
while (p->previous != (Image *) NULL)
p=p->previous;
for (count=1; p->next != (Image *) NULL; count++)
p=p->next;
if (count > 1)
(void) fprintf(file," Scene: %lu of %lu\n",image->scene,count);
else if (image->scene != 0)
(void) fprintf(file," Scene: %lu\n",image->scene);
(void) fprintf(file," Compression: %s\n",
CompressionTypeToString(image->compression));
/*
Display formatted image attributes. This must happen before we access
any pseudo attributes like EXIF since doing so causes real attributes
to be created and we would get duplicates in the output.
*/
attribute=GetImageAttribute(image,(char *) NULL);
{
for ( ; attribute != (const ImageAttribute *) NULL; attribute=attribute->next)
{
if (LocaleNCompare("EXIF",attribute->key,4) != 0)
{
(void) fprintf(file," %c", toupper((int)attribute->key[0]));
if (strlen(attribute->key) > 1)
(void) fprintf(file,"%.1024s",attribute->key+1);
(void) fprintf(file,": ");
(void) fprintf(file,"%s\n",attribute->value);
}
}
}
if((profile=GetImageProfile(image,"ICM",&profile_length)) != 0)
(void) fprintf(file," Profile-color: %lu bytes\n",(unsigned long)
profile_length);
if((profile=GetImageProfile(image,"IPTC",&profile_length)) != 0)
{
char
*tag,
*text;
size_t
length;
/*
Describe IPTC data.
*/
(void) fprintf(file," Profile-iptc: %lu bytes\n",(unsigned long)
profile_length);
for (i=0; i < (long) profile_length; )
{
if (profile[i] != 0x1c)
{
i++;
continue;
}
i++; /* skip file separator */
i++; /* skip record number */
switch (profile[i])
{
case 5:
tag=(char *) "Image Name";
break;
case 7:
tag=(char *) "Edit Status";
break;
case 10:
tag=(char *) "Priority";
break;
case 15:
tag=(char *) "Category";
break;
case 20:
tag=(char *) "Supplemental Category";
break;
case 22:
tag=(char *) "Fixture Identifier";
break;
case 25:
tag=(char *) "Keyword";
break;
case 30:
tag=(char *) "Release Date";
break;
case 35:
tag=(char *) "Release Time";
break;
case 40:
tag=(char *) "Special Instructions";
break;
case 45:
tag=(char *) "Reference Service";
break;
case 47:
tag=(char *) "Reference Date";
break;
case 50:
tag=(char *) "Reference Number";
break;
case 55:
tag=(char *) "Created Date";
break;
case 60:
tag=(char *) "Created Time";
break;
case 65:
tag=(char *) "Originating Program";
break;
case 70:
tag=(char *) "Program Version";
break;
case 75:
tag=(char *) "Object Cycle";
break;
case 80:
tag=(char *) "Byline";
break;
case 85:
tag=(char *) "Byline Title";
break;
case 90:
tag=(char *) "City";
break;
case 95:
tag=(char *) "Province State";
break;
case 100:
tag=(char *) "Country Code";
break;
case 101:
tag=(char *) "Country";
break;
case 103:
tag=(char *) "Original Transmission Reference";
break;
case 105:
tag=(char *) "Headline";
break;
case 110:
tag=(char *) "Credit";
break;
case 115:
tag=(char *) "Source";
break;
case 116:
tag=(char *) "Copyright String";
break;
case 120:
tag=(char *) "Caption";
break;
case 121:
tag=(char *) "Local Caption";
break;
case 122:
tag=(char *) "Caption Writer";
break;
case 200:
tag=(char *) "Custom Field 1";
break;
case 201:
tag=(char *) "Custom Field 2";
break;
case 202:
tag=(char *) "Custom Field 3";
break;
case 203:
tag=(char *) "Custom Field 4";
break;
case 204:
tag=(char *) "Custom Field 5";
break;
case 205:
tag=(char *) "Custom Field 6";
break;
case 206:
tag=(char *) "Custom Field 7";
break;
case 207:
tag=(char *) "Custom Field 8";
break;
case 208:
tag=(char *) "Custom Field 9";
break;
case 209:
tag=(char *) "Custom Field 10";
break;
case 210:
tag=(char *) "Custom Field 11";
break;
case 211:
tag=(char *) "Custom Field 12";
break;
case 212:
tag=(char *) "Custom Field 13";
break;
case 213:
tag=(char *) "Custom Field 14";
break;
case 214:
tag=(char *) "Custom Field 15";
break;
case 215:
tag=(char *) "Custom Field 16";
break;
case 216:
tag=(char *) "Custom Field 17";
break;
case 217:
tag=(char *) "Custom Field 18";
break;
case 218:
tag=(char *) "Custom Field 19";
break;
case 219:
tag=(char *) "Custom Field 20";
break;
default:
tag=(char *) "unknown";
break;
}
i++;
(void) fprintf(file," %.1024s:\n",tag);
length=profile[i++] << 8;
length|=profile[i++];
text=MagickAllocateMemory(char *,length+1);
if (text != (char *) NULL)
{
char
**textlist;
register long
j;
(void) strncpy(text,(char *) profile+i,length);
text[length]='\0';
textlist=StringToList(text);
if (textlist != (char **) NULL)
{
for (j=0; textlist[j] != (char *) NULL; j++)
{
(void) fprintf(file," %s\n",textlist[j]);
MagickFreeMemory(textlist[j]);
}
MagickFreeMemory(textlist);
}
MagickFreeMemory(text);
}
i+=length;
}
}
