]> granicus.if.org Git - imagemagick/blobdiff - MagickCore/composite.c
(no commit message)
[imagemagick] / MagickCore / composite.c
index b597a65f8b3ebb675a5668220dedf56d4c507f63..d604805b56d59fac1bc8555ffd647bd6b901edd0 100644 (file)
 %                     MagickCore Image Composite Methods                      %
 %                                                                             %
 %                              Software Design                                %
-%                                John Cristy                                  %
+%                                   Cristy                                    %
 %                                 July 1992                                   %
 %                                                                             %
 %                                                                             %
-%  Copyright 1999-2012 ImageMagick Studio LLC, a non-profit organization      %
+%  Copyright 1999-2014 ImageMagick Studio LLC, a non-profit organization      %
 %  dedicated to making software imaging solutions freely available.           %
 %                                                                             %
 %  You may not use this file except in compliance with the License.  You may  %
@@ -45,6 +45,7 @@
 #include "MagickCore/cache.h"
 #include "MagickCore/cache-private.h"
 #include "MagickCore/cache-view.h"
+#include "MagickCore/channel.h"
 #include "MagickCore/client.h"
 #include "MagickCore/color.h"
 #include "MagickCore/color-private.h"
        operations.
 */
 
-static inline double MagickMin(const double x,const double y)
+static inline MagickRealType MagickMin(const MagickRealType x,
+  const MagickRealType y)
 {
   if (x < y)
     return(x);
   return(y);
 }
 
-static inline double MagickMax(const double x,const double y)
+static inline MagickRealType MagickMax(const MagickRealType x,
+  const MagickRealType y)
 {
   if (x > y)
     return(x);
   return(y);
 }
 
