Add a bunch of failing tests for Unitful plotting

test/REQUIRE

@@ -1,3 +1,4 @@
 RDatasets
 Cairo
 Fontconfig
+Unitful

test/testscripts/unitful_basic.jl

@@ -0,0 +1,19 @@
+using Unitful, Gadfly, DataFrames
+
+a = -9.81u"m/s^2"
+t = (1:0.5:10)u"s"
+v = a * t
+v2 = v + 2u"m/s"
+h = 0.5 * a * t.^2
+df = DataFrame(time=t, velocity=v, position=h,
+               unitlesst=ustrip.(t), unitlessv=ustrip.(v), unitlessh=ustrip.(h),
+               position2=reverse(h))
+
+# test basics of point/line plots with Unitful
+gridstack(reshape([plot(df, x=:time, y=:velocity),
+                   plot(df, x=:position, y=:time),
+                   plot(df, x=:unitlesst, y=:velocity),
+                   plot(df, x=:time, y=:unitlessh),
+                   plot(df, x=:time, y=:position, Geom.line),
+                   plot(df, layer(x=:time, y=:position, Geom.line),
+                        layer(x=:time, y=:position2))], (3,2)))

test/testscripts/unitful_color.jl

@@ -0,0 +1,15 @@
+using Unitful, Gadfly, DataFrames
+
+a = -9.81u"m/s^2"
+t = (1:0.5:10)u"s"
+v = a * t
+v2 = v + 2u"m/s"
+h = 0.5 * a * t.^2
+df = DataFrame(time=t, velocity=v, position=h,
+               unitlesst=ustrip.(t), unitlessv=ustrip.(v), unitlessh=ustrip.(h),
+               position2=reverse(h))
+
+# Test that it's possible to categorize by Unitful quantities
+vstack(plot(df, x=:time, y=:position, color=:velocity, Geom.point),
+       plot(df, x=:time, y=:velocity, color=:position, Geom.line),
+       plot(df, x=:time, y=:position, color=:unitlessv))

test/testscripts/unitful_geoms.jl

@@ -0,0 +1,26 @@
+using Unitful, Gadfly, DataFrames
+
+a = -9.81u"m/s^2"
+t = (1:0.5:10)u"s"
+v = a * t
+v2 = v + 2u"m/s"
+h = 0.5 * a * t.^2
+df = DataFrame(time=t, velocity=v, position=h,
+               unitlesst=ustrip.(t), unitlessv=ustrip.(v), unitlessh=ustrip.(h),
+               position2=reverse(h))
+
+# Test various geometries with Unitful
+p1 = plot(df, x=:time, y=:velocity, Geom.bar)
+p2 = plot(df, x=:time, y=:position, Geom.point,
+          intercept=[-80u"m"], slope=[10u"m/s"], Geom.abline)
+p3 = plot(df, x=:time, y=:velocity, Geom.line,
+          yintercept=[-20u"m/s", -40u"m/s"], Geom.hline)
+# Currently explicitly stated that it loess and lm require arrays of plain numbers
+#p4 = plot(df, x=:time, y=:position, Geom.point,
+#          Geom.smooth(method=:loess, smoothing=0.2))
+p4 = plot(df, x=:position, y=:velocity, Geom.path)
+p5 = plot(df, x=:time, y=:position, Geom.step,
+          xintercept=[3u"s", 8u"s"], Geom.vline)
+p6 = plot(df, layer(x=:time, y=:position, Geom.line),
+          layer(x=:time, y=:position2, Geom.bar))
+gridstack(reshape([p1,p2,p3,p4,p5,p6], (2,3)))