Some minor changes that apparently workaround Julia bugs that forces recompilation.

src/common_options.jl

@@ -97,7 +97,7 @@ end
 function GMTsyntax_opt(d::Dict, cmd::String, del::Bool=true)::Tuple{String, String}
 	o::String = ""
 	if haskey(d, :compact)
-		o = d[:compact]
+		o = d[:compact]::String
 		(o == "") && (delete!(d, :compact); return cmd, "")
 		cmd = (o[1] == ' ') ? cmd * o : cmd * " " * o
 		del && delete!(d, :compact)
@@ -239,22 +239,16 @@ function merge_R_and_xyzlims(d::Dict, opt_R::String)::String
 	xlim, ylim, zlim = parse_lims(d, 'x'), parse_lims(d, 'y'), parse_lims(d, 'z')
 =#
 	function fetch_xyz_lims(d::Dict, symbs)::String
-		return ((val = find_in_dict(d, symbs, false)[1]) !== nothing) ? @sprintf("%.15g/%.15g", val[1], val[2]) : ""
+		((val = find_in_dict(d, symbs, false)[1]) === nothing) && return ""
+		f1::Float64, f2::Float64 = Float64(val[1]), Float64(val[2])		# Avoid sending Anys to printf. The pure horror acording to JET
+		return @sprintf("%.15g/%.15g", f1, f2)
 	end
 
 	xlim::String = fetch_xyz_lims(d, [:xlim :xlims :xlimits])
 	ylim::String = fetch_xyz_lims(d, [:ylim :ylims :ylimits])
 	zlim::String = fetch_xyz_lims(d, [:zlim :zlims :zlimits])
 
-#=
-	xlim::String = ((val = find_in_dict(d, [:xlim :xlims :xlimits], false)[1]) !== nothing) ?
-	                sprintf("%.15g/%.15g", val[1], val[2]) : ""
-	ylim::String = ((val = find_in_dict(d, [:ylim :ylims :ylimits], false)[1]) !== nothing) ?
-	                sprintf("%.15g/%.15g", val[1], val[2]) : ""
-	zlim::String = ((val = find_in_dict(d, [:zlim :zlims :zlimits], false)[1]) !== nothing) ?
-	                sprintf("%.15g/%.15g", val[1], val[2]) : ""
-=#
-	(xlim == "" && ylim == "" && zlim == "") && return opt_R
+	(xlim === "" && ylim === "" && zlim === "") && return opt_R
 
 	function clear_xyzlims(d::Dict, xlim, ylim, zlim)
 		# When calling this fun, if they exist they were used too so must remove them
@@ -1452,7 +1446,7 @@ function parse_UXY(cmd::String, d::Dict, aliases, opt::Char)::String
 	end
 	# If -X|Y assume it's a new fig so plot any legend that may be trailing around.
 	# This may screw but also screws if we don't do it. 
-	(symb in [:X :xshift :x_offset :Y :yshift :y_offset]) && digests_legend_bag(d)
+	(symb in [:X :xshift :x_offset :Y :yshift :y_offset]) && digests_legend_bag(d)	# FORCES RECOMPILE
 	return cmd
 end
 
@@ -2926,8 +2920,8 @@ function data_type(val)
 end
 
 # ---------------------------------------------------------------------------------------------------
-axis(nt::NamedTuple, D::Dict=Dict(); x::Bool=false, y::Bool=false, z::Bool=false, secondary::Bool=false) = axis(D; x=x, y=y, z=z, secondary=secondary, nt...)
-function axis(D::Dict=Dict(); x::Bool=false, y::Bool=false, z::Bool=false, secondary::Bool=false, kwargs...)::Tuple{String, Vector{Bool}}
+axis(nt::NamedTuple, D::Dict=Dict{Symbol,Any}(); x::Bool=false, y::Bool=false, z::Bool=false, secondary::Bool=false) = axis(D; x=x, y=y, z=z, secondary=secondary, nt...)
+function axis(D::Dict=Dict{Symbol,Any}(); x::Bool=false, y::Bool=false, z::Bool=false, secondary::Bool=false, kwargs...)::Tuple{String, Vector{Bool}}
 	d = KW(kwargs)			# These kwargs always come from the fields of a NamedTuple 
 	axis(D, x, y, z, secondary, d)
 end
