ruby-gr 0.0.19 → 0.0.24

data/lib/gr/ffi.rb

@@ -174,7 +174,7 @@ module GR
     try_extern 'void gr_shadepoints(int, double *, double *, int, int, int)'
     try_extern 'void gr_shadelines(int, double *, double *, int, int, int)'
     try_extern 'void gr_panzoom(double, double, double, double, double *, double *, double *, double *)'
-    # try_extern 'int gr_findboundary(int, double *, double *, double, double (*)(double, double), int, int *)'
+    try_extern 'int gr_findboundary(int, double *, double *, double, double (*)(double, double), int, int *)'
     try_extern 'void gr_setresamplemethod(unsigned int flag)'
     try_extern 'void gr_inqresamplemethod(unsigned int *flag)'
     try_extern 'void gr_path(int, double *, double *, const char *)'
@@ -202,5 +202,6 @@ module GR
     try_extern 'void gr_inqtext3d(double, double, double, char *, int axis, double *, double *)'
     try_extern 'void gr_settextencoding(int)'
     try_extern 'void gr_inqtextencoding(int *)'
+    # gr_setcallback(char *(*)(const char *));
   end
 end

data/lib/gr/plot.rb

@@ -38,59 +38,95 @@ module GR
 
     # Keyword options conform to GR.jl.
     KW_ARGS = %i[accelerate algorithm alpha ax backgroundcolor barwidth baseline
-                 clabels clear clim color colormap crange figsize horizontal
-                 isovalue kind label labels levels location nbins ratio rotation
-                 scale size spec subplot tilt title update xaxis xflip xform
-                 xlabel xlim xlog xrange yaxis yflip ylabel ylim ylog zflip
-                 yrange viewport vp where window zaxis zlabel zlim zlog zrange].freeze
+                 clabels clear clim color colormap crange figsize font grid
+                 horizontal isovalue kind label labels levels location nbins
+                 ratio rotation scale size spec subplot tilt title update xaxis
+                 xflip xform xlabel xlim xlog xrange xticks yaxis yflip ylabel
+                 ylim ylog zflip yrange yticks viewport vp where window zaxis
+                 zlabel zlim zlog zrange zticks].freeze
+
+    FONTS = {
+      times_roman: 101,
+      times_italic: 102,
+      times_bold: 103,
+      times_bolditalic: 104,
+      helvetica_regular: 105,
+      helvetica_oblique: 106,
+      helvetica_bold: 107,
+      helvetica_boldoblique: 108,
+      courier_regular: 109,
+      courier_oblique: 110,
+      courier_bold: 111,
+      courier_boldoblique: 112,
+      symbol: 113,
+      bookman_light: 114,
+      bookman_lightitalic: 115,
+      bookman_demi: 116,
+      bookman_demiitalic: 117,
+      newcenturyschlbk_roman: 118,
+      newcenturyschlbk_italic: 119,
+      newcenturyschlbk_bold: 120,
+      newcenturyschlbk_bolditalic: 121,
+      avantgarde_book: 122,
+      avantgarde_bookoblique: 123,
+      avantgarde_demi: 124,
+      avantgarde_demioblique: 125,
+      palatino_roman: 126,
+      palatino_italic: 127,
+      palatino_bold: 128,
+      palatino_bolditalic: 129,
+      zapfchancery_mediumitalic: 130,
+      zapfdingbats: 131,
+      cmuserif_math: 232, # original: cmuserif-math
+      dejavusans: 233,
+      pingfangsc: 234
+    }.freeze
 
     @last_plot = nil
     class << self
       attr_accessor :last_plot
     end
 
-    def initialize(*args)
-      @kvs = if args[-1].is_a? Hash
-               args.pop
-             else
-               {}
-             end
+    attr_accessor :args, :kvs, :scheme
+
+    def initialize(*raw_args)
+      @kvs = raw_args.last.is_a?(Hash) ? raw_args.pop : {}
+      @args = plot_args(raw_args) # method name is the same as Julia/Python
+
       # Check keyword options.
-      @kvs.each_key do |k|
-        warn "Unknown keyword: #{k}" unless KW_ARGS.include? k
-      end
+      kvs.each_key { |k| warn "Unknown keyword: #{k}" unless KW_ARGS.include? k }
 
       # label(singular form) is a original keyword arg which GR.jl does not have.