-static inline double ConvertHueToRGB(double m1,
-  double m2,double hue)
+static inline MagickRealType ConvertHueToRGB(MagickRealType m1,
+  MagickRealType m2,MagickRealType hue)
 {
   if (hue < 0.0)
     hue+=1.0;
@@ -207,25 +210,25 @@ static inline double ConvertHueToRGB(double m1,
   return(m1);
 }
 
-static void HCLComposite(const double hue,const double chroma,const double luma,
-  double *red,double *green,double *blue)
+static void HCLComposite(const MagickRealType hue,const MagickRealType chroma,
+  const MagickRealType luma,MagickRealType *red,MagickRealType *green,
+  MagickRealType *blue)
 {
-  double
+  MagickRealType
     b,
     c,
     g,
     h,
     m,
     r,
-    x,
-    z;
+    x;
 
   /*
     Convert HCL to RGB colorspace.
   */
-  assert(red != (double *) NULL);
-  assert(green != (double *) NULL);
-  assert(blue != (double *) NULL);
+  assert(red != (MagickRealType *) NULL);
+  assert(green != (MagickRealType *) NULL);
+  assert(blue != (MagickRealType *) NULL);
   h=6.0*hue;
   c=chroma;
   x=c*(1.0-fabs(fmod(h,2.0)-1.0));
@@ -268,31 +271,16 @@ static void HCLComposite(const double hue,const double chroma,const double luma,
                 b=x;
               }
   m=luma-(0.298839*r+0.586811*g+0.114350*b);
-  /*
-    Choose saturation strategy to clip it into the RGB cube; hue and luma are
-    preserved and chroma may be changed.
-  */
-  z=1.0;
-  if (m < 0.0)
-    {
-      z=luma/(luma-m);
-      m=0.0;
-    }
-  else
-    if (m+c > 1.0)
-      {
-        z=(1.0-luma)/(m+c-luma);
-        m=1.0-z*c;
-      }
-  *red=QuantumRange*(z*r+m);
-  *green=QuantumRange*(z*g+m);
-  *blue=QuantumRange*(z*b+m);
+  *red=QuantumRange*(r+m);
+  *green=QuantumRange*(g+m);
+  *blue=QuantumRange*(b+m);
 }
 
-static void CompositeHCL(const double red,const double green,const double blue,
-  double *hue,double *chroma,double *luma)
+static void CompositeHCL(const MagickRealType red,const MagickRealType green,
+  const MagickRealType blue,MagickRealType *hue,MagickRealType *chroma,
+  MagickRealType *luma)
 {
-  double
+  MagickRealType
     b,
     c,
     g,
@@ -303,20 +291,20 @@ static void CompositeHCL(const double red,const double green,const double blue,
   /*
     Convert RGB to HCL colorspace.
   */
-  assert(hue != (double *) NULL);
-  assert(chroma != (double *) NULL);
-  assert(luma != (double *) NULL);
+  assert(hue != (MagickRealType *) NULL);
+  assert(chroma != (MagickRealType *) NULL);
+  assert(luma != (MagickRealType *) NULL);
   r=red;
   g=green;
   b=blue;
   max=MagickMax(r,MagickMax(g,b));
-  c=max-(double) MagickMin(r,MagickMin(g,b));
+  c=max-(MagickRealType) MagickMin(r,MagickMin(g,b));
   h=0.0;
   if (c == 0)
     h=0.0;
   else
     if (red == max)
-      h=fmod(6.0+(g-b)/c,6.0);
+      h=fmod((g-b)/c+6.0,6.0);
     else
       if (green == max)
         h=((b-r)/c)+2.0;
@@ -356,7 +344,7 @@ static MagickBooleanType CompositeOverImage(Image *image,
   image_view=AcquireAuthenticCacheView(image,exception);
 #if defined(MAGICKCORE_OPENMP_SUPPORT)
   #pragma omp parallel for schedule(static,4) shared(progress,status) \
-    dynamic_number_threads(image,image->columns,image->rows,1)
+    magick_threads(composite_image,image,image->rows,1)
 #endif
   for (y=0; y < (ssize_t) image->rows; y++)
   {
@@ -411,10 +399,12 @@ static MagickBooleanType CompositeOverImage(Image *image,
     for (x=0; x < (ssize_t) image->columns; x++)
     {
       double
+        gamma;
+
+      MagickRealType
         alpha,
         Da,
         Dc,
-        gamma,
         Sa,
         Sc;
 
@@ -442,7 +432,7 @@ static MagickBooleanType CompositeOverImage(Image *image,
               Sc: source color.
               Dc: destination color.
           */
-          if (GetPixelMask(image,q) != 0)
+          if (GetPixelReadMask(image,q) == 0)
             {
               q+=GetPixelChannels(image);
               continue;
@@ -451,16 +441,10 @@ static MagickBooleanType CompositeOverImage(Image *image,
             source,exception);
           for (i=0; i < (ssize_t) GetPixelChannels(image); i++)
           {
-            PixelChannel
-              channel;
-
-            PixelTrait
-              composite_traits,
-              traits;
-
-            channel=GetPixelChannelChannel(image,i);
-            traits=GetPixelChannelTraits(image,channel);
-            composite_traits=GetPixelChannelTraits(composite_image,channel);
+            PixelChannel channel=GetPixelChannelChannel(image,i);
+            PixelTrait traits=GetPixelChannelTraits(image,channel);
+            PixelTrait composite_traits=GetPixelChannelTraits(composite_image,
+              channel);
             if ((traits == UndefinedPixelTrait) ||
                 (composite_traits == UndefinedPixelTrait))
               continue;
@@ -474,7 +458,7 @@ static MagickBooleanType CompositeOverImage(Image *image,
           Sa:  normalized source alpha.
           Da:  normalized destination alpha.
       */
-      if (GetPixelMask(composite_image,p) != 0)
+      if (GetPixelReadMask(composite_image,p) == 0)
         {
           p+=GetPixelChannels(composite_image);
           channels=GetPixelChannels(composite_image);
@@ -488,16 +472,10 @@ static MagickBooleanType CompositeOverImage(Image *image,
       alpha=Sa*(-Da)+Sa+Da;
       for (i=0; i < (ssize_t) GetPixelChannels(image); i++)
       {
-        PixelChannel
-          channel;
-
-        PixelTrait
-          composite_traits,
-          traits;
-
-        channel=GetPixelChannelChannel(image,i);
-        traits=GetPixelChannelTraits(image,channel);
-        composite_traits=GetPixelChannelTraits(composite_image,channel);
+        PixelChannel channel=GetPixelChannelChannel(image,i);
+        PixelTrait traits=GetPixelChannelTraits(image,channel);
+        PixelTrait composite_traits=GetPixelChannelTraits(composite_image,
+          channel);
         if ((traits == UndefinedPixelTrait) ||
             (composite_traits == UndefinedPixelTrait))
           continue;
@@ -521,9 +499,9 @@ static MagickBooleanType CompositeOverImage(Image *image,
           Sc: source color.
           Dc: destination color.
         */
-        Sc=(double) GetPixelChannel(composite_image,channel,p);
-        Dc=(double) q[i];
-        gamma=MagickEpsilonReciprocal(alpha);
+        Sc=(MagickRealType) GetPixelChannel(composite_image,channel,p);
+        Dc=(MagickRealType) q[i];
+        gamma=PerceptibleReciprocal(alpha);
         q[i]=ClampToQuantum(gamma*(Sa*Sc-Sa*Da*Dc+Da*Dc));
       }
       p+=GetPixelChannels(composite_image);
@@ -577,7 +555,7 @@ MagickExport MagickBooleanType CompositeImage(Image *image,
   MagickOffsetType
     progress;
 
-  double
+  MagickRealType
     amount,
     destination_dissolve,
     midpoint,
@@ -604,8 +582,11 @@ MagickExport MagickBooleanType CompositeImage(Image *image,
   if (composite_image == (const Image *) NULL)
     return(MagickFalse);
   if (IsGrayColorspace(image->colorspace) != MagickFalse)
-    (void) SetImageColorspace(image,RGBColorspace,exception);
+    (void) SetImageColorspace(image,sRGBColorspace,exception);
   (void) SetImageColorspace(composite_image,image->colorspace,exception);
+  if ((image->alpha_trait == BlendPixelTrait) &&
+      (composite_image->alpha_trait != BlendPixelTrait))
+    (void) SetImageAlphaChannel(composite_image,SetAlphaChannel,exception);
   if ((compose == OverCompositeOp) || (compose == SrcOverCompositeOp))
     {
       status=CompositeOverImage(image,composite_image,clip_to_self,x_offset,
@@ -635,7 +616,7 @@ MagickExport MagickBooleanType CompositeImage(Image *image,
       image_view=AcquireAuthenticCacheView(image,exception);
 #if defined(MAGICKCORE_OPENMP_SUPPORT)
       #pragma omp parallel for schedule(static,4) shared(status) \
-        dynamic_number_threads(image,image->columns,image->rows,1)
+        magick_threads(composite_image,image,composite_image->rows,1)
 #endif
       for (y=0; y < (ssize_t) composite_image->rows; y++)
       {
@@ -667,7 +648,7 @@ MagickExport MagickBooleanType CompositeImage(Image *image,
           register ssize_t
             i;
 
-          if (GetPixelMask(composite_image,p) != 0)
+          if (GetPixelReadMask(composite_image,p) == 0)
             {
               p+=GetPixelChannels(composite_image);
               q+=GetPixelChannels(image);
@@ -675,16 +656,10 @@ MagickExport MagickBooleanType CompositeImage(Image *image,
             }
           for (i=0; i < (ssize_t) GetPixelChannels(composite_image); i++)
           {
-            PixelChannel
-              channel;
-
-            PixelTrait
-              composite_traits,
-              traits;
-
-            channel=GetPixelChannelChannel(composite_image,i);
-            composite_traits=GetPixelChannelTraits(composite_image,channel);
-            traits=GetPixelChannelTraits(image,channel);
+            PixelChannel channel=GetPixelChannelChannel(composite_image,i);
+            PixelTrait composite_traits=GetPixelChannelTraits(composite_image,
+              channel);
+            PixelTrait traits=GetPixelChannelTraits(image,channel);
             if ((traits == UndefinedPixelTrait) ||
                 (composite_traits == UndefinedPixelTrait))
               continue;
@@ -736,15 +711,15 @@ MagickExport MagickBooleanType CompositeImage(Image *image,
       const char
         *value;
 
-      PixelInfo
-        pixel;
-
-      double
+      MagickRealType
         angle_range,
         angle_start,
         height,
         width;
 
+      PixelInfo
+        pixel;
+
       ResampleFilter
         *resample_filter;
 
@@ -769,13 +744,13 @@ MagickExport MagickBooleanType CompositeImage(Image *image,
       */
       SetGeometryInfo(&geometry_info);
       flags=NoValue;
-      value=GetImageArtifact(composite_image,"compose:args");
-      if (value != (char *) NULL)
+      value=GetImageArtifact(image,"compose:args");
+      if (value != (const char *) NULL)
         flags=ParseGeometry(value,&geometry_info);
-      if ((flags & WidthValue) == 0 ) {
-          (void) ThrowMagickException(exception,GetMagickModule(),
-               OptionWarning,"InvalidSetting","'%s' '%s'",
-               "compose:args",value);
+      if ((flags & WidthValue) == 0)
+        {
+          (void) ThrowMagickException(exception,GetMagickModule(),OptionWarning,
+            "InvalidSetting","'%s' '%s'","compose:args",value);
           composite_image=DestroyImage(composite_image);
           destination_image=DestroyImage(destination_image);
           return(MagickFalse);
@@ -788,8 +763,9 @@ MagickExport MagickBooleanType CompositeImage(Image *image,
       width=height=geometry_info.rho*2.0;
       if ((flags & HeightValue) != 0 )
         height=geometry_info.sigma*2.0;
-
-      /* default the unrotated ellipse width and height axis vectors */
+      /*
+        Default the unrotated ellipse width and height axis vectors.
+      */
       blur.x1=width;
       blur.x2=0.0;
       blur.y1=0.0;
@@ -797,7 +773,7 @@ MagickExport MagickBooleanType CompositeImage(Image *image,
       /* rotate vectors if a rotation angle is given */
       if ((flags & XValue) != 0 )
         {
-          double
+          MagickRealType
             angle;
 
           angle=DegreesToRadians(geometry_info.xi);
@@ -863,7 +839,7 @@ MagickExport MagickBooleanType CompositeImage(Image *image,
             }
           if (fabs(angle_range) > MagickEpsilon)
             {
-              double
+              MagickRealType
                 angle;
 
               angle=angle_start+angle_range*QuantumScale*
@@ -875,10 +851,10 @@ MagickExport MagickBooleanType CompositeImage(Image *image,
             }
 #if 0
           if ( x == 10 && y == 60 ) {
-            fprintf(stderr, "blur.x=%lf,%lf, blur.y=%lf,%lf\n",
-                blur.x1, blur.x2, blur.y1, blur.y2);
-            fprintf(stderr, "scaled by=%lf,%lf\n",
-                QuantumScale*GetPixelRed(p), QuantumScale*GetPixelGreen(p));
+            (void) fprintf(stderr, "blur.x=%lf,%lf, blur.y=%lf,%lf\n",blur.x1,
+              blur.x2,blur.y1, blur.y2);
+            (void) fprintf(stderr, "scaled by=%lf,%lf\n",QuantumScale*
+              GetPixelRed(p),QuantumScale*GetPixelGreen(p));
 #endif
           ScaleResampleFilter(resample_filter,
             blur.x1*QuantumScale*GetPixelRed(composite_image,p),
@@ -916,7 +892,7 @@ MagickExport MagickBooleanType CompositeImage(Image *image,
       PixelInfo
         pixel;
 
-      double
+      MagickRealType
         horizontal_scale,
         vertical_scale;
 
@@ -945,14 +921,14 @@ MagickExport MagickBooleanType CompositeImage(Image *image,
         {
           if ((flags & AspectValue) == 0)
             {
-              horizontal_scale=(double) (composite_image->columns-1.0)/
+              horizontal_scale=(MagickRealType) (composite_image->columns-1.0)/
                 2.0;
-              vertical_scale=(double) (composite_image->rows-1.0)/2.0;
+              vertical_scale=(MagickRealType) (composite_image->rows-1.0)/2.0;
             }
           else
             {
-              horizontal_scale=(double) (image->columns-1.0)/2.0;
-              vertical_scale=(double) (image->rows-1.0)/2.0;
+              horizontal_scale=(MagickRealType) (image->columns-1.0)/2.0;
+              vertical_scale=(MagickRealType) (image->rows-1.0)/2.0;
             }
         }
       else
@@ -983,29 +959,30 @@ MagickExport MagickBooleanType CompositeImage(Image *image,
            default = center of overlay image
            arg flag '!' = locations/percentage relative to background image
       */
-      center.x=(double) x_offset;
-      center.y=(double) y_offset;
+      center.x=(MagickRealType) x_offset;
+      center.y=(MagickRealType) y_offset;
       if (compose == DistortCompositeOp)
         {
           if ((flags & XValue) == 0)
             if ((flags & AspectValue) == 0)
-              center.x=(double) x_offset+(composite_image->columns-1)/
-                2.0;
+              center.x=(MagickRealType) (x_offset+(composite_image->columns-1)/
+                2.0);
             else
-              center.x=((double) image->columns-1)/2.0;
+              center.x=(MagickRealType) ((image->columns-1)/2);
           else
             if ((flags & AspectValue) == 0)
-              center.x=(double) x_offset+geometry_info.xi;
+              center.x=(MagickRealType) x_offset+geometry_info.xi;
             else
               center.x=geometry_info.xi;
           if ((flags & YValue) == 0)
             if ((flags & AspectValue) == 0)
-              center.y=(double) y_offset+(composite_image->rows-1)/2.0;
+              center.y=(MagickRealType) (y_offset+(composite_image->rows-1)/
+                2.0);
             else
-              center.y=((double) image->rows-1)/2.0;
+              center.y=(MagickRealType) ((image->rows-1)/2);
           else
             if ((flags & AspectValue) == 0)
-              center.y=(double) y_offset+geometry_info.psi;
+              center.y=(MagickRealType) y_offset+geometry_info.psi;
             else
               center.y=geometry_info.psi;
         }
@@ -1049,12 +1026,12 @@ MagickExport MagickBooleanType CompositeImage(Image *image,
           /*
             Displace the offset.
           */
-          offset.x=(horizontal_scale*(GetPixelRed(composite_image,p)-
-            (((double) QuantumRange+1.0)/2.0)))/(((double)
+          offset.x=(double) (horizontal_scale*(GetPixelRed(composite_image,p)-
+            (((MagickRealType) QuantumRange+1.0)/2.0)))/(((MagickRealType)
             QuantumRange+1.0)/2.0)+center.x+((compose == DisplaceCompositeOp) ?
             x : 0);
-          offset.y=(vertical_scale*(GetPixelGreen(composite_image,p)-
-            (((double) QuantumRange+1.0)/2.0)))/(((double)
+          offset.y=(double) (vertical_scale*(GetPixelGreen(composite_image,p)-
+            (((MagickRealType) QuantumRange+1.0)/2.0)))/(((MagickRealType)
             QuantumRange+1.0)/2.0)+center.y+((compose == DisplaceCompositeOp) ?
             y : 0);
           (void) InterpolatePixelInfo(image,image_view,
@@ -1063,7 +1040,7 @@ MagickExport MagickBooleanType CompositeImage(Image *image,
           /*
             Mask with the 'invalid pixel mask' in alpha channel.
           */
-          pixel.alpha=(double) QuantumRange*(1.0-(1.0-QuantumScale*
+          pixel.alpha=(MagickRealType) QuantumRange*(1.0-(1.0-QuantumScale*
             pixel.alpha)*(1.0-QuantumScale*GetPixelAlpha(composite_image,p)));
           SetPixelInfoPixel(destination_image,&pixel,q);
           p+=GetPixelChannels(composite_image);
@@ -1197,19 +1174,19 @@ MagickExport MagickBooleanType CompositeImage(Image *image,
   */
   status=MagickTrue;
   progress=0;
-  midpoint=((double) QuantumRange+1.0)/2;
+  midpoint=((MagickRealType) QuantumRange+1.0)/2;
   composite_view=AcquireVirtualCacheView(composite_image,exception);
   image_view=AcquireAuthenticCacheView(image,exception);
 #if defined(MAGICKCORE_OPENMP_SUPPORT)
   #pragma omp parallel for schedule(static,4) shared(progress,status) \
-    dynamic_number_threads(image,image->columns,image->rows,1)
+    magick_threads(composite_image,image,image->rows,1)
 #endif
   for (y=0; y < (ssize_t) image->rows; y++)
   {
     const Quantum
       *pixels;
 
-    double
+    MagickRealType
       blue,
       luma,
       green,
@@ -1271,11 +1248,13 @@ MagickExport MagickBooleanType CompositeImage(Image *image,
     for (x=0; x < (ssize_t) image->columns; x++)
     {
       double
+        gamma;
+
+      MagickRealType
         alpha,
         Da,
         Dc,
         Dca,
-        gamma,
         Sa,
         Sc,
         Sca;
@@ -1309,26 +1288,20 @@ MagickExport MagickBooleanType CompositeImage(Image *image,
           */
           (void) GetOneVirtualPixel(composite_image,x-x_offset,y-y_offset,
             source,exception);
-          if (GetPixelMask(image,q) != 0)
+          if (GetPixelReadMask(image,q) == 0)
             {
               q+=GetPixelChannels(image);
               continue;
             }
           for (i=0; i < (ssize_t) GetPixelChannels(image); i++)
           {
-            double
+            MagickRealType
               pixel;
 
-            PixelChannel
-              channel;
-
-            PixelTrait
-              composite_traits,
-              traits;
-
-            channel=GetPixelChannelChannel(image,i);
-            traits=GetPixelChannelTraits(image,channel);
-            composite_traits=GetPixelChannelTraits(composite_image,channel);
+            PixelChannel channel=GetPixelChannelChannel(image,i);
+            PixelTrait traits=GetPixelChannelTraits(image,channel);
+            PixelTrait composite_traits=GetPixelChannelTraits(composite_image,
+              channel);
             if ((traits == UndefinedPixelTrait) ||
                 (composite_traits == UndefinedPixelTrait))
               continue;
@@ -1345,9 +1318,9 @@ MagickExport MagickBooleanType CompositeImage(Image *image,
               case SrcInCompositeOp:
               case SrcOutCompositeOp:
               {
-                pixel=(double) q[i];
+                pixel=(MagickRealType) q[i];
                 if (channel == AlphaPixelChannel)
-                  pixel=(double) TransparentAlpha;
+                  pixel=(MagickRealType) TransparentAlpha;
                 break;
               }
               case ClearCompositeOp:
@@ -1357,7 +1330,7 @@ MagickExport MagickBooleanType CompositeImage(Image *image,
               {
                 if (channel == AlphaPixelChannel)
                   {
-                    pixel=(double) TransparentAlpha;
+                    pixel=(MagickRealType) TransparentAlpha;
                     break;
                   }
                 pixel=0.0;
@@ -1372,12 +1345,12 @@ MagickExport MagickBooleanType CompositeImage(Image *image,
                       source);
                     break;
                   }
-                pixel=(double) source[channel];
+                pixel=(MagickRealType) source[channel];
                 break;
               }
               default:
               {
-                pixel=(double) source[channel];
+                pixel=(MagickRealType) source[channel];
                 break;
               }
             }
@@ -1481,7 +1454,7 @@ MagickExport MagickBooleanType CompositeImage(Image *image,
           break;
         }
       }
-      if (GetPixelMask(image,p) != 0)
+      if (GetPixelReadMask(image,q) == 0)
         {
           p+=GetPixelChannels(composite_image);
           q+=GetPixelChannels(image);
@@ -1504,22 +1477,14 @@ MagickExport MagickBooleanType CompositeImage(Image *image,
       }
       for (i=0; i < (ssize_t) GetPixelChannels(image); i++)
       {
-        double
+        MagickRealType
+          pixel,
           sans;
 
-        double
-          pixel;
-
-        PixelChannel
-          channel;
-
-        PixelTrait
-          composite_traits,
-          traits;
-
-        channel=GetPixelChannelChannel(image,i);
-        traits=GetPixelChannelTraits(image,channel);
-        composite_traits=GetPixelChannelTraits(composite_image,channel);
+        PixelChannel channel=GetPixelChannelChannel(image,i);
+        PixelTrait traits=GetPixelChannelTraits(image,channel);
+        PixelTrait composite_traits=GetPixelChannelTraits(composite_image,
+          channel);
         if (traits == UndefinedPixelTrait)
           continue;
         if ((compose != IntensityCompositeOp) &&
@@ -1529,8 +1494,8 @@ MagickExport MagickBooleanType CompositeImage(Image *image,
           Sc: source color.
           Dc: destination color.
         */
-        Sc=(double) GetPixelChannel(composite_image,channel,p);
-        Dc=(double) q[i];
+        Sc=(MagickRealType) GetPixelChannel(composite_image,channel,p);
+        Dc=(MagickRealType) q[i];
         if ((traits & CopyPixelTrait) != 0)
           {
             if (channel != AlphaPixelChannel)
@@ -1571,23 +1536,23 @@ MagickExport MagickBooleanType CompositeImage(Image *image,
                 MagickBooleanType
                   equivalent;
 
-                if (Da > ((double) QuantumRange/2.0))
+                if (Da > ((MagickRealType) QuantumRange/2.0))
                   {
-                    pixel=(double) TransparentAlpha;
+                    pixel=(MagickRealType) TransparentAlpha;
                     break;
                   }
                 equivalent=IsFuzzyEquivalencePixel(composite_image,p,image,q);
                 if (equivalent != MagickFalse)
                   {
-                    pixel=(double) TransparentAlpha;
+                    pixel=(MagickRealType) TransparentAlpha;
                     break;
                   }
-                pixel=(double) OpaqueAlpha;
+                pixel=(MagickRealType) OpaqueAlpha;
                 break;
               }
               case ClearCompositeOp:
               {
-                pixel=(double) TransparentAlpha;
+                pixel=(MagickRealType) TransparentAlpha;
                 break;
               }
               case ColorizeCompositeOp:
@@ -1616,7 +1581,7 @@ MagickExport MagickBooleanType CompositeImage(Image *image,
               case CopyAlphaCompositeOp:
               {
                 pixel=QuantumRange*Sa;
-                if (composite_image->alpha_trait != BlendPixelTrait)
+                if (composite_image->alpha_trait == BlendPixelTrait)
                   pixel=GetPixelIntensity(composite_image,p);
                 break;
               }
@@ -1685,7 +1650,7 @@ MagickExport MagickBooleanType CompositeImage(Image *image,
           default:
             break;
         }
-        gamma=MagickEpsilonReciprocal(alpha);
+        gamma=PerceptibleReciprocal(alpha);
         pixel=Dc;
         switch (compose)
         {
@@ -1756,12 +1721,7 @@ MagickExport MagickBooleanType CompositeImage(Image *image,
           }
           case ColorDodgeCompositeOp:
           {
-            if ((fabs(Sca-Sa) < MagickEpsilon) && (fabs(Dca) < MagickEpsilon))
-              {
-                pixel=QuantumRange*gamma*(Sca*(1.0-Da)+Dca*(1.0-Sa));
-                break;
-              }
-            if (fabs(Sca-Sa) < MagickEpsilon)
+            if ((Sca*Da+Dca*Sa) >= Sa*Da)
               {
                 pixel=QuantumRange*gamma*(Sa*Da+Sca*(1.0-Da)+Dca*(1.0-Sa));
                 break;
@@ -1805,28 +1765,29 @@ MagickExport MagickBooleanType CompositeImage(Image *image,
           case CopyBlackCompositeOp:
           {
             if (channel == BlackPixelChannel)
-              pixel=(double) GetPixelBlack(composite_image,p);
+              pixel=(MagickRealType) (QuantumRange-
+                GetPixelBlack(composite_image,p));
             break;
           }
           case CopyBlueCompositeOp:
           case CopyYellowCompositeOp:
           {
             if (channel == BluePixelChannel)
-              pixel=(double) GetPixelBlue(composite_image,p);
+              pixel=(MagickRealType) GetPixelBlue(composite_image,p);
             break;
           }
           case CopyGreenCompositeOp:
           case CopyMagentaCompositeOp:
           {
             if (channel == GreenPixelChannel)
-              pixel=(double) GetPixelGreen(composite_image,p);
+              pixel=(MagickRealType) GetPixelGreen(composite_image,p);
             break;
           }
           case CopyRedCompositeOp:
           case CopyCyanCompositeOp:
           {
             if (channel == RedPixelChannel)
-              pixel=(double) GetPixelRed(composite_image,p);
+              pixel=(MagickRealType) GetPixelRed(composite_image,p);
             break;
           }
           case DarkenCompositeOp:
@@ -2116,16 +2077,16 @@ MagickExport MagickBooleanType CompositeImage(Image *image,
           {
             pixel=Sc+Dc;
             if (pixel > QuantumRange)
-              pixel-=(QuantumRange+1.0);
-            pixel=gamma*(pixel*Sa*Da+Sa*Sc*(1.0-Da)+Da*Dc*(1.0-Sa));
+              pixel-=QuantumRange;
+            pixel=gamma*(Sa*Da*pixel+Sa*Sc*(1.0-Da)+Da*Dc*(1.0-Sa));
             break;
           }
           case ModulusSubtractCompositeOp:
           {
             pixel=Sc-Dc;
             if (pixel < 0.0)
-              pixel+=(QuantumRange+1.0);
-            pixel=gamma*(pixel*Sa*Da+Sa*Sc*(1.0-Da)+Da*Dc*(1.0-Sa));
+              pixel+=QuantumRange;
+            pixel=gamma*(Sa*Da*pixel+Sa*Sc*(1.0-Da)+Da*Dc*(1.0-Sa));
             break;
           }
           case MultiplyCompositeOp:
@@ -2262,11 +2223,11 @@ MagickExport MagickBooleanType CompositeImage(Image *image,
           }
           case ThresholdCompositeOp:
           {
-            double
+            MagickRealType
               delta;
 
             delta=Sc-Dc;
-            if ((double) fabs((double) (2.0*delta)) < threshold)
+            if ((MagickRealType) fabs((double) (2.0*delta)) < threshold)
               {
                 pixel=gamma*Dc;
                 break;
@@ -2338,8 +2299,6 @@ MagickExport MagickBooleanType CompositeImage(Image *image,
     destination_image=DestroyImage(destination_image);
   else
     composite_image=DestroyImage(composite_image);
-  if (status != MagickFalse)
-    (void) ClampImage(image,exception);
   return(status);
 }
 \f
@@ -2403,16 +2362,13 @@ MagickExport MagickBooleanType TextureImage(Image *image,const Image *texture,
     exception);
   status=MagickTrue;
   if ((image->compose != CopyCompositeOp) &&
-      ((image->compose != OverCompositeOp) || (image->alpha_trait == BlendPixelTrait) ||
+      ((image->compose != OverCompositeOp) ||
+       (image->alpha_trait == BlendPixelTrait) ||
        (texture_image->alpha_trait == BlendPixelTrait)))
     {
       /*
         Tile texture onto the image background.
       */
-#if defined(MAGICKCORE_OPENMP_SUPPORT)
-      #pragma omp parallel for schedule(static) shared(status) \
-        dynamic_number_threads(image,image->columns,image->rows,1)
-#endif
       for (y=0; y < (ssize_t) image->rows; y+=(ssize_t) texture_image->rows)
       {
         register ssize_t
@@ -2439,9 +2395,6 @@ MagickExport MagickBooleanType TextureImage(Image *image,const Image *texture,
             MagickBooleanType
               proceed;
 
-#if defined(MAGICKCORE_OPENMP_SUPPORT)
-           #pragma omp critical (MagickCore_TextureImage)
-#endif
             proceed=SetImageProgress(image,TextureImageTag,(MagickOffsetType)
               y,image->rows);
             if (proceed == MagickFalse)
@@ -2460,8 +2413,8 @@ MagickExport MagickBooleanType TextureImage(Image *image,const Image *texture,
   texture_view=AcquireVirtualCacheView(texture_image,exception);
   image_view=AcquireAuthenticCacheView(image,exception);
 #if defined(MAGICKCORE_OPENMP_SUPPORT)
-  #pragma omp parallel for schedule(static) shared(status) \
-    dynamic_number_threads(image,image->columns,image->rows,1)
+  #pragma omp parallel for schedule(static,4) shared(status) \
+    magick_threads(texture_image,image,1,1)
 #endif
   for (y=0; y < (ssize_t) image->rows; y++)
   {
@@ -2506,7 +2459,7 @@ MagickExport MagickBooleanType TextureImage(Image *image,const Image *texture,
         register ssize_t
           i;
 
-        if (GetPixelMask(image,p) != 0)
+        if (GetPixelReadMask(image,q) == 0)
           {
             p+=GetPixelChannels(texture_image);
             q+=GetPixelChannels(image);
@@ -2514,16 +2467,10 @@ MagickExport MagickBooleanType TextureImage(Image *image,const Image *texture,
           }
         for (i=0; i < (ssize_t) GetPixelChannels(texture_image); i++)
         {
-          PixelChannel
-            channel;
-
-          PixelTrait
-            texture_traits,
-            traits;
-
-          channel=GetPixelChannelChannel(texture_image,i);
-          texture_traits=GetPixelChannelTraits(texture_image,channel);
-          traits=GetPixelChannelTraits(image,channel);
+          PixelChannel channel=GetPixelChannelChannel(texture_image,i);
+          PixelTrait traits=GetPixelChannelTraits(image,channel);
+          PixelTrait texture_traits=GetPixelChannelTraits(texture_image,
+            channel);
           if ((traits == UndefinedPixelTrait) ||
               (texture_traits == UndefinedPixelTrait))
             continue;
@@ -2541,9 +2488,6 @@ MagickExport MagickBooleanType TextureImage(Image *image,const Image *texture,
         MagickBooleanType
           proceed;
 
-#if defined(MAGICKCORE_OPENMP_SUPPORT)
-        #pragma omp critical (MagickCore_TextureImage)
-#endif
         proceed=SetImageProgress(image,TextureImageTag,(MagickOffsetType) y,
           image->rows);
         if (proceed == MagickFalse)