text: /*
  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  %                                                                             %
  %                                                                             %
  %                                                                             %
  %                        DDDD   RRRR    AAA   W   W                           %
  %                        D   D  R   R  A   A  W   W                           %
  %                        D   D  RRRR   AAAAA  W W W                           %
  %                        D   D  R RN   A   A  WW WW                           %
  %                        DDDD   R  R   A   A  W   W                           %
  %                                                                             %
  %                                                                             %
  %                     MagickCore Image Drawing Methods                        %
  %                                                                             %
  %                                                                             %
  %                              Software Design                                %
  %                                   Cristy                                    %
  %                                 July 1998                                   %
  %                                                                             %
  %                                                                             %
  %  Copyright 1999-2020 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  %
  %  obtain a copy of the License at                                            %
  %                                                                             %
  %    https://imagemagick.org/script/license.php                               %
  %                                                                             %
  %  Unless required by applicable law or agreed to in writing, software        %
  %  distributed under the License is distributed on an "AS IS" BASIS,          %
  %  WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.   %
  %  See the License for the specific language governing permissions and        %
  %  limitations under the License.                                             %
  %                                                                             %
  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  %
  % Bill Radcliffe of Corbis (www.corbis.com) contributed the polygon
  % rendering code based on Paul Heckbert's "Concave Polygon Scan Conversion",
  % Graphics Gems, 1990.  Leonard Rosenthal and David Harr of Appligent
  % (www.appligent.com) contributed the dash pattern, linecap stroking
  % algorithm, and minor rendering improvements.
  %
  */
  
  /*
    Include declarations.
  */
  #include "MagickCore/studio.h"
  #include "MagickCore/annotate.h"
  #include "MagickCore/artifact.h"
  #include "MagickCore/blob.h"
  #include "MagickCore/cache.h"
  #include "MagickCore/cache-private.h"
  #include "MagickCore/cache-view.h"
  #include "MagickCore/channel.h"
  #include "MagickCore/color.h"
  #include "MagickCore/colorspace-private.h"
  #include "MagickCore/composite.h"
  #include "MagickCore/composite-private.h"
  #include "MagickCore/constitute.h"
  #include "MagickCore/draw.h"
  #include "MagickCore/draw-private.h"
  #include "MagickCore/enhance.h"
  #include "MagickCore/exception.h"
  #include "MagickCore/exception-private.h"
  #include "MagickCore/gem.h"
  #include "MagickCore/geometry.h"
  #include "MagickCore/image-private.h"
  #include "MagickCore/list.h"
  #include "MagickCore/log.h"
  #include "MagickCore/memory-private.h"
  #include "MagickCore/monitor.h"
  #include "MagickCore/monitor-private.h"
  #include "MagickCore/option.h"
  #include "MagickCore/paint.h"
  #include "MagickCore/pixel-accessor.h"
  #include "MagickCore/pixel-private.h"
  #include "MagickCore/property.h"
  #include "MagickCore/resample.h"
  #include "MagickCore/resample-private.h"
  #include "MagickCore/resource_.h"
  #include "MagickCore/splay-tree.h"
  #include "MagickCore/string_.h"
  #include "MagickCore/string-private.h"
  #include "MagickCore/thread-private.h"
  #include "MagickCore/token.h"
  #include "MagickCore/transform-private.h"
  #include "MagickCore/utility.h"
  
  /*
    Define declarations.
  */
  #define BezierQuantum  200
  #define PrimitiveExtentPad  2048
  #define MaxBezierCoordinates  67108864
  #define ThrowPointExpectedException(token,exception) \
  { \
    (void) ThrowMagickException(exception,GetMagickModule(),DrawError, \
      "NonconformingDrawingPrimitiveDefinition","`%s'",token); \
    status=MagickFalse; \
    break; \
  }
  
  /*
    Typedef declarations.
  */
  typedef struct _EdgeInfo
  {
    SegmentInfo
      bounds;
  
    double
      scanline;
  
    PointInfo
      *points;
  
    size_t
      number_points;
  
    ssize_t
      direction;
  
    MagickBooleanType
      ghostline;
  
    size_t
      highwater;
  } EdgeInfo;
  
  typedef struct _ElementInfo
  {
    double
      cx,
      cy,
      major,
      minor,
      angle;
  } ElementInfo;
  
  typedef struct _MVGInfo
  {
    PrimitiveInfo
      **primitive_info;
  
    size_t
      *extent;
  
    ssize_t
      offset;
  
    PointInfo
      point;
  
    ExceptionInfo
      *exception;
  } MVGInfo;
  
  typedef struct _PolygonInfo
  {
    EdgeInfo
      *edges;
  
    size_t
      number_edges;
  } PolygonInfo;
  
  typedef enum
  {
    MoveToCode,
    OpenCode,
    GhostlineCode,
    LineToCode,
    EndCode
  } PathInfoCode;
  
  typedef struct _PathInfo
  {
    PointInfo
      point;
  
    PathInfoCode
      code;
  } PathInfo;
  
  /*
    Forward declarations.
  */
  static Image
    *DrawClippingMask(Image *,const DrawInfo *,const char *,const char *,
      ExceptionInfo *);
  
  static MagickBooleanType
    DrawStrokePolygon(Image *,const DrawInfo *,const PrimitiveInfo *,
      ExceptionInfo *),
    RenderMVGContent(Image *,const DrawInfo *,const size_t,ExceptionInfo *),
    TraceArc(MVGInfo *,const PointInfo,const PointInfo,const PointInfo),
    TraceArcPath(MVGInfo *,const PointInfo,const PointInfo,const PointInfo,
      const double,const MagickBooleanType,const MagickBooleanType),
    TraceBezier(MVGInfo *,const size_t),
    TraceCircle(MVGInfo *,const PointInfo,const PointInfo),
    TraceEllipse(MVGInfo *,const PointInfo,const PointInfo,const PointInfo),
    TraceLine(PrimitiveInfo *,const PointInfo,const PointInfo),
    TraceRectangle(PrimitiveInfo *,const PointInfo,const PointInfo),
    TraceRoundRectangle(MVGInfo *,const PointInfo,const PointInfo,PointInfo),
    TraceSquareLinecap(PrimitiveInfo *,const size_t,const double);
  
  static PrimitiveInfo
    *TraceStrokePolygon(const Image *,const DrawInfo *,const PrimitiveInfo *);
  
  static ssize_t
    TracePath(MVGInfo *,const char *,ExceptionInfo *);
  
  /*
  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  %                                                                             %
  %                                                                             %
  %                                                                             %
  %   A c q u i r e D r a w I n f o                                             %
  %                                                                             %
  %                                                                             %
  %                                                                             %
  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  %
  %  AcquireDrawInfo() returns a DrawInfo structure properly initialized.
  %
  %  The format of the AcquireDrawInfo method is:
  %
  %      DrawInfo *AcquireDrawInfo(void)
  %
  */
  MagickExport DrawInfo *AcquireDrawInfo(void)
  {
    DrawInfo
      *draw_info;
  
    draw_info=(DrawInfo *) AcquireCriticalMemory(sizeof(*draw_info));
    GetDrawInfo((ImageInfo *) NULL,draw_info);
    return(draw_info);
  }
  
  /*
  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  %                                                                             %
  %                                                                             %
  %                                                                             %
  %   C l o n e D r a w I n f o                                                 %
  %                                                                             %
  %                                                                             %
  %                                                                             %
  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  %
  %  CloneDrawInfo() makes a copy of the given draw_info structure.  If NULL
  %  is specified, a new DrawInfo structure is created initialized to default
  %  values.
  %
  %  The format of the CloneDrawInfo method is:
  %
  %      DrawInfo *CloneDrawInfo(const ImageInfo *image_info,
  %        const DrawInfo *draw_info)
  %
  %  A description of each parameter follows:
  %
  %    o image_info: the image info.
  %
  %    o draw_info: the draw info.
  %
  */
  MagickExport DrawInfo *CloneDrawInfo(const ImageInfo *image_info,
    const DrawInfo *draw_info)
  {
    DrawInfo
      *clone_info;
  
    ExceptionInfo
      *exception;
  
    clone_info=(DrawInfo *) AcquireCriticalMemory(sizeof(*clone_info));
    GetDrawInfo(image_info,clone_info);
    if (draw_info == (DrawInfo *) NULL)
      return(clone_info);
    exception=AcquireExceptionInfo();
    if (draw_info->id != (char *) NULL)
      (void) CloneString(&clone_info->id,draw_info->id);
    if (draw_info->primitive != (char *) NULL)
      (void) CloneString(&clone_info->primitive,draw_info->primitive);
    if (draw_info->geometry != (char *) NULL)
      (void) CloneString(&clone_info->geometry,draw_info->geometry);
    clone_info->compliance=draw_info->compliance;
    clone_info->viewbox=draw_info->viewbox;
    clone_info->affine=draw_info->affine;
    clone_info->gravity=draw_info->gravity;
    clone_info->fill=draw_info->fill;
    clone_info->stroke=draw_info->stroke;
    clone_info->stroke_width=draw_info->stroke_width;
    if (draw_info->fill_pattern != (Image *) NULL)
      clone_info->fill_pattern=CloneImage(draw_info->fill_pattern,0,0,MagickTrue,
        exception);
    if (draw_info->stroke_pattern != (Image *) NULL)
      clone_info->stroke_pattern=CloneImage(draw_info->stroke_pattern,0,0,
        MagickTrue,exception);
    clone_info->stroke_antialias=draw_info->stroke_antialias;
    clone_info->text_antialias=draw_info->text_antialias;
    clone_info->fill_rule=draw_info->fill_rule;
    clone_info->linecap=draw_info->linecap;
    clone_info->linejoin=draw_info->linejoin;
    clone_info->miterlimit=draw_info->miterlimit;
    clone_info->dash_offset=draw_info->dash_offset;
    clone_info->decorate=draw_info->decorate;
    clone_info->compose=draw_info->compose;
    if (draw_info->text != (char *) NULL)
      (void) CloneString(&clone_info->text,draw_info->text);
    if (draw_info->font != (char *) NULL)
      (void) CloneString(&clone_info->font,draw_info->font);
    if (draw_info->metrics != (char *) NULL)
      (void) CloneString(&clone_info->metrics,draw_info->metrics);
    if (draw_info->family != (char *) NULL)
      (void) CloneString(&clone_info->family,draw_info->family);
    clone_info->style=draw_info->style;
    clone_info->stretch=draw_info->stretch;
    clone_info->weight=draw_info->weight;
    if (draw_info->encoding != (char *) NULL)
      (void) CloneString(&clone_info->encoding,draw_info->encoding);
    clone_info->pointsize=draw_info->pointsize;
    clone_info->kerning=draw_info->kerning;
    clone_info->interline_spacing=draw_info->interline_spacing;
    clone_info->interword_spacing=draw_info->interword_spacing;
    clone_info->direction=draw_info->direction;
    if (draw_info->density != (char *) NULL)
      (void) CloneString(&clone_info->density,draw_info->density);
    clone_info->align=draw_info->align;
    clone_info->undercolor=draw_info->undercolor;
    clone_info->border_color=draw_info->border_color;
    if (draw_info->server_name != (char *) NULL)
      (void) CloneString(&clone_info->server_name,draw_info->server_name);
    if (draw_info->dash_pattern != (double *) NULL)
      {
        register ssize_t
          x;
  
        for (x=0; fabs(draw_info->dash_pattern[x]) >= MagickEpsilon; x++) ;
        clone_info->dash_pattern=(double *) AcquireQuantumMemory((size_t) (2*x+2),
          sizeof(*clone_info->dash_pattern));
        if (clone_info->dash_pattern == (double *) NULL)
          ThrowFatalException(ResourceLimitFatalError,
            "UnableToAllocateDashPattern");
        (void) memset(clone_info->dash_pattern,0,(size_t) (2*x+2)*
          sizeof(*clone_info->dash_pattern));
        (void) memcpy(clone_info->dash_pattern,draw_info->dash_pattern,(size_t)
          (x+1)*sizeof(*clone_info->dash_pattern));
      }
    clone_info->gradient=draw_info->gradient;
    if (draw_info->gradient.stops != (StopInfo *) NULL)
      {
        size_t
          number_stops;
  
        number_stops=clone_info->gradient.number_stops;
        clone_info->gradient.stops=(StopInfo *) AcquireQuantumMemory((size_t)
          number_stops,sizeof(*clone_info->gradient.stops));
        if (clone_info->gradient.stops == (StopInfo *) NULL)
          ThrowFatalException(ResourceLimitFatalError,
            "UnableToAllocateDashPattern");
        (void) memcpy(clone_info->gradient.stops,draw_info->gradient.stops,
          (size_t) number_stops*sizeof(*clone_info->gradient.stops));
      }
    clone_info->bounds=draw_info->bounds;
    clone_info->fill_alpha=draw_info->fill_alpha;
    clone_info->stroke_alpha=draw_info->stroke_alpha;
    clone_info->element_reference=draw_info->element_reference;
    clone_info->clip_path=draw_info->clip_path;
    clone_info->clip_units=draw_info->clip_units;
    if (draw_info->clip_mask != (char *) NULL)
      (void) CloneString(&clone_info->clip_mask,draw_info->clip_mask);
    if (draw_info->clipping_mask != (Image *) NULL)
      clone_info->clipping_mask=CloneImage(draw_info->clipping_mask,0,0,
        MagickTrue,exception);
    if (draw_info->composite_mask != (Image *) NULL)
      clone_info->composite_mask=CloneImage(draw_info->composite_mask,0,0,
        MagickTrue,exception);
    clone_info->render=draw_info->render;
    clone_info->debug=IsEventLogging();
    exception=DestroyExceptionInfo(exception);
    return(clone_info);
  }
  
  /*
  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  %                                                                             %
  %                                                                             %
  %                                                                             %
  +   C o n v e r t P a t h T o P o l y g o n                                   %
  %                                                                             %
  %                                                                             %
  %                                                                             %
  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  %
  %  ConvertPathToPolygon() converts a path to the more efficient sorted
  %  rendering form.
  %
  %  The format of the ConvertPathToPolygon method is:
  %
  %      PolygonInfo *ConvertPathToPolygon(const PathInfo *path_info)
  %
  %  A description of each parameter follows:
  %
  %    o Method ConvertPathToPolygon returns the path in a more efficient sorted
  %      rendering form of type PolygonInfo.
  %
  %    o draw_info: Specifies a pointer to an DrawInfo structure.
  %
  %    o path_info: Specifies a pointer to an PathInfo structure.
  %
  %
  */
  
  #if defined(__cplusplus) || defined(c_plusplus)
  extern "C" {
  #endif
  
  static int DrawCompareEdges(const void *p_edge,const void *q_edge)
  {
  #define DrawCompareEdge(p,q) \
  { \
    if (((p)-(q)) < 0.0) \
      return(-1); \
    if (((p)-(q)) > 0.0) \
      return(1); \
  }
  
    register const PointInfo
      *p,
      *q;
  
    /*
      Edge sorting for right-handed coordinate system.
    */
    p=((const EdgeInfo *) p_edge)->points;
    q=((const EdgeInfo *) q_edge)->points;
    DrawCompareEdge(p[0].y,q[0].y);
    DrawCompareEdge(p[0].x,q[0].x);
    DrawCompareEdge((p[1].x-p[0].x)*(q[1].y-q[0].y),(p[1].y-p[0].y)*
      (q[1].x-q[0].x));
    DrawCompareEdge(p[1].y,q[1].y);
    DrawCompareEdge(p[1].x,q[1].x);
    return(0);
  }
  
  #if defined(__cplusplus) || defined(c_plusplus)
  }
  #endif
  
  static void LogPolygonInfo(const PolygonInfo *polygon_info)
  {
    register EdgeInfo
      *p;
  
    register ssize_t
      i,
      j;
  
    (void) LogMagickEvent(DrawEvent,GetMagickModule(),"    begin active-edge");
    p=polygon_info->edges;
    for (i=0; i < (ssize_t) polygon_info->number_edges; i++)
    {
      (void) LogMagickEvent(DrawEvent,GetMagickModule(),"      edge %.20g:",
        (double) i);
      (void) LogMagickEvent(DrawEvent,GetMagickModule(),"      direction: %s",
        p->direction != MagickFalse ? "down" : "up");
      (void) LogMagickEvent(DrawEvent,GetMagickModule(),"      ghostline: %s",
        p->ghostline != MagickFalse ? "transparent" : "opaque");
      (void) LogMagickEvent(DrawEvent,GetMagickModule(),
        "      bounds: %g,%g - %g,%g",p->bounds.x1,p->bounds.y1,
        p->bounds.x2,p->bounds.y2);
      for (j=0; j < (ssize_t) p->number_points; j++)
        (void) LogMagickEvent(DrawEvent,GetMagickModule(),"        %g,%g",
          p->points[j].x,p->points[j].y);
      p++;
    }
    (void) LogMagickEvent(DrawEvent,GetMagickModule(),"    end active-edge");
  }
  
  static void ReversePoints(PointInfo *points,const size_t number_points)
  {
    PointInfo
      point;
  
    register ssize_t
      i;
  
    for (i=0; i < (ssize_t) (number_points >> 1); i++)
    {
      point=points[i];
      points[i]=points[number_points-(i+1)];
      points[number_points-(i+1)]=point;
    }
  }
  
  static PolygonInfo *ConvertPathToPolygon(const PathInfo *path_info)
  {
    long
      direction,
      next_direction;
  
    PointInfo
      point,
      *points;
  
    PolygonInfo
      *polygon_info;
  
    SegmentInfo
      bounds;
  
    register ssize_t
      i,
      n;
  
    MagickBooleanType
      ghostline;
  
    size_t
      edge,
      number_edges,
      number_points;
  
    /*
      Convert a path to the more efficient sorted rendering form.
    */
    polygon_info=(PolygonInfo *) AcquireMagickMemory(sizeof(*polygon_info));
    if (polygon_info == (PolygonInfo *) NULL)
      return((PolygonInfo *) NULL);
    number_edges=16;
    polygon_info->edges=(EdgeInfo *) AcquireQuantumMemory(number_edges,
      sizeof(*polygon_info->edges));
    if (polygon_info->edges == (EdgeInfo *) NULL)
      return((PolygonInfo *) NULL);
    (void) memset(polygon_info->edges,0,number_edges*
      sizeof(*polygon_info->edges));
    direction=0;
    edge=0;
    ghostline=MagickFalse;
    n=0;
    number_points=0;
    points=(PointInfo *) NULL;
    (void) memset(&point,0,sizeof(point));
    (void) memset(&bounds,0,sizeof(bounds));
    polygon_info->edges[edge].number_points=(size_t) n;
    polygon_info->edges[edge].scanline=0.0;
    polygon_info->edges[edge].highwater=0;
    polygon_info->edges[edge].ghostline=ghostline;
    polygon_info->edges[edge].direction=(ssize_t) direction;
    polygon_info->edges[edge].points=points;
    polygon_info->edges[edge].bounds=bounds;
    polygon_info->number_edges=0;
    for (i=0; path_info[i].code != EndCode; i++)
    {
      if ((path_info[i].code == MoveToCode) || (path_info[i].code == OpenCode) ||
          (path_info[i].code == GhostlineCode))
        {
          /*
            Move to.
          */
          if ((points != (PointInfo *) NULL) && (n >= 2))
            {
              if (edge == number_edges)
                {
                  number_edges<<=1;
                  polygon_info->edges=(EdgeInfo *) ResizeQuantumMemory(
                    polygon_info->edges,(size_t) number_edges,
                    sizeof(*polygon_info->edges));
                  if (polygon_info->edges == (EdgeInfo *) NULL)
                    return((PolygonInfo *) NULL);
                }
              polygon_info->edges[edge].number_points=(size_t) n;
              polygon_info->edges[edge].scanline=(-1.0);
              polygon_info->edges[edge].highwater=0;
              polygon_info->edges[edge].ghostline=ghostline;
              polygon_info->edges[edge].direction=(ssize_t) (direction > 0);
              if (direction < 0)
                ReversePoints(points,(size_t) n);
              polygon_info->edges[edge].points=points;
              polygon_info->edges[edge].bounds=bounds;
              polygon_info->edges[edge].bounds.y1=points[0].y;
              polygon_info->edges[edge].bounds.y2=points[n-1].y;
              points=(PointInfo *) NULL;
              ghostline=MagickFalse;
              edge++;
            }
          if (points == (PointInfo *) NULL)
            {
              number_points=16;
              points=(PointInfo *) AcquireQuantumMemory((size_t) number_points,
                sizeof(*points));
              if (points == (PointInfo *) NULL)
                return((PolygonInfo *) NULL);
            }
          ghostline=path_info[i].code == GhostlineCode ? MagickTrue : MagickFalse;
          point=path_info[i].point;
          points[0]=point;
          bounds.x1=point.x;
          bounds.x2=point.x;
          direction=0;
          n=1;
          continue;
        }
      /*
        Line to.
      */
      next_direction=((path_info[i].point.y > point.y) ||
        ((fabs(path_info[i].point.y-point.y) < MagickEpsilon) &&
         (path_info[i].point.x > point.x))) ? 1 : -1;
      if ((points != (PointInfo *) NULL) && (direction != 0) &&
          (direction != next_direction))
        {
          /*
            New edge.
          */
          point=points[n-1];
          if (edge == number_edges)
            {
              number_edges<<=1;
              polygon_info->edges=(EdgeInfo *) ResizeQuantumMemory(
                polygon_info->edges,(size_t) number_edges,
                sizeof(*polygon_info->edges));
              if (polygon_info->edges == (EdgeInfo *) NULL)
                return((PolygonInfo *) NULL);
            }
          polygon_info->edges[edge].number_points=(size_t) n;
          polygon_info->edges[edge].scanline=(-1.0);
          polygon_info->edges[edge].highwater=0;
          polygon_info->edges[edge].ghostline=ghostline;
          polygon_info->edges[edge].direction=(ssize_t) (direction > 0);
          if (direction < 0)
            ReversePoints(points,(size_t) n);
          polygon_info->edges[edge].points=points;
          polygon_info->edges[edge].bounds=bounds;
          polygon_info->edges[edge].bounds.y1=points[0].y;
          polygon_info->edges[edge].bounds.y2=points[n-1].y;
          number_points=16;
          points=(PointInfo *) AcquireQuantumMemory((size_t) number_points,
            sizeof(*points));
          if (points == (PointInfo *) NULL)
            return((PolygonInfo *) NULL);
          n=1;
          ghostline=MagickFalse;
          points[0]=point;
          bounds.x1=point.x;
          bounds.x2=point.x;
          edge++;
        }
      direction=next_direction;
      if (points == (PointInfo *) NULL)
        continue;
      if (n == (ssize_t) number_points)
        {
          number_points<<=1;
          points=(PointInfo *) ResizeQuantumMemory(points,(size_t) number_points,
            sizeof(*points));
          if (points == (PointInfo *) NULL)
            return((PolygonInfo *) NULL);
        }
      point=path_info[i].point;
      points[n]=point;
      if (point.x < bounds.x1)
        bounds.x1=point.x;
      if (point.x > bounds.x2)
        bounds.x2=point.x;
      n++;
    }
    if (points != (PointInfo *) NULL)
      {
        if (n < 2)
          points=(PointInfo *) RelinquishMagickMemory(points);
        else
          {
            if (edge == number_edges)
              {
                number_edges<<=1;
                polygon_info->edges=(EdgeInfo *) ResizeQuantumMemory(
                  polygon_info->edges,(size_t) number_edges,
                  sizeof(*polygon_info->edges));
                if (polygon_info->edges == (EdgeInfo *) NULL)
                  return((PolygonInfo *) NULL);
              }
            polygon_info->edges[edge].number_points=(size_t) n;
            polygon_info->edges[edge].scanline=(-1.0);
            polygon_info->edges[edge].highwater=0;
            polygon_info->edges[edge].ghostline=ghostline;
            polygon_info->edges[edge].direction=(ssize_t) (direction > 0);
            if (direction < 0)
              ReversePoints(points,(size_t) n);
            polygon_info->edges[edge].points=points;
            polygon_info->edges[edge].bounds=bounds;
            polygon_info->edges[edge].bounds.y1=points[0].y;
            polygon_info->edges[edge].bounds.y2=points[n-1].y;
            ghostline=MagickFalse;
            edge++;
          }
      }
    polygon_info->number_edges=edge;
    qsort(polygon_info->edges,(size_t) polygon_info->number_edges,
      sizeof(*polygon_info->edges),DrawCompareEdges);
    if (IsEventLogging() != MagickFalse)
      LogPolygonInfo(polygon_info);
    return(polygon_info);
  }
  
  /*
  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  %                                                                             %
  %                                                                             %
  %                                                                             %
  +   C o n v e r t P r i m i t i v e T o P a t h                               %
  %                                                                             %
  %                                                                             %
  %                                                                             %
  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  %
  %  ConvertPrimitiveToPath() converts a PrimitiveInfo structure into a vector
  %  path structure.
  %
  %  The format of the ConvertPrimitiveToPath method is:
  %
  %      PathInfo *ConvertPrimitiveToPath(const DrawInfo *draw_info,
  %        const PrimitiveInfo *primitive_info)
  %
  %  A description of each parameter follows:
  %
  %    o Method ConvertPrimitiveToPath returns a vector path structure of type
  %      PathInfo.
  %
  %    o draw_info: a structure of type DrawInfo.
  %
  %    o primitive_info: Specifies a pointer to an PrimitiveInfo structure.
  %
  %
  */
  
  static void LogPathInfo(const PathInfo *path_info)
  {
    register const PathInfo
      *p;
  
    (void) LogMagickEvent(DrawEvent,GetMagickModule(),"    begin vector-path");
    for (p=path_info; p->code != EndCode; p++)
      (void) LogMagickEvent(DrawEvent,GetMagickModule(),
        "      %g,%g %s",p->point.x,p->point.y,p->code == GhostlineCode ?
        "moveto ghostline" : p->code == OpenCode ? "moveto open" :
        p->code == MoveToCode ? "moveto" : p->code == LineToCode ? "lineto" :
        "?");
    (void) LogMagickEvent(DrawEvent,GetMagickModule(),"    end vector-path");
  }
  
  static PathInfo *ConvertPrimitiveToPath(const PrimitiveInfo *primitive_info)
  {
    MagickBooleanType
      closed_subpath;
  
    PathInfo
      *path_info;
  
    PathInfoCode
      code;
  
    PointInfo
      p,
      q;
  
    register ssize_t
      i,
      n;
  
    ssize_t
      coordinates,
      start;
  
    /*
      Converts a PrimitiveInfo structure into a vector path structure.
    */
    switch (primitive_info->primitive)
    {
      case AlphaPrimitive:
      case ColorPrimitive:
      case ImagePrimitive:
      case PointPrimitive:
      case TextPrimitive:
        return((PathInfo *) NULL);
      default:
        break;
    }
    for (i=0; primitive_info[i].primitive != UndefinedPrimitive; i++) ;
    path_info=(PathInfo *) AcquireQuantumMemory((size_t) (3UL*i+1UL),
      sizeof(*path_info));
    if (path_info == (PathInfo *) NULL)
      return((PathInfo *) NULL);
    coordinates=0;
    closed_subpath=MagickFalse;
    n=0;
    p.x=(-1.0);
    p.y=(-1.0);
    q.x=(-1.0);
    q.y=(-1.0);
    start=0;
    for (i=0; primitive_info[i].primitive != UndefinedPrimitive; i++)
    {
      code=LineToCode;
      if (coordinates <= 0)
        {
          /*
            New subpath.
          */
          coordinates=(ssize_t) primitive_info[i].coordinates;
          p=primitive_info[i].point;
          start=n;
          code=MoveToCode;
          closed_subpath=primitive_info[i].closed_subpath;
        }
      coordinates--;
      if ((code == MoveToCode) || (coordinates <= 0) ||
          (fabs(q.x-primitive_info[i].point.x) >= MagickEpsilon) ||
          (fabs(q.y-primitive_info[i].point.y) >= MagickEpsilon))
        {
          /*
            Eliminate duplicate points.
          */
          path_info[n].code=code;
          path_info[n].point=primitive_info[i].point;
          q=primitive_info[i].point;
          n++;
        }
      if (coordinates > 0)
        continue;  /* next point in current subpath */
      if (closed_subpath != MagickFalse)
        {
          closed_subpath=MagickFalse;
          continue;
        }
      /*
        Mark the p point as open if the subpath is not closed.
      */
      path_info[start].code=OpenCode;
      path_info[n].code=GhostlineCode;
      path_info[n].point=primitive_info[i].point;
      n++;
      path_info[n].code=LineToCode;
      path_info[n].point=p;
      n++;
    }
    path_info[n].code=EndCode;
    path_info[n].point.x=0.0;
    path_info[n].point.y=0.0;
    if (IsEventLogging() != MagickFalse)
      LogPathInfo(path_info);
    path_info=(PathInfo *) ResizeQuantumMemory(path_info,(size_t) (n+1),
      sizeof(*path_info));
    return(path_info);
  }
  
  /*
  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  %                                                                             %
  %                                                                             %
  %                                                                             %
  %   D e s t r o y D r a w I n f o                                             %
  %                                                                             %
  %                                                                             %
  %                                                                             %
  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  %
  %  DestroyDrawInfo() deallocates memory associated with an DrawInfo structure.
  %
  %  The format of the DestroyDrawInfo method is:
  %
  %      DrawInfo *DestroyDrawInfo(DrawInfo *draw_info)
  %
  %  A description of each parameter follows:
  %
  %    o draw_info: the draw info.
  %
  */
  MagickExport DrawInfo *DestroyDrawInfo(DrawInfo *draw_info)
  {
    assert(draw_info != (DrawInfo *) NULL);
    if (draw_info->debug != MagickFalse)
      (void) LogMagickEvent(TraceEvent,GetMagickModule(),"...");
    assert(draw_info->signature == MagickCoreSignature);
    if (draw_info->id != (char *) NULL)
      draw_info->id=DestroyString(draw_info->id);
    if (draw_info->primitive != (char *) NULL)
      draw_info->primitive=DestroyString(draw_info->primitive);
    if (draw_info->text != (char *) NULL)
      draw_info->text=DestroyString(draw_info->text);
    if (draw_info->geometry != (char *) NULL)
      draw_info->geometry=DestroyString(draw_info->geometry);
    if (draw_info->fill_pattern != (Image *) NULL)
      draw_info->fill_pattern=DestroyImage(draw_info->fill_pattern);
    if (draw_info->stroke_pattern != (Image *) NULL)
      draw_info->stroke_pattern=DestroyImage(draw_info->stroke_pattern);
    if (draw_info->font != (char *) NULL)
      draw_info->font=DestroyString(draw_info->font);
    if (draw_info->metrics != (char *) NULL)
      draw_info->metrics=DestroyString(draw_info->metrics);
    if (draw_info->family != (char *) NULL)
      draw_info->family=DestroyString(draw_info->family);
    if (draw_info->encoding != (char *) NULL)
      draw_info->encoding=DestroyString(draw_info->encoding);
    if (draw_info->density != (char *) NULL)
      draw_info->density=DestroyString(draw_info->density);
    if (draw_info->server_name != (char *) NULL)
      draw_info->server_name=(char *)
       RelinquishMagickMemory(draw_info->server_name);
    if (draw_info->dash_pattern != (double *) NULL)
      draw_info->dash_pattern=(double *) RelinquishMagickMemory(
        draw_info->dash_pattern);
    if (draw_info->gradient.stops != (StopInfo *) NULL)
      draw_info->gradient.stops=(StopInfo *) RelinquishMagickMemory(
        draw_info->gradient.stops);
    if (draw_info->clip_mask != (char *) NULL)
      draw_info->clip_mask=DestroyString(draw_info->clip_mask);
    if (draw_info->clipping_mask != (Image *) NULL)
      draw_info->clipping_mask=DestroyImage(draw_info->clipping_mask);
    if (draw_info->composite_mask != (Image *) NULL)
      draw_info->composite_mask=DestroyImage(draw_info->composite_mask);
    draw_info->signature=(~MagickCoreSignature);
    draw_info=(DrawInfo *) RelinquishMagickMemory(draw_info);
    return(draw_info);
  }
  
  /*
  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  %                                                                             %
  %                                                                             %
  %                                                                             %
  +   D e s t r o y E d g e                                                     %
  %                                                                             %
  %                                                                             %
  %                                                                             %
  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  %
  %  DestroyEdge() destroys the specified polygon edge.
  %
  %  The format of the DestroyEdge method is:
  %
  %      ssize_t DestroyEdge(PolygonInfo *polygon_info,const int edge)
  %
  %  A description of each parameter follows:
  %
  %    o polygon_info: Specifies a pointer to an PolygonInfo structure.
  %
  %    o edge: the polygon edge number to destroy.
  %
  */
  static size_t DestroyEdge(PolygonInfo *polygon_info,
    const size_t edge)
  {
    assert(edge < polygon_info->number_edges);
    polygon_info->edges[edge].points=(PointInfo *) RelinquishMagickMemory(
      polygon_info->edges[edge].points);
    polygon_info->number_edges--;
    if (edge < polygon_info->number_edges)
      (void) memmove(polygon_info->edges+edge,polygon_info->edges+edge+1,
        (size_t) (polygon_info->number_edges-edge)*sizeof(*polygon_info->edges));
    return(polygon_info->number_edges);
  }
  
  /*
  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  %                                                                             %
  %                                                                             %
  %                                                                             %
  +   D e s t r o y P o l y g o n I n f o                                       %
  %                                                                             %
  %                                                                             %
  %                                                                             %
  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  %
  %  DestroyPolygonInfo() destroys the PolygonInfo data structure.
  %
  %  The format of the DestroyPolygonInfo method is:
  %
  %      PolygonInfo *DestroyPolygonInfo(PolygonInfo *polygon_info)
  %
  %  A description of each parameter follows:
  %
  %    o polygon_info: Specifies a pointer to an PolygonInfo structure.
  %
  */
  static PolygonInfo *DestroyPolygonInfo(PolygonInfo *polygon_info)
  {
    register ssize_t
      i;
  
    if (polygon_info->edges != (EdgeInfo *) NULL)
      {
        for (i=0; i < (ssize_t) polygon_info->number_edges; i++)
          polygon_info->edges[i].points=(PointInfo *)
            RelinquishMagickMemory(polygon_info->edges[i].points);
        polygon_info->edges=(EdgeInfo *) RelinquishMagickMemory(
          polygon_info->edges);
      }
    return((PolygonInfo *) RelinquishMagickMemory(polygon_info));
  }
  
  /*
  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  %                                                                             %
  %                                                                             %
  %                                                                             %
  %     D r a w A f f i n e I m a g e                                           %
  %                                                                             %
  %                                                                             %
  %                                                                             %
  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  %
  %  DrawAffineImage() composites the source over the destination image as
  %  dictated by the affine transform.
  %
  %  The format of the DrawAffineImage method is:
  %
  %      MagickBooleanType DrawAffineImage(Image *image,const Image *source,
  %        const AffineMatrix *affine,ExceptionInfo *exception)
  %
  %  A description of each parameter follows:
  %
  %    o image: the image.
  %
  %    o source: the source image.
  %
  %    o affine: the affine transform.
  %
  %    o exception: return any errors or warnings in this structure.
  %
  */
  
  static SegmentInfo AffineEdge(const Image *image,const AffineMatrix *affine,
    const double y,const SegmentInfo *edge)
  {
    double
      intercept,
      z;
  
    register double
      x;
  
    SegmentInfo
      inverse_edge;
  
    /*
      Determine left and right edges.
    */
    inverse_edge.x1=edge->x1;
    inverse_edge.y1=edge->y1;
    inverse_edge.x2=edge->x2;
    inverse_edge.y2=edge->y2;
    z=affine->ry*y+affine->tx;
    if (affine->sx >= MagickEpsilon)
      {
        intercept=(-z/affine->sx);
        x=intercept;
        if (x > inverse_edge.x1)
          inverse_edge.x1=x;
        intercept=(-z+(double) image->columns)/affine->sx;
        x=intercept;
        if (x < inverse_edge.x2)
          inverse_edge.x2=x;
      }
    else
      if (affine->sx < -MagickEpsilon)
        {
          intercept=(-z+(double) image->columns)/affine->sx;
          x=intercept;
          if (x > inverse_edge.x1)
            inverse_edge.x1=x;
          intercept=(-z/affine->sx);
          x=intercept;
          if (x < inverse_edge.x2)
            inverse_edge.x2=x;
        }
      else
        if ((z < 0.0) || ((size_t) floor(z+0.5) >= image->columns))
          {
            inverse_edge.x2=edge->x1;
            return(inverse_edge);
          }
    /*
      Determine top and bottom edges.
    */
    z=affine->sy*y+affine->ty;
    if (affine->rx >= MagickEpsilon)
      {
        intercept=(-z/affine->rx);
        x=intercept;
        if (x > inverse_edge.x1)
          inverse_edge.x1=x;
        intercept=(-z+(double) image->rows)/affine->rx;
        x=intercept;
        if (x < inverse_edge.x2)
          inverse_edge.x2=x;
      }
    else
      if (affine->rx < -MagickEpsilon)
        {
          intercept=(-z+(double) image->rows)/affine->rx;
          x=intercept;
          if (x > inverse_edge.x1)
            inverse_edge.x1=x;
          intercept=(-z/affine->rx);
          x=intercept;
          if (x < inverse_edge.x2)
            inverse_edge.x2=x;
        }
      else
        if ((z < 0.0) || ((size_t) floor(z+0.5) >= image->rows))
          {
            inverse_edge.x2=edge->x2;
            return(inverse_edge);
          }
    return(inverse_edge);
  }
  
  static AffineMatrix InverseAffineMatrix(const AffineMatrix *affine)
  {
    AffineMatrix
      inverse_affine;
  
    double
      determinant;
  
    determinant=PerceptibleReciprocal(affine->sx*affine->sy-affine->rx*
      affine->ry);
    inverse_affine.sx=determinant*affine->sy;
    inverse_affine.rx=determinant*(-affine->rx);
    inverse_affine.ry=determinant*(-affine->ry);
    inverse_affine.sy=determinant*affine->sx;
    inverse_affine.tx=(-affine->tx)*inverse_affine.sx-affine->ty*
      inverse_affine.ry;
    inverse_affine.ty=(-affine->tx)*inverse_affine.rx-affine->ty*
      inverse_affine.sy;
    return(inverse_affine);
  }
  
  MagickExport MagickBooleanType DrawAffineImage(Image *image,
    const Image *source,const AffineMatrix *affine,ExceptionInfo *exception)
  {
    AffineMatrix
      inverse_affine;
  
    CacheView
      *image_view,
      *source_view;
  
    MagickBooleanType
      status;
  
    PixelInfo
      zero;
  
    PointInfo
      extent[4],
      min,
      max;
  
    register ssize_t
      i;
  
    SegmentInfo
      edge;
  
    ssize_t
      start,
      stop,
      y;
  
    /*
      Determine bounding box.
    */
    assert(image != (Image *) NULL);
    assert(image->signature == MagickCoreSignature);
    if (image->debug != MagickFalse)
      (void) LogMagickEvent(TraceEvent,GetMagickModule(),"%s",image->filename);
    assert(source != (const Image *) NULL);
    assert(source->signature == MagickCoreSignature);
    assert(affine != (AffineMatrix *) NULL);
    extent[0].x=0.0;
    extent[0].y=0.0;
    extent[1].x=(double) source->columns-1.0;
    extent[1].y=0.0;
    extent[2].x=(double) source->columns-1.0;
    extent[2].y=(double) source->rows-1.0;
    extent[3].x=0.0;
    extent[3].y=(double) source->rows-1.0;
    for (i=0; i < 4; i++)
    {
      PointInfo
        point;
  
      point=extent[i];
      extent[i].x=point.x*affine->sx+point.y*affine->ry+affine->tx;
      extent[i].y=point.x*affine->rx+point.y*affine->sy+affine->ty;
    }
    min=extent[0];
    max=extent[0];
    for (i=1; i < 4; i++)
    {
      if (min.x > extent[i].x)
        min.x=extent[i].x;
      if (min.y > extent[i].y)
        min.y=extent[i].y;
      if (max.x < extent[i].x)
        max.x=extent[i].x;
      if (max.y < extent[i].y)
        max.y=extent[i].y;
    }
    /*
      Affine transform image.
    */
    if (SetImageStorageClass(image,DirectClass,exception) == MagickFalse)
      return(MagickFalse);
    status=MagickTrue;
    edge.x1=MagickMax(min.x,0.0);
    edge.y1=MagickMax(min.y,0.0);
    edge.x2=MagickMin(max.x,(double) image->columns-1.0);
    edge.y2=MagickMin(max.y,(double) image->rows-1.0);
    inverse_affine=InverseAffineMatrix(affine);
    GetPixelInfo(image,&zero);
    start=(ssize_t) ceil(edge.y1-0.5);
    stop=(ssize_t) floor(edge.y2+0.5);
    source_view=AcquireVirtualCacheView(source,exception);
    image_view=AcquireAuthenticCacheView(image,exception);
  #if defined(MAGICKCORE_OPENMP_SUPPORT)
    #pragma omp parallel for schedule(static) shared(status) \
      magick_number_threads(source,image,stop-start,1)
  #endif
    for (y=start; y <= stop; y++)
    {
      PixelInfo
        composite,
        pixel;
  
      PointInfo
        point;
  
      register ssize_t
        x;
  
      register Quantum
        *magick_restrict q;
  
      SegmentInfo
        inverse_edge;
  
      ssize_t
        x_offset;
  
      inverse_edge=AffineEdge(source,&inverse_affine,(double) y,&edge);
      if (inverse_edge.x2 < inverse_edge.x1)
        continue;
      q=GetCacheViewAuthenticPixels(image_view,(ssize_t) ceil(inverse_edge.x1-
        0.5),y,(size_t) (floor(inverse_edge.x2+0.5)-ceil(inverse_edge.x1-0.5)+1),
        1,exception);
      if (q == (Quantum *) NULL)
        continue;
      pixel=zero;
      composite=zero;
      x_offset=0;
      for (x=(ssize_t) ceil(inverse_edge.x1-0.5); x <= (ssize_t) floor(inverse_edge.x2+0.5); x++)
      {
        point.x=(double) x*inverse_affine.sx+y*inverse_affine.ry+
          inverse_affine.tx;
        point.y=(double) x*inverse_affine.rx+y*inverse_affine.sy+
          inverse_affine.ty;
        status=InterpolatePixelInfo(source,source_view,UndefinedInterpolatePixel,
          point.x,point.y,&pixel,exception);
        if (status == MagickFalse)
          break;
        GetPixelInfoPixel(image,q,&composite);
        CompositePixelInfoOver(&pixel,pixel.alpha,&composite,composite.alpha,
          &composite);
        SetPixelViaPixelInfo(image,&composite,q);
        x_offset++;
        q+=GetPixelChannels(image);
      }
      if (SyncCacheViewAuthenticPixels(image_view,exception) == MagickFalse)
        status=MagickFalse;
    }
    source_view=DestroyCacheView(source_view);
    image_view=DestroyCacheView(image_view);
    return(status);
  }
  
  /*
  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  %                                                                             %
  %                                                                             %
  %                                                                             %
  +   D r a w B o u n d i n g R e c t a n g l e s                               %
  %                                                                             %
  %                                                                             %
  %                                                                             %
  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  %
  %  DrawBoundingRectangles() draws the bounding rectangles on the image.  This
  %  is only useful for developers debugging the rendering algorithm.
  %
  %  The format of the DrawBoundingRectangles method is:
  %
  %      MagickBooleanType DrawBoundingRectangles(Image *image,
  %        const DrawInfo *draw_info,PolygonInfo *polygon_info,
  %        ExceptionInfo *exception)
  %
  %  A description of each parameter follows:
  %
  %    o image: the image.
  %
  %    o draw_info: the draw info.
  %
  %    o polygon_info: Specifies a pointer to a PolygonInfo structure.
  %
  %    o exception: return any errors or warnings in this structure.
  %
  */
  
  static MagickBooleanType DrawBoundingRectangles(Image *image,
    const DrawInfo *draw_info,const PolygonInfo *polygon_info,
    ExceptionInfo *exception)
  {
    double
      mid;
  
    DrawInfo
      *clone_info;
  
    MagickStatusType
      status;
  
    PointInfo
      end,
      resolution,
      start;
  
    PrimitiveInfo
      primitive_info[6];
  
    register ssize_t
      i;
  
    SegmentInfo
      bounds;
  
    ssize_t
      coordinates;
  
    (void) memset(primitive_info,0,sizeof(primitive_info));
    clone_info=CloneDrawInfo((ImageInfo *) NULL,draw_info);
    status=QueryColorCompliance("#000F",AllCompliance,&clone_info->fill,
      exception);
    if (status == MagickFalse)
      {
        clone_info=DestroyDrawInfo(clone_info);
        return(MagickFalse);
      }
    resolution.x=96.0;
    resolution.y=96.0;
    if (clone_info->density != (char *) NULL)
      {
        GeometryInfo
          geometry_info;
  
        MagickStatusType
          flags;
  
        flags=ParseGeometry(clone_info->density,&geometry_info);
        resolution.x=geometry_info.rho;
        resolution.y=geometry_info.sigma;
        if ((flags & SigmaValue) == MagickFalse)
          resolution.y=resolution.x;
      }
    mid=(resolution.x/96.0)*ExpandAffine(&clone_info->affine)*
      clone_info->stroke_width/2.0;
    bounds.x1=0.0;
    bounds.y1=0.0;
    bounds.x2=0.0;
    bounds.y2=0.0;
    if (polygon_info != (PolygonInfo *) NULL)
      {
        bounds=polygon_info->edges[0].bounds;
        for (i=1; i < (ssize_t) polygon_info->number_edges; i++)
        {
          if (polygon_info->edges[i].bounds.x1 < (double) bounds.x1)
            bounds.x1=polygon_info->edges[i].bounds.x1;
          if (polygon_info->edges[i].bounds.y1 < (double) bounds.y1)
            bounds.y1=polygon_info->edges[i].bounds.y1;
          if (polygon_info->edges[i].bounds.x2 > (double) bounds.x2)
            bounds.x2=polygon_info->edges[i].bounds.x2;
          if (polygon_info->edges[i].bounds.y2 > (double) bounds.y2)
            bounds.y2=polygon_info->edges[i].bounds.y2;
        }
        bounds.x1-=mid;
        bounds.x1=bounds.x1 < 0.0 ? 0.0 : bounds.x1 >= (double)
          image->columns ? (double) image->columns-1 : bounds.x1;
        bounds.y1-=mid;
        bounds.y1=bounds.y1 < 0.0 ? 0.0 : bounds.y1 >= (double)
          image->rows ? (double) image->rows-1 : bounds.y1;
        bounds.x2+=mid;
        bounds.x2=bounds.x2 < 0.0 ? 0.0 : bounds.x2 >= (double)
          image->columns ? (double) image->columns-1 : bounds.x2;
        bounds.y2+=mid;
        bounds.y2=bounds.y2 < 0.0 ? 0.0 : bounds.y2 >= (double)
          image->rows ? (double) image->rows-1 : bounds.y2;
        for (i=0; i < (ssize_t) polygon_info->number_edges; i++)
        {
          if (polygon_info->edges[i].direction != 0)
            status=QueryColorCompliance("#f00",AllCompliance,&clone_info->stroke,
              exception);
          else
            status=QueryColorCompliance("#0f0",AllCompliance,&clone_info->stroke,
              exception);
          if (status == MagickFalse)
            break;
          start.x=(double) (polygon_info->edges[i].bounds.x1-mid);
          start.y=(double) (polygon_info->edges[i].bounds.y1-mid);
          end.x=(double) (polygon_info->edges[i].bounds.x2+mid);
          end.y=(double) (polygon_info->edges[i].bounds.y2+mid);
          primitive_info[0].primitive=RectanglePrimitive;
          status&=TraceRectangle(primitive_info,start,end);
          primitive_info[0].method=ReplaceMethod;
          coordinates=(ssize_t) primitive_info[0].coordinates;
          primitive_info[coordinates].primitive=UndefinedPrimitive;
          status=DrawPrimitive(image,clone_info,primitive_info,exception);
          if (status == MagickFalse)
            break;
        }
        if (i < (ssize_t) polygon_info->number_edges)
          {
            clone_info=DestroyDrawInfo(clone_info);
            return(status == 0 ? MagickFalse : MagickTrue);
          }
      }
    status=QueryColorCompliance("#00f",AllCompliance,&clone_info->stroke,
      exception);
    if (status == MagickFalse)
      {
        clone_info=DestroyDrawInfo(clone_info);
        return(MagickFalse);
      }
    start.x=(double) (bounds.x1-mid);
    start.y=(double) (bounds.y1-mid);
    end.x=(double) (bounds.x2+mid);
    end.y=(double) (bounds.y2+mid);
    primitive_info[0].primitive=RectanglePrimitive;
    status&=TraceRectangle(primitive_info,start,end);
    primitive_info[0].method=ReplaceMethod;
    coordinates=(ssize_t) primitive_info[0].coordinates;
    primitive_info[coordinates].primitive=UndefinedPrimitive;
    status=DrawPrimitive(image,clone_info,primitive_info,exception);
    clone_info=DestroyDrawInfo(clone_info);
    return(status == 0 ? MagickFalse : MagickTrue);
  }
  
  /*
  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  %                                                                             %
  %                                                                             %
  %                                                                             %
  %   D r a w C l i p P a t h                                                   %
  %                                                                             %
  %                                                                             %
  %                                                                             %
  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  %
  %  DrawClipPath() draws the clip path on the image mask.
  %
  %  The format of the DrawClipPath method is:
  %
  %      MagickBooleanType DrawClipPath(Image *image,const DrawInfo *draw_info,
  %        const char *id,ExceptionInfo *exception)
  %
  %  A description of each parameter follows:
  %
  %    o image: the image.
  %
  %    o draw_info: the draw info.
  %
  %    o id: the clip path id.
  %
  %    o exception: return any errors or warnings in this structure.
  %
  */
  MagickExport MagickBooleanType DrawClipPath(Image *image,
    const DrawInfo *draw_info,const char *id,ExceptionInfo *exception)
  {
    const char
      *clip_path;
  
    Image
      *clipping_mask;
  
    MagickBooleanType
      status;
  
    clip_path=GetImageArtifact(image,id);
    if (clip_path == (const char *) NULL)
      return(MagickFalse);
    clipping_mask=DrawClippingMask(image,draw_info,draw_info->clip_mask,clip_path,
      exception);
    if (clipping_mask == (Image *) NULL)
      return(MagickFalse);
    status=SetImageMask(image,WritePixelMask,clipping_mask,exception);
    clipping_mask=DestroyImage(clipping_mask);
    return(status);
  }
  
  /*
  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  %                                                                             %
  %                                                                             %
  %                                                                             %
  %   D r a w C l i p p i n g M a s k                                           %
  %                                                                             %
  %                                                                             %
  %                                                                             %
  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  %
  %  DrawClippingMask() draws the clip path and returns it as an image clipping
  %  mask.
  %
  %  The format of the DrawClippingMask method is:
  %
  %      Image *DrawClippingMask(Image *image,const DrawInfo *draw_info,
  %        const char *id,const char *clip_path,ExceptionInfo *exception)
  %
  %  A description of each parameter follows:
  %
  %    o image: the image.
  %
  %    o draw_info: the draw info.
  %
  %    o id: the clip path id.
  %
  %    o clip_path: the clip path.
  %
  %    o exception: return any errors or warnings in this structure.
  %
  */
  static Image *DrawClippingMask(Image *image,const DrawInfo *draw_info,
    const char *id,const char *clip_path,ExceptionInfo *exception)
  {
    DrawInfo
      *clone_info;
  
    Image
      *clip_mask,
      *separate_mask;
  
    MagickStatusType
      status;
  
    /*
      Draw a clip path.
    */
    assert(image != (Image *) NULL);
    assert(image->signature == MagickCoreSignature);
    if (image->debug != MagickFalse)
      (void) LogMagickEvent(TraceEvent,GetMagickModule(),"%s",image->filename);
    assert(draw_info != (const DrawInfo *) NULL);
    clip_mask=AcquireImage((const ImageInfo *) NULL,exception);
    status=SetImageExtent(clip_mask,image->columns,image->rows,exception);
    if (status == MagickFalse)
      return(DestroyImage(clip_mask));
    status=SetImageMask(clip_mask,WritePixelMask,(Image *) NULL,exception);
    status=QueryColorCompliance("#0000",AllCompliance,
      &clip_mask->background_color,exception);
    clip_mask->background_color.alpha=(MagickRealType) TransparentAlpha;
    clip_mask->background_color.alpha_trait=BlendPixelTrait;
    status=SetImageBackgroundColor(clip_mask,exception);
    if (image->debug != MagickFalse)
      (void) LogMagickEvent(DrawEvent,GetMagickModule(),"\nbegin clip-path %s",
        id);
    clone_info=CloneDrawInfo((ImageInfo *) NULL,draw_info);
    (void) CloneString(&clone_info->primitive,clip_path);
    status=QueryColorCompliance("#ffffff",AllCompliance,&clone_info->fill,
      exception);
    if (clone_info->clip_mask != (char *) NULL)
      clone_info->clip_mask=DestroyString(clone_info->clip_mask);
    status=QueryColorCompliance("#00000000",AllCompliance,&clone_info->stroke,
      exception);
    clone_info->stroke_width=0.0;
    clone_info->alpha=OpaqueAlpha;
    clone_info->clip_path=MagickTrue;
    status=RenderMVGContent(clip_mask,clone_info,0,exception);
    clone_info=DestroyDrawInfo(clone_info);
    separate_mask=SeparateImage(clip_mask,AlphaChannel,exception);
    if (separate_mask != (Image *) NULL)
      {
        clip_mask=DestroyImage(clip_mask);
        clip_mask=separate_mask;
        status=NegateImage(clip_mask,MagickFalse,exception);
        if (status == MagickFalse)
          clip_mask=DestroyImage(clip_mask);
      }
    if (status == MagickFalse)
      clip_mask=DestroyImage(clip_mask);
    if (image->debug != MagickFalse)
      (void) LogMagickEvent(DrawEvent,GetMagickModule(),"end clip-path");
    return(clip_mask);
  }
  
  /*
  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  %                                                                             %
  %                                                                             %
  %                                                                             %
  %   D r a w C o m p o s i t e M a s k                                         %
  %                                                                             %
  %                                                                             %
  %                                                                             %
  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  %
  %  DrawCompositeMask() draws the mask path and returns it as an image mask.
  %
  %  The format of the DrawCompositeMask method is:
  %
  %      Image *DrawCompositeMask(Image *image,const DrawInfo *draw_info,
  %        const char *id,const char *mask_path,ExceptionInfo *exception)
  %
  %  A description of each parameter follows:
  %
  %    o image: the image.
  %
  %    o draw_info: the draw info.
  %
  %    o id: the mask path id.
  %
  %    o mask_path: the mask path.
  %
  %    o exception: return any errors or warnings in this structure.
  %
  */
  static Image *DrawCompositeMask(Image *image,const DrawInfo *draw_info,
    const char *id,const char *mask_path,ExceptionInfo *exception)
  {
    Image
      *composite_mask,
      *separate_mask;
  
    DrawInfo
      *clone_info;
  
    MagickStatusType
      status;
  
    /*
      Draw a mask path.
    */
    assert(image != (Image *) NULL);
    assert(image->signature == MagickCoreSignature);
    if (image->debug != MagickFalse)
      (void) LogMagickEvent(TraceEvent,GetMagickModule(),"%s",image->filename);
    assert(draw_info != (const DrawInfo *) NULL);
    composite_mask=AcquireImage((const ImageInfo *) NULL,exception);
    status=SetImageExtent(composite_mask,image->columns,image->rows,exception);
    if (status == MagickFalse)
      return(DestroyImage(composite_mask));
    status=SetImageMask(composite_mask,CompositePixelMask,(Image *) NULL,
      exception);
    status=QueryColorCompliance("#0000",AllCompliance,
      &composite_mask->background_color,exception);
    composite_mask->background_color.alpha=(MagickRealType) TransparentAlpha;
    composite_mask->background_color.alpha_trait=BlendPixelTrait;
    (void) SetImageBackgroundColor(composite_mask,exception);
    if (image->debug != MagickFalse)
      (void) LogMagickEvent(DrawEvent,GetMagickModule(),"\nbegin mask-path %s",
        id);
    clone_info=CloneDrawInfo((ImageInfo *) NULL,draw_info);
    (void) CloneString(&clone_info->primitive,mask_path);
    status=QueryColorCompliance("#ffffff",AllCompliance,&clone_info->fill,
      exception);
    status=QueryColorCompliance("#00000000",AllCompliance,&clone_info->stroke,
      exception);
    clone_info->stroke_width=0.0;
    clone_info->alpha=OpaqueAlpha;
    status=RenderMVGContent(composite_mask,clone_info,0,exception);
    clone_info=DestroyDrawInfo(clone_info);
    separate_mask=SeparateImage(composite_mask,AlphaChannel,exception);
    if (separate_mask != (Image *) NULL)
      {
        composite_mask=DestroyImage(composite_mask);
        composite_mask=separate_mask;
        status=NegateImage(composite_mask,MagickFalse,exception);
        if (status == MagickFalse)
          composite_mask=DestroyImage(composite_mask);
      }
    if (status == MagickFalse)
      composite_mask=DestroyImage(composite_mask);
    if (image->debug != MagickFalse)
      (void) LogMagickEvent(DrawEvent,GetMagickModule(),"end mask-path");
    return(composite_mask);
  }
  
  /*
  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  %                                                                             %
  %                                                                             %
  %                                                                             %
  +   D r a w D a s h P o l y g o n                                             %
  %                                                                             %
  %                                                                             %
  %                                                                             %
  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  %
  %  DrawDashPolygon() draws a dashed polygon (line, rectangle, ellipse) on the
  %  image while respecting the dash offset and dash pattern attributes.
  %
  %  The format of the DrawDashPolygon method is:
  %
  %      MagickBooleanType DrawDashPolygon(const DrawInfo *draw_info,
  %        const PrimitiveInfo *primitive_info,Image *image,
  %        ExceptionInfo *exception)
  %
  %  A description of each parameter follows:
  %
  %    o draw_info: the draw info.
  %
  %    o primitive_info: Specifies a pointer to a PrimitiveInfo structure.
  %
  %    o image: the image.
  %
  %    o exception: return any errors or warnings in this structure.
  %
  */
  static MagickBooleanType DrawDashPolygon(const DrawInfo *draw_info,
    const PrimitiveInfo *primitive_info,Image *image,ExceptionInfo *exception)
  {
    double
      length,
      maximum_length,
      offset,
      scale,
      total_length;
  
    DrawInfo
      *clone_info;
  
    MagickStatusType
      status;
  
    PrimitiveInfo
      *dash_polygon;
  
    register double
      dx,
      dy;
  
    register ssize_t
      i;
  
    size_t
      number_vertices;
  
    ssize_t
      j,
      n;
  
    assert(draw_info != (const DrawInfo *) NULL);
    if (image->debug != MagickFalse)
      (void) LogMagickEvent(DrawEvent,GetMagickModule(),"    begin draw-dash");
    for (i=0; primitive_info[i].primitive != UndefinedPrimitive; i++) ;
    number_vertices=(size_t) i;
    dash_polygon=(PrimitiveInfo *) AcquireQuantumMemory((size_t)
      (2UL*number_vertices+32UL),sizeof(*dash_polygon));
    if (dash_polygon == (PrimitiveInfo *) NULL)
      return(MagickFalse);
    (void) memset(dash_polygon,0,(2UL*number_vertices+32UL)*
      sizeof(*dash_polygon));
    clone_info=CloneDrawInfo((ImageInfo *) NULL,draw_info);
    clone_info->miterlimit=0;
    dash_polygon[0]=primitive_info[0];
    scale=ExpandAffine(&draw_info->affine);
    length=scale*draw_info->dash_pattern[0];
    offset=fabs(draw_info->dash_offset) >= MagickEpsilon ?
      scale*draw_info->dash_offset : 0.0;
    j=1;
    for (n=0; offset > 0.0; j=0)
    {
      if (draw_info->dash_pattern[n] <= 0.0)
        break;
      length=scale*(draw_info->dash_pattern[n]+(n == 0 ? -0.5 : 0.5));
      if (offset > length)
        {
          offset-=length;
          n++;
          length=scale*draw_info->dash_pattern[n];
          continue;
        }
      if (offset < length)
        {
          length-=offset;
          offset=0.0;
          break;
        }
      offset=0.0;
      n++;
    }
    status=MagickTrue;
    maximum_length=0.0;
    total_length=0.0;
    for (i=1; (i < (ssize_t) number_vertices) && (length >= 0.0); i++)
    {
      dx=primitive_info[i].point.x-primitive_info[i-1].point.x;
      dy=primitive_info[i].point.y-primitive_info[i-1].point.y;
      maximum_length=hypot(dx,dy);
      if (maximum_length > (MaxBezierCoordinates >> 2))
        break;
      if (fabs(length) < MagickEpsilon)
        {
          if (fabs(draw_info->dash_pattern[n]) >= MagickEpsilon)
            n++;
          if (fabs(draw_info->dash_pattern[n]) < MagickEpsilon)
            n=0;
          length=scale*draw_info->dash_pattern[n];
        }
      for (total_length=0.0; (length >= 0.0) && (maximum_length >= (total_length+length)); )
      {
        total_length+=length;
        if ((n & 0x01) != 0)
          {
            dash_polygon[0]=primitive_info[0];
            dash_polygon[0].point.x=(double) (primitive_info[i-1].point.x+dx*
              total_length*PerceptibleReciprocal(maximum_length));
            dash_polygon[0].point.y=(double) (primitive_info[i-1].point.y+dy*
              total_length*PerceptibleReciprocal(maximum_length));
            j=1;
          }
        else
          {
            if ((j+1) > (ssize_t) number_vertices)
              break;
            dash_polygon[j]=primitive_info[i-1];
            dash_polygon[j].point.x=(double) (primitive_info[i-1].point.x+dx*
              total_length*PerceptibleReciprocal(maximum_length));
            dash_polygon[j].point.y=(double) (primitive_info[i-1].point.y+dy*
              total_length*PerceptibleReciprocal(maximum_length));
            dash_polygon[j].coordinates=1;
            j++;
            dash_polygon[0].coordinates=(size_t) j;
            dash_polygon[j].primitive=UndefinedPrimitive;
            status&=DrawStrokePolygon(image,clone_info,dash_polygon,exception);
            if (status == MagickFalse)
              break;
          }
        if (fabs(draw_info->dash_pattern[n]) >= MagickEpsilon)
          n++;
        if (fabs(draw_info->dash_pattern[n]) < MagickEpsilon)
          n=0;
        length=scale*draw_info->dash_pattern[n];
      }
      length-=(maximum_length-total_length);
      if ((n & 0x01) != 0)
        continue;
      dash_polygon[j]=primitive_info[i];
      dash_polygon[j].coordinates=1;
      j++;
    }
    if ((status != MagickFalse) && (total_length < maximum_length) &&
        ((n & 0x01) == 0) && (j > 1))
      {
        dash_polygon[j]=primitive_info[i-1];
        dash_polygon[j].point.x+=MagickEpsilon;
        dash_polygon[j].point.y+=MagickEpsilon;
        dash_polygon[j].coordinates=1;
        j++;
        dash_polygon[0].coordinates=(size_t) j;
        dash_polygon[j].primitive=UndefinedPrimitive;
        status&=DrawStrokePolygon(image,clone_info,dash_polygon,exception);
      }
    dash_polygon=(PrimitiveInfo *) RelinquishMagickMemory(dash_polygon);
    clone_info=DestroyDrawInfo(clone_info);
    if (image->debug != MagickFalse)
      (void) LogMagickEvent(DrawEvent,GetMagickModule(),"    end draw-dash");
    return(status != 0 ? MagickTrue : MagickFalse);
  }
  
  /*
  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  %                                                                             %
  %                                                                             %
  %                                                                             %
  %     D r a w G r a d i e n t I m a g e                                       %
  %                                                                             %
  %                                                                             %
  %                                                                             %
  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  %
  %  DrawGradientImage() draws a linear gradient on the image.
  %
  %  The format of the DrawGradientImage method is:
  %
  %      MagickBooleanType DrawGradientImage(Image *image,
  %        const DrawInfo *draw_info,ExceptionInfo *exception)
  %
  %  A description of each parameter follows:
  %
  %    o image: the image.
  %
  %    o draw_info: the draw info.
  %
  %    o exception: return any errors or warnings in this structure.
  %
  */
  
  static inline double GetStopColorOffset(const GradientInfo *gradient,
    const ssize_t x,const ssize_t y)
  {
    switch (gradient->type)
    {
      case UndefinedGradient:
      case LinearGradient:
      {
        double
          gamma,
          length,
          offset,
          scale;
  
        PointInfo
          p,
          q;
  
        const SegmentInfo
          *gradient_vector;
  
        gradient_vector=(&gradient->gradient_vector);
        p.x=gradient_vector->x2-gradient_vector->x1;
        p.y=gradient_vector->y2-gradient_vector->y1;
        q.x=(double) x-gradient_vector->x1;
        q.y=(double) y-gradient_vector->y1;
        length=sqrt(q.x*q.x+q.y*q.y);
        gamma=sqrt(p.x*p.x+p.y*p.y)*length;
        gamma=PerceptibleReciprocal(gamma);
        scale=p.x*q.x+p.y*q.y;
        offset=gamma*scale*length;
        return(offset);
      }
      case RadialGradient:
      {
        PointInfo
          v;
  
        if (gradient->spread == RepeatSpread)
          {
            v.x=(double) x-gradient->center.x;
            v.y=(double) y-gradient->center.y;
            return(sqrt(v.x*v.x+v.y*v.y));
          }
        v.x=(double) (((x-gradient->center.x)*cos(DegreesToRadians(
          gradient->angle)))+((y-gradient->center.y)*sin(DegreesToRadians(
          gradient->angle))))*PerceptibleReciprocal(gradient->radii.x);
        v.y=(double) (((x-gradient->center.x)*sin(DegreesToRadians(
          gradient->angle)))-((y-gradient->center.y)*cos(DegreesToRadians(
          gradient->angle))))*PerceptibleReciprocal(gradient->radii.y);
        return(sqrt(v.x*v.x+v.y*v.y));
      }
    }
    return(0.0);
  }
  
  static int StopInfoCompare(const void *x,const void *y)
  {
    StopInfo
      *stop_1,
      *stop_2;
  
    stop_1=(StopInfo *) x;
    stop_2=(StopInfo *) y;
    if (stop_1->offset > stop_2->offset)
      return(1);
    if (fabs(stop_1->offset-stop_2->offset) <= MagickEpsilon)
      return(0);
    return(-1);
  }
  
  MagickExport MagickBooleanType DrawGradientImage(Image *image,
    const DrawInfo *draw_info,ExceptionInfo *exception)
  {
    CacheView
      *image_view;
  
    const GradientInfo
      *gradient;
  
    const SegmentInfo
      *gradient_vector;
  
    double
      length;
  
    MagickBooleanType
      status;
  
    PixelInfo
      zero;
  
    PointInfo
      point;
  
    RectangleInfo
      bounding_box;
  
    ssize_t
      y;
  
    /*
      Draw linear or radial gradient on image.
    */
    assert(image != (Image *) NULL);
    assert(image->signature == MagickCoreSignature);
    if (image->debug != MagickFalse)
      (void) LogMagickEvent(TraceEvent,GetMagickModule(),"%s",image->filename);
    assert(draw_info != (const DrawInfo *) NULL);
    gradient=(&draw_info->gradient);
    qsort(gradient->stops,gradient->number_stops,sizeof(StopInfo),
      StopInfoCompare);
    gradient_vector=(&gradient->gradient_vector);
    point.x=gradient_vector->x2-gradient_vector->x1;
    point.y=gradient_vector->y2-gradient_vector->y1;
    length=sqrt(point.x*point.x+point.y*point.y);
    bounding_box=gradient->bounding_box;
    status=MagickTrue;
    GetPixelInfo(image,&zero);
    image_view=AcquireAuthenticCacheView(image,exception);
  #if defined(MAGICKCORE_OPENMP_SUPPORT)
    #pragma omp parallel for schedule(static) shared(status) \
      magick_number_threads(image,image,bounding_box.height-bounding_box.y,1)
  #endif
    for (y=bounding_box.y; y < (ssize_t) bounding_box.height; y++)
    {
      double
        alpha,
        offset;
  
      PixelInfo
        composite,
        pixel;
  
      register Quantum
        *magick_restrict q;
  
      register ssize_t
        i,
        x;
  
      ssize_t
        j;
  
      if (status == MagickFalse)
        continue;
      q=GetCacheViewAuthenticPixels(image_view,0,y,image->columns,1,exception);
      if (q == (Quantum *) NULL)
        {
          status=MagickFalse;
          continue;
        }
      pixel=zero;
      composite=zero;
      offset=GetStopColorOffset(gradient,0,y);
      if (gradient->type != RadialGradient)
        offset*=PerceptibleReciprocal(length);
      for (x=bounding_box.x; x < (ssize_t) bounding_box.width; x++)
      {
        GetPixelInfoPixel(image,q,&pixel);
        switch (gradient->spread)
        {
          case UndefinedSpread:
          case PadSpread:
          {
            if ((x != (ssize_t) ceil(gradient_vector->x1-0.5)) ||
                (y != (ssize_t) ceil(gradient_vector->y1-0.5)))
              {
                offset=GetStopColorOffset(gradient,x,y);
                if (gradient->type != RadialGradient)
                  offset*=PerceptibleReciprocal(length);
              }
            for (i=0; i < (ssize_t) gradient->number_stops; i++)
              if (offset < gradient->stops[i].offset)
                break;
            if ((offset < 0.0) || (i == 0))
              composite=gradient->stops[0].color;
            else
              if ((offset > 1.0) || (i == (ssize_t) gradient->number_stops))
                composite=gradient->stops[gradient->number_stops-1].color;
              else
                {
                  j=i;
                  i--;
                  alpha=(offset-gradient->stops[i].offset)/
                    (gradient->stops[j].offset-gradient->stops[i].offset);
                  CompositePixelInfoBlend(&gradient->stops[i].color,1.0-alpha,
                    &gradient->stops[j].color,alpha,&composite);
                }
            break;
          }
          case ReflectSpread:
          {
            if ((x != (ssize_t) ceil(gradient_vector->x1-0.5)) ||
                (y != (ssize_t) ceil(gradient_vector->y1-0.5)))
              {
                offset=GetStopColorOffset(gradient,x,y);
                if (gradient->type != RadialGradient)
                  offset*=PerceptibleReciprocal(length);
              }
            if (offset < 0.0)
              offset=(-offset);
            if ((ssize_t) fmod(offset,2.0) == 0)
              offset=fmod(offset,1.0);
            else
              offset=1.0-fmod(offset,1.0);
            for (i=0; i < (ssize_t) gradient->number_stops; i++)
              if (offset < gradient->stops[i].offset)
                break;
            if (i == 0)
              composite=gradient->stops[0].color;
            else
              if (i == (ssize_t) gradient->number_stops)
                composite=gradient->stops[gradient->number_stops-1].color;
              else
                {
                  j=i;
                  i--;
                  alpha=(offset-gradient->stops[i].offset)/
                    (gradient->stops[j].offset-gradient->stops[i].offset);
                  CompositePixelInfoBlend(&gradient->stops[i].color,1.0-alpha,
                    &gradient->stops[j].color,alpha,&composite);
                }
            break;
          }
          case RepeatSpread:
          {
            double
              repeat;
  
            MagickBooleanType
              antialias;
  
            antialias=MagickFalse;
            repeat=0.0;
            if ((x != (ssize_t) ceil(gradient_vector->x1-0.5)) ||
                (y != (ssize_t) ceil(gradient_vector->y1-0.5)))
              {
                offset=GetStopColorOffset(gradient,x,y);
                if (gradient->type == LinearGradient)
                  {
                    repeat=fmod(offset,length);
                    if (repeat < 0.0)
                      repeat=length-fmod(-repeat,length);
                    else
                      repeat=fmod(offset,length);
                    antialias=(repeat < length) && ((repeat+1.0) > length) ?
                      MagickTrue : MagickFalse;
                    offset=PerceptibleReciprocal(length)*repeat;
                  }
                else
                  {
                    repeat=fmod(offset,gradient->radius);
                    if (repeat < 0.0)
                      repeat=gradient->radius-fmod(-repeat,gradient->radius);
                    else
                      repeat=fmod(offset,gradient->radius);
                    antialias=repeat+1.0 > gradient->radius ? MagickTrue :
                      MagickFalse;
                    offset=repeat/gradient->radius;
                  }
              }
            for (i=0; i < (ssize_t) gradient->number_stops; i++)
              if (offset < gradient->stops[i].offset)
                break;
            if (i == 0)
              composite=gradient->stops[0].color;
            else
              if (i == (ssize_t) gradient->number_stops)
                composite=gradient->stops[gradient->number_stops-1].color;
              else
                {
                  j=i;
                  i--;
                  alpha=(offset-gradient->stops[i].offset)/
                    (gradient->stops[j].offset-gradient->stops[i].offset);
                  if (antialias != MagickFalse)
                    {
                      if (gradient->type == LinearGradient)
                        alpha=length-repeat;
                      else
                        alpha=gradient->radius-repeat;
                      i=0;
                      j=(ssize_t) gradient->number_stops-1L;
                    }
                  CompositePixelInfoBlend(&gradient->stops[i].color,1.0-alpha,
                    &gradient->stops[j].color,alpha,&composite);
                }
            break;
          }
        }
        CompositePixelInfoOver(&composite,composite.alpha,&pixel,pixel.alpha,
          &pixel);
        SetPixelViaPixelInfo(image,&pixel,q);
        q+=GetPixelChannels(image);
      }
      if (SyncCacheViewAuthenticPixels(image_view,exception) == MagickFalse)
        status=MagickFalse;
    }
    image_view=DestroyCacheView(image_view);
    return(status);
  }
  
  /*
  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  %                                                                             %
  %                                                                             %
  %                                                                             %
  %   D r a w I m a g e                                                         %
  %                                                                             %
  %                                                                             %
  %                                                                             %
  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  %
  %  DrawImage() draws a graphic primitive on your image.  The primitive
  %  may be represented as a string or filename.  Precede the filename with an
  %  "at" sign (@) and the contents of the file are drawn on the image.  You
  %  can affect how text is drawn by setting one or more members of the draw
  %  info structure.
  %
  %  The format of the DrawImage method is:
  %
  %      MagickBooleanType DrawImage(Image *image,const DrawInfo *draw_info,
  %        ExceptionInfo *exception)
  %
  %  A description of each parameter follows:
  %
  %    o image: the image.
  %
  %    o draw_info: the draw info.
  %
  %    o exception: return any errors or warnings in this structure.
  %
  */
  
  static MagickBooleanType CheckPrimitiveExtent(MVGInfo *mvg_info,
    const size_t pad)
  {
    double
      extent;
  
    size_t
      quantum;
  
    /*
      Check if there is enough storage for drawing pimitives.
    */
    extent=(double) mvg_info->offset+pad+PrimitiveExtentPad;
    quantum=sizeof(**mvg_info->primitive_info);
    if (((extent*quantum) < (double) SSIZE_MAX) &&
        ((extent*quantum) < (double) GetMaxMemoryRequest()))
      {
        if (extent <= (double) *mvg_info->extent)
          return(MagickTrue);
        *mvg_info->primitive_info=(PrimitiveInfo *) ResizeQuantumMemory(
          *mvg_info->primitive_info,(size_t) extent,quantum);
        if (*mvg_info->primitive_info != (PrimitiveInfo *) NULL)
          {
            register ssize_t
              i;
  
            *mvg_info->extent=(size_t) extent;
            for (i=mvg_info->offset+1; i < (ssize_t) extent; i++)
              (*mvg_info->primitive_info)[i].primitive=UndefinedPrimitive;
            return(MagickTrue);
          }
      }
    /*
      Reallocation failed, allocate a primitive to facilitate unwinding.
    */
    (void) ThrowMagickException(mvg_info->exception,GetMagickModule(),
      ResourceLimitError,"MemoryAllocationFailed","`%s'","");
    if (*mvg_info->primitive_info != (PrimitiveInfo *) NULL)
      *mvg_info->primitive_info=(PrimitiveInfo *) RelinquishMagickMemory(
        *mvg_info->primitive_info);
    *mvg_info->primitive_info=(PrimitiveInfo *) AcquireCriticalMemory(
      PrimitiveExtentPad*quantum);
    (void) memset(*mvg_info->primitive_info,0,PrimitiveExtentPad*quantum);
    *mvg_info->extent=1;
    return(MagickFalse);
  }
  
  static inline double GetDrawValue(const char *magick_restrict string,
    char **magick_restrict sentinal)
  {
    char
      **magick_restrict q;
  
    double
      value;
  
    q=sentinal;
    value=InterpretLocaleValue(string,q);
    if ((IsNaN(value) != 0) || (value < -(SSIZE_MAX-512.0)) ||
        (value > (SSIZE_MAX-512.0)))
      return(0.0);
    sentinal=q;
    return(value);
  }
  
  static int MVGMacroCompare(const void *target,const void *source)
  {
    const char
      *p,
      *q;
  
    p=(const char *) target;
    q=(const char *) source;
    return(strcmp(p,q));
  }
  
  static SplayTreeInfo *GetMVGMacros(const char *primitive)
  {
    char
      *macro,
      *token;
  
    const char
      *q;
  
    size_t
      extent;
  
    SplayTreeInfo
      *macros;
  
    /*
      Scan graphic primitives for definitions and classes.
    */
    if (primitive == (const char *) NULL)
      return((SplayTreeInfo *) NULL);
    macros=NewSplayTree(MVGMacroCompare,RelinquishMagickMemory,
      RelinquishMagickMemory);
    macro=AcquireString(primitive);
    token=AcquireString(primitive);
    extent=strlen(token)+MagickPathExtent;
    for (q=primitive; *q != '\0'; )
    {
      if (GetNextToken(q,&q,extent,token) < 1)
        break;
      if (*token == '\0')
        break;
      if (LocaleCompare("push",token) == 0)
        {
          register const char
            *end,
            *start;
  
          (void) GetNextToken(q,&q,extent,token);
          if (*q == '"')
            {
              char
                name[MagickPathExtent];
  
              const char
                *p;
  
              ssize_t
                n;
  
              /*
                Named macro (e.g. push graphic-context "wheel").
              */
              (void) GetNextToken(q,&q,extent,token);
              start=q;
              end=q;
              (void) CopyMagickString(name,token,MagickPathExtent);
              n=1;
              for (p=q; *p != '\0'; )
              {
                if (GetNextToken(p,&p,extent,token) < 1)
                  break;
                if (*token == '\0')
                  break;
                if (LocaleCompare(token,"pop") == 0)
                  {
                    end=p-strlen(token)-1;
                    n--;
                  }
                if (LocaleCompare(token,"push") == 0)
                  n++;
                if ((n == 0) && (end > start))
                  {
                    /*
                      Extract macro.
                    */
                    (void) GetNextToken(p,&p,extent,token);
                    (void) CopyMagickString(macro,start,(size_t) (end-start));
                    (void) AddValueToSplayTree(macros,ConstantString(name),
                      ConstantString(macro));
                    break;
                  }
              }
            }
        }
    }
    token=DestroyString(token);
    macro=DestroyString(macro);
    return(macros);
  }
  
  static inline MagickBooleanType IsPoint(const char *point)
  {
    char
      *p;
  
    double
      value;
  
    value=GetDrawValue(point,&p);
    return((fabs(value) < MagickEpsilon) && (p == point) ? MagickFalse :
      MagickTrue);
  }
  
  static inline MagickBooleanType TracePoint(PrimitiveInfo *primitive_info,
    const PointInfo point)
  {
    primitive_info->coordinates=1;
    primitive_info->closed_subpath=MagickFalse;
    primitive_info->point=point;
    return(MagickTrue);
  }
  
  static MagickBooleanType RenderMVGContent(Image *image,
    const DrawInfo *draw_info,const size_t depth,ExceptionInfo *exception)
  {
  #define RenderImageTag  "Render/Image"
  
    AffineMatrix
      affine,
      current;
  
    char
      keyword[MagickPathExtent],
      geometry[MagickPathExtent],
      *next_token,
      pattern[MagickPathExtent],
      *primitive,
      *token;
  
    const char
      *q;
  
    double
      angle,
      coordinates,
      cursor,
      factor,
      primitive_extent;
  
    DrawInfo
      *clone_info,
      **graphic_context;
  
    MagickBooleanType
      proceed;
  
    MagickStatusType
      status;
  
    MVGInfo
      mvg_info;
  
    PointInfo
      point;
  
    PrimitiveInfo
      *primitive_info;
  
    PrimitiveType
      primitive_type;
  
    register const char
      *p;
  
    register ssize_t
      i,
      x;
  
    SegmentInfo
      bounds;
  
    size_t
      extent,
      number_points,
      number_stops;
  
    SplayTreeInfo
      *macros;
  
    ssize_t
      defsDepth,
      j,
      k,
      n,
      symbolDepth;
  
    StopInfo
      *stops;
  
    TypeMetric
      metrics;
  
    assert(image != (Image *) NULL);
    assert(image->signature == MagickCoreSignature);
    if (image->debug != MagickFalse)
      (void) LogMagickEvent(TraceEvent,GetMagickModule(),"%s",image->filename);
    assert(draw_info != (DrawInfo *) NULL);
    assert(draw_info->signature == MagickCoreSignature);
    if (image->debug != MagickFalse)
      (void) LogMagickEvent(TraceEvent,GetMagickModule(),"...");
    if (depth > MagickMaxRecursionDepth)
      ThrowBinaryException(DrawError,"VectorGraphicsNestedTooDeeply",
        image->filename);
    if ((draw_info->primitive == (char *) NULL) ||
        (*draw_info->primitive == '\0'))
      return(MagickFalse);
    if (image->debug != MagickFalse)
      (void) LogMagickEvent(DrawEvent,GetMagickModule(),"begin draw-image");
    if (SetImageStorageClass(image,DirectClass,exception) == MagickFalse)
      return(MagickFalse);
    if (image->alpha_trait == UndefinedPixelTrait)
      {
        status=SetImageAlphaChannel(image,OpaqueAlphaChannel,exception);
        if (status == MagickFalse)
          return(MagickFalse);
      }
    if ((*draw_info->primitive == '@') && (strlen(draw_info->primitive) > 1) &&
        (*(draw_info->primitive+1) != '-') && (depth == 0))
      primitive=FileToString(draw_info->primitive+1,~0UL,exception);
    else
      primitive=AcquireString(draw_info->primitive);
    if (primitive == (char *) NULL)
      return(MagickFalse);
    primitive_extent=(double) strlen(primitive);
    (void) SetImageArtifact(image,"mvg:vector-graphics",primitive);
    n=0;
    number_stops=0;
    stops=(StopInfo *) NULL;
    /*
      Allocate primitive info memory.
    */
    graphic_context=(DrawInfo **) AcquireMagickMemory(sizeof(*graphic_context));
    if (graphic_context == (DrawInfo **) NULL)
      {
        primitive=DestroyString(primitive);
        ThrowBinaryException(ResourceLimitError,"MemoryAllocationFailed",
          image->filename);
      }
    number_points=PrimitiveExtentPad;
    primitive_info=(PrimitiveInfo *) AcquireQuantumMemory((size_t) number_points,
      sizeof(*primitive_info));
    if (primitive_info == (PrimitiveInfo *) NULL)
      {
        primitive=DestroyString(primitive);
        for ( ; n >= 0; n--)
          graphic_context[n]=DestroyDrawInfo(graphic_context[n]);
        graphic_context=(DrawInfo **) RelinquishMagickMemory(graphic_context);
        ThrowBinaryException(ResourceLimitError,"MemoryAllocationFailed",
          image->filename);
      }
    (void) memset(primitive_info,0,(size_t) number_points*
      sizeof(*primitive_info));
    (void) memset(&mvg_info,0,sizeof(mvg_info));
    mvg_info.primitive_info=(&primitive_info);
    mvg_info.extent=(&number_points);
    mvg_info.exception=exception;
    graphic_context[n]=CloneDrawInfo((ImageInfo *) NULL,draw_info);
    graphic_context[n]->viewbox=image->page;
    if ((image->page.width == 0) || (image->page.height == 0))
      {
        graphic_context[n]->viewbox.width=image->columns;
        graphic_context[n]->viewbox.height=image->rows;
      }
    token=AcquireString(primitive);
    extent=strlen(token)+MagickPathExtent;
    defsDepth=0;
    symbolDepth=0;
    cursor=0.0;
    macros=GetMVGMacros(primitive);
    status=MagickTrue;
    for (q=primitive; *q != '\0'; )
    {
      /*
        Interpret graphic primitive.
      */
      if (GetNextToken(q,&q,MagickPathExtent,keyword) < 1)
        break;
      if (*keyword == '\0')
        break;
      if (*keyword == '#')
        {
          /*
            Comment.
          */
          while ((*q != '\n') && (*q != '\0'))
            q++;
          continue;
        }
      p=q-strlen(keyword)-1;
      primitive_type=UndefinedPrimitive;
      current=graphic_context[n]->affine;
      GetAffineMatrix(&affine);
      *token='\0';
      switch (*keyword)
      {
        case ';':
          break;
        case 'a':
        case 'A':
        {
          if (LocaleCompare("affine",keyword) == 0)
            {
              (void) GetNextToken(q,&q,extent,token);
              affine.sx=GetDrawValue(token,&next_token);
              if (token == next_token)
                ThrowPointExpectedException(token,exception);
              (void) GetNextToken(q,&q,extent,token);
              if (*token == ',')
                (void) GetNextToken(q,&q,extent,token);
              affine.rx=GetDrawValue(token,&next_token);
              if (token == next_token)
                ThrowPointExpectedException(token,exception);
              (void) GetNextToken(q,&q,extent,token);
              if (*token == ',')
                (void) GetNextToken(q,&q,extent,token);
              affine.ry=GetDrawValue(token,&next_token);
              if (token == next_token)
                ThrowPointExpectedException(token,exception);
              (void) GetNextToken(q,&q,extent,token);
              if (*token == ',')
                (void) GetNextToken(q,&q,extent,token);
              affine.sy=GetDrawValue(token,&next_token);
              if (token == next_token)
                ThrowPointExpectedException(token,exception);
              (void) GetNextToken(q,&q,extent,token);
              if (*token == ',')
                (void) GetNextToken(q,&q,extent,token);
              affine.tx=GetDrawValue(token,&next_token);
              if (token == next_token)
                ThrowPointExpectedException(token,exception);
              (void) GetNextToken(q,&q,extent,token);
              if (*token == ',')
                (void) GetNextToken(q,&q,extent,token);
              affine.ty=GetDrawValue(token,&next_token);
              if (token == next_token)
                ThrowPointExpectedException(token,exception);
              break;
            }
          if (LocaleCompare("alpha",keyword) == 0)
            {
              primitive_type=AlphaPrimitive;
              break;
            }
          if (LocaleCompare("arc",keyword) == 0)
            {
              primitive_type=ArcPrimitive;
              break;
            }
          status=MagickFalse;
          break;
        }
        case 'b':
        case 'B':
        {
          if (LocaleCompare("bezier",keyword) == 0)
            {
              primitive_type=BezierPrimitive;
              break;
            }
          if (LocaleCompare("border-color",keyword) == 0)
            {
              (void) GetNextToken(q,&q,extent,token);
              status&=QueryColorCompliance(token,AllCompliance,
                &graphic_context[n]->border_color,exception);
              break;
            }
          status=MagickFalse;
          break;
        }
        case 'c':
        case 'C':
        {
          if (LocaleCompare("class",keyword) == 0)
            {
              const char
                *mvg_class;
  
              (void) GetNextToken(q,&q,extent,token);
              if (*token == '\0')
                {
                  status=MagickFalse;
                  break;
                }
              if (LocaleCompare(token,graphic_context[n]->id) == 0)
                break;
              mvg_class=(const char *) GetValueFromSplayTree(macros,token);
              if ((mvg_class != (const char *) NULL) && (p > primitive))
                {
                  char
                    *elements;
  
                  ssize_t
                    offset;
  
                  /*
                    Inject class elements in stream.
                  */
                  offset=(ssize_t) (p-primitive);
                  elements=AcquireString(primitive);
                  elements[offset]='\0';
                  (void) ConcatenateString(&elements,mvg_class);
                  (void) ConcatenateString(&elements,"\n");
                  (void) ConcatenateString(&elements,q);
                  primitive=DestroyString(primitive);
                  primitive=elements;
                  q=primitive+offset;
                }
              break;
            }
          if (LocaleCompare("clip-path",keyword) == 0)
            {
              const char
                *clip_path;
  
              /*
                Take a node from within the MVG document, and duplicate it here.
              */
              (void) GetNextToken(q,&q,extent,token);
              if (*token == '\0')
                {
                  status=MagickFalse;
                  break;
                }
              (void) CloneString(&graphic_context[n]->clip_mask,token);
              clip_path=(const char *) GetValueFromSplayTree(macros,token);
              if (clip_path != (const char *) NULL)
                {
                  if (graphic_context[n]->clipping_mask != (Image *) NULL)
                    graphic_context[n]->clipping_mask=
                      DestroyImage(graphic_context[n]->clipping_mask);
                  graphic_context[n]->clipping_mask=DrawClippingMask(image,
                    graphic_context[n],token,clip_path,exception);
                  if (graphic_context[n]->compliance != SVGCompliance)
                    {
                      clip_path=(const char *) GetValueFromSplayTree(macros,
                        graphic_context[n]->clip_mask);
                      if (clip_path != (const char *) NULL)
                        (void) SetImageArtifact(image,
                          graphic_context[n]->clip_mask,clip_path);
                      status&=DrawClipPath(image,graphic_context[n],
                        graphic_context[n]->clip_mask,exception);
                    }
                }
              break;
            }
          if (LocaleCompare("clip-rule",keyword) == 0)
            {
              ssize_t
                fill_rule;
  
              (void) GetNextToken(q,&q,extent,token);
              fill_rule=ParseCommandOption(MagickFillRuleOptions,MagickFalse,
                token);
              if (fill_rule == -1)
                {
                  status=MagickFalse;
                  break;
                }
              graphic_context[n]->fill_rule=(FillRule) fill_rule;
              break;
            }
          if (LocaleCompare("clip-units",keyword) == 0)
            {
              ssize_t
                clip_units;
  
              (void) GetNextToken(q,&q,extent,token);
              clip_units=ParseCommandOption(MagickClipPathOptions,MagickFalse,
                token);
              if (clip_units == -1)
                {
                  status=MagickFalse;
                  break;
                }
              graphic_context[n]->clip_units=(ClipPathUnits) clip_units;
              if (clip_units == ObjectBoundingBox)
                {
                  GetAffineMatrix(&current);
                  affine.sx=draw_info->bounds.x2;
                  affine.sy=draw_info->bounds.y2;
                  affine.tx=draw_info->bounds.x1;
                  affine.ty=draw_info->bounds.y1;
                  break;
                }
              break;
            }
          if (LocaleCompare("circle",keyword) == 0)
            {
              primitive_type=CirclePrimitive;
              break;
            }
          if (LocaleCompare("color",keyword) == 0)
            {
              primitive_type=ColorPrimitive;
              break;
            }
          if (LocaleCompare("compliance",keyword) == 0)
            {
              /*
                MVG compliance associates a clipping mask with an image; SVG
                compliance associates a clipping mask with a graphics context.
              */
              (void) GetNextToken(q,&q,extent,token);
              graphic_context[n]->compliance=(ComplianceType) ParseCommandOption(
                MagickComplianceOptions,MagickFalse,token);
              break;
            }
          status=MagickFalse;
          break;
        }
        case 'd':
        case 'D':
        {
          if (LocaleCompare("decorate",keyword) == 0)
            {
              ssize_t
                decorate;
  
              (void) GetNextToken(q,&q,extent,token);
              decorate=ParseCommandOption(MagickDecorateOptions,MagickFalse,
                token);
              if (decorate == -1)
                {
                  status=MagickFalse;
                  break;
                }
              graphic_context[n]->decorate=(DecorationType) decorate;
              break;
            }
          if (LocaleCompare("density",keyword) == 0)
            {
              (void) GetNextToken(q,&q,extent,token);
              (void) CloneString(&graphic_context[n]->density,token);
              break;
            }
          if (LocaleCompare("direction",keyword) == 0)
            {
              ssize_t
                direction;
  
              (void) GetNextToken(q,&q,extent,token);
              direction=ParseCommandOption(MagickDirectionOptions,MagickFalse,
                token);
              if (direction == -1)
                status=MagickFalse;
              else
                graphic_context[n]->direction=(DirectionType) direction;
              break;
            }
          status=MagickFalse;
          break;
        }
        case 'e':
        case 'E':
        {
          if (LocaleCompare("ellipse",keyword) == 0)
            {
              primitive_type=EllipsePrimitive;
              break;
            }
          if (LocaleCompare("encoding",keyword) == 0)
            {
              (void) GetNextToken(q,&q,extent,token);
              (void) CloneString(&graphic_context[n]->encoding,token);
              break;
            }
          status=MagickFalse;
          break;
        }
        case 'f':
        case 'F':
        {
          if (LocaleCompare("fill",keyword) == 0)
            {
              (void) GetNextToken(q,&q,extent,token);
              if (graphic_context[n]->clip_path != MagickFalse)
                break;
              (void) FormatLocaleString(pattern,MagickPathExtent,"%s",token);
              if (GetImageArtifact(image,pattern) != (const char *) NULL)
                (void) DrawPatternPath(image,draw_info,token,
                  &graphic_context[n]->fill_pattern,exception);
              else
                {
                  status&=QueryColorCompliance(token,AllCompliance,
                    &graphic_context[n]->fill,exception);
                  if (graphic_context[n]->fill_alpha != OpaqueAlpha)
                    graphic_context[n]->fill.alpha=graphic_context[n]->fill_alpha;
                }
              break;
            }
          if (LocaleCompare("fill-opacity",keyword) == 0)
            {
              double
                opacity;
  
              (void) GetNextToken(q,&q,extent,token);
              if (graphic_context[n]->clip_path != MagickFalse)
                break;
              factor=strchr(token,'%') != (char *) NULL ? 0.01 : 1.0;
              opacity=MagickMin(MagickMax(factor*
                GetDrawValue(token,&next_token),0.0),1.0);
              if (token == next_token)
                ThrowPointExpectedException(token,exception);
              if (graphic_context[n]->compliance == SVGCompliance)
                graphic_context[n]->fill_alpha*=opacity;
              else
                graphic_context[n]->fill_alpha=QuantumRange*opacity;
              if (graphic_context[n]->fill.alpha != TransparentAlpha)
                graphic_context[n]->fill.alpha=graphic_context[n]->fill_alpha;
              else
                graphic_context[n]->fill.alpha=(MagickRealType)
                  ClampToQuantum(QuantumRange*(1.0-opacity));
              break;
            }
          if (LocaleCompare("fill-rule",keyword) == 0)
            {
              ssize_t
                fill_rule;
  
              (void) GetNextToken(q,&q,extent,token);
              fill_rule=ParseCommandOption(MagickFillRuleOptions,MagickFalse,
                token);
              if (fill_rule == -1)
                {
                  status=MagickFalse;
                  break;
                }
              graphic_context[n]->fill_rule=(FillRule) fill_rule;
              break;
            }
          if (LocaleCompare("font",keyword) == 0)
            {
              (void) GetNextToken(q,&q,extent,token);
              (void) CloneString(&graphic_context[n]->font,token);
              if (LocaleCompare("none",token) == 0)
                graphic_context[n]->font=(char *) RelinquishMagickMemory(
                  graphic_context[n]->font);
              break;
            }
          if (LocaleCompare("font-family",keyword) == 0)
            {
              (void) GetNextToken(q,&q,extent,token);
              (void) CloneString(&graphic_context[n]->family,token);
              break;
            }
          if (LocaleCompare("font-size",keyword) == 0)
            {
              (void) GetNextToken(q,&q,extent,token);
              graphic_context[n]->pointsize=GetDrawValue(token,&next_token);
              if (token == next_token)
                ThrowPointExpectedException(token,exception);
              break;
            }
          if (LocaleCompare("font-stretch",keyword) == 0)
            {
              ssize_t
                stretch;
  
              (void) GetNextToken(q,&q,extent,token);
              stretch=ParseCommandOption(MagickStretchOptions,MagickFalse,token);
              if (stretch == -1)
                {
                  status=MagickFalse;
                  break;
                }
              graphic_context[n]->stretch=(StretchType) stretch;
              break;
            }
          if (LocaleCompare("font-style",keyword) == 0)
            {
              ssize_t
                style;
  
              (void) GetNextToken(q,&q,extent,token);
              style=ParseCommandOption(MagickStyleOptions,MagickFalse,token);
              if (style == -1)
                {
                  status=MagickFalse;
                  break;
                }
              graphic_context[n]->style=(StyleType) style;
              break;
            }
          if (LocaleCompare("font-weight",keyword) == 0)
            {
              ssize_t
                weight;
  
              (void) GetNextToken(q,&q,extent,token);
              weight=ParseCommandOption(MagickWeightOptions,MagickFalse,token);
              if (weight == -1)
                weight=(ssize_t) StringToUnsignedLong(token);
              graphic_context[n]->weight=(size_t) weight;
              break;
            }
          status=MagickFalse;
          break;
        }
        case 'g':
        case 'G':
        {
          if (LocaleCompare("gradient-units",keyword) == 0)
            {
              (void) GetNextToken(q,&q,extent,token);
              break;
            }
          if (LocaleCompare("gravity",keyword) == 0)
            {
              ssize_t
                gravity;
  
              (void) GetNextToken(q,&q,extent,token);
              gravity=ParseCommandOption(MagickGravityOptions,MagickFalse,token);
              if (gravity == -1)
                {
                  status=MagickFalse;
                  break;
                }
              graphic_context[n]->gravity=(GravityType) gravity;
              break;
            }
          status=MagickFalse;
          break;
        }
        case 'i':
        case 'I':
        {
          if (LocaleCompare("image",keyword) == 0)
            {
              ssize_t
                compose;
  
              primitive_type=ImagePrimitive;
              (void) GetNextToken(q,&q,extent,token);
              compose=ParseCommandOption(MagickComposeOptions,MagickFalse,token);
              if (compose == -1)
                {
                  status=MagickFalse;
                  break;
                }
              graphic_context[n]->compose=(CompositeOperator) compose;
              break;
            }
          if (LocaleCompare("interline-spacing",keyword) == 0)
            {
              (void) GetNextToken(q,&q,extent,token);
              graphic_context[n]->interline_spacing=GetDrawValue(token,
                &next_token);
              if (token == next_token)
                ThrowPointExpectedException(token,exception);
              break;
            }
          if (LocaleCompare("interword-spacing",keyword) == 0)
            {
              (void) GetNextToken(q,&q,extent,token);
              graphic_context[n]->interword_spacing=GetDrawValue(token,
                &next_token);
              if (token == next_token)
                ThrowPointExpectedException(token,exception);
              break;
            }
          status=MagickFalse;
          break;
        }
        case 'k':
        case 'K':
        {
          if (LocaleCompare("kerning",keyword) == 0)
            {
              (void) GetNextToken(q,&q,extent,token);
              graphic_context[n]->kerning=GetDrawValue(token,&next_token);
              if (token == next_token)
                ThrowPointExpectedException(token,exception);
              break;
            }
          status=MagickFalse;
          break;
        }
        case 'l':
        case 'L':
        {
          if (LocaleCompare("letter-spacing",keyword) == 0)
            {
              (void) GetNextToken(q,&q,extent,token);
              if (IsPoint(token) == MagickFalse)
                break;
              clone_info=CloneDrawInfo((ImageInfo *) NULL,graphic_context[n]);
              clone_info->text=AcquireString(" ");
              status&=GetTypeMetrics(image,clone_info,&metrics,exception);
              graphic_context[n]->kerning=metrics.width*
                GetDrawValue(token,&next_token);
              clone_info=DestroyDrawInfo(clone_info);
              if (token == next_token)
                ThrowPointExpectedException(token,exception);
              break;
            }
          if (LocaleCompare("line",keyword) == 0)
            {
              primitive_type=LinePrimitive;
              break;
            }
          status=MagickFalse;
          break;
        }
        case 'm':
        case 'M':
        {
          if (LocaleCompare("mask",keyword) == 0)
            {
              const char
                *mask_path;
  
              /*
                Take a node from within the MVG document, and duplicate it here.
              */
              (void) GetNextToken(q,&q,extent,token);
              mask_path=(const char *) GetValueFromSplayTree(macros,token);
              if (mask_path != (const char *) NULL)
                {
                  if (graphic_context[n]->composite_mask != (Image *) NULL)
                    graphic_context[n]->composite_mask=
                      DestroyImage(graphic_context[n]->composite_mask);
                  graphic_context[n]->composite_mask=DrawCompositeMask(image,
                    graphic_context[n],token,mask_path,exception);
                  if (graphic_context[n]->compliance != SVGCompliance)
                    status=SetImageMask(image,CompositePixelMask,
                      graphic_context[n]->composite_mask,exception);
                }
              break;
            }
          status=MagickFalse;
          break;
        }
        case 'o':
        case 'O':
        {
          if (LocaleCompare("offset",keyword) == 0)
            {
              (void) GetNextToken(q,&q,extent,token);
              break;
            }
          if (LocaleCompare("opacity",keyword) == 0)
            {
              double
                opacity;
  
              (void) GetNextToken(q,&q,extent,token);
              if (graphic_context[n]->clip_path != MagickFalse)
                break;
              factor=strchr(token,'%') != (char *) NULL ? 0.01 : 1.0;
              opacity=MagickMin(MagickMax(factor*
                GetDrawValue(token,&next_token),0.0),1.0);
              if (token == next_token)
                ThrowPointExpectedException(token,exception);
              if (graphic_context[n]->compliance == SVGCompliance)
                {
                  graphic_context[n]->fill_alpha*=opacity;
                  graphic_context[n]->stroke_alpha*=opacity;
                }
              else
                {
                  graphic_context[n]->fill_alpha=QuantumRange*opacity;
                  graphic_context[n]->stroke_alpha=QuantumRange*opacity;
                }
              break;
            }
          status=MagickFalse;
          break;
        }
        case 'p':
        case 'P':
        {
          if (LocaleCompare("path",keyword) == 0)
            {
              primitive_type=PathPrimitive;
              break;
            }
          if (LocaleCompare("point",keyword) == 0)
            {
              primitive_type=PointPrimitive;
              break;
            }
          if (LocaleCompare("polyline",keyword) == 0)
            {
              primitive_type=PolylinePrimitive;
              break;
            }
          if (LocaleCompare("polygon",keyword) == 0)
            {
              primitive_type=PolygonPrimitive;
              break;
            }
          if (LocaleCompare("pop",keyword) == 0)
            {
              if (GetNextToken(q,&q,extent,token) < 1)
                break;
              if (LocaleCompare("class",token) == 0)
                break;
              if (LocaleCompare("clip-path",token) == 0)
                break;
              if (LocaleCompare("defs",token) == 0)
                {
                  defsDepth--;
                  graphic_context[n]->render=defsDepth > 0 ? MagickFalse :
                    MagickTrue;
                  break;
                }
              if (LocaleCompare("gradient",token) == 0)
                break;
              if (LocaleCompare("graphic-context",token) == 0)
                {
                  if (n <= 0)
                    {
                      (void) ThrowMagickException(exception,GetMagickModule(),
                        DrawError,"UnbalancedGraphicContextPushPop","`%s'",token);
                      status=MagickFalse;
                      n=0;
                      break;
                    }
                  if ((graphic_context[n]->clip_mask != (char *) NULL) &&
                      (graphic_context[n]->compliance != SVGCompliance))
                    if (LocaleCompare(graphic_context[n]->clip_mask,
                        graphic_context[n-1]->clip_mask) != 0)
                      status=SetImageMask(image,WritePixelMask,(Image *) NULL,
                        exception);
                  graphic_context[n]=DestroyDrawInfo(graphic_context[n]);
                  n--;
                  break;
                }
              if (LocaleCompare("mask",token) == 0)
                break;
              if (LocaleCompare("pattern",token) == 0)
                break;
              if (LocaleCompare("symbol",token) == 0)
                {
                  symbolDepth--;
                  graphic_context[n]->render=symbolDepth > 0 ? MagickFalse :
                    MagickTrue;
                  break;
                }
              status=MagickFalse;
              break;
            }
          if (LocaleCompare("push",keyword) == 0)
            {
              if (GetNextToken(q,&q,extent,token) < 1)
                break;
              if (LocaleCompare("class",token) == 0)
                {
                  /*
                    Class context.
                  */
                  for (p=q; *q != '\0'; )
                  {
                    if (GetNextToken(q,&q,extent,token) < 1)
                      break;
                    if (LocaleCompare(token,"pop") != 0)
                      continue;
                    (void) GetNextToken(q,(const char **) NULL,extent,token);
                    if (LocaleCompare(token,"class") != 0)
                      continue;
                    break;
                  }
                  (void) GetNextToken(q,&q,extent,token);
                  break;
                }
              if (LocaleCompare("clip-path",token) == 0)
                {
                  (void) GetNextToken(q,&q,extent,token);
                  for (p=q; *q != '\0'; )
                  {
                  	if (GetNextToken(q,&q,extent,token) < 1)
                      break;
                  	if (LocaleCompare(token,"pop") != 0)
                  		continue;
                  	(void) GetNextToken(q,(const char **) NULL,extent,token);
                  	if (LocaleCompare(token,"clip-path") != 0)
                  		continue;
                  	break;
                  }
                  if ((q == (char *) NULL) || (p == (char *) NULL) || ((q-4) < p))
                    {
                      status=MagickFalse;
                      break;
                    }
                  (void) GetNextToken(q,&q,extent,token);
                  break;
                }
              if (LocaleCompare("defs",token) == 0)
                {
                  defsDepth++;
                  graphic_context[n]->render=defsDepth > 0 ? MagickFalse :
                    MagickTrue;
                  break;
                }
              if (LocaleCompare("gradient",token) == 0)
                {
                  char
                    key[2*MagickPathExtent],
                    name[MagickPathExtent],
                    type[MagickPathExtent];
  
                  SegmentInfo
                    segment;
  
                  (void) GetNextToken(q,&q,extent,token);
                  (void) CopyMagickString(name,token,MagickPathExtent);
                  (void) GetNextToken(q,&q,extent,token);
                  (void) CopyMagickString(type,token,MagickPathExtent);
                  (void) GetNextToken(q,&q,extent,token);
                  segment.x1=GetDrawValue(token,&next_token);
                  if (token == next_token)
                    ThrowPointExpectedException(token,exception);
                  (void) GetNextToken(q,&q,extent,token);
                  if (*token == ',')
                    (void) GetNextToken(q,&q,extent,token);
                  segment.y1=GetDrawValue(token,&next_token);
                  if (token == next_token)
                    ThrowPointExpectedException(token,exception);
                  (void) GetNextToken(q,&q,extent,token);
                  if (*token == ',')
                    (void) GetNextToken(q,&q,extent,token);
                  segment.x2=GetDrawValue(token,&next_token);
                  if (token == next_token)
                    ThrowPointExpectedException(token,exception);
                  (void) GetNextToken(q,&q,extent,token);
                  if (*token == ',')
                    (void) GetNextToken(q,&q,extent,token);
                  segment.y2=GetDrawValue(token,&next_token);
                  if (token == next_token)
                    ThrowPointExpectedException(token,exception);
                  if (LocaleCompare(type,"radial") == 0)
                    {
                      (void) GetNextToken(q,&q,extent,token);
                      if (*token == ',')
                        (void) GetNextToken(q,&q,extent,token);
                    }
                  for (p=q; *q != '\0'; )
                  {
                    if (GetNextToken(q,&q,extent,token) < 1)
                      break;
                    if (LocaleCompare(token,"pop") != 0)
                      continue;
                    (void) GetNextToken(q,(const char **) NULL,extent,token);
                    if (LocaleCompare(token,"gradient") != 0)
                      continue;
                    break;
                  }
                  if ((q == (char *) NULL) || (p == (char *) NULL) || ((q-4) < p))
                    {
                      status=MagickFalse;
                      break;
                    }
                  (void) CopyMagickString(token,p,(size_t) (q-p-4+1));
                  bounds.x1=graphic_context[n]->affine.sx*segment.x1+
                    graphic_context[n]->affine.ry*segment.y1+
                    graphic_context[n]->affine.tx;
                  bounds.y1=graphic_context[n]->affine.rx*segment.x1+
                    graphic_context[n]->affine.sy*segment.y1+
                    graphic_context[n]->affine.ty;
                  bounds.x2=graphic_context[n]->affine.sx*segment.x2+
                    graphic_context[n]->affine.ry*segment.y2+
                    graphic_context[n]->affine.tx;
                  bounds.y2=graphic_context[n]->affine.rx*segment.x2+
                    graphic_context[n]->affine.sy*segment.y2+
                    graphic_context[n]->affine.ty;
                  (void) FormatLocaleString(key,MagickPathExtent,"%s",name);
                  (void) SetImageArtifact(image,key,token);
                  (void) FormatLocaleString(key,MagickPathExtent,"%s-type",name);
                  (void) SetImageArtifact(image,key,type);
                  (void) FormatLocaleString(key,MagickPathExtent,"%s-geometry",
                    name);
                  (void) FormatLocaleString(geometry,MagickPathExtent,
                    "%gx%g%+.15g%+.15g",
                    MagickMax(fabs(bounds.x2-bounds.x1+1.0),1.0),
                    MagickMax(fabs(bounds.y2-bounds.y1+1.0),1.0),
                    bounds.x1,bounds.y1);
                  (void) SetImageArtifact(image,key,geometry);
                  (void) GetNextToken(q,&q,extent,token);
                  break;
                }
              if (LocaleCompare("graphic-context",token) == 0)
                {
                  n++;
                  graphic_context=(DrawInfo **) ResizeQuantumMemory(
                    graphic_context,(size_t) (n+1),sizeof(*graphic_context));
                  if (graphic_context == (DrawInfo **) NULL)
                    {
                      (void) ThrowMagickException(exception,GetMagickModule(),
                        ResourceLimitError,"MemoryAllocationFailed","`%s'",
                        image->filename);
                      break;
                    }
                  graphic_context[n]=CloneDrawInfo((ImageInfo *) NULL,
                    graphic_context[n-1]);
                  if (*q == '"')
                    {
                      (void) GetNextToken(q,&q,extent,token);
                      (void) CloneString(&graphic_context[n]->id,token);
                    }
                  break;
                }
              if (LocaleCompare("mask",token) == 0)
                {
                  (void) GetNextToken(q,&q,extent,token);
                  break;
                }
              if (LocaleCompare("pattern",token) == 0)
                {
                  char
                    key[2*MagickPathExtent],
                    name[MagickPathExtent];
  
                  RectangleInfo
                    bounds;
  
                  (void) GetNextToken(q,&q,extent,token);
                  (void) CopyMagickString(name,token,MagickPathExtent);
                  (void) GetNextToken(q,&q,extent,token);
                  bounds.x=(ssize_t) ceil(GetDrawValue(token,&next_token)-0.5);
                  if (token == next_token)
                    ThrowPointExpectedException(token,exception);
                  (void) GetNextToken(q,&q,extent,token);
                  if (*token == ',')
                    (void) GetNextToken(q,&q,extent,token);
                  bounds.y=(ssize_t) ceil(GetDrawValue(token,&next_token)-0.5);
                  if (token == next_token)
                    ThrowPointExpectedException(token,exception);
                  (void) GetNextToken(q,&q,extent,token);
                  if (*token == ',')
                    (void) GetNextToken(q,&q,extent,token);
                  bounds.width=(size_t) floor(GetDrawValue(token,&next_token)+
                    0.5);
                  if (token == next_token)
                    ThrowPointExpectedException(token,exception);
                  (void) GetNextToken(q,&q,extent,token);
                  if (*token == ',')
                    (void) GetNextToken(q,&q,extent,token);
                  bounds.height=(size_t) floor(GetDrawValue(token,&next_token)+
                    0.5);
                  if (token == next_token)
                    ThrowPointExpectedException(token,exception);
                  for (p=q; *q != '\0'; )
                  {
                    if (GetNextToken(q,&q,extent,token) < 1)
                      break;
                    if (LocaleCompare(token,"pop") != 0)
                      continue;
                    (void) GetNextToken(q,(const char **) NULL,extent,token);
                    if (LocaleCompare(token,"pattern") != 0)
                      continue;
                    break;
                  }
                  if ((q == (char *) NULL) || (p == (char *) NULL) || ((q-4) < p))
                    {
                      status=MagickFalse;
                      break;
                    }
                  (void) CopyMagickString(token,p,(size_t) (q-p-4+1));
                  (void) FormatLocaleString(key,MagickPathExtent,"%s",name);
                  (void) SetImageArtifact(image,key,token);
                  (void) FormatLocaleString(key,MagickPathExtent,"%s-geometry",
                    name);
                  (void) FormatLocaleString(geometry,MagickPathExtent,
                    "%.20gx%.20g%+.20g%+.20g",(double) bounds.width,(double)
                    bounds.height,(double) bounds.x,(double) bounds.y);
                  (void) SetImageArtifact(image,key,geometry);
                  (void) GetNextToken(q,&q,extent,token);
                  break;
                }
              if (LocaleCompare("symbol",token) == 0)
                {
                  symbolDepth++;
                  graphic_context[n]->render=symbolDepth > 0 ? MagickFalse :
                    MagickTrue;
                  break;
                }
              status=MagickFalse;
              break;
            }
          status=MagickFalse;
          break;
        }
        case 'r':
        case 'R':
        {
          if (LocaleCompare("rectangle",keyword) == 0)
            {
              primitive_type=RectanglePrimitive;
              break;
            }
          if (LocaleCompare("rotate",keyword) == 0)
            {
              (void) GetNextToken(q,&q,extent,token);
              angle=GetDrawValue(token,&next_token);
              if (token == next_token)
                ThrowPointExpectedException(token,exception);
              affine.sx=cos(DegreesToRadians(fmod((double) angle,360.0)));
              affine.rx=sin(DegreesToRadians(fmod((double) angle,360.0)));
              affine.ry=(-sin(DegreesToRadians(fmod((double) angle,360.0))));
              affine.sy=cos(DegreesToRadians(fmod((double) angle,360.0)));
              break;
            }
          if (LocaleCompare("roundRectangle",keyword) == 0)
            {
              primitive_type=RoundRectanglePrimitive;
              break;
            }
          status=MagickFalse;
          break;
        }
        case 's':
        case 'S':
        {
          if (LocaleCompare("scale",keyword) == 0)
            {
              (void) GetNextToken(q,&q,extent,token);
              affine.sx=GetDrawValue(token,&next_token);
              if (token == next_token)
                ThrowPointExpectedException(token,exception);
              (void) GetNextToken(q,&q,extent,token);
              if (*token == ',')
                (void) GetNextToken(q,&q,extent,token);
              affine.sy=GetDrawValue(token,&next_token);
              if (token == next_token)
                ThrowPointExpectedException(token,exception);
              break;
            }
          if (LocaleCompare("skewX",keyword) == 0)
            {
              (void) GetNextToken(q,&q,extent,token);
              angle=GetDrawValue(token,&next_token);
              if (token == next_token)
                ThrowPointExpectedException(token,exception);
              affine.ry=sin(DegreesToRadians(angle));
              break;
            }
          if (LocaleCompare("skewY",keyword) == 0)
            {
              (void) GetNextToken(q,&q,extent,token);
              angle=GetDrawValue(token,&next_token);
              if (token == next_token)
                ThrowPointExpectedException(token,exception);
              affine.rx=(-tan(DegreesToRadians(angle)/2.0));
              break;
            }
          if (LocaleCompare("stop-color",keyword) == 0)
            {
              PixelInfo
                stop_color;
  
              number_stops++;
              if (number_stops == 1)
                stops=(StopInfo *) AcquireQuantumMemory(2,sizeof(*stops));
              else
                if (number_stops > 2)
                  stops=(StopInfo *) ResizeQuantumMemory(stops,number_stops,
                    sizeof(*stops));
              if (stops == (StopInfo *) NULL)
                {
                  (void) ThrowMagickException(exception,GetMagickModule(),
                    ResourceLimitError,"MemoryAllocationFailed","`%s'",
                    image->filename);
                  break;
                }
              (void) GetNextToken(q,&q,extent,token);
              status&=QueryColorCompliance(token,AllCompliance,&stop_color,
                exception);
              stops[number_stops-1].color=stop_color;
              (void) GetNextToken(q,&q,extent,token);
              factor=strchr(token,'%') != (char *) NULL ? 0.01 : 1.0;
              stops[number_stops-1].offset=factor*GetDrawValue(token,
                &next_token);
              if (token == next_token)
                ThrowPointExpectedException(token,exception);
              break;
            }
          if (LocaleCompare("stroke",keyword) == 0)
            {
              (void) GetNextToken(q,&q,extent,token);
              if (graphic_context[n]->clip_path != MagickFalse)
                break;
              (void) FormatLocaleString(pattern,MagickPathExtent,"%s",token);
              if (GetImageArtifact(image,pattern) != (const char *) NULL)
                (void) DrawPatternPath(image,draw_info,token,
                  &graphic_context[n]->stroke_pattern,exception);
              else
                {
                  status&=QueryColorCompliance(token,AllCompliance,
                    &graphic_context[n]->stroke,exception);
                  if (graphic_context[n]->stroke_alpha != OpaqueAlpha)
                    graphic_context[n]->stroke.alpha=
                      graphic_context[n]->stroke_alpha;
                 }
              break;
            }
          if (LocaleCompare("stroke-antialias",keyword) == 0)
            {
              (void) GetNextToken(q,&q,extent,token);
              graphic_context[n]->stroke_antialias=StringToLong(token) != 0 ?
                MagickTrue : MagickFalse;
              break;
            }
          if (LocaleCompare("stroke-dasharray",keyword) == 0)
            {
              if (graphic_context[n]->dash_pattern != (double *) NULL)
                graphic_context[n]->dash_pattern=(double *)
                  RelinquishMagickMemory(graphic_context[n]->dash_pattern);
              if (IsPoint(q) != MagickFalse)
                {
                  const char
                    *r;
  
                  r=q;
                  (void) GetNextToken(r,&r,extent,token);
                  if (*token == ',')
                    (void) GetNextToken(r,&r,extent,token);
                  for (x=0; IsPoint(token) != MagickFalse; x++)
                  {
                    (void) GetNextToken(r,&r,extent,token);
                    if (*token == ',')
                      (void) GetNextToken(r,&r,extent,token);
                  }
                  graphic_context[n]->dash_pattern=(double *)
                    AcquireQuantumMemory((size_t) (2*x+2),
                    sizeof(*graphic_context[n]->dash_pattern));
                  if (graphic_context[n]->dash_pattern == (double *) NULL)
                    {
                      (void) ThrowMagickException(exception,GetMagickModule(),
                        ResourceLimitError,"MemoryAllocationFailed","`%s'",
                        image->filename);
                      status=MagickFalse;
                      break;
                    }
                  (void) memset(graphic_context[n]->dash_pattern,0,(size_t)
                    (2*x+2)*sizeof(*graphic_context[n]->dash_pattern));
                  for (j=0; j < x; j++)
                  {
                    (void) GetNextToken(q,&q,extent,token);
                    if (*token == ',')
                      (void) GetNextToken(q,&q,extent,token);
                    graphic_context[n]->dash_pattern[j]=GetDrawValue(token,
                      &next_token);
                    if (token == next_token)
                      ThrowPointExpectedException(token,exception);
                    if (graphic_context[n]->dash_pattern[j] < 0.0)
                      status=MagickFalse;
                  }
                  if ((x & 0x01) != 0)
                    for ( ; j < (2*x); j++)
                      graphic_context[n]->dash_pattern[j]=
                        graphic_context[n]->dash_pattern[j-x];
                  graphic_context[n]->dash_pattern[j]=0.0;
                  break;
                }
              (void) GetNextToken(q,&q,extent,token);
              break;
            }
          if (LocaleCompare("stroke-dashoffset",keyword) == 0)
            {
              (void) GetNextToken(q,&q,extent,token);
              graphic_context[n]->dash_offset=GetDrawValue(token,&next_token);
              if (token == next_token)
                ThrowPointExpectedException(token,exception);
              break;
            }
          if (LocaleCompare("stroke-linecap",keyword) == 0)
            {
              ssize_t
                linecap;
  
              (void) GetNextToken(q,&q,extent,token);
              linecap=ParseCommandOption(MagickLineCapOptions,MagickFalse,token);
              if (linecap == -1)
                {
                  status=MagickFalse;
                  break;
                }
              graphic_context[n]->linecap=(LineCap) linecap;
              break;
            }
          if (LocaleCompare("stroke-linejoin",keyword) == 0)
            {
              ssize_t
                linejoin;
  
              (void) GetNextToken(q,&q,extent,token);
              linejoin=ParseCommandOption(MagickLineJoinOptions,MagickFalse,
                token);
              if (linejoin == -1)
                {
                  status=MagickFalse;
                  break;
                }
              graphic_context[n]->linejoin=(LineJoin) linejoin;
              break;
            }
          if (LocaleCompare("stroke-miterlimit",keyword) == 0)
            {
              (void) GetNextToken(q,&q,extent,token);
              graphic_context[n]->miterlimit=StringToUnsignedLong(token);
              break;
            }
          if (LocaleCompare("stroke-opacity",keyword) == 0)
            {
              double
                opacity;
  
              (void) GetNextToken(q,&q,extent,token);
              if (graphic_context[n]->clip_path != MagickFalse)
                break;
              factor=strchr(token,'%') != (char *) NULL ? 0.01 : 1.0;
              opacity=MagickMin(MagickMax(factor*
                GetDrawValue(token,&next_token),0.0),1.0);
              if (token == next_token)
                ThrowPointExpectedException(token,exception);
              if (graphic_context[n]->compliance == SVGCompliance)
                graphic_context[n]->stroke_alpha*=opacity;
              else
                graphic_context[n]->stroke_alpha=QuantumRange*opacity;
              if (graphic_context[n]->stroke.alpha != TransparentAlpha)
                graphic_context[n]->stroke.alpha=graphic_context[n]->stroke_alpha;
              else
                graphic_context[n]->stroke.alpha=(MagickRealType)
                  ClampToQuantum(QuantumRange*(1.0-opacity));
              break;
            }
          if (LocaleCompare("stroke-width",keyword) == 0)
            {
              (void) GetNextToken(q,&q,extent,token);
              if (graphic_context[n]->clip_path != MagickFalse)
                break;
              graphic_context[n]->stroke_width=GetDrawValue(token,&next_token);
              if (token == next_token)
                ThrowPointExpectedException(token,exception);
              break;
            }
          status=MagickFalse;
          break;
        }
        case 't':
        case 'T':
        {
          if (LocaleCompare("text",keyword) == 0)
            {
              primitive_type=TextPrimitive;
              cursor=0.0;
              break;
            }
          if (LocaleCompare("text-align",keyword) == 0)
            {
              ssize_t
                align;
  
              (void) GetNextToken(q,&q,extent,token);
              align=ParseCommandOption(MagickAlignOptions,MagickFalse,token);
              if (align == -1)
                {
                  status=MagickFalse;
                  break;
                }
              graphic_context[n]->align=(AlignType) align;
              break;
            }
          if (LocaleCompare("text-anchor",keyword) == 0)
            {
              ssize_t
                align;
  
              (void) GetNextToken(q,&q,extent,token);
              align=ParseCommandOption(MagickAlignOptions,MagickFalse,token);
              if (align == -1)
                {
                  status=MagickFalse;
                  break;
                }
              graphic_context[n]->align=(AlignType) align;
              break;
            }
          if (LocaleCompare("text-antialias",keyword) == 0)
            {
              (void) GetNextToken(q,&q,extent,token);
              graphic_context[n]->text_antialias=StringToLong(token) != 0 ?
                MagickTrue : MagickFalse;
              break;
            }
          if (LocaleCompare("text-undercolor",keyword) == 0)
            {
              (void) GetNextToken(q,&q,extent,token);
              status&=QueryColorCompliance(token,AllCompliance,
                &graphic_context[n]->undercolor,exception);
              break;
            }
          if (LocaleCompare("translate",keyword) == 0)
            {
              (void) GetNextToken(q,&q,extent,token);
              affine.tx=GetDrawValue(token,&next_token);
              if (token == next_token)
                ThrowPointExpectedException(token,exception);
              (void) GetNextToken(q,&q,extent,token);
              if (*token == ',')
                (void) GetNextToken(q,&q,extent,token);
              affine.ty=GetDrawValue(token,&next_token);
              if (token == next_token)
                ThrowPointExpectedException(token,exception);
              cursor=0.0;
              break;
            }
          status=MagickFalse;
          break;
        }
        case 'u':
        case 'U':
        {
          if (LocaleCompare("use",keyword) == 0)
            {
              const char
                *use;
  
              /*
                Get a macro from the MVG document, and "use" it here.
              */
              (void) GetNextToken(q,&q,extent,token);
              use=(const char *) GetValueFromSplayTree(macros,token);
              if (use != (const char *) NULL)
                {
                  clone_info=CloneDrawInfo((ImageInfo *) NULL,graphic_context[n]);
                  (void) CloneString(&clone_info->primitive,use);
                  status=RenderMVGContent(image,clone_info,depth+1,exception);
                  clone_info=DestroyDrawInfo(clone_info);
                }
              break;
            }
          status=MagickFalse;
          break;
        }
        case 'v':
        case 'V':
        {
          if (LocaleCompare("viewbox",keyword) == 0)
            {
              (void) GetNextToken(q,&q,extent,token);
              graphic_context[n]->viewbox.x=(ssize_t) ceil(GetDrawValue(token,
                &next_token)-0.5);
              if (token == next_token)
                ThrowPointExpectedException(token,exception);
              (void) GetNextToken(q,&q,extent,token);
              if (*token == ',')
                (void) GetNextToken(q,&q,extent,token);
              graphic_context[n]->viewbox.y=(ssize_t) ceil(GetDrawValue(token,
                &next_token)-0.5);
              if (token == next_token)
                ThrowPointExpectedException(token,exception);
              (void) GetNextToken(q,&q,extent,token);
              if (*token == ',')
                (void) GetNextToken(q,&q,extent,token);
              graphic_context[n]->viewbox.width=(size_t) floor(GetDrawValue(
                token,&next_token)+0.5);
              if (token == next_token)
                ThrowPointExpectedException(token,exception);
              (void) GetNextToken(q,&q,extent,token);
              if (*token == ',')
                (void) GetNextToken(q,&q,extent,token);
              graphic_context[n]->viewbox.height=(size_t) floor(GetDrawValue(
                token,&next_token)+0.5);
              if (token == next_token)
                ThrowPointExpectedException(token,exception);
              break;
            }
          status=MagickFalse;
          break;
        }
        case 'w':
        case 'W':
        {
          if (LocaleCompare("word-spacing",keyword) == 0)
            {
              (void) GetNextToken(q,&q,extent,token);
              graphic_context[n]->interword_spacing=GetDrawValue(token,
                &next_token);
              if (token == next_token)
                ThrowPointExpectedException(token,exception);
              break;
            }
          status=MagickFalse;
          break;
        }
        default:
        {
          status=MagickFalse;
          break;
        }
      }
      if (status == MagickFalse)
        break;
      if ((fabs(affine.sx-1.0) >= MagickEpsilon) ||
          (fabs(affine.rx) >= MagickEpsilon) || (fabs(affine.ry) >= MagickEpsilon) ||
          (fabs(affine.sy-1.0) >= MagickEpsilon) ||
          (fabs(affine.tx) >= MagickEpsilon) || (fabs(affine.ty) >= MagickEpsilon))
        {
          graphic_context[n]->affine.sx=current.sx*affine.sx+current.ry*affine.rx;
          graphic_context[n]->affine.rx=current.rx*affine.sx+current.sy*affine.rx;
          graphic_context[n]->affine.ry=current.sx*affine.ry+current.ry*affine.sy;
          graphic_context[n]->affine.sy=current.rx*affine.ry+current.sy*affine.sy;
          graphic_context[n]->affine.tx=current.sx*affine.tx+current.ry*affine.ty+
            current.tx;
          graphic_context[n]->affine.ty=current.rx*affine.tx+current.sy*affine.ty+
            current.ty;
        }
      if (primitive_type == UndefinedPrimitive)
        {
          if (*q == '\0')
            {
              if (number_stops > 1)
                {
                  GradientType
                    type;
  
                type=LinearGradient;
                if (draw_info->gradient.type == RadialGradient)
                  type=RadialGradient;
                (void) GradientImage(image,type,PadSpread,stops,number_stops,
                  exception);
               }
             if (number_stops > 0)
               stops=(StopInfo *) RelinquishMagickMemory(stops);
            }
          if ((image->debug != MagickFalse) && (q > p))
            (void) LogMagickEvent(DrawEvent,GetMagickModule(),"  %.*s",(int)
              (q-p-1),p);
          continue;
        }
      /*
        Parse the primitive attributes.
      */
      for (i=0; primitive_info[i].primitive != UndefinedPrimitive; i++)
        if ((primitive_info[i].primitive == TextPrimitive) ||
            (primitive_info[i].primitive == ImagePrimitive))
          if (primitive_info[i].text != (char *) NULL)
            primitive_info[i].text=DestroyString(primitive_info[i].text);
      i=0;
      mvg_info.offset=i;
      j=0;
      primitive_info[0].point.x=0.0;
      primitive_info[0].point.y=0.0;
      primitive_info[0].coordinates=0;
      primitive_info[0].method=FloodfillMethod;
      primitive_info[0].closed_subpath=MagickFalse;
      for (x=0; *q != '\0'; x++)
      {
        /*
          Define points.
        */
        if (IsPoint(q) == MagickFalse)
          break;
        (void) GetNextToken(q,&q,extent,token);
        point.x=GetDrawValue(token,&next_token);
        if (token == next_token)
          ThrowPointExpectedException(token,exception);
        (void) GetNextToken(q,&q,extent,token);
        if (*token == ',')
          (void) GetNextToken(q,&q,extent,token);
        point.y=GetDrawValue(token,&next_token);
        if (token == next_token)
          ThrowPointExpectedException(token,exception);
        (void) GetNextToken(q,(const char **) NULL,extent,token);
        if (*token == ',')
          (void) GetNextToken(q,&q,extent,token);
        primitive_info[i].primitive=primitive_type;
        primitive_info[i].point=point;
        primitive_info[i].coordinates=0;
        primitive_info[i].method=FloodfillMethod;
        primitive_info[i].closed_subpath=MagickFalse;
        i++;
        mvg_info.offset=i;
        if (i < (ssize_t) number_points)
          continue;
        status&=CheckPrimitiveExtent(&mvg_info,number_points);
      }
      if (status == MagickFalse)
        break;
      if ((primitive_info[j].primitive == TextPrimitive) ||
          (primitive_info[j].primitive == ImagePrimitive))
        if (primitive_info[j].text != (char *) NULL)
          primitive_info[j].text=DestroyString(primitive_info[j].text);
      primitive_info[j].primitive=primitive_type;
      primitive_info[j].coordinates=(size_t) x;
      primitive_info[j].method=FloodfillMethod;
      primitive_info[j].closed_subpath=MagickFalse;
      /*
        Circumscribe primitive within a circle.
      */
      bounds.x1=primitive_info[j].point.x;
      bounds.y1=primitive_info[j].point.y;
      bounds.x2=primitive_info[j].point.x;
      bounds.y2=primitive_info[j].point.y;
      for (k=1; k < (ssize_t) primitive_info[j].coordinates; k++)
      {
        point=primitive_info[j+k].point;
        if (point.x < bounds.x1)
          bounds.x1=point.x;
        if (point.y < bounds.y1)
          bounds.y1=point.y;
        if (point.x > bounds.x2)
          bounds.x2=point.x;
        if (point.y > bounds.y2)
          bounds.y2=point.y;
      }
      /*
        Speculate how many points our primitive might consume.
      */
      coordinates=(double) primitive_info[j].coordinates;
      switch (primitive_type)
      {
        case RectanglePrimitive:
        {
          coordinates*=5.0;
          break;
        }
        case RoundRectanglePrimitive:
        {
          double
            alpha,
            beta,
            radius;
  
          alpha=bounds.x2-bounds.x1;
          beta=bounds.y2-bounds.y1;
          radius=hypot((double) alpha,(double) beta);
          coordinates*=5.0;
          coordinates+=2.0*((size_t) ceil((double) MagickPI*radius))+6.0*
            BezierQuantum+360.0;
          break;
        }
        case BezierPrimitive:
        {
          coordinates=(double) (BezierQuantum*primitive_info[j].coordinates);
          if (primitive_info[j].coordinates > (107*BezierQuantum))
            {
              (void) ThrowMagickException(exception,GetMagickModule(),DrawError,
                "TooManyBezierCoordinates","`%s'",token);
              status=MagickFalse;
              break;
            }
          break;
        }
        case PathPrimitive:
        {
          char
            *s,
            *t;
  
          (void) GetNextToken(q,&q,extent,token);
          coordinates=1.0;
          t=token;
          for (s=token; *s != '\0'; s=t)
          {
            double
              value;
  
            value=GetDrawValue(s,&t);
            (void) value;
            if (s == t)
              {
                t++;
                continue;
              }
            coordinates++;
          }
          for (s=token; *s != '\0'; s++)
            if (strspn(s,"AaCcQqSsTt") != 0)
              coordinates+=(20.0*BezierQuantum)+360.0;
          break;
        }
        case CirclePrimitive:
        case ArcPrimitive:
        case EllipsePrimitive:
        {
          double
            alpha,
            beta,
            radius;
  
          alpha=bounds.x2-bounds.x1;
          beta=bounds.y2-bounds.y1;
          radius=hypot(alpha,beta);
          coordinates=2.0*(ceil(MagickPI*radius))+6.0*BezierQuantum+360.0;
          if (coordinates > (107*BezierQuantum))
            {
              (void) ThrowMagickException(exception,GetMagickModule(),DrawError,
                "TooManyBezierCoordinates","`%s'",token);
              status=MagickFalse;
              break;
            }
          break;
        }
        default:
          break;
      }
      if (status == MagickFalse)
        break;
      if (((size_t) (i+coordinates)) >= number_points)
        {
          /*
            Resize based on speculative points required by primitive.
          */
          number_points+=coordinates+1;
          if (number_points < (size_t) coordinates)
            {
              (void) ThrowMagickException(exception,GetMagickModule(),
                ResourceLimitError,"MemoryAllocationFailed","`%s'",
                image->filename);
              break;
            }
          mvg_info.offset=i;
          status&=CheckPrimitiveExtent(&mvg_info,number_points);
        }
      status&=CheckPrimitiveExtent(&mvg_info,PrimitiveExtentPad);
      if (status == MagickFalse)
        break;
      mvg_info.offset=j;
      switch (primitive_type)
      {
        case PointPrimitive:
        default:
        {
          if (primitive_info[j].coordinates != 1)
            {
              status=MagickFalse;
              break;
            }
          status&=TracePoint(primitive_info+j,primitive_info[j].point);
          i=(ssize_t) (j+primitive_info[j].coordinates);
          break;
        }
        case LinePrimitive:
        {
          double
            dx,
            dy,
            maximum_length;
  
          if (primitive_info[j].coordinates != 2)
            {
              status=MagickFalse;
              break;
            }
          dx=primitive_info[i].point.x-primitive_info[i-1].point.x;
          dy=primitive_info[i].point.y-primitive_info[i-1].point.y;
          maximum_length=hypot(dx,dy);
          if (maximum_length > (MaxBezierCoordinates/100.0))
            ThrowPointExpectedException(keyword,exception);
          status&=TraceLine(primitive_info+j,primitive_info[j].point,
            primitive_info[j+1].point);
          i=(ssize_t) (j+primitive_info[j].coordinates);
          break;
        }
        case RectanglePrimitive:
        {
          if (primitive_info[j].coordinates != 2)
            {
              status=MagickFalse;
              break;
            }
          status&=TraceRectangle(primitive_info+j,primitive_info[j].point,
            primitive_info[j+1].point);
          i=(ssize_t) (j+primitive_info[j].coordinates);
          break;
        }
        case RoundRectanglePrimitive:
        {
          if (primitive_info[j].coordinates != 3)
            {
              status=MagickFalse;
              break;
            }
          if ((primitive_info[j+2].point.x < 0.0) ||
              (primitive_info[j+2].point.y < 0.0))
            {
              status=MagickFalse;
              break;
            }
          if ((primitive_info[j+1].point.x-primitive_info[j].point.x) < 0.0)
            {
              status=MagickFalse;
              break;
            }
          if ((primitive_info[j+1].point.y-primitive_info[j].point.y) < 0.0)
            {
              status=MagickFalse;
              break;
            }
          status&=TraceRoundRectangle(&mvg_info,primitive_info[j].point,
            primitive_info[j+1].point,primitive_info[j+2].point);
          i=(ssize_t) (j+primitive_info[j].coordinates);
          break;
        }
        case ArcPrimitive:
        {
          if (primitive_info[j].coordinates != 3)
            {
              status=MagickFalse;
              break;
            }
          status&=TraceArc(&mvg_info,primitive_info[j].point,
            primitive_info[j+1].point,primitive_info[j+2].point);
          i=(ssize_t) (j+primitive_info[j].coordinates);
          break;
        }
        case EllipsePrimitive:
        {
          if (primitive_info[j].coordinates != 3)
            {
              status=MagickFalse;
              break;
            }
          if ((primitive_info[j+1].point.x < 0.0) ||
              (primitive_info[j+1].point.y < 0.0))
            {
              status=MagickFalse;
              break;
            }
          status&=TraceEllipse(&mvg_info,primitive_info[j].point,
            primitive_info[j+1].point,primitive_info[j+2].point);
          i=(ssize_t) (j+primitive_info[j].coordinates);
          break;
        }
        case CirclePrimitive:
        {
          if (primitive_info[j].coordinates != 2)
            {
              status=MagickFalse;
              break;
            }
          status&=TraceCircle(&mvg_info,primitive_info[j].point,
            primitive_info[j+1].point);
          i=(ssize_t) (j+primitive_info[j].coordinates);
          break;
        }
        case PolylinePrimitive:
        {
          if (primitive_info[j].coordinates < 1)
            {
              status=MagickFalse;
              break;
            }
          break;
        }
        case PolygonPrimitive:
        {
          if (primitive_info[j].coordinates < 3)
            {
              status=MagickFalse;
              break;
            }
          primitive_info[i]=primitive_info[j];
          primitive_info[i].coordinates=0;
          primitive_info[j].coordinates++;
          primitive_info[j].closed_subpath=MagickTrue;
          i++;
          break;
        }
        case BezierPrimitive:
        {
          if (primitive_info[j].coordinates < 3)
            {
              status=MagickFalse;
              break;
            }
          status&=TraceBezier(&mvg_info,primitive_info[j].coordinates);
          i=(ssize_t) (j+primitive_info[j].coordinates);
          break;
        }
        case PathPrimitive:
        {
          coordinates=(double) TracePath(&mvg_info,token,exception);
          if (coordinates < 0.0)
            {
              status=MagickFalse;
              break;
            }
          i=(ssize_t) (j+coordinates);
          break;
        }
        case AlphaPrimitive:
        case ColorPrimitive:
        {
          ssize_t
            method;
  
          if (primitive_info[j].coordinates != 1)
            {
              status=MagickFalse;
              break;
            }
          (void) GetNextToken(q,&q,extent,token);
          method=ParseCommandOption(MagickMethodOptions,MagickFalse,token);
          if (method == -1)
            {
              status=MagickFalse;
              break;
            }
          primitive_info[j].method=(PaintMethod) method;
          break;
        }
        case TextPrimitive:
        {
          char
            geometry[MagickPathExtent];
  
          if (primitive_info[j].coordinates != 1)
            {
              status=MagickFalse;
              break;
            }
          if (*token != ',')
            (void) GetNextToken(q,&q,extent,token);
          (void) CloneString(&primitive_info[j].text,token);
          /*
            Compute text cursor offset.
          */
          clone_info=CloneDrawInfo((ImageInfo *) NULL,graphic_context[n]);
          if ((fabs(mvg_info.point.x-primitive_info->point.x) < MagickEpsilon) &&
              (fabs(mvg_info.point.y-primitive_info->point.y) < MagickEpsilon))
            {
              mvg_info.point=primitive_info->point;
              primitive_info->point.x+=cursor;
            }
          else
            {
              mvg_info.point=primitive_info->point;
              cursor=0.0;
            }
          (void) FormatLocaleString(geometry,MagickPathExtent,"%+f%+f",
            primitive_info->point.x,primitive_info->point.y);
          clone_info->render=MagickFalse;
          clone_info->text=AcquireString(token);
          status&=GetTypeMetrics(image,clone_info,&metrics,exception);
          clone_info=DestroyDrawInfo(clone_info);
          cursor+=metrics.width;
          if (graphic_context[n]->compliance != SVGCompliance)
            cursor=0.0;
          break;
        }
        case ImagePrimitive:
        {
          if (primitive_info[j].coordinates != 2)
            {
              status=MagickFalse;
              break;
            }
          (void) GetNextToken(q,&q,extent,token);
          (void) CloneString(&primitive_info[j].text,token);
          break;
        }
      }
      mvg_info.offset=i;
      if (status == 0)
        break;
      primitive_info[i].primitive=UndefinedPrimitive;
      if ((image->debug != MagickFalse) && (q > p))
        (void) LogMagickEvent(DrawEvent,GetMagickModule(),"  %.*s",(int) (q-p-1),
          p);
      /*
        Sanity check.
      */
      status&=CheckPrimitiveExtent(&mvg_info,
        ExpandAffine(&graphic_context[n]->affine));
      if (status == 0)
        break;
      status&=CheckPrimitiveExtent(&mvg_info,graphic_context[n]->stroke_width);
      if (status == 0)
        break;
      if (i == 0)
        continue;
      /*
        Transform points.
      */
      for (i=0; primitive_info[i].primitive != UndefinedPrimitive; i++)
      {
        point=primitive_info[i].point;
        primitive_info[i].point.x=graphic_context[n]->affine.sx*point.x+
          graphic_context[n]->affine.ry*point.y+graphic_context[n]->affine.tx;
        primitive_info[i].point.y=graphic_context[n]->affine.rx*point.x+
          graphic_context[n]->affine.sy*point.y+graphic_context[n]->affine.ty;
        point=primitive_info[i].point;
        if (point.x < graphic_context[n]->bounds.x1)
          graphic_context[n]->bounds.x1=point.x;
        if (point.y < graphic_context[n]->bounds.y1)
          graphic_context[n]->bounds.y1=point.y;
        if (point.x > graphic_context[n]->bounds.x2)
          graphic_context[n]->bounds.x2=point.x;
        if (point.y > graphic_context[n]->bounds.y2)
          graphic_context[n]->bounds.y2=point.y;
        if (primitive_info[i].primitive == ImagePrimitive)
          break;
        if (i >= (ssize_t) number_points)
          ThrowFatalException(ResourceLimitFatalError,"MemoryAllocationFailed");
      }
      if (graphic_context[n]->render != MagickFalse)
        {
          if ((n != 0) && (graphic_context[n]->compliance != SVGCompliance) &&
              (graphic_context[n]->clip_mask != (char *) NULL) &&
              (LocaleCompare(graphic_context[n]->clip_mask,
               graphic_context[n-1]->clip_mask) != 0))
            {
              const char
                *clip_path;
  
              clip_path=(const char *) GetValueFromSplayTree(macros,
                graphic_context[n]->clip_mask);
              if (clip_path != (const char *) NULL)
                (void) SetImageArtifact(image,graphic_context[n]->clip_mask,
                  clip_path);
              status&=DrawClipPath(image,graphic_context[n],
                graphic_context[n]->clip_mask,exception);
            }
          status&=DrawPrimitive(image,graphic_context[n],primitive_info,
            exception);
        }
      proceed=SetImageProgress(image,RenderImageTag,q-primitive,(MagickSizeType)
        primitive_extent);
      if (proceed == MagickFalse)
        break;
      if (status == 0)
        break;
    }
    if (image->debug != MagickFalse)
      (void) LogMagickEvent(DrawEvent,GetMagickModule(),"end draw-image");
    /*
      Relinquish resources.
    */
    macros=DestroySplayTree(macros);
    token=DestroyString(token);
    if (primitive_info != (PrimitiveInfo *) NULL)
      {
        for (i=0; primitive_info[i].primitive != UndefinedPrimitive; i++)
          if ((primitive_info[i].primitive == TextPrimitive) ||
              (primitive_info[i].primitive == ImagePrimitive))
            if (primitive_info[i].text != (char *) NULL)
              primitive_info[i].text=DestroyString(primitive_info[i].text);
        primitive_info=(PrimitiveInfo *) RelinquishMagickMemory(primitive_info);
      }
    primitive=DestroyString(primitive);
    if (stops != (StopInfo *) NULL)
      stops=(StopInfo *) RelinquishMagickMemory(stops);
    for ( ; n >= 0; n--)
      graphic_context[n]=DestroyDrawInfo(graphic_context[n]);
    graphic_context=(DrawInfo **) RelinquishMagickMemory(graphic_context);
    if (status == MagickFalse)
      ThrowBinaryException(DrawError,"NonconformingDrawingPrimitiveDefinition",
        keyword);
    return(status != 0 ? MagickTrue : MagickFalse);
  }
  
  MagickExport MagickBooleanType DrawImage(Image *image,const DrawInfo *draw_info,
    ExceptionInfo *exception)
  {
    return(RenderMVGContent(image,draw_info,0,exception));
  }
  
  /*
  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  %                                                                             %
  %                                                                             %
  %                                                                             %
  %   D r a w P a t t e r n P a t h                                             %
  %                                                                             %
  %                                                                             %
  %                                                                             %
  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  %
  %  DrawPatternPath() draws a pattern.
  %
  %  The format of the DrawPatternPath method is:
  %
  %      MagickBooleanType DrawPatternPath(Image *image,const DrawInfo *draw_info,
  %        const char *name,Image **pattern,ExceptionInfo *exception)
  %
  %  A description of each parameter follows:
  %
  %    o image: the image.
  %
  %    o draw_info: the draw info.
  %
  %    o name: the pattern name.
  %
  %    o image: the image.
  %
  %    o exception: return any errors or warnings in this structure.
  %
  */
  MagickExport MagickBooleanType DrawPatternPath(Image *image,
    const DrawInfo *draw_info,const char *name,Image **pattern,
    ExceptionInfo *exception)
  {
    char
      property[MagickPathExtent];
  
    const char
      *geometry,
      *path,
      *type;
  
    DrawInfo
      *clone_info;
  
    ImageInfo
      *image_info;
  
    MagickBooleanType
      status;
  
    assert(image != (Image *) NULL);
    assert(image->signature == MagickCoreSignature);
    if (image->debug != MagickFalse)
      (void) LogMagickEvent(TraceEvent,GetMagickModule(),"%s",image->filename);
    assert(draw_info != (const DrawInfo *) NULL);
    assert(name != (const char *) NULL);
    (void) FormatLocaleString(property,MagickPathExtent,"%s",name);
    path=GetImageArtifact(image,property);
    if (path == (const char *) NULL)
      return(MagickFalse);
    (void) FormatLocaleString(property,MagickPathExtent,"%s-geometry",name);
    geometry=GetImageArtifact(image,property);
    if (geometry == (const char *) NULL)
      return(MagickFalse);
    if ((*pattern) != (Image *) NULL)
      *pattern=DestroyImage(*pattern);
    image_info=AcquireImageInfo();
    image_info->size=AcquireString(geometry);
    *pattern=AcquireImage(image_info,exception);
    image_info=DestroyImageInfo(image_info);
    (void) QueryColorCompliance("#00000000",AllCompliance,
      &(*pattern)->background_color,exception);
    (void) SetImageBackgroundColor(*pattern,exception);
    if (image->debug != MagickFalse)
      (void) LogMagickEvent(DrawEvent,GetMagickModule(),
        "begin pattern-path %s %s",name,geometry);
    clone_info=CloneDrawInfo((ImageInfo *) NULL,draw_info);
    if (clone_info->fill_pattern != (Image *) NULL)
      clone_info->fill_pattern=DestroyImage(clone_info->fill_pattern);
    if (clone_info->stroke_pattern != (Image *) NULL)
      clone_info->stroke_pattern=DestroyImage(clone_info->stroke_pattern);
    (void) FormatLocaleString(property,MagickPathExtent,"%s-type",name);
    type=GetImageArtifact(image,property);
    if (type != (const char *) NULL)
      clone_info->gradient.type=(GradientType) ParseCommandOption(
        MagickGradientOptions,MagickFalse,type);
    (void) CloneString(&clone_info->primitive,path);
    status=RenderMVGContent(*pattern,clone_info,0,exception);
    clone_info=DestroyDrawInfo(clone_info);
    if (image->debug != MagickFalse)
      (void) LogMagickEvent(DrawEvent,GetMagickModule(),"end pattern-path");
    return(status);
  }
  
  /*
  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  %                                                                             %
  %                                                                             %
  %                                                                             %
  +   D r a w P o l y g o n P r i m i t i v e                                   %
  %                                                                             %
  %                                                                             %
  %                                                                             %
  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  %
  %  DrawPolygonPrimitive() draws a polygon on the image.
  %
  %  The format of the DrawPolygonPrimitive method is:
  %
  %      MagickBooleanType DrawPolygonPrimitive(Image *image,
  %        const DrawInfo *draw_info,const PrimitiveInfo *primitive_info,
  %        ExceptionInfo *exception)
  %
  %  A description of each parameter follows:
  %
  %    o image: the image.
  %
  %    o draw_info: the draw info.
  %
  %    o primitive_info: Specifies a pointer to a PrimitiveInfo structure.
  %
  %    o exception: return any errors or warnings in this structure.
  %
  */
  
  static PolygonInfo **DestroyPolygonThreadSet(PolygonInfo **polygon_info)
  {
    register ssize_t
      i;
  
    assert(polygon_info != (PolygonInfo **) NULL);
    for (i=0; i < (ssize_t) GetMagickResourceLimit(ThreadResource); i++)
      if (polygon_info[i] != (PolygonInfo *) NULL)
        polygon_info[i]=DestroyPolygonInfo(polygon_info[i]);
    polygon_info=(PolygonInfo **) RelinquishMagickMemory(polygon_info);
    return(polygon_info);
  }
  
  static PolygonInfo **AcquirePolygonThreadSet(
    const PrimitiveInfo *primitive_info)
  {
    PathInfo
      *magick_restrict path_info;
  
    PolygonInfo
      **polygon_info;
  
    register ssize_t
      i;
  
    size_t
      number_threads;
  
    number_threads=(size_t) GetMagickResourceLimit(ThreadResource);
    polygon_info=(PolygonInfo **) AcquireQuantumMemory(number_threads,
      sizeof(*polygon_info));
    if (polygon_info == (PolygonInfo **) NULL)
      return((PolygonInfo **) NULL);
    (void) memset(polygon_info,0,number_threads*sizeof(*polygon_info));
    path_info=ConvertPrimitiveToPath(primitive_info);
    if (path_info == (PathInfo *) NULL)
      return(DestroyPolygonThreadSet(polygon_info));
    for (i=0; i < (ssize_t) number_threads; i++)
    {
      polygon_info[i]=ConvertPathToPolygon(path_info);
      if (polygon_info[i] == (PolygonInfo *) NULL)
        return(DestroyPolygonThreadSet(polygon_info));
    }
    path_info=(PathInfo *) RelinquishMagickMemory(path_info);
    return(polygon_info);
  }
  
  static double GetFillAlpha(PolygonInfo *polygon_info,const double mid,
    const MagickBooleanType fill,const FillRule fill_rule,const ssize_t x,
    const ssize_t y,double *stroke_alpha)
  {
    double
      alpha,
      beta,
      distance,
      subpath_alpha;
  
    PointInfo
      delta;
  
    register const PointInfo
      *q;
  
    register EdgeInfo
      *p;
  
    register ssize_t
      i;
  
    ssize_t
      j,
      winding_number;
  
    /*
      Compute fill & stroke opacity for this (x,y) point.
    */
    *stroke_alpha=0.0;
    subpath_alpha=0.0;
    p=polygon_info->edges;
    for (j=0; j < (ssize_t) polygon_info->number_edges; j++, p++)
    {
      if ((double) y <= (p->bounds.y1-mid-0.5))
        break;
      if ((double) y > (p->bounds.y2+mid+0.5))
        {
          (void) DestroyEdge(polygon_info,(size_t) j);
          continue;
        }
      if (((double) x <= (p->bounds.x1-mid-0.5)) ||
          ((double) x > (p->bounds.x2+mid+0.5)))
        continue;
      i=(ssize_t) MagickMax((double) p->highwater,1.0);
      for ( ; i < (ssize_t) p->number_points; i++)
      {
        if ((double) y <= (p->points[i-1].y-mid-0.5))
          break;
        if ((double) y > (p->points[i].y+mid+0.5))
          continue;
        if (p->scanline != (double) y)
          {
            p->scanline=(double) y;
            p->highwater=(size_t) i;
          }
        /*
          Compute distance between a point and an edge.
        */
        q=p->points+i-1;
        delta.x=(q+1)->x-q->x;
        delta.y=(q+1)->y-q->y;
        beta=delta.x*(x-q->x)+delta.y*(y-q->y);
        if (beta <= 0.0)
          {
            delta.x=(double) x-q->x;
            delta.y=(double) y-q->y;
            distance=delta.x*delta.x+delta.y*delta.y;
          }
        else
          {
            alpha=delta.x*delta.x+delta.y*delta.y;
            if (beta >= alpha)
              {
                delta.x=(double) x-(q+1)->x;
                delta.y=(double) y-(q+1)->y;
                distance=delta.x*delta.x+delta.y*delta.y;
              }
            else
              {
                alpha=PerceptibleReciprocal(alpha);
                beta=delta.x*(y-q->y)-delta.y*(x-q->x)+MagickEpsilon;
                distance=alpha*beta*beta;
              }
          }
        /*
          Compute stroke & subpath opacity.
        */
        beta=0.0;
        if (p->ghostline == MagickFalse)
          {
            alpha=mid+0.5;
            if ((*stroke_alpha < 1.0) &&
                (distance <= ((alpha+0.25)*(alpha+0.25))))
              {
                alpha=mid-0.5;
                if (distance <= ((alpha+0.25)*(alpha+0.25)))
                  *stroke_alpha=1.0;
                else
                  {
                    beta=1.0;
                    if (fabs(distance-1.0) >= MagickEpsilon)
                      beta=sqrt((double) distance);
                    alpha=beta-mid-0.5;
                    if (*stroke_alpha < ((alpha-0.25)*(alpha-0.25)))
                      *stroke_alpha=(alpha-0.25)*(alpha-0.25);
                  }
              }
          }
        if ((fill == MagickFalse) || (distance > 1.0) || (subpath_alpha >= 1.0))
          continue;
        if (distance <= 0.0)
          {
            subpath_alpha=1.0;
            continue;
          }
        if (distance > 1.0)
          continue;
        if (fabs(beta) < MagickEpsilon)
          {
            beta=1.0;
            if (fabs(distance-1.0) >= MagickEpsilon)
              beta=sqrt(distance);
          }
        alpha=beta-1.0;
        if (subpath_alpha < (alpha*alpha))
          subpath_alpha=alpha*alpha;
      }
    }
    /*
      Compute fill opacity.
    */
    if (fill == MagickFalse)
      return(0.0);
    if (subpath_alpha >= 1.0)
      return(1.0);
    /*
      Determine winding number.
    */
    winding_number=0;
    p=polygon_info->edges;
    for (j=0; j < (ssize_t) polygon_info->number_edges; j++, p++)
    {
      if ((double) y <= p->bounds.y1)
        break;
      if (((double) y > p->bounds.y2) || ((double) x <= p->bounds.x1))
        continue;
      if ((double) x > p->bounds.x2)
        {
          winding_number+=p->direction ? 1 : -1;
          continue;
        }
      i=(ssize_t) MagickMax((double) p->highwater,1.0);
      for ( ; i < (ssize_t) (p->number_points-1); i++)
        if ((double) y <= p->points[i].y)
          break;
      q=p->points+i-1;
      if ((((q+1)->x-q->x)*(y-q->y)) <= (((q+1)->y-q->y)*(x-q->x)))
        winding_number+=p->direction ? 1 : -1;
    }
    if (fill_rule != NonZeroRule)
      {
        if ((MagickAbsoluteValue(winding_number) & 0x01) != 0)
          return(1.0);
      }
    else
      if (MagickAbsoluteValue(winding_number) != 0)
        return(1.0);
    return(subpath_alpha);
  }
  
  static MagickBooleanType DrawPolygonPrimitive(Image *image,
    const DrawInfo *draw_info,const PrimitiveInfo *primitive_info,
    ExceptionInfo *exception)
  {
    CacheView
      *image_view;
  
    MagickBooleanType
      fill,
      status;
  
    double
      mid;
  
    PolygonInfo
      **magick_restrict polygon_info;
  
    register EdgeInfo
      *p;
  
    register ssize_t
      i;
  
    SegmentInfo
      bounds;
  
    ssize_t
      start_y,
      stop_y,
      y;
  
    assert(image != (Image *) NULL);
    assert(image->signature == MagickCoreSignature);
    if (image->debug != MagickFalse)
      (void) LogMagickEvent(TraceEvent,GetMagickModule(),"%s",image->filename);
    assert(draw_info != (DrawInfo *) NULL);
    assert(draw_info->signature == MagickCoreSignature);
    assert(primitive_info != (PrimitiveInfo *) NULL);
    if (primitive_info->coordinates <= 1)
      return(MagickTrue);
    /*
      Compute bounding box.
    */
    polygon_info=AcquirePolygonThreadSet(primitive_info);
    if (polygon_info == (PolygonInfo **) NULL)
      return(MagickFalse);
    if (image->debug != MagickFalse)
      (void) LogMagickEvent(DrawEvent,GetMagickModule(),"    begin draw-polygon");
    fill=(primitive_info->method == FillToBorderMethod) ||
      (primitive_info->method == FloodfillMethod) ? MagickTrue : MagickFalse;
    mid=ExpandAffine(&draw_info->affine)*draw_info->stroke_width/2.0;
    bounds=polygon_info[0]->edges[0].bounds;
    for (i=1; i < (ssize_t) polygon_info[0]->number_edges; i++)
    {
      p=polygon_info[0]->edges+i;
      if (p->bounds.x1 < bounds.x1)
        bounds.x1=p->bounds.x1;
      if (p->bounds.y1 < bounds.y1)
        bounds.y1=p->bounds.y1;
      if (p->bounds.x2 > bounds.x2)
        bounds.x2=p->bounds.x2;
      if (p->bounds.y2 > bounds.y2)
        bounds.y2=p->bounds.y2;
    }
    bounds.x1-=(mid+1.0);
    bounds.y1-=(mid+1.0);
    bounds.x2+=(mid+1.0);
    bounds.y2+=(mid+1.0);
    if ((bounds.x1 >= (double) image->columns) ||
        (bounds.y1 >= (double) image->rows) ||
        (bounds.x2 <= 0.0) || (bounds.y2 <= 0.0))
      {
        polygon_info=DestroyPolygonThreadSet(polygon_info);
        return(MagickTrue);  /* virtual polygon */
      }
    bounds.x1=bounds.x1 < 0.0 ? 0.0 : bounds.x1 >= (double) image->columns-1.0 ?
      (double) image->columns-1.0 : bounds.x1;
    bounds.y1=bounds.y1 < 0.0 ? 0.0 : bounds.y1 >= (double) image->rows-1.0 ?
      (double) image->rows-1.0 : bounds.y1;
    bounds.x2=bounds.x2 < 0.0 ? 0.0 : bounds.x2 >= (double) image->columns-1.0 ?
      (double) image->columns-1.0 : bounds.x2;
    bounds.y2=bounds.y2 < 0.0 ? 0.0 : bounds.y2 >= (double) image->rows-1.0 ?
      (double) image->rows-1.0 : bounds.y2;
    status=MagickTrue;
    image_view=AcquireAuthenticCacheView(image,exception);
    if ((primitive_info->coordinates == 1) ||
        (polygon_info[0]->number_edges == 0))
      {
        /*
          Draw point.
        */
        start_y=(ssize_t) ceil(bounds.y1-0.5);
        stop_y=(ssize_t) floor(bounds.y2+0.5);
  #if defined(MAGICKCORE_OPENMP_SUPPORT)
        #pragma omp parallel for schedule(static) shared(status) \
          magick_number_threads(image,image,stop_y-start_y+1,1)
  #endif
        for (y=start_y; y <= stop_y; y++)
        {
          MagickBooleanType
            sync;
  
          PixelInfo
            pixel;
  
          register ssize_t
            x;
  
          register Quantum
            *magick_restrict q;
  
          ssize_t
            start_x,
            stop_x;
  
          if (status == MagickFalse)
            continue;
          start_x=(ssize_t) ceil(bounds.x1-0.5);
          stop_x=(ssize_t) floor(bounds.x2+0.5);
          x=start_x;
          q=GetCacheViewAuthenticPixels(image_view,x,y,(size_t) (stop_x-x+1),1,
            exception);
          if (q == (Quantum *) NULL)
            {
              status=MagickFalse;
              continue;
            }
          GetPixelInfo(image,&pixel);
          for ( ; x <= stop_x; x++)
          {
            if ((x == (ssize_t) ceil(primitive_info->point.x-0.5)) &&
                (y == (ssize_t) ceil(primitive_info->point.y-0.5)))
              {
                GetFillColor(draw_info,x-start_x,y-start_y,&pixel,exception);
                SetPixelViaPixelInfo(image,&pixel,q);
              }
            q+=GetPixelChannels(image);
          }
          sync=SyncCacheViewAuthenticPixels(image_view,exception);
          if (sync == MagickFalse)
            status=MagickFalse;
        }
        image_view=DestroyCacheView(image_view);
        polygon_info=DestroyPolygonThreadSet(polygon_info);
        if (image->debug != MagickFalse)
          (void) LogMagickEvent(DrawEvent,GetMagickModule(),
            "    end draw-polygon");
        return(status);
      }
    /*
      Draw polygon or line.
    */
    start_y=(ssize_t) ceil(bounds.y1-0.5);
    stop_y=(ssize_t) floor(bounds.y2+0.5);
  #if defined(MAGICKCORE_OPENMP_SUPPORT)
    #pragma omp parallel for schedule(static) shared(status) \
      magick_number_threads(image,image,stop_y-start_y+1,1)
  #endif
    for (y=start_y; y <= stop_y; y++)
    {
      const int
        id = GetOpenMPThreadId();
  
      register Quantum
        *magick_restrict q;
  
      register ssize_t
        x;
  
      ssize_t
        start_x,
        stop_x;
  
      if (status == MagickFalse)
        continue;
      start_x=(ssize_t) ceil(bounds.x1-0.5);
      stop_x=(ssize_t) floor(bounds.x2+0.5);
      q=GetCacheViewAuthenticPixels(image_view,start_x,y,(size_t) (stop_x-start_x+
        1),1,exception);
      if (q == (Quantum *) NULL)
        {
          status=MagickFalse;
          continue;
        }
      for (x=start_x; x <= stop_x; x++)
      {
        double
          fill_alpha,
          stroke_alpha;
  
        PixelInfo
          fill_color,
          stroke_color;
  
        /*
          Fill and/or stroke.
        */
        fill_alpha=GetFillAlpha(polygon_info[id],mid,fill,draw_info->fill_rule,
          x,y,&stroke_alpha);
        if (draw_info->stroke_antialias == MagickFalse)
          {
            fill_alpha=fill_alpha > 0.5 ? 1.0 : 0.0;
            stroke_alpha=stroke_alpha > 0.5 ? 1.0 : 0.0;
          }
        GetFillColor(draw_info,x-start_x,y-start_y,&fill_color,exception);
        CompositePixelOver(image,&fill_color,fill_alpha*fill_color.alpha,q,
          (double) GetPixelAlpha(image,q),q);
        GetStrokeColor(draw_info,x-start_x,y-start_y,&stroke_color,exception);
        CompositePixelOver(image,&stroke_color,stroke_alpha*stroke_color.alpha,q,
          (double) GetPixelAlpha(image,q),q);
        q+=GetPixelChannels(image);
      }
      if (SyncCacheViewAuthenticPixels(image_view,exception) == MagickFalse)
        status=MagickFalse;
    }
    image_view=DestroyCacheView(image_view);
    polygon_info=DestroyPolygonThreadSet(polygon_info);
    if (image->debug != MagickFalse)
      (void) LogMagickEvent(DrawEvent,GetMagickModule(),"    end draw-polygon");
    return(status);
  }
  
  /*
  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  %                                                                             %
  %                                                                             %
  %                                                                             %
  %   D r a w P r i m i t i v e                                                 %
  %                                                                             %
  %                                                                             %
  %                                                                             %
  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  %
  %  DrawPrimitive() draws a primitive (line, rectangle, ellipse) on the image.
  %
  %  The format of the DrawPrimitive method is:
  %
  %      MagickBooleanType DrawPrimitive(Image *image,const DrawInfo *draw_info,
  %        PrimitiveInfo *primitive_info,ExceptionInfo *exception)
  %
  %  A description of each parameter follows:
  %
  %    o image: the image.
  %
  %    o draw_info: the draw info.
  %
  %    o primitive_info: Specifies a pointer to a PrimitiveInfo structure.
  %
  %    o exception: return any errors or warnings in this structure.
  %
  */
  static void LogPrimitiveInfo(const PrimitiveInfo *primitive_info)
  {
    const char
      *methods[] =
      {
        "point",
        "replace",
        "floodfill",
        "filltoborder",
        "reset",
        "?"
      };
  
    PointInfo
      p,
      point,
      q;
  
    register ssize_t
      i,
      x;
  
    ssize_t
      coordinates,
      y;
  
    x=(ssize_t) ceil(primitive_info->point.x-0.5);
    y=(ssize_t) ceil(primitive_info->point.y-0.5);
    switch (primitive_info->primitive)
    {
      case AlphaPrimitive:
      {
        (void) LogMagickEvent(DrawEvent,GetMagickModule(),
          "AlphaPrimitive %.20g,%.20g %s",(double) x,(double) y,
          methods[primitive_info->method]);
        return;
      }
      case ColorPrimitive:
      {
        (void) LogMagickEvent(DrawEvent,GetMagickModule(),
          "ColorPrimitive %.20g,%.20g %s",(double) x,(double) y,
          methods[primitive_info->method]);
        return;
      }
      case ImagePrimitive:
      {
        (void) LogMagickEvent(DrawEvent,GetMagickModule(),
          "ImagePrimitive %.20g,%.20g",(double) x,(double) y);
        return;
      }
      case PointPrimitive:
      {
        (void) LogMagickEvent(DrawEvent,GetMagickModule(),
          "PointPrimitive %.20g,%.20g %s",(double) x,(double) y,
          methods[primitive_info->method]);
        return;
      }
      case TextPrimitive:
      {
        (void) LogMagickEvent(DrawEvent,GetMagickModule(),
          "TextPrimitive %.20g,%.20g",(double) x,(double) y);
        return;
      }
      default:
        break;
    }
    coordinates=0;
    p=primitive_info[0].point;
    q.x=(-1.0);
    q.y=(-1.0);
    for (i=0; primitive_info[i].primitive != UndefinedPrimitive; i++)
    {
      point=primitive_info[i].point;
      if (coordinates <= 0)
        {
          coordinates=(ssize_t) primitive_info[i].coordinates;
          (void) LogMagickEvent(DrawEvent,GetMagickModule(),
            "    begin open (%.20g)",(double) coordinates);
          p=point;
        }
      point=primitive_info[i].point;
      if ((fabs(q.x-point.x) >= MagickEpsilon) ||
          (fabs(q.y-point.y) >= MagickEpsilon))
        (void) LogMagickEvent(DrawEvent,GetMagickModule(),
          "      %.20g: %.18g,%.18g",(double) coordinates,point.x,point.y);
      else
        (void) LogMagickEvent(DrawEvent,GetMagickModule(),
          "      %.20g: %g %g (duplicate)",(double) coordinates,point.x,point.y);
      q=point;
      coordinates--;
      if (coordinates > 0)
        continue;
      if ((fabs(p.x-point.x) >= MagickEpsilon) ||
          (fabs(p.y-point.y) >= MagickEpsilon))
        (void) LogMagickEvent(DrawEvent,GetMagickModule(),"    end last (%.20g)",
          (double) coordinates);
      else
        (void) LogMagickEvent(DrawEvent,GetMagickModule(),"    end open (%.20g)",
          (double) coordinates);
    }
  }
  
  MagickExport MagickBooleanType DrawPrimitive(Image *image,
    const DrawInfo *draw_info,const PrimitiveInfo *primitive_info,
    ExceptionInfo *exception)
  {
    CacheView
      *image_view;
  
    MagickStatusType
      status;
  
    register ssize_t
      i,
      x;
  
    ssize_t
      y;
  
    if (image->debug != MagickFalse)
      {
        (void) LogMagickEvent(DrawEvent,GetMagickModule(),
          "  begin draw-primitive");
        (void) LogMagickEvent(DrawEvent,GetMagickModule(),
          "    affine: %g,%g,%g,%g,%g,%g",draw_info->affine.sx,
          draw_info->affine.rx,draw_info->affine.ry,draw_info->affine.sy,
          draw_info->affine.tx,draw_info->affine.ty);
      }
    status=MagickTrue;
    if ((IsGrayColorspace(image->colorspace) != MagickFalse) &&
        ((IsPixelInfoGray(&draw_info->fill) == MagickFalse) ||
         (IsPixelInfoGray(&draw_info->stroke) == MagickFalse)))
      status&=SetImageColorspace(image,sRGBColorspace,exception);
    if (draw_info->compliance == SVGCompliance)
      {
        status&=SetImageMask(image,WritePixelMask,draw_info->clipping_mask,
          exception);
        status&=SetImageMask(image,CompositePixelMask,draw_info->composite_mask,
          exception);
      }
    x=(ssize_t) ceil(primitive_info->point.x-0.5);
    y=(ssize_t) ceil(primitive_info->point.y-0.5);
    image_view=AcquireAuthenticCacheView(image,exception);
    switch (primitive_info->primitive)
    {
      case AlphaPrimitive:
      {
        if (image->alpha_trait == UndefinedPixelTrait)
          status&=SetImageAlphaChannel(image,OpaqueAlphaChannel,exception);
        switch (primitive_info->method)
        {
          case PointMethod:
          default:
          {
            PixelInfo
              pixel;
  
            register Quantum
              *q;
  
            q=GetCacheViewAuthenticPixels(image_view,x,y,1,1,exception);
            if (q == (Quantum *) NULL)
              break;
            GetFillColor(draw_info,x,y,&pixel,exception);
            SetPixelAlpha(image,ClampToQuantum(pixel.alpha),q);
            status&=SyncCacheViewAuthenticPixels(image_view,exception);
            break;
          }
          case ReplaceMethod:
          {
            PixelInfo
              pixel,
              target;
  
            status&=GetOneCacheViewVirtualPixelInfo(image_view,x,y,&target,
              exception);
            GetPixelInfo(image,&pixel);
            for (y=0; y < (ssize_t) image->rows; y++)
            {
              register Quantum
                *magick_restrict q;
  
              q=GetCacheViewAuthenticPixels(image_view,0,y,image->columns,1,
                exception);
              if (q == (Quantum *) NULL)
                break;
              for (x=0; x < (ssize_t) image->columns; x++)
              {
                GetPixelInfoPixel(image,q,&pixel);
                if (IsFuzzyEquivalencePixelInfo(&pixel,&target) == MagickFalse)
                  {
                    q+=GetPixelChannels(image);
                    continue;
                  }
                GetFillColor(draw_info,x,y,&pixel,exception);
                SetPixelAlpha(image,ClampToQuantum(pixel.alpha),q);
                q+=GetPixelChannels(image);
              }
              status&=SyncCacheViewAuthenticPixels(image_view,exception);
              if (status == MagickFalse)
                break;
            }
            break;
          }
          case FloodfillMethod:
          case FillToBorderMethod:
          {
            ChannelType
              channel_mask;
  
            PixelInfo
              target;
  
            status&=GetOneVirtualPixelInfo(image,TileVirtualPixelMethod,x,y,
              &target,exception);
            if (primitive_info->method == FillToBorderMethod)
              {
                target.red=(double) draw_info->border_color.red;
                target.green=(double) draw_info->border_color.green;
                target.blue=(double) draw_info->border_color.blue;
              }
            channel_mask=SetImageChannelMask(image,AlphaChannel);
            status&=FloodfillPaintImage(image,draw_info,&target,x,y,
              primitive_info->method == FloodfillMethod ? MagickFalse :
              MagickTrue,exception);
            (void) SetImageChannelMask(image,channel_mask);
            break;
          }
          case ResetMethod:
          {
            PixelInfo
              pixel;
  
            for (y=0; y < (ssize_t) image->rows; y++)
            {
              register Quantum
                *magick_restrict q;
  
              q=GetCacheViewAuthenticPixels(image_view,0,y,image->columns,1,
                exception);
              if (q == (Quantum *) NULL)
                break;
              for (x=0; x < (ssize_t) image->columns; x++)
              {
                GetFillColor(draw_info,x,y,&pixel,exception);
                SetPixelAlpha(image,ClampToQuantum(pixel.alpha),q);
                q+=GetPixelChannels(image);
              }
              status&=SyncCacheViewAuthenticPixels(image_view,exception);
              if (status == MagickFalse)
                break;
            }
            break;
          }
        }
        break;
      }
      case ColorPrimitive:
      {
        switch (primitive_info->method)
        {
          case PointMethod:
          default:
          {
            PixelInfo
              pixel;
  
            register Quantum
              *q;
  
            q=GetCacheViewAuthenticPixels(image_view,x,y,1,1,exception);
            if (q == (Quantum *) NULL)
              break;
            GetPixelInfo(image,&pixel);
            GetFillColor(draw_info,x,y,&pixel,exception);
            SetPixelViaPixelInfo(image,&pixel,q);
            status&=SyncCacheViewAuthenticPixels(image_view,exception);
            break;
          }
          case ReplaceMethod:
          {
            PixelInfo
              pixel,
              target;
  
            status&=GetOneCacheViewVirtualPixelInfo(image_view,x,y,&target,
              exception);
            for (y=0; y < (ssize_t) image->rows; y++)
            {
              register Quantum
                *magick_restrict q;
  
              q=GetCacheViewAuthenticPixels(image_view,0,y,image->columns,1,
                exception);
              if (q == (Quantum *) NULL)
                break;
              for (x=0; x < (ssize_t) image->columns; x++)
              {
                GetPixelInfoPixel(image,q,&pixel);
                if (IsFuzzyEquivalencePixelInfo(&pixel,&target) == MagickFalse)
                  {
                    q+=GetPixelChannels(image);
                    continue;
                  }
                GetFillColor(draw_info,x,y,&pixel,exception);
                SetPixelViaPixelInfo(image,&pixel,q);
                q+=GetPixelChannels(image);
              }
              status&=SyncCacheViewAuthenticPixels(image_view,exception);
              if (status == MagickFalse)
                break;
            }
            break;
          }
          case FloodfillMethod:
          case FillToBorderMethod:
          {
            PixelInfo
              target;
  
            status&=GetOneVirtualPixelInfo(image,TileVirtualPixelMethod,x,y,
              &target,exception);
            if (primitive_info->method == FillToBorderMethod)
              {
                target.red=(double) draw_info->border_color.red;
                target.green=(double) draw_info->border_color.green;
                target.blue=(double) draw_info->border_color.blue;
              }
            status&=FloodfillPaintImage(image,draw_info,&target,x,y,
              primitive_info->method == FloodfillMethod ? MagickFalse :
              MagickTrue,exception);
            break;
          }
          case ResetMethod:
          {
            PixelInfo
              pixel;
  
            GetPixelInfo(image,&pixel);
            for (y=0; y < (ssize_t) image->rows; y++)
            {
              register Quantum
                *magick_restrict q;
  
              q=GetCacheViewAuthenticPixels(image_view,0,y,image->columns,1,
                exception);
              if (q == (Quantum *) NULL)
                break;
              for (x=0; x < (ssize_t) image->columns; x++)
              {
                GetFillColor(draw_info,x,y,&pixel,exception);
                SetPixelViaPixelInfo(image,&pixel,q);
                q+=GetPixelChannels(image);
              }
              status&=SyncCacheViewAuthenticPixels(image_view,exception);
              if (status == MagickFalse)
                break;
            }
            break;
          }
        }
        break;
      }
      case ImagePrimitive:
      {
        AffineMatrix
          affine;
  
        char
          composite_geometry[MagickPathExtent];
  
        Image
          *composite_image,
          *composite_images;
  
        ImageInfo
          *clone_info;
  
        RectangleInfo
          geometry;
  
        ssize_t
          x1,
          y1;
  
        if (primitive_info->text == (char *) NULL)
          break;
        clone_info=AcquireImageInfo();
        composite_images=(Image *) NULL;
        if (LocaleNCompare(primitive_info->text,"data:",5) == 0)
          composite_images=ReadInlineImage(clone_info,primitive_info->text,
            exception);
        else
          if (*primitive_info->text != '\0')
            {
              (void) CopyMagickString(clone_info->filename,primitive_info->text,
                MagickPathExtent);
              status&=SetImageInfo(clone_info,0,exception);
              if ((LocaleNCompare(clone_info->magick,"http",4) == 0) ||
                  (LocaleCompare(clone_info->magick,"mpri") == 0))
                (void) CopyMagickString(clone_info->filename,primitive_info->text,
                  MagickPathExtent);
              composite_images=ReadImage(clone_info,exception);
            }
        clone_info=DestroyImageInfo(clone_info);
        if (composite_images == (Image *) NULL)
          {
            status=MagickFalse;
            break;
          }
        composite_image=RemoveFirstImageFromList(&composite_images);
        composite_images=DestroyImageList(composite_images);
        (void) SetImageProgressMonitor(composite_image,(MagickProgressMonitor)
          NULL,(void *) NULL);
        x1=(ssize_t) ceil(primitive_info[1].point.x-0.5);
        y1=(ssize_t) ceil(primitive_info[1].point.y-0.5);
        if (((x1 != 0L) && (x1 != (ssize_t) composite_image->columns)) ||
            ((y1 != 0L) && (y1 != (ssize_t) composite_image->rows)))
          {
            /*
              Resize image.
            */
            (void) FormatLocaleString(composite_geometry,MagickPathExtent,
              "%gx%g!",primitive_info[1].point.x,primitive_info[1].point.y);
            composite_image->filter=image->filter;
            status&=TransformImage(&composite_image,(char *) NULL,
              composite_geometry,exception);
          }
        if (composite_image->alpha_trait == UndefinedPixelTrait)
          status&=SetImageAlphaChannel(composite_image,OpaqueAlphaChannel,
            exception);
        if (draw_info->alpha != OpaqueAlpha)
          status&=SetImageAlpha(composite_image,draw_info->alpha,exception);
        SetGeometry(image,&geometry);
        image->gravity=draw_info->gravity;
        geometry.x=x;
        geometry.y=y;
        (void) FormatLocaleString(composite_geometry,MagickPathExtent,
          "%.20gx%.20g%+.20g%+.20g",(double) composite_image->columns,(double)
          composite_image->rows,(double) geometry.x,(double) geometry.y);
        (void) ParseGravityGeometry(image,composite_geometry,&geometry,exception);
        affine=draw_info->affine;
        affine.tx=(double) geometry.x;
        affine.ty=(double) geometry.y;
        composite_image->interpolate=image->interpolate;
        if ((draw_info->compose == OverCompositeOp) ||
            (draw_info->compose == SrcOverCompositeOp))
          status&=DrawAffineImage(image,composite_image,&affine,exception);
        else
          status&=CompositeImage(image,composite_image,draw_info->compose,
            MagickTrue,geometry.x,geometry.y,exception);
        composite_image=DestroyImage(composite_image);
        break;
      }
      case PointPrimitive:
      {
        PixelInfo
          fill_color;
  
        register Quantum
          *q;
  
        if ((y < 0) || (y >= (ssize_t) image->rows))
          break;
        if ((x < 0) || (x >= (ssize_t) image->columns))
          break;
        q=GetCacheViewAuthenticPixels(image_view,x,y,1,1,exception);
        if (q == (Quantum *) NULL)
          break;
        GetFillColor(draw_info,x,y,&fill_color,exception);
        CompositePixelOver(image,&fill_color,(double) fill_color.alpha,q,(double)
          GetPixelAlpha(image,q),q);
        status&=SyncCacheViewAuthenticPixels(image_view,exception);
        break;
      }
      case TextPrimitive:
      {
        char
          geometry[MagickPathExtent];
  
        DrawInfo
          *clone_info;
  
        if (primitive_info->text == (char *) NULL)
          break;
        clone_info=CloneDrawInfo((ImageInfo *) NULL,draw_info);
        (void) CloneString(&clone_info->text,primitive_info->text);
        (void) FormatLocaleString(geometry,MagickPathExtent,"%+f%+f",
          primitive_info->point.x,primitive_info->point.y);
        (void) CloneString(&clone_info->geometry,geometry);
        status&=AnnotateImage(image,clone_info,exception);
        clone_info=DestroyDrawInfo(clone_info);
        break;
      }
      default:
      {
        double
          mid,
          scale;
  
        DrawInfo
          *clone_info;
  
        if (IsEventLogging() != MagickFalse)
          LogPrimitiveInfo(primitive_info);
        scale=ExpandAffine(&draw_info->affine);
        if ((draw_info->dash_pattern != (double *) NULL) &&
            (fabs(draw_info->dash_pattern[0]) >= MagickEpsilon) &&
            (fabs(scale*draw_info->stroke_width) >= MagickEpsilon) &&
            (draw_info->stroke.alpha != (Quantum) TransparentAlpha))
          {
            /*
              Draw dash polygon.
            */
            clone_info=CloneDrawInfo((ImageInfo *) NULL,draw_info);
            clone_info->stroke_width=0.0;
            clone_info->stroke.alpha=(MagickRealType) TransparentAlpha;
            status&=DrawPolygonPrimitive(image,clone_info,primitive_info,
              exception);
            clone_info=DestroyDrawInfo(clone_info);
            if (status != MagickFalse)
              status&=DrawDashPolygon(draw_info,primitive_info,image,exception);
            break;
          }
        mid=ExpandAffine(&draw_info->affine)*draw_info->stroke_width/2.0;
        if ((mid > 1.0) &&
            ((draw_info->stroke.alpha != (Quantum) TransparentAlpha) ||
             (draw_info->stroke_pattern != (Image *) NULL)))
          {
            double
              x,
              y;
  
            MagickBooleanType
              closed_path;
  
            /*
              Draw strokes while respecting line cap/join attributes.
            */
            closed_path=primitive_info[0].closed_subpath;
            i=(ssize_t) primitive_info[0].coordinates;
            x=fabs(primitive_info[i-1].point.x-primitive_info[0].point.x);
            y=fabs(primitive_info[i-1].point.y-primitive_info[0].point.y);
            if ((x < MagickEpsilon) && (y < MagickEpsilon))
              closed_path=MagickTrue;
            if ((((draw_info->linecap == RoundCap) ||
                  (closed_path != MagickFalse)) &&
                 (draw_info->linejoin == RoundJoin)) ||
                 (primitive_info[i].primitive != UndefinedPrimitive))
              {
                status&=DrawPolygonPrimitive(image,draw_info,primitive_info,
                  exception);
                break;
              }
            clone_info=CloneDrawInfo((ImageInfo *) NULL,draw_info);
            clone_info->stroke_width=0.0;
            clone_info->stroke.alpha=(MagickRealType) TransparentAlpha;
            status&=DrawPolygonPrimitive(image,clone_info,primitive_info,
              exception);
            clone_info=DestroyDrawInfo(clone_info);
            if (status != MagickFalse)
              status&=DrawStrokePolygon(image,draw_info,primitive_info,exception);
            break;
          }
        status&=DrawPolygonPrimitive(image,draw_info,primitive_info,exception);
        break;
      }
    }
    image_view=DestroyCacheView(image_view);
    if (draw_info->compliance == SVGCompliance)
      {
        status&=SetImageMask(image,WritePixelMask,(Image *) NULL,exception);
        status&=SetImageMask(image,CompositePixelMask,(Image *) NULL,exception);
      }
    if (image->debug != MagickFalse)
      (void) LogMagickEvent(DrawEvent,GetMagickModule(),"  end draw-primitive");
    return(status != 0 ? MagickTrue : MagickFalse);
  }
  
  /*
  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  %                                                                             %
  %                                                                             %
  %                                                                             %
  +   D r a w S t r o k e P o l y g o n                                         %
  %                                                                             %
  %                                                                             %
  %                                                                             %
  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  %
  %  DrawStrokePolygon() draws a stroked polygon (line, rectangle, ellipse) on
  %  the image while respecting the line cap and join attributes.
  %
  %  The format of the DrawStrokePolygon method is:
  %
  %      MagickBooleanType DrawStrokePolygon(Image *image,
  %        const DrawInfo *draw_info,const PrimitiveInfo *primitive_info)
  %
  %  A description of each parameter follows:
  %
  %    o image: the image.
  %
  %    o draw_info: the draw info.
  %
  %    o primitive_info: Specifies a pointer to a PrimitiveInfo structure.
  %
  %
  */
  
  static MagickBooleanType DrawRoundLinecap(Image *image,
    const DrawInfo *draw_info,const PrimitiveInfo *primitive_info,
    ExceptionInfo *exception)
  {
    PrimitiveInfo
      linecap[5];
  
    register ssize_t
      i;
  
    for (i=0; i < 4; i++)
      linecap[i]=(*primitive_info);
    linecap[0].coordinates=4;
    linecap[1].point.x+=2.0*MagickEpsilon;
    linecap[2].point.x+=2.0*MagickEpsilon;
    linecap[2].point.y+=2.0*MagickEpsilon;
    linecap[3].point.y+=2.0*MagickEpsilon;
    linecap[4].primitive=UndefinedPrimitive;
    return(DrawPolygonPrimitive(image,draw_info,linecap,exception));
  }
  
  static MagickBooleanType DrawStrokePolygon(Image *image,
    const DrawInfo *draw_info,const PrimitiveInfo *primitive_info,
    ExceptionInfo *exception)
  {
    DrawInfo
      *clone_info;
  
    MagickBooleanType
      closed_path;
  
    MagickStatusType
      status;
  
    PrimitiveInfo
      *stroke_polygon;
  
    register const PrimitiveInfo
      *p,
      *q;
  
    /*
      Draw stroked polygon.
    */
    if (image->debug != MagickFalse)
      (void) LogMagickEvent(DrawEvent,GetMagickModule(),
        "    begin draw-stroke-polygon");
    clone_info=CloneDrawInfo((ImageInfo *) NULL,draw_info);
    clone_info->fill=draw_info->stroke;
    if (clone_info->fill_pattern != (Image *) NULL)
      clone_info->fill_pattern=DestroyImage(clone_info->fill_pattern);
    if (clone_info->stroke_pattern != (Image *) NULL)
      clone_info->fill_pattern=CloneImage(clone_info->stroke_pattern,0,0,
        MagickTrue,exception);
    clone_info->stroke.alpha=(MagickRealType) TransparentAlpha;
    clone_info->stroke_width=0.0;
    clone_info->fill_rule=NonZeroRule;
    status=MagickTrue;
    for (p=primitive_info; p->primitive != UndefinedPrimitive; p+=p->coordinates)
    {
      if (p->coordinates == 1)
        continue;
      stroke_polygon=TraceStrokePolygon(image,draw_info,p);
      if (stroke_polygon == (PrimitiveInfo *) NULL)
        {
          status=0;
          break;
        }
      status&=DrawPolygonPrimitive(image,clone_info,stroke_polygon,exception);
      stroke_polygon=(PrimitiveInfo *) RelinquishMagickMemory(stroke_polygon);
      if (status == 0)
        break;
      q=p+p->coordinates-1;
      closed_path=p->closed_subpath;
      if ((draw_info->linecap == RoundCap) && (closed_path == MagickFalse))
        {
          status&=DrawRoundLinecap(image,draw_info,p,exception);
          status&=DrawRoundLinecap(image,draw_info,q,exception);
        }
    }
    clone_info=DestroyDrawInfo(clone_info);
    if (image->debug != MagickFalse)
      (void) LogMagickEvent(DrawEvent,GetMagickModule(),
        "    end draw-stroke-polygon");
    return(status != 0 ? MagickTrue : MagickFalse);
  }
  
  /*
  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  %                                                                             %
  %                                                                             %
  %                                                                             %
  %   G e t A f f i n e M a t r i x                                             %
  %                                                                             %
  %                                                                             %
  %                                                                             %
  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  %
  %  GetAffineMatrix() returns an AffineMatrix initialized to the identity
  %  matrix.
  %
  %  The format of the GetAffineMatrix method is:
  %
  %      void GetAffineMatrix(AffineMatrix *affine_matrix)
  %
  %  A description of each parameter follows:
  %
  %    o affine_matrix: the affine matrix.
  %
  */
  MagickExport void GetAffineMatrix(AffineMatrix *affine_matrix)
  {
    (void) LogMagickEvent(TraceEvent,GetMagickModule(),"...");
    assert(affine_matrix != (AffineMatrix *) NULL);
    (void) memset(affine_matrix,0,sizeof(*affine_matrix));
    affine_matrix->sx=1.0;
    affine_matrix->sy=1.0;
  }
  
  /*
  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  %                                                                             %
  %                                                                             %
  %                                                                             %
  +   G e t D r a w I n f o                                                     %
  %                                                                             %
  %                                                                             %
  %                                                                             %
  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  %
  %  GetDrawInfo() initializes draw_info to default values from image_info.
  %
  %  The format of the GetDrawInfo method is:
  %
  %      void GetDrawInfo(const ImageInfo *image_info,DrawInfo *draw_info)
  %
  %  A description of each parameter follows:
  %
  %    o image_info: the image info..
  %
  %    o draw_info: the draw info.
  %
  */
  MagickExport void GetDrawInfo(const ImageInfo *image_info,DrawInfo *draw_info)
  {
    char
      *next_token;
  
    const char
      *option;
  
    ExceptionInfo
      *exception;
  
    ImageInfo
      *clone_info;
  
    /*
      Initialize draw attributes.
    */
    (void) LogMagickEvent(TraceEvent,GetMagickModule(),"...");
    assert(draw_info != (DrawInfo *) NULL);
    (void) memset(draw_info,0,sizeof(*draw_info));
    clone_info=CloneImageInfo(image_info);
    GetAffineMatrix(&draw_info->affine);
    exception=AcquireExceptionInfo();
    (void) QueryColorCompliance("#000F",AllCompliance,&draw_info->fill,
      exception);
    (void) QueryColorCompliance("#FFF0",AllCompliance,&draw_info->stroke,
      exception);
    draw_info->stroke_antialias=clone_info->antialias;
    draw_info->stroke_width=1.0;
    draw_info->fill_rule=EvenOddRule;
    draw_info->alpha=OpaqueAlpha;
    draw_info->fill_alpha=OpaqueAlpha;
    draw_info->stroke_alpha=OpaqueAlpha;
    draw_info->linecap=ButtCap;
    draw_info->linejoin=MiterJoin;
    draw_info->miterlimit=10;
    draw_info->decorate=NoDecoration;
    draw_info->pointsize=12.0;
    draw_info->undercolor.alpha=(MagickRealType) TransparentAlpha;
    draw_info->compose=OverCompositeOp;
    draw_info->render=MagickTrue;
    draw_info->clip_path=MagickFalse;
    draw_info->debug=IsEventLogging();
    if (clone_info->font != (char *) NULL)
      draw_info->font=AcquireString(clone_info->font);
    if (clone_info->density != (char *) NULL)
      draw_info->density=AcquireString(clone_info->density);
    draw_info->text_antialias=clone_info->antialias;
    if (fabs(clone_info->pointsize) >= MagickEpsilon)
      draw_info->pointsize=clone_info->pointsize;
    draw_info->border_color=clone_info->border_color;
    if (clone_info->server_name != (char *) NULL)
      draw_info->server_name=AcquireString(clone_info->server_name);
    option=GetImageOption(clone_info,"direction");
    if (option != (const char *) NULL)
      draw_info->direction=(DirectionType) ParseCommandOption(
        MagickDirectionOptions,MagickFalse,option);
    else
      draw_info->direction=UndefinedDirection;
    option=GetImageOption(clone_info,"encoding");
    if (option != (const char *) NULL)
      (void) CloneString(&draw_info->encoding,option);
    option=GetImageOption(clone_info,"family");
    if (option != (const char *) NULL)
      (void) CloneString(&draw_info->family,option);
    option=GetImageOption(clone_info,"fill");
    if (option != (const char *) NULL)
      (void) QueryColorCompliance(option,AllCompliance,&draw_info->fill,
        exception);
    option=GetImageOption(clone_info,"gravity");
    if (option != (const char *) NULL)
      draw_info->gravity=(GravityType) ParseCommandOption(MagickGravityOptions,
        MagickFalse,option);
    option=GetImageOption(clone_info,"interline-spacing");
    if (option != (const char *) NULL)
      draw_info->interline_spacing=GetDrawValue(option,&next_token);
    option=GetImageOption(clone_info,"interword-spacing");
    if (option != (const char *) NULL)
      draw_info->interword_spacing=GetDrawValue(option,&next_token);
    option=GetImageOption(clone_info,"kerning");
    if (option != (const char *) NULL)
      draw_info->kerning=GetDrawValue(option,&next_token);
    option=GetImageOption(clone_info,"stroke");
    if (option != (const char *) NULL)
      (void) QueryColorCompliance(option,AllCompliance,&draw_info->stroke,
        exception);
    option=GetImageOption(clone_info,"strokewidth");
    if (option != (const char *) NULL)
      draw_info->stroke_width=GetDrawValue(option,&next_token);
    option=GetImageOption(clone_info,"style");
    if (option != (const char *) NULL)
      draw_info->style=(StyleType) ParseCommandOption(MagickStyleOptions,
        MagickFalse,option);
    option=GetImageOption(clone_info,"undercolor");
    if (option != (const char *) NULL)
      (void) QueryColorCompliance(option,AllCompliance,&draw_info->undercolor,
        exception);
    option=GetImageOption(clone_info,"weight");
    if (option != (const char *) NULL)
      {
        ssize_t
          weight;
  
        weight=ParseCommandOption(MagickWeightOptions,MagickFalse,option);
        if (weight == -1)
          weight=(ssize_t) StringToUnsignedLong(option);
        draw_info->weight=(size_t) weight;
      }
    exception=DestroyExceptionInfo(exception);
    draw_info->signature=MagickCoreSignature;
    clone_info=DestroyImageInfo(clone_info);
  }
  
  /*
  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  %                                                                             %
  %                                                                             %
  %                                                                             %
  +   P e r m u t a t e                                                         %
  %                                                                             %
  %                                                                             %
  %                                                                             %
  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  %
  %  Permutate() returns the permuation of the (n,k).
  %
  %  The format of the Permutate method is:
  %
  %      void Permutate(ssize_t n,ssize_t k)
  %
  %  A description of each parameter follows:
  %
  %    o n:
  %
  %    o k:
  %
  %
  */
  static inline double Permutate(const ssize_t n,const ssize_t k)
  {
    double
      r;
  
    register ssize_t
      i;
  
    r=1.0;
    for (i=k+1; i <= n; i++)
      r*=i;
    for (i=1; i <= (n-k); i++)
      r/=i;
    return(r);
  }
  
  /*
  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  %                                                                             %
  %                                                                             %
  %                                                                             %
  +   T r a c e P r i m i t i v e                                               %
  %                                                                             %
  %                                                                             %
  %                                                                             %
  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  %
  %  TracePrimitive is a collection of methods for generating graphic
  %  primitives such as arcs, ellipses, paths, etc.
  %
  */
  
  static MagickBooleanType TraceArc(MVGInfo *mvg_info,const PointInfo start,
    const PointInfo end,const PointInfo degrees)
  {
    PointInfo
      center,
      radius;
  
    center.x=0.5*(end.x+start.x);
    center.y=0.5*(end.y+start.y);
    radius.x=fabs(center.x-start.x);
    radius.y=fabs(center.y-start.y);
    return(TraceEllipse(mvg_info,center,radius,degrees));
  }
  
  static MagickBooleanType TraceArcPath(MVGInfo *mvg_info,const PointInfo start,
    const PointInfo end,const PointInfo arc,const double angle,
    const MagickBooleanType large_arc,const MagickBooleanType sweep)
  {
    double
      alpha,
      beta,
      delta,
      factor,
      gamma,
      theta;
  
    MagickStatusType
      status;
  
    PointInfo
      center,
      points[3],
      radii;
  
    register double
      cosine,
      sine;
  
    PrimitiveInfo
      *primitive_info;
  
    register PrimitiveInfo
      *p;
  
    register ssize_t
      i;
  
    size_t
      arc_segments;
  
    ssize_t
      offset;
  
    offset=mvg_info->offset;
    primitive_info=(*mvg_info->primitive_info)+mvg_info->offset;
    primitive_info->coordinates=0;
    if ((fabs(start.x-end.x) < MagickEpsilon) &&
        (fabs(start.y-end.y) < MagickEpsilon))
      return(TracePoint(primitive_info,end));
    radii.x=fabs(arc.x);
    radii.y=fabs(arc.y);
    if ((radii.x < MagickEpsilon) || (radii.y < MagickEpsilon))
      return(TraceLine(primitive_info,start,end));
    cosine=cos(DegreesToRadians(fmod((double) angle,360.0)));
    sine=sin(DegreesToRadians(fmod((double) angle,360.0)));
    center.x=(double) (cosine*(end.x-start.x)/2+sine*(end.y-start.y)/2);
    center.y=(double) (cosine*(end.y-start.y)/2-sine*(end.x-start.x)/2);
    delta=(center.x*center.x)/(radii.x*radii.x)+(center.y*center.y)/
      (radii.y*radii.y);
    if (delta < MagickEpsilon)
      return(TraceLine(primitive_info,start,end));
    if (delta > 1.0)
      {
        radii.x*=sqrt((double) delta);
        radii.y*=sqrt((double) delta);
      }
    points[0].x=(double) (cosine*start.x/radii.x+sine*start.y/radii.x);
    points[0].y=(double) (cosine*start.y/radii.y-sine*start.x/radii.y);
    points[1].x=(double) (cosine*end.x/radii.x+sine*end.y/radii.x);
    points[1].y=(double) (cosine*end.y/radii.y-sine*end.x/radii.y);
    alpha=points[1].x-points[0].x;
    beta=points[1].y-points[0].y;
    if (fabs(alpha*alpha+beta*beta) < MagickEpsilon)
      return(TraceLine(primitive_info,start,end));
    factor=PerceptibleReciprocal(alpha*alpha+beta*beta)-0.25;
    if (factor <= 0.0)
      factor=0.0;
    else
      {
        factor=sqrt((double) factor);
        if (sweep == large_arc)
          factor=(-factor);
      }
    center.x=(double) ((points[0].x+points[1].x)/2-factor*beta);
    center.y=(double) ((points[0].y+points[1].y)/2+factor*alpha);
    alpha=atan2(points[0].y-center.y,points[0].x-center.x);
    theta=atan2(points[1].y-center.y,points[1].x-center.x)-alpha;
    if ((theta < 0.0) && (sweep != MagickFalse))
      theta+=2.0*MagickPI;
    else
      if ((theta > 0.0) && (sweep == MagickFalse))
        theta-=2.0*MagickPI;
    arc_segments=(size_t) ceil(fabs((double) (theta/(0.5*MagickPI+
      MagickEpsilon))));
    status=MagickTrue;
    p=primitive_info;
    for (i=0; i < (ssize_t) arc_segments; i++)
    {
      beta=0.5*((alpha+(i+1)*theta/arc_segments)-(alpha+i*theta/arc_segments));
      gamma=(8.0/3.0)*sin(fmod((double) (0.5*beta),DegreesToRadians(360.0)))*
        sin(fmod((double) (0.5*beta),DegreesToRadians(360.0)))/
        sin(fmod((double) beta,DegreesToRadians(360.0)));
      points[0].x=(double) (center.x+cos(fmod((double) (alpha+(double) i*theta/
        arc_segments),DegreesToRadians(360.0)))-gamma*sin(fmod((double) (alpha+
        (double) i*theta/arc_segments),DegreesToRadians(360.0))));
      points[0].y=(double) (center.y+sin(fmod((double) (alpha+(double) i*theta/
        arc_segments),DegreesToRadians(360.0)))+gamma*cos(fmod((double) (alpha+
        (double) i*theta/arc_segments),DegreesToRadians(360.0))));
      points[2].x=(double) (center.x+cos(fmod((double) (alpha+(double) (i+1)*
        theta/arc_segments),DegreesToRadians(360.0))));
      points[2].y=(double) (center.y+sin(fmod((double) (alpha+(double) (i+1)*
        theta/arc_segments),DegreesToRadians(360.0))));
      points[1].x=(double) (points[2].x+gamma*sin(fmod((double) (alpha+(double)
        (i+1)*theta/arc_segments),DegreesToRadians(360.0))));
      points[1].y=(double) (points[2].y-gamma*cos(fmod((double) (alpha+(double)
        (i+1)*theta/arc_segments),DegreesToRadians(360.0))));
      p->point.x=(p == primitive_info) ? start.x : (p-1)->point.x;
      p->point.y=(p == primitive_info) ? start.y : (p-1)->point.y;
      (p+1)->point.x=(double) (cosine*radii.x*points[0].x-sine*radii.y*
        points[0].y);
      (p+1)->point.y=(double) (sine*radii.x*points[0].x+cosine*radii.y*
        points[0].y);
      (p+2)->point.x=(double) (cosine*radii.x*points[1].x-sine*radii.y*
        points[1].y);
      (p+2)->point.y=(double) (sine*radii.x*points[1].x+cosine*radii.y*
        points[1].y);
      (p+3)->point.x=(double) (cosine*radii.x*points[2].x-sine*radii.y*
        points[2].y);
      (p+3)->point.y=(double) (sine*radii.x*points[2].x+cosine*radii.y*
        points[2].y);
      if (i == (ssize_t) (arc_segments-1))
        (p+3)->point=end;
      status&=TraceBezier(mvg_info,4);
      if (status == 0)
        break;
      p=(*mvg_info->primitive_info)+mvg_info->offset;
      mvg_info->offset+=p->coordinates;
      p+=p->coordinates;
    }
    if (status == 0)
      return(MagickFalse);
    mvg_info->offset=offset;
    primitive_info=(*mvg_info->primitive_info)+mvg_info->offset;
    primitive_info->coordinates=(size_t) (p-primitive_info);
    primitive_info->closed_subpath=MagickFalse;
    for (i=0; i < (ssize_t) primitive_info->coordinates; i++)
    {
      p->primitive=primitive_info->primitive;
      p--;
    }
    return(MagickTrue);
  }
  
  static MagickBooleanType TraceBezier(MVGInfo *mvg_info,
    const size_t number_coordinates)
  {
    double
      alpha,
      *coefficients,
      weight;
  
    PointInfo
      end,
      point,
      *points;
  
    PrimitiveInfo
      *primitive_info;
  
    register PrimitiveInfo
      *p;
  
    register ssize_t
      i,
      j;
  
    size_t
      control_points,
      quantum;
  
    /*
      Allocate coefficients.
    */
    primitive_info=(*mvg_info->primitive_info)+mvg_info->offset;
    quantum=number_coordinates;
    for (i=0; i < (ssize_t) number_coordinates; i++)
    {
      for (j=i+1; j < (ssize_t) number_coordinates; j++)
      {
        alpha=fabs(primitive_info[j].point.x-primitive_info[i].point.x);
        if (alpha > (double) SSIZE_MAX)
          {
            (void) ThrowMagickException(mvg_info->exception,GetMagickModule(),
              ResourceLimitError,"MemoryAllocationFailed","`%s'","");
            return(MagickFalse);
          }
        if (alpha > (double) quantum)
          quantum=(size_t) alpha;
        alpha=fabs(primitive_info[j].point.y-primitive_info[i].point.y);
        if (alpha > (double) SSIZE_MAX)
          {
            (void) ThrowMagickException(mvg_info->exception,GetMagickModule(),
              ResourceLimitError,"MemoryAllocationFailed","`%s'","");
            return(MagickFalse);
          }
        if (alpha > (double) quantum)
          quantum=(size_t) alpha;
      }
    }
    primitive_info=(*mvg_info->primitive_info)+mvg_info->offset;
    quantum=MagickMin(quantum/number_coordinates,BezierQuantum);
    coefficients=(double *) AcquireQuantumMemory(number_coordinates,
      sizeof(*coefficients));
    points=(PointInfo *) AcquireQuantumMemory(quantum,number_coordinates*
      sizeof(*points));
    if ((coefficients == (double *) NULL) || (points == (PointInfo *) NULL))
      {
        if (points != (PointInfo *) NULL)
          points=(PointInfo *) RelinquishMagickMemory(points);
        if (coefficients != (double *) NULL)
          coefficients=(double *) RelinquishMagickMemory(coefficients);
        (void) ThrowMagickException(mvg_info->exception,GetMagickModule(),
          ResourceLimitError,"MemoryAllocationFailed","`%s'","");
        return(MagickFalse);
      }
    control_points=quantum*number_coordinates;
    if (CheckPrimitiveExtent(mvg_info,control_points+1) == MagickFalse)
      {
        points=(PointInfo *) RelinquishMagickMemory(points);
        coefficients=(double *) RelinquishMagickMemory(coefficients);
        return(MagickFalse);
      }
    primitive_info=(*mvg_info->primitive_info)+mvg_info->offset;
    /*
      Compute bezier points.
    */
    end=primitive_info[number_coordinates-1].point;
    for (i=0; i < (ssize_t) number_coordinates; i++)
      coefficients[i]=Permutate((ssize_t) number_coordinates-1,i);
    weight=0.0;
    for (i=0; i < (ssize_t) control_points; i++)
    {
      p=primitive_info;
      point.x=0.0;
      point.y=0.0;
      alpha=pow((double) (1.0-weight),(double) number_coordinates-1.0);
      for (j=0; j < (ssize_t) number_coordinates; j++)
      {
        point.x+=alpha*coefficients[j]*p->point.x;
        point.y+=alpha*coefficients[j]*p->point.y;
        alpha*=weight/(1.0-weight);
        p++;
      }
      points[i]=point;
      weight+=1.0/control_points;
    }
    /*
      Bezier curves are just short segmented polys.
    */
    p=primitive_info;
    for (i=0; i < (ssize_t) control_points; i++)
    {
      if (TracePoint(p,points[i]) == MagickFalse)
        {
          points=(PointInfo *) RelinquishMagickMemory(points);
          coefficients=(double *) RelinquishMagickMemory(coefficients);
          return(MagickFalse);
        }
      p+=p->coordinates;
    }
    if (TracePoint(p,end) == MagickFalse)
      {
        points=(PointInfo *) RelinquishMagickMemory(points);
        coefficients=(double *) RelinquishMagickMemory(coefficients);
        return(MagickFalse);
      }
    p+=p->coordinates;
    primitive_info->coordinates=(size_t) (p-primitive_info);
    primitive_info->closed_subpath=MagickFalse;
    for (i=0; i < (ssize_t) primitive_info->coordinates; i++)
    {
      p->primitive=primitive_info->primitive;
      p--;
    }
    points=(PointInfo *) RelinquishMagickMemory(points);
    coefficients=(double *) RelinquishMagickMemory(coefficients);
    return(MagickTrue);
  }
  
  static MagickBooleanType TraceCircle(MVGInfo *mvg_info,const PointInfo start,
    const PointInfo end)
  {
    double
      alpha,
      beta,
      radius;
  
    PointInfo
      offset,
      degrees;
  
    alpha=end.x-start.x;
    beta=end.y-start.y;
    radius=hypot((double) alpha,(double) beta);
    offset.x=(double) radius;
    offset.y=(double) radius;
    degrees.x=0.0;
    degrees.y=360.0;
    return(TraceEllipse(mvg_info,start,offset,degrees));
  }
  
  static MagickBooleanType TraceEllipse(MVGInfo *mvg_info,const PointInfo center,
    const PointInfo radii,const PointInfo arc)
  {
    double
      coordinates,
      delta,
      step,
      x,
      y;
  
    PointInfo
      angle,
      point;
  
    PrimitiveInfo
      *primitive_info;
  
    register PrimitiveInfo
      *p;
  
    register ssize_t
      i;
  
    /*
      Ellipses are just short segmented polys.
    */
    primitive_info=(*mvg_info->primitive_info)+mvg_info->offset;
    primitive_info->coordinates=0;
    if ((fabs(radii.x) < MagickEpsilon) || (fabs(radii.y) < MagickEpsilon))
      return(MagickTrue);
    delta=2.0*PerceptibleReciprocal(MagickMax(radii.x,radii.y));
    step=MagickPI/8.0;
    if ((delta >= 0.0) && (delta < (MagickPI/8.0)))
      step=MagickPI/4.0/(MagickPI*PerceptibleReciprocal(delta)/2.0);
    angle.x=DegreesToRadians(arc.x);
    y=arc.y;
    while (y < arc.x)
      y+=360.0;
    angle.y=DegreesToRadians(y);
    coordinates=ceil((angle.y-angle.x)/step+1.0);
    if (coordinates > (double) SSIZE_MAX)
      {
        (void) ThrowMagickException(mvg_info->exception,GetMagickModule(),
          ResourceLimitError,"MemoryAllocationFailed","`%s'","");
        return(MagickFalse);
      }
    if (CheckPrimitiveExtent(mvg_info,(size_t) coordinates) == MagickFalse)
      return(MagickFalse);
    primitive_info=(*mvg_info->primitive_info)+mvg_info->offset;
    for (p=primitive_info; angle.x < angle.y; angle.x+=step)
    {
      point.x=cos(fmod(angle.x,DegreesToRadians(360.0)))*radii.x+center.x;
      point.y=sin(fmod(angle.x,DegreesToRadians(360.0)))*radii.y+center.y;
      if (TracePoint(p,point) == MagickFalse)
        return(MagickFalse);
      p+=p->coordinates;
    }
    point.x=cos(fmod(angle.y,DegreesToRadians(360.0)))*radii.x+center.x;
    point.y=sin(fmod(angle.y,DegreesToRadians(360.0)))*radii.y+center.y;
    if (TracePoint(p,point) == MagickFalse)
      return(MagickFalse);
    p+=p->coordinates;
    primitive_info->coordinates=(size_t) (p-primitive_info);
    primitive_info->closed_subpath=MagickFalse;
    x=fabs(primitive_info[0].point.x-
      primitive_info[primitive_info->coordinates-1].point.x);
    y=fabs(primitive_info[0].point.y-
      primitive_info[primitive_info->coordinates-1].point.y);
    if ((x < MagickEpsilon) && (y < MagickEpsilon))
      primitive_info->closed_subpath=MagickTrue;
    for (i=0; i < (ssize_t) primitive_info->coordinates; i++)
    {
      p->primitive=primitive_info->primitive;
      p--;
    }
    return(MagickTrue);
  }
  
  static MagickBooleanType TraceLine(PrimitiveInfo *primitive_info,
    const PointInfo start,const PointInfo end)
  {
    if (TracePoint(primitive_info,start) == MagickFalse)
      return(MagickFalse);
    if ((fabs(start.x-end.x) < MagickEpsilon) &&
        (fabs(start.y-end.y) < MagickEpsilon))
      {
        primitive_info->primitive=PointPrimitive;
        primitive_info->coordinates=1;
        return(MagickTrue);
      }
    if (TracePoint(primitive_info+1,end) == MagickFalse)
      return(MagickFalse);
    (primitive_info+1)->primitive=primitive_info->primitive;
    primitive_info->coordinates=2;
    primitive_info->closed_subpath=MagickFalse;
    return(MagickTrue);
  }
  
  static ssize_t TracePath(MVGInfo *mvg_info,const char *path,
    ExceptionInfo *exception)
  {
    char
      *next_token,
      token[MagickPathExtent];
  
    const char
      *p;
  
    double
      x,
      y;
  
    int
      attribute,
      last_attribute;
  
    MagickBooleanType
      status;
  
    PointInfo
      end = {0.0, 0.0},
      points[4] = { {0.0, 0.0}, {0.0, 0.0}, {0.0, 0.0}, {0.0, 0.0} },
      point = {0.0, 0.0},
      start = {0.0, 0.0};
  
    PrimitiveInfo
      *primitive_info;
  
    PrimitiveType
      primitive_type;
  
    register PrimitiveInfo
      *q;
  
    register ssize_t
      i;
  
    size_t
      number_coordinates,
      z_count;
  
    ssize_t
      subpath_offset;
  
    subpath_offset=mvg_info->offset;
    primitive_info=(*mvg_info->primitive_info)+mvg_info->offset;
    status=MagickTrue;
    attribute=0;
    number_coordinates=0;
    z_count=0;
    primitive_type=primitive_info->primitive;
    q=primitive_info;
    for (p=path; *p != '\0'; )
    {
      if (status == MagickFalse)
        break;
      while (isspace((int) ((unsigned char) *p)) != 0)
        p++;
      if (*p == '\0')
        break;
      last_attribute=attribute;
      attribute=(int) (*p++);
      switch (attribute)
      {
        case 'a':
        case 'A':
        {
          double
            angle = 0.0;
  
          MagickBooleanType
            large_arc = MagickFalse,
            sweep = MagickFalse;
  
          PointInfo
            arc = {0.0, 0.0};
  
          /*
            Elliptical arc.
          */
          do
          {
            (void) GetNextToken(p,&p,MagickPathExtent,token);
            if (*token == ',')
              (void) GetNextToken(p,&p,MagickPathExtent,token);
            arc.x=GetDrawValue(token,&next_token);
            if (token == next_token)
              ThrowPointExpectedException(token,exception);
            (void) GetNextToken(p,&p,MagickPathExtent,token);
            if (*token == ',')
              (void) GetNextToken(p,&p,MagickPathExtent,token);
            arc.y=GetDrawValue(token,&next_token);
            if (token == next_token)
              ThrowPointExpectedException(token,exception);
            (void) GetNextToken(p,&p,MagickPathExtent,token);
            if (*token == ',')
              (void) GetNextToken(p,&p,MagickPathExtent,token);
            angle=GetDrawValue(token,&next_token);
            if (token == next_token)
              ThrowPointExpectedException(token,exception);
            (void) GetNextToken(p,&p,MagickPathExtent,token);
            if (*token == ',')
              (void) GetNextToken(p,&p,MagickPathExtent,token);
            large_arc=StringToLong(token) != 0 ? MagickTrue : MagickFalse;
            (void) GetNextToken(p,&p,MagickPathExtent,token);
            if (*token == ',')
              (void) GetNextToken(p,&p,MagickPathExtent,token);
            sweep=StringToLong(token) != 0 ? MagickTrue : MagickFalse;
            if (*token == ',')
              (void) GetNextToken(p,&p,MagickPathExtent,token);
            (void) GetNextToken(p,&p,MagickPathExtent,token);
            if (*token == ',')
              (void) GetNextToken(p,&p,MagickPathExtent,token);
            x=GetDrawValue(token,&next_token);
            if (token == next_token)
              ThrowPointExpectedException(token,exception);
            (void) GetNextToken(p,&p,MagickPathExtent,token);
            if (*token == ',')
              (void) GetNextToken(p,&p,MagickPathExtent,token);
            y=GetDrawValue(token,&next_token);
            if (token == next_token)
              ThrowPointExpectedException(token,exception);
            end.x=(double) (attribute == (int) 'A' ? x : point.x+x);
            end.y=(double) (attribute == (int) 'A' ? y : point.y+y);
            if (TraceArcPath(mvg_info,point,end,arc,angle,large_arc,sweep) == MagickFalse)
              return(-1);
            q=(*mvg_info->primitive_info)+mvg_info->offset;
            mvg_info->offset+=q->coordinates;
            q+=q->coordinates;
            point=end;
            while (isspace((int) ((unsigned char) *p)) != 0)
              p++;
            if (*p == ',')
              p++;
          } while (IsPoint(p) != MagickFalse);
          break;
        }
        case 'c':
        case 'C':
        {
          /*
            Cubic Bézier curve.
          */
          do
          {
            points[0]=point;
            for (i=1; i < 4; i++)
            {
              (void) GetNextToken(p,&p,MagickPathExtent,token);
              if (*token == ',')
                (void) GetNextToken(p,&p,MagickPathExtent,token);
              x=GetDrawValue(token,&next_token);
              if (token == next_token)
                ThrowPointExpectedException(token,exception);
              (void) GetNextToken(p,&p,MagickPathExtent,token);
              if (*token == ',')
                (void) GetNextToken(p,&p,MagickPathExtent,token);
              y=GetDrawValue(token,&next_token);
              if (token == next_token)
                ThrowPointExpectedException(token,exception);
              end.x=(double) (attribute == (int) 'C' ? x : point.x+x);
              end.y=(double) (attribute == (int) 'C' ? y : point.y+y);
              points[i]=end;
            }
            for (i=0; i < 4; i++)
              (q+i)->point=points[i];
            if (TraceBezier(mvg_info,4) == MagickFalse)
              return(-1);
            q=(*mvg_info->primitive_info)+mvg_info->offset;
            mvg_info->offset+=q->coordinates;
            q+=q->coordinates;
            point=end;
            while (isspace((int) ((unsigned char) *p)) != 0)
              p++;
            if (*p == ',')
              p++;
          } while (IsPoint(p) != MagickFalse);
          break;
        }
        case 'H':
        case 'h':
        {
          do
          {
            (void) GetNextToken(p,&p,MagickPathExtent,token);
            if (*token == ',')
              (void) GetNextToken(p,&p,MagickPathExtent,token);
            x=GetDrawValue(token,&next_token);
            if (token == next_token)
              ThrowPointExpectedException(token,exception);
            point.x=(double) (attribute == (int) 'H' ? x: point.x+x);
            if (CheckPrimitiveExtent(mvg_info,PrimitiveExtentPad) == MagickFalse)
              return(-1);
            q=(*mvg_info->primitive_info)+mvg_info->offset;
            if (TracePoint(q,point) == MagickFalse)
              return(-1);
            mvg_info->offset+=q->coordinates;
            q+=q->coordinates;
            while (isspace((int) ((unsigned char) *p)) != 0)
              p++;
            if (*p == ',')
              p++;
          } while (IsPoint(p) != MagickFalse);
          break;
        }
        case 'l':
        case 'L':
        {
          /*
            Line to.
          */
          do
          {
            (void) GetNextToken(p,&p,MagickPathExtent,token);
            if (*token == ',')
              (void) GetNextToken(p,&p,MagickPathExtent,token);
            x=GetDrawValue(token,&next_token);
            if (token == next_token)
              ThrowPointExpectedException(token,exception);
            (void) GetNextToken(p,&p,MagickPathExtent,token);
            if (*token == ',')
              (void) GetNextToken(p,&p,MagickPathExtent,token);
            y=GetDrawValue(token,&next_token);
            if (token == next_token)
              ThrowPointExpectedException(token,exception);
            point.x=(double) (attribute == (int) 'L' ? x : point.x+x);
            point.y=(double) (attribute == (int) 'L' ? y : point.y+y);
            if (CheckPrimitiveExtent(mvg_info,PrimitiveExtentPad) == MagickFalse)
              return(-1);
            q=(*mvg_info->primitive_info)+mvg_info->offset;
            if (TracePoint(q,point) == MagickFalse)
              return(-1);
            mvg_info->offset+=q->coordinates;
            q+=q->coordinates;
            while (isspace((int) ((unsigned char) *p)) != 0)
              p++;
            if (*p == ',')
              p++;
          } while (IsPoint(p) != MagickFalse);
          break;
        }
        case 'M':
        case 'm':
        {
          /*
            Move to.
          */
          if (mvg_info->offset != subpath_offset)
            {
              primitive_info=(*mvg_info->primitive_info)+subpath_offset;
              primitive_info->coordinates=(size_t) (q-primitive_info);
              number_coordinates+=primitive_info->coordinates;
              primitive_info=q;
              subpath_offset=mvg_info->offset;
            }
          i=0;
          do
          {
            (void) GetNextToken(p,&p,MagickPathExtent,token);
            if (*token == ',')
              (void) GetNextToken(p,&p,MagickPathExtent,token);
            x=GetDrawValue(token,&next_token);
            if (token == next_token)
              ThrowPointExpectedException(token,exception);
            (void) GetNextToken(p,&p,MagickPathExtent,token);
            if (*token == ',')
              (void) GetNextToken(p,&p,MagickPathExtent,token);
            y=GetDrawValue(token,&next_token);
            if (token == next_token)
              ThrowPointExpectedException(token,exception);
            point.x=(double) (attribute == (int) 'M' ? x : point.x+x);
            point.y=(double) (attribute == (int) 'M' ? y : point.y+y);
            if (i == 0)
              start=point;
            i++;
            if (CheckPrimitiveExtent(mvg_info,PrimitiveExtentPad) == MagickFalse)
              return(-1);
            q=(*mvg_info->primitive_info)+mvg_info->offset;
            if (TracePoint(q,point) == MagickFalse)
              return(-1);
            mvg_info->offset+=q->coordinates;
            q+=q->coordinates;
            while (isspace((int) ((unsigned char) *p)) != 0)
              p++;
            if (*p == ',')
              p++;
          } while (IsPoint(p) != MagickFalse);
          break;
        }
        case 'q':
        case 'Q':
        {
          /*
            Quadratic Bézier curve.
          */
          do
          {
            points[0]=point;
            for (i=1; i < 3; i++)
            {
              (void) GetNextToken(p,&p,MagickPathExtent,token);
              if (*token == ',')
                (void) GetNextToken(p,&p,MagickPathExtent,token);
              x=GetDrawValue(token,&next_token);
              if (token == next_token)
                ThrowPointExpectedException(token,exception);
              (void) GetNextToken(p,&p,MagickPathExtent,token);
              if (*token == ',')
                (void) GetNextToken(p,&p,MagickPathExtent,token);
              y=GetDrawValue(token,&next_token);
              if (token == next_token)
                ThrowPointExpectedException(token,exception);
              if (*p == ',')
                p++;
              end.x=(double) (attribute == (int) 'Q' ? x : point.x+x);
              end.y=(double) (attribute == (int) 'Q' ? y : point.y+y);
              points[i]=end;
            }
            for (i=0; i < 3; i++)
              (q+i)->point=points[i];
            if (TraceBezier(mvg_info,3) == MagickFalse)
              return(-1);
            q=(*mvg_info->primitive_info)+mvg_info->offset;
            mvg_info->offset+=q->coordinates;
            q+=q->coordinates;
            point=end;
            while (isspace((int) ((unsigned char) *p)) != 0)
              p++;
            if (*p == ',')
              p++;
          } while (IsPoint(p) != MagickFalse);
          break;
        }
        case 's':
        case 'S':
        {
          /*
            Cubic Bézier curve.
          */
          do
          {
            points[0]=points[3];
            points[1].x=2.0*points[3].x-points[2].x;
            points[1].y=2.0*points[3].y-points[2].y;
            for (i=2; i < 4; i++)
            {
              (void) GetNextToken(p,&p,MagickPathExtent,token);
              if (*token == ',')
                (void) GetNextToken(p,&p,MagickPathExtent,token);
              x=GetDrawValue(token,&next_token);
              if (token == next_token)
                ThrowPointExpectedException(token,exception);
              (void) GetNextToken(p,&p,MagickPathExtent,token);
              if (*token == ',')
                (void) GetNextToken(p,&p,MagickPathExtent,token);
              y=GetDrawValue(token,&next_token);
              if (token == next_token)
                ThrowPointExpectedException(token,exception);
              if (*p == ',')
                p++;
              end.x=(double) (attribute == (int) 'S' ? x : point.x+x);
              end.y=(double) (attribute == (int) 'S' ? y : point.y+y);
              points[i]=end;
            }
            if (strchr("CcSs",last_attribute) == (char *) NULL)
              {
                points[0]=point;
                points[1]=point;
              }
            for (i=0; i < 4; i++)
              (q+i)->point=points[i];
            if (TraceBezier(mvg_info,4) == MagickFalse)
              return(-1);
            q=(*mvg_info->primitive_info)+mvg_info->offset;
            mvg_info->offset+=q->coordinates;
            q+=q->coordinates;
            point=end;
            last_attribute=attribute;
            while (isspace((int) ((unsigned char) *p)) != 0)
              p++;
            if (*p == ',')
              p++;
          } while (IsPoint(p) != MagickFalse);
          break;
        }
        case 't':
        case 'T':
        {
          /*
            Quadratic Bézier curve.
          */
          do
          {
            points[0]=points[2];
            points[1].x=2.0*points[2].x-points[1].x;
            points[1].y=2.0*points[2].y-points[1].y;
            for (i=2; i < 3; i++)
            {
              (void) GetNextToken(p,&p,MagickPathExtent,token);
              if (*token == ',')
                (void) GetNextToken(p,&p,MagickPathExtent,token);
              x=GetDrawValue(token,&next_token);
              if (token == next_token)
                ThrowPointExpectedException(token,exception);
              (void) GetNextToken(p,&p,MagickPathExtent,token);
              if (*token == ',')
                (void) GetNextToken(p,&p,MagickPathExtent,token);
              y=GetDrawValue(token,&next_token);
              if (token == next_token)
                ThrowPointExpectedException(token,exception);
              end.x=(double) (attribute == (int) 'T' ? x : point.x+x);
              end.y=(double) (attribute == (int) 'T' ? y : point.y+y);
              points[i]=end;
            }
            if (status == MagickFalse)
              break;
            if (strchr("QqTt",last_attribute) == (char *) NULL)
              {
                points[0]=point;
                points[1]=point;
              }
            for (i=0; i < 3; i++)
              (q+i)->point=points[i];
            if (TraceBezier(mvg_info,3) == MagickFalse)
              return(-1);
            q=(*mvg_info->primitive_info)+mvg_info->offset;
            mvg_info->offset+=q->coordinates;
            q+=q->coordinates;
            point=end;
            last_attribute=attribute;
            while (isspace((int) ((unsigned char) *p)) != 0)
              p++;
            if (*p == ',')
              p++;
          } while (IsPoint(p) != MagickFalse);
          break;
        }
        case 'v':
        case 'V':
        {
          /*
            Line to.
          */
          do
          {
            (void) GetNextToken(p,&p,MagickPathExtent,token);
            if (*token == ',')
              (void) GetNextToken(p,&p,MagickPathExtent,token);
            y=GetDrawValue(token,&next_token);
            if (token == next_token)
              ThrowPointExpectedException(token,exception);
            point.y=(double) (attribute == (int) 'V' ? y : point.y+y);
            if (CheckPrimitiveExtent(mvg_info,PrimitiveExtentPad) == MagickFalse)
              return(-1);
            q=(*mvg_info->primitive_info)+mvg_info->offset;
            if (TracePoint(q,point) == MagickFalse)
              return(-1);
            mvg_info->offset+=q->coordinates;
            q+=q->coordinates;
            while (isspace((int) ((unsigned char) *p)) != 0)
              p++;
            if (*p == ',')
              p++;
          } while (IsPoint(p) != MagickFalse);
          break;
        }
        case 'z':
        case 'Z':
        {
          /*
            Close path.
          */
          point=start;
          if (CheckPrimitiveExtent(mvg_info,PrimitiveExtentPad) == MagickFalse)
            return(-1);
          q=(*mvg_info->primitive_info)+mvg_info->offset;
          if (TracePoint(q,point) == MagickFalse)
            return(-1);
          mvg_info->offset+=q->coordinates;
          q+=q->coordinates;
          primitive_info=(*mvg_info->primitive_info)+subpath_offset;
          primitive_info->coordinates=(size_t) (q-primitive_info);
          primitive_info->closed_subpath=MagickTrue;
          number_coordinates+=primitive_info->coordinates;
          primitive_info=q;
          subpath_offset=mvg_info->offset;
          z_count++;
          break;
        }
        default:
        {
          ThrowPointExpectedException(token,exception);
          break;
        }
      }
    }
    if (status == MagickFalse)
      return(-1);
    primitive_info=(*mvg_info->primitive_info)+subpath_offset;
    primitive_info->coordinates=(size_t) (q-primitive_info);
    number_coordinates+=primitive_info->coordinates;
    for (i=0; i < (ssize_t) number_coordinates; i++)
    {
      q--;
      q->primitive=primitive_type;
      if (z_count > 1)
        q->method=FillToBorderMethod;
    }
    q=primitive_info;
    return((ssize_t) number_coordinates);
  }
  
  static MagickBooleanType TraceRectangle(PrimitiveInfo *primitive_info,
    const PointInfo start,const PointInfo end)
  {
    PointInfo
      point;
  
    register PrimitiveInfo
      *p;
  
    register ssize_t
      i;
  
    p=primitive_info;
    if (TracePoint(p,start) == MagickFalse)
      return(MagickFalse);
    p+=p->coordinates;
    point.x=start.x;
    point.y=end.y;
    if (TracePoint(p,point) == MagickFalse)
      return(MagickFalse);
    p+=p->coordinates;
    if (TracePoint(p,end) == MagickFalse)
      return(MagickFalse);
    p+=p->coordinates;
    point.x=end.x;
    point.y=start.y;
    if (TracePoint(p,point) == MagickFalse)
      return(MagickFalse);
    p+=p->coordinates;
    if (TracePoint(p,start) == MagickFalse)
      return(MagickFalse);
    p+=p->coordinates;
    primitive_info->coordinates=(size_t) (p-primitive_info);
    primitive_info->closed_subpath=MagickTrue;
    for (i=0; i < (ssize_t) primitive_info->coordinates; i++)
    {
      p->primitive=primitive_info->primitive;
      p--;
    }
    return(MagickTrue);
  }
  
  static MagickBooleanType TraceRoundRectangle(MVGInfo *mvg_info,
    const PointInfo start,const PointInfo end,PointInfo arc)
  {
    PointInfo
      degrees,
      point,
      segment;
  
    PrimitiveInfo
      *primitive_info;
  
    register PrimitiveInfo
      *p;
  
    register ssize_t
      i;
  
    ssize_t
      offset;
  
    offset=mvg_info->offset;
    segment.x=fabs(end.x-start.x);
    segment.y=fabs(end.y-start.y);
    if ((segment.x < MagickEpsilon) || (segment.y < MagickEpsilon))
      {
        (*mvg_info->primitive_info+mvg_info->offset)->coordinates=0;
        return(MagickTrue);
      }
    if (arc.x > (0.5*segment.x))
      arc.x=0.5*segment.x;
    if (arc.y > (0.5*segment.y))
      arc.y=0.5*segment.y;
    point.x=start.x+segment.x-arc.x;
    point.y=start.y+arc.y;
    degrees.x=270.0;
    degrees.y=360.0;
    if (TraceEllipse(mvg_info,point,arc,degrees) == MagickFalse)
      return(MagickFalse);
    p=(*mvg_info->primitive_info)+mvg_info->offset;
    mvg_info->offset+=p->coordinates;
    point.x=start.x+segment.x-arc.x;
    point.y=start.y+segment.y-arc.y;
    degrees.x=0.0;
    degrees.y=90.0;
    if (TraceEllipse(mvg_info,point,arc,degrees) == MagickFalse)
      return(MagickFalse);
    p=(*mvg_info->primitive_info)+mvg_info->offset;
    mvg_info->offset+=p->coordinates;
    point.x=start.x+arc.x;
    point.y=start.y+segment.y-arc.y;
    degrees.x=90.0;
    degrees.y=180.0;
    if (TraceEllipse(mvg_info,point,arc,degrees) == MagickFalse)
      return(MagickFalse);
    p=(*mvg_info->primitive_info)+mvg_info->offset;
    mvg_info->offset+=p->coordinates;
    point.x=start.x+arc.x;
    point.y=start.y+arc.y;
    degrees.x=180.0;
    degrees.y=270.0;
    if (TraceEllipse(mvg_info,point,arc,degrees) == MagickFalse)
      return(MagickFalse);
    p=(*mvg_info->primitive_info)+mvg_info->offset;
    mvg_info->offset+=p->coordinates;
    if (CheckPrimitiveExtent(mvg_info,PrimitiveExtentPad) == MagickFalse)
      return(MagickFalse);
    p=(*mvg_info->primitive_info)+mvg_info->offset;
    if (TracePoint(p,(*mvg_info->primitive_info+offset)->point) == MagickFalse)
      return(MagickFalse);
    p+=p->coordinates;
    mvg_info->offset=offset;
    primitive_info=(*mvg_info->primitive_info)+offset;
    primitive_info->coordinates=(size_t) (p-primitive_info);
    primitive_info->closed_subpath=MagickTrue;
    for (i=0; i < (ssize_t) primitive_info->coordinates; i++)
    {
      p->primitive=primitive_info->primitive;
      p--;
    }
    return(MagickTrue);
  }
  
  static MagickBooleanType TraceSquareLinecap(PrimitiveInfo *primitive_info,
    const size_t number_vertices,const double offset)
  {
    double
      distance;
  
    register double
      dx,
      dy;
  
    register ssize_t
      i;
  
    ssize_t
      j;
  
    dx=0.0;
    dy=0.0;
    for (i=1; i < (ssize_t) number_vertices; i++)
    {
      dx=primitive_info[0].point.x-primitive_info[i].point.x;
      dy=primitive_info[0].point.y-primitive_info[i].point.y;
      if ((fabs((double) dx) >= MagickEpsilon) ||
          (fabs((double) dy) >= MagickEpsilon))
        break;
    }
    if (i == (ssize_t) number_vertices)
      i=(ssize_t) number_vertices-1L;
    distance=hypot((double) dx,(double) dy);
    primitive_info[0].point.x=(double) (primitive_info[i].point.x+
      dx*(distance+offset)/distance);
    primitive_info[0].point.y=(double) (primitive_info[i].point.y+
      dy*(distance+offset)/distance);
    for (j=(ssize_t) number_vertices-2; j >= 0;  j--)
    {
      dx=primitive_info[number_vertices-1].point.x-primitive_info[j].point.x;
      dy=primitive_info[number_vertices-1].point.y-primitive_info[j].point.y;
      if ((fabs((double) dx) >= MagickEpsilon) ||
          (fabs((double) dy) >= MagickEpsilon))
        break;
    }
    distance=hypot((double) dx,(double) dy);
    primitive_info[number_vertices-1].point.x=(double) (primitive_info[j].point.x+
      dx*(distance+offset)/distance);
    primitive_info[number_vertices-1].point.y=(double) (primitive_info[j].point.y+
      dy*(distance+offset)/distance);
    return(MagickTrue);
  }
  
  static PrimitiveInfo *TraceStrokePolygon(const Image *image,
    const DrawInfo *draw_info,const PrimitiveInfo *primitive_info)
  {
  #define MaxStrokePad  (6*BezierQuantum+360)
  #define CheckPathExtent(pad_p,pad_q) \
  {   \
    if ((pad_p) > MaxBezierCoordinates) \
      stroke_p=(PointInfo *) RelinquishMagickMemory(stroke_p); \
    else \
      if ((ssize_t) (p+(pad_p)) >= (ssize_t) extent_p) \
        { \
          if (~extent_p < (pad_p)) \
            stroke_p=(PointInfo *) RelinquishMagickMemory(stroke_p); \
          else \
            { \
              extent_p+=(pad_p); \
              stroke_p=(PointInfo *) ResizeQuantumMemory(stroke_p,extent_p+ \
                MaxStrokePad,sizeof(*stroke_p)); \
            } \
        } \
    if ((pad_q) > MaxBezierCoordinates) \
      stroke_q=(PointInfo *) RelinquishMagickMemory(stroke_q); \
    else \
      if ((ssize_t) (q+(pad_q)) >= (ssize_t) extent_q) \
        { \
          if (~extent_q < (pad_q)) \
            stroke_q=(PointInfo *) RelinquishMagickMemory(stroke_q); \
          else \
            { \
              extent_q+=(pad_q); \
              stroke_q=(PointInfo *) ResizeQuantumMemory(stroke_q,extent_q+ \
                MaxStrokePad,sizeof(*stroke_q)); \
            } \
        } \
    if ((stroke_p == (PointInfo *) NULL) || (stroke_q == (PointInfo *) NULL)) \
      { \
        if (stroke_p != (PointInfo *) NULL) \
          stroke_p=(PointInfo *) RelinquishMagickMemory(stroke_p); \
        if (stroke_q != (PointInfo *) NULL) \
          stroke_q=(PointInfo *) RelinquishMagickMemory(stroke_q); \
        polygon_primitive=(PrimitiveInfo *) \
          RelinquishMagickMemory(polygon_primitive); \
        return((PrimitiveInfo *) NULL); \
      } \
  }
  
    typedef struct _StrokeSegment
    {
      double
        p,
        q;
    } StrokeSegment;
  
    double
      delta_theta,
      dot_product,
      mid,
      miterlimit;
  
    MagickBooleanType
      closed_path;
  
    PointInfo
      box_p[5],
      box_q[5],
      center,
      offset,
      *stroke_p,
      *stroke_q;
  
    PrimitiveInfo
      *polygon_primitive,
      *stroke_polygon;
  
    register ssize_t
      i;
  
    size_t
      arc_segments,
      extent_p,
      extent_q,
      number_vertices;
  
    ssize_t
      j,
      n,
      p,
      q;
  
    StrokeSegment
      dx = {0.0, 0.0},
      dy = {0.0, 0.0},
      inverse_slope = {0.0, 0.0},
      slope = {0.0, 0.0},
      theta = {0.0, 0.0};
  
    /*
      Allocate paths.
    */
    number_vertices=primitive_info->coordinates;
    polygon_primitive=(PrimitiveInfo *) AcquireQuantumMemory((size_t)
      number_vertices+2UL,sizeof(*polygon_primitive));
    if (polygon_primitive == (PrimitiveInfo *) NULL)
      return((PrimitiveInfo *) NULL);
    (void) memcpy(polygon_primitive,primitive_info,(size_t) number_vertices*
      sizeof(*polygon_primitive));
    offset.x=primitive_info[number_vertices-1].point.x-primitive_info[0].point.x;
    offset.y=primitive_info[number_vertices-1].point.y-primitive_info[0].point.y;
    closed_path=(fabs(offset.x) < MagickEpsilon) &&
      (fabs(offset.y) < MagickEpsilon) ?  MagickTrue : MagickFalse;
    if (((draw_info->linejoin == RoundJoin) ||
         (draw_info->linejoin == MiterJoin)) && (closed_path != MagickFalse))
      {
        polygon_primitive[number_vertices]=primitive_info[1];
        number_vertices++;
      }
    polygon_primitive[number_vertices].primitive=UndefinedPrimitive;
    /*
      Compute the slope for the first line segment, p.
    */
    dx.p=0.0;
    dy.p=0.0;
    for (n=1; n < (ssize_t) number_vertices; n++)
    {
      dx.p=polygon_primitive[n].point.x-polygon_primitive[0].point.x;
      dy.p=polygon_primitive[n].point.y-polygon_primitive[0].point.y;
      if ((fabs(dx.p) >= MagickEpsilon) || (fabs(dy.p) >= MagickEpsilon))
        break;
    }
    if (n == (ssize_t) number_vertices)
      {
        if ((draw_info->linecap != RoundCap) || (closed_path != MagickFalse))
          {
            /*
              Zero length subpath.
            */
            stroke_polygon=(PrimitiveInfo *) AcquireCriticalMemory(
              sizeof(*stroke_polygon));
            stroke_polygon[0]=polygon_primitive[0];
            stroke_polygon[0].coordinates=0;
            polygon_primitive=(PrimitiveInfo *) RelinquishMagickMemory(
              polygon_primitive);
            return(stroke_polygon);
          }
        n=(ssize_t) number_vertices-1L;
      }
    extent_p=2*number_vertices;
    extent_q=2*number_vertices;
    stroke_p=(PointInfo *) AcquireQuantumMemory((size_t) extent_p+MaxStrokePad,
      sizeof(*stroke_p));
    stroke_q=(PointInfo *) AcquireQuantumMemory((size_t) extent_q+MaxStrokePad,
      sizeof(*stroke_q));
    if ((stroke_p == (PointInfo *) NULL) || (stroke_q == (PointInfo *) NULL))
      {
        if (stroke_p != (PointInfo *) NULL)
          stroke_p=(PointInfo *) RelinquishMagickMemory(stroke_p);
        if (stroke_q != (PointInfo *) NULL)
          stroke_q=(PointInfo *) RelinquishMagickMemory(stroke_q);
        polygon_primitive=(PrimitiveInfo *)
          RelinquishMagickMemory(polygon_primitive);
        return((PrimitiveInfo *) NULL);
      }
    slope.p=0.0;
    inverse_slope.p=0.0;
    if (fabs(dx.p) < MagickEpsilon)
      {
        if (dx.p >= 0.0)
          slope.p=dy.p < 0.0 ? -1.0/MagickEpsilon : 1.0/MagickEpsilon;
        else
          slope.p=dy.p < 0.0 ? 1.0/MagickEpsilon : -1.0/MagickEpsilon;
      }
    else
      if (fabs(dy.p) < MagickEpsilon)
        {
          if (dy.p >= 0.0)
            inverse_slope.p=dx.p < 0.0 ? -1.0/MagickEpsilon : 1.0/MagickEpsilon;
          else
            inverse_slope.p=dx.p < 0.0 ? 1.0/MagickEpsilon : -1.0/MagickEpsilon;
        }
      else
        {
          slope.p=dy.p/dx.p;
          inverse_slope.p=(-1.0/slope.p);
        }
    mid=ExpandAffine(&draw_info->affine)*draw_info->stroke_width/2.0;
    miterlimit=(double) (draw_info->miterlimit*draw_info->miterlimit*mid*mid);
    if ((draw_info->linecap == SquareCap) && (closed_path == MagickFalse))
      (void) TraceSquareLinecap(polygon_primitive,number_vertices,mid);
    offset.x=sqrt((double) (mid*mid/(inverse_slope.p*inverse_slope.p+1.0)));
    offset.y=(double) (offset.x*inverse_slope.p);
    if ((dy.p*offset.x-dx.p*offset.y) > 0.0)
      {
        box_p[0].x=polygon_primitive[0].point.x-offset.x;
        box_p[0].y=polygon_primitive[0].point.y-offset.x*inverse_slope.p;
        box_p[1].x=polygon_primitive[n].point.x-offset.x;
        box_p[1].y=polygon_primitive[n].point.y-offset.x*inverse_slope.p;
        box_q[0].x=polygon_primitive[0].point.x+offset.x;
        box_q[0].y=polygon_primitive[0].point.y+offset.x*inverse_slope.p;
        box_q[1].x=polygon_primitive[n].point.x+offset.x;
        box_q[1].y=polygon_primitive[n].point.y+offset.x*inverse_slope.p;
      }
    else
      {
        box_p[0].x=polygon_primitive[0].point.x+offset.x;
        box_p[0].y=polygon_primitive[0].point.y+offset.y;
        box_p[1].x=polygon_primitive[n].point.x+offset.x;
        box_p[1].y=polygon_primitive[n].point.y+offset.y;
        box_q[0].x=polygon_primitive[0].point.x-offset.x;
        box_q[0].y=polygon_primitive[0].point.y-offset.y;
        box_q[1].x=polygon_primitive[n].point.x-offset.x;
        box_q[1].y=polygon_primitive[n].point.y-offset.y;
      }
    /*
      Create strokes for the line join attribute: bevel, miter, round.
    */
    p=0;
    q=0;
    stroke_q[p++]=box_q[0];
    stroke_p[q++]=box_p[0];
    for (i=(ssize_t) n+1; i < (ssize_t) number_vertices; i++)
    {
      /*
        Compute the slope for this line segment, q.
      */
      dx.q=polygon_primitive[i].point.x-polygon_primitive[n].point.x;
      dy.q=polygon_primitive[i].point.y-polygon_primitive[n].point.y;
      dot_product=dx.q*dx.q+dy.q*dy.q;
      if (dot_product < 0.25)
        continue;
      slope.q=0.0;
      inverse_slope.q=0.0;
      if (fabs(dx.q) < MagickEpsilon)
        {
          if (dx.q >= 0.0)
            slope.q=dy.q < 0.0 ? -1.0/MagickEpsilon : 1.0/MagickEpsilon;
          else
            slope.q=dy.q < 0.0 ? 1.0/MagickEpsilon : -1.0/MagickEpsilon;
        }
      else
        if (fabs(dy.q) < MagickEpsilon)
          {
            if (dy.q >= 0.0)
              inverse_slope.q=dx.q < 0.0 ? -1.0/MagickEpsilon : 1.0/MagickEpsilon;
            else
              inverse_slope.q=dx.q < 0.0 ? 1.0/MagickEpsilon : -1.0/MagickEpsilon;
          }
        else
          {
            slope.q=dy.q/dx.q;
            inverse_slope.q=(-1.0/slope.q);
          }
      offset.x=sqrt((double) (mid*mid/(inverse_slope.q*inverse_slope.q+1.0)));
      offset.y=(double) (offset.x*inverse_slope.q);
      dot_product=dy.q*offset.x-dx.q*offset.y;
      if (dot_product > 0.0)
        {
          box_p[2].x=polygon_primitive[n].point.x-offset.x;
          box_p[2].y=polygon_primitive[n].point.y-offset.y;
          box_p[3].x=polygon_primitive[i].point.x-offset.x;
          box_p[3].y=polygon_primitive[i].point.y-offset.y;
          box_q[2].x=polygon_primitive[n].point.x+offset.x;
          box_q[2].y=polygon_primitive[n].point.y+offset.y;
          box_q[3].x=polygon_primitive[i].point.x+offset.x;
          box_q[3].y=polygon_primitive[i].point.y+offset.y;
        }
      else
        {
          box_p[2].x=polygon_primitive[n].point.x+offset.x;
          box_p[2].y=polygon_primitive[n].point.y+offset.y;
          box_p[3].x=polygon_primitive[i].point.x+offset.x;
          box_p[3].y=polygon_primitive[i].point.y+offset.y;
          box_q[2].x=polygon_primitive[n].point.x-offset.x;
          box_q[2].y=polygon_primitive[n].point.y-offset.y;
          box_q[3].x=polygon_primitive[i].point.x-offset.x;
          box_q[3].y=polygon_primitive[i].point.y-offset.y;
        }
      if (fabs((double) (slope.p-slope.q)) < MagickEpsilon)
        {
          box_p[4]=box_p[1];
          box_q[4]=box_q[1];
        }
      else
        {
          box_p[4].x=(double) ((slope.p*box_p[0].x-box_p[0].y-slope.q*box_p[3].x+
            box_p[3].y)/(slope.p-slope.q));
          box_p[4].y=(double) (slope.p*(box_p[4].x-box_p[0].x)+box_p[0].y);
          box_q[4].x=(double) ((slope.p*box_q[0].x-box_q[0].y-slope.q*box_q[3].x+
            box_q[3].y)/(slope.p-slope.q));
          box_q[4].y=(double) (slope.p*(box_q[4].x-box_q[0].x)+box_q[0].y);
        }
      CheckPathExtent(MaxStrokePad,MaxStrokePad);
      dot_product=dx.q*dy.p-dx.p*dy.q;
      if (dot_product <= 0.0)
        switch (draw_info->linejoin)
        {
          case BevelJoin:
          {
            stroke_q[q++]=box_q[1];
            stroke_q[q++]=box_q[2];
            dot_product=(box_q[4].x-box_p[4].x)*(box_q[4].x-box_p[4].x)+
              (box_q[4].y-box_p[4].y)*(box_q[4].y-box_p[4].y);
            if (dot_product <= miterlimit)
              stroke_p[p++]=box_p[4];
            else
              {
                stroke_p[p++]=box_p[1];
                stroke_p[p++]=box_p[2];
              }
            break;
          }
          case MiterJoin:
          {
            dot_product=(box_q[4].x-box_p[4].x)*(box_q[4].x-box_p[4].x)+
              (box_q[4].y-box_p[4].y)*(box_q[4].y-box_p[4].y);
            if (dot_product <= miterlimit)
              {
                stroke_q[q++]=box_q[4];
                stroke_p[p++]=box_p[4];
              }
            else
              {
                stroke_q[q++]=box_q[1];
                stroke_q[q++]=box_q[2];
                stroke_p[p++]=box_p[1];
                stroke_p[p++]=box_p[2];
              }
            break;
          }
          case RoundJoin:
          {
            dot_product=(box_q[4].x-box_p[4].x)*(box_q[4].x-box_p[4].x)+
              (box_q[4].y-box_p[4].y)*(box_q[4].y-box_p[4].y);
            if (dot_product <= miterlimit)
              stroke_p[p++]=box_p[4];
            else
              {
                stroke_p[p++]=box_p[1];
                stroke_p[p++]=box_p[2];
              }
            center=polygon_primitive[n].point;
            theta.p=atan2(box_q[1].y-center.y,box_q[1].x-center.x);
            theta.q=atan2(box_q[2].y-center.y,box_q[2].x-center.x);
            if (theta.q < theta.p)
              theta.q+=2.0*MagickPI;
            arc_segments=(size_t) ceil((double) ((theta.q-theta.p)/
              (2.0*sqrt((double) (1.0/mid)))));
            CheckPathExtent(MaxStrokePad,arc_segments+MaxStrokePad);
            stroke_q[q].x=box_q[1].x;
            stroke_q[q].y=box_q[1].y;
            q++;
            for (j=1; j < (ssize_t) arc_segments; j++)
            {
              delta_theta=(double) (j*(theta.q-theta.p)/arc_segments);
              stroke_q[q].x=(double) (center.x+mid*cos(fmod((double)
                (theta.p+delta_theta),DegreesToRadians(360.0))));
              stroke_q[q].y=(double) (center.y+mid*sin(fmod((double)
                (theta.p+delta_theta),DegreesToRadians(360.0))));
              q++;
            }
            stroke_q[q++]=box_q[2];
            break;
          }
          default:
            break;
        }
      else
        switch (draw_info->linejoin)
        {
          case BevelJoin:
          {
            stroke_p[p++]=box_p[1];
            stroke_p[p++]=box_p[2];
            dot_product=(box_q[4].x-box_p[4].x)*(box_q[4].x-box_p[4].x)+
              (box_q[4].y-box_p[4].y)*(box_q[4].y-box_p[4].y);
            if (dot_product <= miterlimit)
              stroke_q[q++]=box_q[4];
            else
              {
                stroke_q[q++]=box_q[1];
                stroke_q[q++]=box_q[2];
              }
            break;
          }
          case MiterJoin:
          {
            dot_product=(box_q[4].x-box_p[4].x)*(box_q[4].x-box_p[4].x)+
              (box_q[4].y-box_p[4].y)*(box_q[4].y-box_p[4].y);
            if (dot_product <= miterlimit)
              {
                stroke_q[q++]=box_q[4];
                stroke_p[p++]=box_p[4];
              }
            else
              {
                stroke_q[q++]=box_q[1];
                stroke_q[q++]=box_q[2];
                stroke_p[p++]=box_p[1];
                stroke_p[p++]=box_p[2];
              }
            break;
          }
          case RoundJoin:
          {
            dot_product=(box_q[4].x-box_p[4].x)*(box_q[4].x-box_p[4].x)+
              (box_q[4].y-box_p[4].y)*(box_q[4].y-box_p[4].y);
            if (dot_product <= miterlimit)
              stroke_q[q++]=box_q[4];
            else
              {
                stroke_q[q++]=box_q[1];
                stroke_q[q++]=box_q[2];
              }
            center=polygon_primitive[n].point;
            theta.p=atan2(box_p[1].y-center.y,box_p[1].x-center.x);
            theta.q=atan2(box_p[2].y-center.y,box_p[2].x-center.x);
            if (theta.p < theta.q)
              theta.p+=2.0*MagickPI;
            arc_segments=(size_t) ceil((double) ((theta.p-theta.q)/
              (2.0*sqrt((double) (1.0/mid)))));
            CheckPathExtent(arc_segments+MaxStrokePad,MaxStrokePad);
            stroke_p[p++]=box_p[1];
            for (j=1; j < (ssize_t) arc_segments; j++)
            {
              delta_theta=(double) (j*(theta.q-theta.p)/arc_segments);
              stroke_p[p].x=(double) (center.x+mid*cos(fmod((double)
                (theta.p+delta_theta),DegreesToRadians(360.0))));
              stroke_p[p].y=(double) (center.y+mid*sin(fmod((double)
                (theta.p+delta_theta),DegreesToRadians(360.0))));
              p++;
            }
            stroke_p[p++]=box_p[2];
            break;
          }
          default:
            break;
        }
      slope.p=slope.q;
      inverse_slope.p=inverse_slope.q;
      box_p[0]=box_p[2];
      box_p[1]=box_p[3];
      box_q[0]=box_q[2];
      box_q[1]=box_q[3];
      dx.p=dx.q;
      dy.p=dy.q;
      n=i;
    }
    stroke_p[p++]=box_p[1];
    stroke_q[q++]=box_q[1];
    /*
      Trace stroked polygon.
    */
    stroke_polygon=(PrimitiveInfo *) AcquireQuantumMemory((size_t)
      (p+q+2UL*closed_path+2UL),sizeof(*stroke_polygon));
    if (stroke_polygon != (PrimitiveInfo *) NULL)
      {
        for (i=0; i < (ssize_t) p; i++)
        {
          stroke_polygon[i]=polygon_primitive[0];
          stroke_polygon[i].point=stroke_p[i];
        }
        if (closed_path != MagickFalse)
          {
            stroke_polygon[i]=polygon_primitive[0];
            stroke_polygon[i].point=stroke_polygon[0].point;
            i++;
          }
        for ( ; i < (ssize_t) (p+q+closed_path); i++)
        {
          stroke_polygon[i]=polygon_primitive[0];
          stroke_polygon[i].point=stroke_q[p+q+closed_path-(i+1)];
        }
        if (closed_path != MagickFalse)
          {
            stroke_polygon[i]=polygon_primitive[0];
            stroke_polygon[i].point=stroke_polygon[p+closed_path].point;
            i++;
          }
        stroke_polygon[i]=polygon_primitive[0];
        stroke_polygon[i].point=stroke_polygon[0].point;
        i++;
        stroke_polygon[i].primitive=UndefinedPrimitive;
        stroke_polygon[0].coordinates=(size_t) (p+q+2*closed_path+1);
      }
    stroke_p=(PointInfo *) RelinquishMagickMemory(stroke_p);
    stroke_q=(PointInfo *) RelinquishMagickMemory(stroke_q);
    polygon_primitive=(PrimitiveInfo *) RelinquishMagickMemory(polygon_primitive);
    return(stroke_polygon);
  }
  