-      @kvs[:labels] = [@kvs[:label]] if @kvs[:label] && @kvs[:labels].nil?
-
-      @args = plot_args(args) # method name is the same as Julia/Python
-      @kvs[:size] ||= [600, 450]
-      @kvs[:ax] = false if @kvs[:ax].nil?
-      @kvs[:subplot] ||= [0, 1, 0, 1]
-      @kvs[:clear] = true if @kvs[:clear].nil?
-      @kvs[:update] = true if @kvs[:update].nil?
-      @scheme = 0
+      kvs[:labels] ||= [kvs[:label]] if kvs.has_key? :label
+
+      # Don't use ||= here, because we need to tell `false` from `nil`
+      kvs[:size]    = [600, 450]   unless kvs.has_key? :size
+      kvs[:ax]      = false        unless kvs.has_key? :ax
+      kvs[:subplot] = [0, 1, 0, 1] unless kvs.has_key? :subplot
+      kvs[:clear]   = true         unless kvs.has_key? :clear
+      kvs[:update]  = true         unless kvs.has_key? :update
+
+      @scheme     = 0
       @background = 0xffffff
-      @handle = nil
+      # @handle     = nil           # This variable will be used in gr_meta
+
       self.class.last_plot = self
     end
-    attr_accessor :args, :kvs, :scheme
 
     def set_viewport(kind, subplot)
       mwidth, mheight, width, height = GR.inqdspsize
-      if kvs[:figsize]
-        w = 0.0254 * width * kvs[:figsize][0] / mwidth
-        h = 0.0254 * height * kvs[:figsize][1] / mheight
-      else
-        dpi = width / mwidth * 0.0254
-        if dpi > 200
-          w, h = kvs[:size].map { |x| x * dpi / 100 }
-        else
-          w, h = kvs[:size]
-        end
-      end
-      viewport = [0, 0, 0, 0]
+      dpi = width / mwidth * 0.0254
+      w, h = if kvs[:figsize]
+               [(0.0254 * width  * kvs[:figsize][0] / mwidth),
+                (0.0254 * height * kvs[:figsize][1] / mheight)]
+             elsif dpi > 200
+               kvs[:size].map { |i| i * dpi / 100 }
+             else
+               kvs[:size]
+             end
+
       vp = subplot.clone
+
       if w > h
         ratio = h / w.to_f
         msize = mwidth * w / width
@@ -106,6 +142,7 @@ module GR
         vp[0] *= ratio
         vp[1] *= ratio
       end
+
       if %i[wireframe surface plot3 scatter3 trisurf volume].include?(kind)
         extent = [vp[1] - vp[0], vp[3] - vp[2]].min
         vp1 = 0.5 * (vp[0] + vp[1] - extent)
@@ -115,10 +152,12 @@ module GR
       else
         vp1, vp2, vp3, vp4 = vp
       end
-      viewport[0] = vp1 + 0.125 * (vp2 - vp1)
-      viewport[1] = vp1 + 0.925 * (vp2 - vp1)
-      viewport[2] = vp3 + 0.125 * (vp4 - vp3)
-      viewport[3] = vp3 + 0.925 * (vp4 - vp3)
+
+      viewport = [vp1 + 0.125 * (vp2 - vp1),
+                  vp1 + 0.925 * (vp2 - vp1),
+                  vp3 + 0.125 * (vp4 - vp3),
+                  vp3 + 0.925 * (vp4 - vp3)]
+
       if %i[contour contourf hexbin heatmap nonuniformheatmap polarheatmap
             surface trisurf volume].include?(kind)
         viewport[1] -= 0.1
@@ -132,11 +171,11 @@ module GR
         end
       end
 
-      GR.setviewport(viewport[0], viewport[1], viewport[2], viewport[3])
+      GR.setviewport(*viewport)
 
       kvs[:viewport] = viewport
-      kvs[:vp] = vp
-      kvs[:ratio] = ratio
+      kvs[:vp]       = vp
+      kvs[:ratio]    = ratio
 
       if kvs[:backgroundcolor]
         GR.savestate
@@ -166,14 +205,15 @@ module GR
     def set_window(kind)
       scale = 0
       unless %i[polar polarhist polarheatmap].include?(kind)
