diff --git a/MagickCore/compare.c b/MagickCore/compare.c
index 7d6c090cdc..92e7bffdb5 100644
--- a/MagickCore/compare.c
+++ b/MagickCore/compare.c
@@ -3674,10 +3674,6 @@ static Image *PhaseSimilarityImage(const Image *image,const Image *reconstruct,
   if (status == MagickFalse)
     ThrowPhaseSimilarityException();
   status=SIMMaximaImage(phase_image,&maxima,offset,exception);
-  if (status == MagickFalse)
-    ThrowPhaseSimilarityException();
-  status=SIMMultiplyImage(phase_image,QuantumRange,
-    (const ChannelStatistics *) NULL,exception);
   if (status == MagickFalse)
     ThrowPhaseSimilarityException();
   *similarity_metric=QuantumScale*maxima;
diff --git a/MagickWand/compare.c b/MagickWand/compare.c
index 6e3c952710..94ccc773bc 100644
--- a/MagickWand/compare.c
+++ b/MagickWand/compare.c
@@ -201,7 +201,7 @@ WandExport MagickBooleanType CompareImagesCommand(ImageInfo *image_info,
   int argc,char **argv,char **metadata,ExceptionInfo *exception)
 {
 #define CompareEpsilon  (1.0e-06)
-#define DefaultDissimilarityThreshold  0.31830988618379067154
+#define DefaultDissimilarityThreshold  (1.0/MagickPI)
 #define DefaultSimilarityThreshold  (-1.0)
 #define DestroyCompare() \
 { \