@@ -3091,7 +3085,7 @@ function axis(D::Dict, x::Bool, y::Bool, z::Bool, secondary::Bool, d::Dict)::Tup
 		end
 	end
 	consume_used(d, [:corners, :internal, :cube, :noframe, :pole, :title, :subtitle, :seclabel, :label, :xlabel, :ylabel, :zlabel, :Yhlabel, :alabel, :blabel, :clabel, :annot, :annot_unit, :ticks, :ticks_unit, :prefix, :suffix, :grid, :slanted, :custom, :customticks, :pi, :scale, :phase_add, :phase_sub])
-	(length(d) > 0) && println("Warning: the following sub-options were not consumed in 'frame' => ", keys(d))
+	(length(d) > 0) && println("Warning: the following sub-options were not consumed in 'frame' => ", "$(keys(d))") # WTF ???. keys(d) alone FORCES RECOMPILATION OF THE ENTIRE plot() PROGRAM
 	# ----------------------------------------------------
 
 	return opt, [have_Baxes, (opt_Bframe != "")]

src/gmt_main.jl

@@ -229,13 +229,11 @@ function gmt(cmd::String, args...)
 		return out[1], out[2]
 	elseif (n_out == 3)
 		return out[1], out[2], out[3]
-	else
-		return out
 	end
 
 end
 
-gmt_GMTdataset(cmd::String, args...)::GMTdataset{Float64,2} = gmt(cmd, args...)
+gmt_GMTdataset(cmd::String, args...)::Union{GMTdataset{Float64,2}, Vector{<:GMTdataset{Float64,2}}} = gmt(cmd, args...)
 gmt_GMTgrid(cmd::String, args...)::GMTgrid{Float32,2} = gmt(cmd, args...)::GMTgrid{Float32,2}
 
 # -----------------------------------------------------------------------------------------------

src/makecpt.jl