-        scale |= GR::OPTION_X_LOG if kvs[:xlog]
-        scale |= GR::OPTION_Y_LOG if kvs[:ylog]
-        scale |= GR::OPTION_Z_LOG if kvs[:zlog]
+        scale |= GR::OPTION_X_LOG  if kvs[:xlog]
+        scale |= GR::OPTION_Y_LOG  if kvs[:ylog]
+        scale |= GR::OPTION_Z_LOG  if kvs[:zlog]
         scale |= GR::OPTION_FLIP_X if kvs[:xflip]
         scale |= GR::OPTION_FLIP_Y if kvs[:yflip]
         scale |= GR::OPTION_FLIP_Z if kvs[:zflip]
       end
       kvs[:scale] = scale
+
       if kvs.has_key?(:panzoom)
         xmin, xmax, ymin, ymax = GR.panzoom(*kvs[:panzoom])
         kvs[:xrange] = [xmin, xmax]
@@ -182,68 +222,81 @@ module GR
         minmax
       end
 
-      major_count = %i[wireframe surface plot3 scatter3 polar polarhist
-                       polarheatmap trisurf volume].include?(kind) ? 2 : 5
+      major_count = if %i[wireframe surface plot3 scatter3 polar polarhist
+                          polarheatmap trisurf volume].include?(kind)
+                      2
+                    else
+                      5
+                    end
+
+      kvs[:xticks] = [kvs[:xticks], major_count] if kvs[:xticks].is_a? Numeric
+      kvs[:yticks] = [kvs[:yticks], major_count] if kvs[:yticks].is_a? Numeric
+      kvs[:zticks] = [kvs[:zticks], major_count] if kvs[:zticks].is_a? Numeric
 
       xmin, xmax = kvs[:xrange]
-      if (scale & GR::OPTION_X_LOG) == 0
-        xmin, xmax = GR.adjustlimits(xmin, xmax) unless kvs.has_key?(:xlim) || kvs.has_key?(:panzoom)
-        if kvs.has_key?(:xticks)
-          xtick, majorx = kvs[:xticks]
-        else
-          majorx = major_count
-          xtick = GR.tick(xmin, xmax) / majorx
-        end
-      else
-        xtick = majorx = 1
+      if %i[heatmap polarheatmap].include?(kind) && kvs.has_key?(:xlim)
+        xmin -= 0.5
+        xmax += 0.5
       end
-      xorg = if (scale & GR::OPTION_FLIP_X) == 0
-               [xmin, xmax]
-             else
-               [xmax, xmin]
-             end
+      xtick, majorx = if (scale & GR::OPTION_X_LOG) == 0
+                        unless %i[heatmap polarheatmap].include?(kind)
+                          unless kvs.has_key?(:xlim)
+                            xmin, xmax = GR.adjustlimits(xmin, xmax) unless kvs[:panzoom]
+                          end
+                        end
+                        if kvs.has_key?(:xticks)
+                          kvs[:xticks]
+                        else
+                          [GR.tick(xmin, xmax) / major_count, major_count]
+                        end
+                      else
+                        [1, 1]
+                      end
+      xorg = (scale & GR::OPTION_FLIP_X) == 0 ? [xmin, xmax] : [xmax, xmin]
       kvs[:xaxis] = xtick, xorg, majorx
 
       ymin, ymax = kvs[:yrange]
-      if kind == :hist && !kvs.has_key?(:ylim)
-        ymin = (scale & GR::OPTION_Y_LOG) == 0 ? 0 : 1
+      if %i[heatmap polarheatmap].include?(kind) && kvs.has_key?(:ylim)
+        ymin -= 0.5
+        ymax += 0.5
       end
-      if (scale & GR::OPTION_Y_LOG) == 0
-        ymin, ymax = GR.adjustlimits(ymin, ymax) unless kvs.has_key?(:ylim) || kvs.has_key?(:panzoom)
-        if kvs.has_key?(:yticks)
-          ytick, majory = kvs[:yticks]
-        else
-          majory = major_count
-          ytick = GR.tick(ymin, ymax) / majory
+      if kind == :hist
+        if kvs[:horizontal] && !kvs.has_key?(:xlim)
+          xmin = (scale & GR::OPTION_X_LOG) == 0 ? 0 : 1
+        elsif !kvs.has_key?(:ylim)
+          ymin = (scale & GR::OPTION_Y_LOG) == 0 ? 0 : 1
         end
-      else
-        ytick = majory = 1
       end
-      yorg = if (scale & GR::OPTION_FLIP_Y) == 0
-               [ymin, ymax]
-             else
-               [ymax, ymin]
-             end
+      ytick, majory = if (scale & GR::OPTION_Y_LOG) == 0
+                        unless %i[heatmap polarheatmap].include?(kind)
+                          unless kvs.has_key?(:ylim)
+                            ymin, ymax = GR.adjustlimits(ymin, ymax) unless kvs[:panzoom]
+                          end
+                        end
+                        if kvs.has_key?(:yticks)
+                          kvs[:yticks]
+                        else
+                          [GR.tick(ymin, ymax) / major_count, major_count]
+                        end
+                      else
+                        [1, 1]
+                      end
+      yorg = (scale & GR::OPTION_FLIP_Y) == 0 ? [ymin, ymax] : [ymax, ymin]
       kvs[:yaxis] = ytick, yorg, majory
 
       if %i[wireframe surface plot3 scatter3 trisurf volume].include?(kind)
         zmin, zmax = kvs[:zrange]
-        if (scale & GR::OPTION_Z_LOG) == 0
-          zmin, zmax = GR.adjustlimits(zmin, zmax) if kvs.has_key?(:zlim)
-          if kvs.has_key?(:zticks)
-            ztick, majorz = kvs[:zticks]
-          else
-            majorz = major_count
-            ztick = GR.tick(zmin, zmax) / majorz
-          end
-        else
-          ztick = majorz = 1
-        end
-        zorg = if (scale & GR::OPTION_FLIP_Z) == 0
-                 [zmin, zmax]
-               else
-                 [zmax, zmin]
-               end
+        ztick, majorz = if (scale & GR::OPTION_Z_LOG) == 0
+                          zmin, zmax = GR.adjustlimits(zmin, zmax) if kvs.has_key?(:zlim)
+                          if kvs.has_key?(:zticks)
+                            kvs[:zticks]
+                          else
+                            [GR.tick(zmin, zmax) / major_count, major_count]
+                          end
+                        else
+                          [1, 1]
+                        end
+        zorg = (scale & GR::OPTION_FLIP_Z) == 0 ? [zmin, zmax] : [zmax, zmin]
         kvs[:zaxis] = ztick, zorg, majorz
       end
 
@@ -255,7 +308,7 @@ module GR
       end
       if %i[wireframe surface plot3 scatter3 trisurf volume].include?(kind)
         rotation = kvs[:rotation] || 40
-        tilt = kvs[:tilt] || 70
+        tilt     = kvs[:tilt]     || 70
         GR.setspace(zmin, zmax, rotation, tilt)
       end
 
@@ -329,29 +382,29 @@ module GR
         GR.axes(xtick, ytick, xorg[1], yorg[1], -majorx, -majory, -ticksize)
       end
 
-      if kvs[:title]
+      if kvs.has_key?(:title)
         GR.savestate
         GR.settextalign(GR::TEXT_HALIGN_CENTER, GR::TEXT_VALIGN_TOP)
-        text(0.5 * (viewport[0] + viewport[1]), vp[3], kvs[:title])
+        text(0.5 * (viewport[0] + viewport[1]), vp[3], kvs[:title].to_s)
         GR.restorestate
       end
       if %i[wireframe surface plot3 scatter3 trisurf volume].include?(kind)
-        xlabel = kvs[:xlabel] || ''
-        ylabel = kvs[:ylabel] || ''
-        zlabel = kvs[:zlabel] || ''
+        xlabel = (kvs[:xlabel] || '').to_s
+        ylabel = (kvs[:ylabel] || '').to_s
+        zlabel = (kvs[:zlabel] || '').to_s
         GR.titles3d(xlabel, ylabel, zlabel)
       else
         if kvs.has_key?(:xlabel)
           GR.savestate
           GR.settextalign(GR::TEXT_HALIGN_CENTER, GR::TEXT_VALIGN_BOTTOM)
-          text(0.5 * (viewport[0] + viewport[1]), vp[2] + 0.5 * charheight, kvs[:xlabel])
+          text(0.5 * (viewport[0] + viewport[1]), vp[2] + 0.5 * charheight, kvs[:xlabel].to_s)
           GR.restorestate
         end
         if kvs.has_key?(:ylabel)
           GR.savestate
           GR.settextalign(GR::TEXT_HALIGN_CENTER, GR::TEXT_VALIGN_TOP)
           GR.setcharup(-1, 0)