@@ -60,7 +60,7 @@ function makecpt(cmd0::String="", arg1=nothing; kwargs...)::Union{String, GMTcpt
 end
 function makecpt(cmd0::String, arg1, d::Dict)::Union{String, GMTcpt}
 
-	cmd, = parse_common_opts(d, "", [:V_params])
+	cmd = parse_V_params(d, "")
 
 	# This deals with the special case of, for example, "makecpt(G, cmap=:gray"). Here, we recieve (makecpt(G,...) is comp)
 	# a (C=nothing, cmap=:gray) and because 'C' is searched first we would loose the cmap=:gray. Solution is to delete :C

src/psxy.jl

@@ -63,7 +63,7 @@ function _common_plot_xyz(cmd0::String, arg1, caller::String, O::Bool, K::Bool,
 		(isa(arg1, GMTdataset) && size(arg1,2) > 1 && !isempty(arg1.colnames)) && (CTRL.XYlabels[1] = arg1.colnames[1]; CTRL.XYlabels[2] = arg1.colnames[2])
 		isa(arg1, Vector{<:GMTdataset}) && !isempty(arg1[1].colnames) && (CTRL.XYlabels[1] = arg1[1].colnames[1]; CTRL.XYlabels[2] = arg1[1].colnames[2])
 		if (is_ternary)  cmd, opt_J = parse_J(d, cmd, default=def_J)
-		else             cmd, opt_B, opt_J, opt_R = parse_BJR(d, cmd, caller, O, def_J)	# FORCES RECOMPILE
+		else             cmd, opt_B, opt_J, opt_R = parse_BJR(d, cmd, caller, O, def_J)
 		end
 		# Current parse_B does not add a default -Baz when 3D AND -J has a projection. More or less fix that.
 		if (is3D && opt_B != "" && !contains(opt_B, " -Bz"))
@@ -194,7 +194,7 @@ function _common_plot_xyz(cmd0::String, arg1, caller::String, O::Bool, K::Bool,
 			else
 				mima = [extrema(last_non_nothing(arg1, arg2, arg3))...]	# Why 'last'?
 			end
-			r = makecpt(@sprintf("-T%f/%f/65+n -Cturbo -Vq", mima[1]-eps(1e10), mima[2]+eps(1e10)))
+			r = makecpt(@sprintf("-T%f/%f/65+n -Cturbo -Vq", mima[1]-eps(1e10), mima[2]+eps(1e10)))::GMTcpt
 		else
 			r = CURRENT_CPT[1]
 		end
@@ -215,14 +215,15 @@ function _common_plot_xyz(cmd0::String, arg1, caller::String, O::Bool, K::Bool,
 	(made_it_vector && opt_S == "") && (cmd *= " -Sv+s")	# No set opt_S because it results in 2 separate commands
 
 	# See if any of the scatter, bar, lines, etc... was the caller and if yes, set sensible defaults.
-	_cmd = set_avatar_defaults(d, cmd, mcc, caller, got_usr_R, opt_B, opt_Gsymb, opt_ML, opt_R, opt_S, opt_W, sub_module, g_bar_fill, opt_Wmarker, opt_UVXY, opt_c, O, arg1)
+	_cmd = set_avatar_defaults(d, cmd, mcc, caller, got_usr_R, opt_B, opt_Gsymb, opt_ML, opt_R, opt_S, opt_W,
+	                           sub_module, g_bar_fill, opt_Wmarker, opt_UVXY, opt_c, O, arg1)	# FORCES RECOMPILE
 
 	(got_Zvars && opt_S == "" && opt_W == "" && !occursin(" -G", _cmd[1])) && (_cmd[1] *= " -W0.5")
 	(opt_W == "" && caller == "feather") && (_cmd[1] *= " -W0.1")		# feathers are normally many so better they are thin
 
 	# Let matrices with more data columns, and for which Color info was NOT set, plot multiple lines at once
 	if (!mcc && sub_module == "bar" && check_bar_group(arg1))	# !mcc because the bar-groups all have mcc = false
-		_cmd[1], arg1, cmd2::String = bar_group(d, _cmd[1], opt_R, g_bar_fill, got_Ebars, got_usr_R, arg1)
+		_cmd[1], arg1, cmd2::String = bar_group(d, _cmd[1], opt_R, g_bar_fill, got_Ebars, got_usr_R, arg1)	# FORCES RECOMPILE
 		(cmd2 != "") && (length(_cmd) == 1 ? (_cmd = [cmd2; _cmd[1]]) :
 			(@warn("Can't plot the connector when 'bar' is already a nested call."); CTRL.pocket_call[3] = nothing))
 	end
@@ -346,7 +347,7 @@ end
 # ---------------------------------------------------------------------------------------------------
 function set_avatar_defaults(d, cmd, mcc, caller, got_usr_R, opt_B, opt_Gsymb, opt_ML, opt_R, opt_S, opt_W, sub_module, g_bar_fill, opt_Wmarker, opt_UVXY, opt_c, O, arg1)::Vector{String}
 	cmd  = check_caller(d, cmd, opt_S, opt_W, sub_module, g_bar_fill, O)
-	(mcc && caller == "bar" && !got_usr_R && opt_R != " -R") && (cmd = recompute_R_4bars!(cmd, opt_R, arg1))	# Often needed
+	(mcc && caller == "bar" && !got_usr_R && opt_R != " -R") && (cmd = recompute_R_4bars!(cmd, opt_R, arg1))	# Often needed  FORCES RECOMPILE
 	_cmd = build_run_cmd(cmd, opt_B, opt_Gsymb, opt_ML, opt_S, opt_W, opt_Wmarker, opt_UVXY, opt_c)
 	return _cmd
 end
@@ -1012,7 +1013,6 @@ function build_run_cmd(cmd, opt_B, opt_Gsymb, opt_ML, opt_S, opt_W, opt_Wmarker,
 			(opt_ML != "") && (cmd2 = cmd2 * opt_ML)			# If we have a symbol outline pen
 			_cmd = [cmd1; cmd2]
 		end
-
 	else
 		_cmd = [cmd * opt_UVXY]
 	end
@@ -1232,13 +1232,14 @@ function recompute_R_4bars!(cmd::String, opt_R::String, arg1)
 	sub_b = ((ind = findfirst("+", opt_S)) !== nothing) ? opt_S[ind[1]:end] : ""	# The +Base modifier
 	(sub_b != "") && (opt_S = opt_S[1:ind[1]-1])# Strip it because we need to (re)find Bar width
 	bw = (isletter(opt_S[end])) ? parse(Float64, opt_S[3:end-1]) : parse(Float64, opt_S[2:end])	# Bar width
-	info = gmt("gmtinfo -C", arg1)
+	info = gmt("gmtinfo -C", arg1)		# FORCES RECOMPILE
 	dx::Float64 = (info.data[2] - info.data[1]) * 0.005 + bw/2;
 	dy::Float64 = (info.data[4] - info.data[3]) * 0.005;
 	info.data[1] -= dx;	info.data[2] += dx;	info.data[4] += dy;
 	info.data = round_wesn(info.data)		# Add a pad if not-tight
 	new_opt_R::String = @sprintf(" -R%.15g/%.15g/%.15g/%.15g", info.data[1], info.data[2], 0, info.data[4])
 	cmd = replace(cmd, opt_R => new_opt_R)
+	return cmd
 end
 
 # ---------------------------------------------------------------------------------------------------
@@ -1762,15 +1763,17 @@ function sort_visible_triangles(Dv::Vector{<:GMTdataset}; del_hidden=false, zfac
 	(del_hidden != 1 && contains(Dv[1].comment[1], "vwall")) && (del_hidden = true)	# If have vwalls, need to del invis
 	if (del_hidden == 1)		# Remove the triangles that are not visible from the normal view_vec
 		bak_view = CURRENT_VIEW[1]	# Save because mapproject will reset it to "" (parsing on a module that has first = true)
-		t = isgeog(Dv) ? mapproject(Dv, J="t$((Dv[1].ds_bbox[1] + Dv[1].ds_bbox[2])/2)/1:1", C=true, F=true) : Dv
+		t = isgeog(Dv) ? mapproject(Dv, J="t$((Dv[1].ds_bbox[1] + Dv[1].ds_bbox[2])/2)/1:1", C=true, F=true) : Dv	# FORCES RECOMPILE
+		#t = isgeog(Dv) ? gmt_GMTdataset("mapproject " * "-Jt$((Dv[1].ds_bbox[1] + Dv[1].ds_bbox[2])/2)/1:1 -C -F") : Dv
 		CURRENT_VIEW[1] = bak_view 
 		view_vec = [sin_az * cosd(elev), cos_az * cosd(elev), sin_el]
 		is_vis = [dot(facenorm(t[k].data, zfact=zfact, normalize=false), view_vec) > 0 for k in eachindex(t)]
 		Dv = Dv[is_vis]
 	end
 
 	# ---------------------- Now sort by distance to the viewer ----------------------
-	Dc = gmtspatial(Dv, Q=true, o="0,1")	# Should this directly in Julia and avoid calling GMT (=> a copy of Dv)
+	#Dc = gmtspatial(Dv, Q=true, o="0,1")	# Should this directly in Julia and avoid calling GMT (=> a copy of Dv)
+	Dc::GMTdataset{Float64,2} = mat2ds(gmt_centroid_area(G_API[1], Dv, Int(isgeog(Dv)), ca=2), geom=wkbPoint)
 	dists = [(Dc.data[1,1] * sin_az + Dc.data[1,2] * cos_az, (Dv[1].bbox[5] + Dv[1].bbox[6]) / 2 * sin_el)]
 	for k = 2:size(Dc, 1)
 		push!(dists, (Dc.data[k,1] * sin_az + Dc.data[k,2] * cos_az, (Dv[k].bbox[5] + Dv[k].bbox[6]) / 2 * sin_el))