-          text(vp[0] + 0.5 * charheight, 0.5 * (viewport[2] + viewport[3]), kvs[:ylabel])
+          text(vp[0] + 0.5 * charheight, 0.5 * (viewport[2] + viewport[3]), kvs[:ylabel].to_s)
           GR.restorestate
         end
       end
@@ -459,7 +512,7 @@ module GR
       if kvs.has_key?(:title)
         GR.savestate
         GR.settextalign(GR::TEXT_HALIGN_CENTER, GR::TEXT_VALIGN_TOP)
-        text(0.5 * (viewport[0] + viewport[1]), vp[3], kvs[:title])
+        text(0.5 * (viewport[0] + viewport[1]), vp[3], kvs[:title].to_s)
         GR.restorestate
       end
       GR.selntran(1)
@@ -561,7 +614,22 @@ module GR
         # Not yet.
       end
 
-      GR.settextfontprec(232, 3)
+      if kvs.has_key?(:font)
+        name = kvs[:font]
+        # 'Cmuserif-Math' => :cmuserif_math
+        sym_name = name.to_s.gsub('-','_').downcase.to_sym 
+        if FONTS.include?(sym_name)
+          font = FONTS[sym_name]
+          GR.settextfontprec(font, font > 200 ? 3 : 0)
+        else
+          warn "Unknown font name: #{name}" # should raise error?
+        end
+      else
+        # The following fonts are the default in GR.jl
+        # Japanese, Chinese, Korean, etc. cannot be displayed.
+
+        # GR.settextfontprec(232, 3) # CM Serif Roman
+      end
 
       set_viewport(kind, kvs[:subplot])
       unless kvs[:ax]
@@ -629,8 +697,7 @@ module GR
           if z || c
             if c
               cmin, cmax = kvs[:crange]
-              c = c.to_a if narray?(c)
-              c.map! { |i| normalize_color(i, cmin, cmax) }
+              c = c.map { |i| normalize_color(i, cmin, cmax) }
               cind = c.map { |i| (1000 + i * 255).round }
             end
             x.length.times do |i|
@@ -655,14 +722,26 @@ module GR
           GR.polymarker(x, y)
 
         when :hist
-          ymin = kvs[:window][2]
-          y.length.times do |i|
-            GR.setfillcolorind(989)
-            GR.setfillintstyle(GR::INTSTYLE_SOLID)
-            GR.fillrect(x[i], x[i + 1], ymin, y[i])
-            GR.setfillcolorind(1)
-            GR.setfillintstyle(GR::INTSTYLE_HOLLOW)
-            GR.fillrect(x[i], x[i + 1], ymin, y[i])
+          if kvs[:horizontal]
+            xmin = kvs[:window][0]
+            x.length.times do |i|
+              GR.setfillcolorind(989)
+              GR.setfillintstyle(GR::INTSTYLE_SOLID)
+              GR.fillrect(xmin, x[i], y[i], y[i + 1])
+              GR.setfillcolorind(1)
+              GR.setfillintstyle(GR::INTSTYLE_HOLLOW)
+              GR.fillrect(xmin, x[i], y[i], y[i + 1])
+            end
+          else
+            ymin = kvs[:window][2]
+            y.length.times do |i|
+              GR.setfillcolorind(989)
+              GR.setfillintstyle(GR::INTSTYLE_SOLID)
+              GR.fillrect(x[i], x[i + 1], ymin, y[i])
+              GR.setfillcolorind(1)
+              GR.setfillintstyle(GR::INTSTYLE_HOLLOW)
+              GR.fillrect(x[i], x[i + 1], ymin, y[i])
+            end
           end
 
         when :polarhist
@@ -719,9 +798,10 @@ module GR
           else
             h = levels
           end
-          if kind == :contour
+          case kind
+          when :contour
             GR._contour_(x, y, h, z, clabels ? 1 : 1000)
-          elsif kind == :contourf
+          when :contourf
             GR._contourf_(x, y, h, z, clabels ? 1 : 0)
           end
           colorbar(0, h.length)
@@ -810,7 +890,7 @@ module GR
           GR.setmarkertype(GR::MARKERTYPE_SOLID_CIRCLE)
           if c
             cmin, cmax = kvs[:crange]
-            c = c.map { |i| normalize_color(i, cmin, cmax) } # NArray -> Array
+            c = c.map { |i| normalize_color(i, cmin, cmax) }
             cind = c.map { |i| (1000 + i * 255).round }
             x.length.times do |i|
               GR.setmarkercolorind(cind[i])
@@ -1016,13 +1096,16 @@ module GR
 
     # Normalize a color c with the range [cmin, cmax]
     #   0 <= normalize_color(c, cmin, cmax) <= 1
+    # Note: narray.map{|i| normalize_color(i)} There's room for speedup.
     def normalize_color(c, cmin, cmax)
+      c = c.to_f # if c is Integer
       c = c.clamp(cmin, cmax) - cmin
       c /= (cmax - cmin) if cmin != cmax
       c
     end
 
     def inqtext(x, y, s)
+      s = s.to_s
       if s.length >= 2 && s[0] == '$' && s[-1] == '$'
         GR.inqmathtex(x, y, s[1..-2])
       elsif s.include?('\\') || s.include?('_') || s.include?('^')
@@ -1033,6 +1116,7 @@ module GR
     end
 
     def text(x, y, s)
+      s = s.to_s
       if s.length >= 2 && s[0] == '$' && s[-1] == '$'
         GR.mathtex(x, y, s[1..-2])
       elsif s.include?('\\') || s.include?('_') || s.include?('^')
@@ -1057,7 +1141,8 @@ module GR
       args.each do |x, y, z, c|
         if x
           if scale & GR::OPTION_X_LOG != 0
-            x.map! { |v| v > 0 ? v : Float::NAN }
+            # duck typing for NArray
+            x = x.map { |v| v > 0 ? v : Float::NAN }
           end
           x0, x1 = x.minmax
           xmin = [x0, xmin].min
@@ -1068,7 +1153,7 @@ module GR
         end
         if y
           if scale & GR::OPTION_Y_LOG != 0
-            y.map! { |v| v > 0 ? v : Float::NAN }
+            y = y.map { |v| v > 0 ? v : Float::NAN }
           end
           y0, y1 = y.minmax
           ymin = [y0, ymin].min
@@ -1079,7 +1164,7 @@ module GR
         end
         if z
           if scale & GR::OPTION_Z_LOG != 0
-            z.map! { |v| v > 0 ? v : Float::NAN }
+            z = z.map { |v| v > 0 ? v : Float::NAN }
           end
           z0, z1 = z.minmax
           zmin = [z0, zmin].min
@@ -1131,20 +1216,12 @@ module GR
       [w, h]
     end
 
-    # NOTE: duplicated definition (GRCommonUtils)
     def equal_length(*args)
-      lengths = args.map(&:length)
-      unless lengths.all? { |l| l == lengths[0] }
-        raise ArgumentError,
-              'Sequences must have same length.'
-      end
-
-      lengths[0]
+      GRCommons::GRCommonUtils.equal_length(*args)
     end
 
-    # NOTE: duplicated definition (GRCommonUtils)
     def narray?(data)
-      defined?(Numo::NArray) && data.is_a?(Numo::NArray)
+      GRCommons::GRCommonUtils.narray?(data)
     end
   end
 
@@ -1293,11 +1370,15 @@ module GR
     end
 
     # (Plot) Draw a histogram.
-    def histogram(x, kv = {})
-      create_plot(:hist, x, kv) do |plt|
+    def histogram(series, kv = {})
+      create_plot(:hist, series, kv) do |plt|
         nbins = plt.kvs[:nbins] || 0
-        x, y = hist(x, nbins)
-        plt.args = [[x, y, nil, nil, '']]
+        x, y = hist(series, nbins)
+        plt.args = if kv[:horizontal]
+                     [[y, x, nil, nil, '']]
+                   else
+                     [[x, y, nil, nil, '']]
+                   end
       end
     end
 
@@ -1323,7 +1404,7 @@ module GR
     end
 
     # Set current subplot index.
-    def subplot(nr, nc, p)
+    def subplot(nr, nc, p, kv = {})
       xmin = 1
       xmax = 0
       ymin = 1
@@ -1342,7 +1423,7 @@ module GR
         # The policy of clearing when p[0]==1 is controversial
         clear: p[0] == 1,
         update: p[-1] == nr * nc
-      }
+      }.merge kv
     end
 
     # (Plot) Save the current figure to a file.