text:
@test supertype(GraphicsArea) == Any

@testset "BarplotGraphics" begin
    @test BarplotGraphics <: GraphicsArea
    @test !(BarplotGraphics <: Canvas)

    @testset "basic API and display" begin
        g = @inferred BarplotGraphics(0:2:10, 30)
        @test @inferred(nrows(g)) === 6
        @test @inferred(ncols(g)) === 30
        @test_throws ArgumentError printrow(stdout, g, 0)
        @test_throws ArgumentError printrow(stdout, g, 7)
        test_ref(
            "references/graphics/bar_printrow.txt",
            @io2str(printrow(IOContext(::IO, :color=>true), g, 3))
        )
        test_ref(
            "references/graphics/bar_print.txt",
            @io2str(print(IOContext(::IO, :color=>true), g))
        )
        test_ref(
            "references/graphics/bar_show.txt",
            @io2str(show(IOContext(::IO, :color=>true), g))
        )
        test_ref(
            "references/graphics/bar_print_nocolor.txt",
            @io2str(print(::IO, g))
        )
        test_ref(
            "references/graphics/bar_show_nocolor.txt",
            @io2str(show(::IO, g))
        )
        @test @inferred(addrow!(g, 3)) === g
        @test @inferred(nrows(g)) === 7
        @test @inferred(ncols(g)) === 30
        @test_throws InexactError addrow!(g, 3.5)
        test_ref(
            "references/graphics/bar_addrow_scalar.txt",
            @io2str(show(IOContext(::IO, :color=>true), g))
        )
        @test @inferred(addrow!(g, [50, 7])) === g
        @test_throws MethodError addrow!(g, [3.5])
        @test @inferred(nrows(g)) === 9
        @test @inferred(ncols(g)) === 30
        test_ref(
            "references/graphics/bar_addrow_vector.txt",
            @io2str(show(IOContext(::IO, :color=>true), g))
        )
    end

    @testset "Constructor parameter" begin
        @test_throws MethodError BarplotGraphics(0:2:10)
        @test_throws Union{ArgumentError,MethodError} BarplotGraphics(Int[], 20)
        @test_throws ArgumentError BarplotGraphics(0:2:10, 30, symb="--")
        g = @inferred BarplotGraphics([0, 0], 20)
        @test @inferred(nrows(g)) === 2
        @test @inferred(ncols(g)) === 20
        test_ref(
            "references/graphics/bar_empty.txt",
            @io2str(show(IOContext(::IO, :color=>true), g))
        )
        g = @inferred BarplotGraphics([1, 1], 20, color=:blue)
        @test @inferred(nrows(g)) === 2
        @test @inferred(ncols(g)) === 20
        test_ref(
            "references/graphics/bar_blue.txt",
            @io2str(show(IOContext(::IO, :color=>true), g))
        )
        g = @inferred BarplotGraphics([0, 1, 10, 100, 1000], 20, log10)
        @test @inferred(nrows(g)) === 5
        @test @inferred(ncols(g)) === 20
        test_ref(
            "references/graphics/bar_log10.txt",
            @io2str(show(IOContext(::IO, :color=>true), g))
        )
        g = @inferred BarplotGraphics([0, 1, 10], 20, symb="#")
        @test @inferred(nrows(g)) === 3
        @test @inferred(ncols(g)) === 20
        test_ref(
            "references/graphics/bar_symb.txt",
            @io2str(show(IOContext(::IO, :color=>true), g))
        )
        g = @inferred BarplotGraphics([0, 1, 9], 1)
        @test @inferred(nrows(g)) === 3
        @test @inferred(ncols(g)) === 10
        test_ref(
            "references/graphics/bar_short.txt",
            @io2str(show(IOContext(::IO, :color=>true), g))
        )
        g = @inferred BarplotGraphics([0, 1, 10000000], 1)
        @test @inferred(nrows(g)) === 3
        @test @inferred(ncols(g)) === 15
        test_ref(
            "references/graphics/bar_short2.txt",
            @io2str(show(IOContext(::IO, :color=>true), g))
        )
    end
end

@testset "BoxplotGraphics" begin
    @test BoxplotGraphics <: GraphicsArea
    @test !(BoxplotGraphics <: Canvas)

    @testset "basic API and display" begin
        g = @inferred BoxplotGraphics([1,2,2,4,5,6], 30)
        @test @inferred(nrows(g)) === 3
        @test @inferred(ncols(g)) === 30
        @test_throws ArgumentError printrow(stdout, g, 0)
        @test_throws ArgumentError printrow(stdout, g, 4)
        test_ref(
            "references/graphics/box_printrow.txt",
            @io2str(printrow(IOContext(::IO, :color=>true), g, 3))
        )
        test_ref(
            "references/graphics/box_print.txt",
            @io2str(print(IOContext(::IO, :color=>true), g))
        )
        test_ref(
            "references/graphics/box_show.txt",
            @io2str(show(IOContext(::IO, :color=>true), g))
        )
        test_ref(
            "references/graphics/box_print_nocolor.txt",
            @io2str(print(::IO, g))
        )
        test_ref(
            "references/graphics/box_show_nocolor.txt",
            @io2str(show(::IO, g))
        )
        @test @inferred(UnicodePlots.addseries!(g, [3.,3,2,4,0,1])) === g
        @test @inferred(nrows(g)) === 6
        @test @inferred(ncols(g)) === 30
        test_ref(
            "references/graphics/box_addseries.txt",
            @io2str(show(IOContext(::IO, :color=>true), g))
        )
        @test @inferred(UnicodePlots.addseries!(g, [2,2,2])) === g
        @test @inferred(nrows(g)) === 9
        @test @inferred(ncols(g)) === 30
        test_ref(
            "references/graphics/box_addseries_small.txt",
            @io2str(show(IOContext(::IO, :color=>true), g))
        )
    end

    @testset "Constructor parameter" begin
        @test_throws MethodError BoxplotGraphics([1,2,2,4,5,6])
        @test_throws Union{ArgumentError,MethodError} BoxplotGraphics(Int[], 30)
        g = @inferred BoxplotGraphics([0, 0], 20)
        @test @inferred(nrows(g)) === 3
        @test @inferred(ncols(g)) === 20
        test_ref(
            "references/graphics/box_empty.txt",
            @io2str(show(IOContext(::IO, :color=>true), g))
        )
        g = @inferred BoxplotGraphics([1, 2], 20, color=:blue)
        @test @inferred(nrows(g)) === 3
        @test @inferred(ncols(g)) === 20
        test_ref(
            "references/graphics/box_blue.txt",
            @io2str(show(IOContext(::IO, :color=>true), g))
        )
        g = @inferred BoxplotGraphics([0, 1, 9], 1)
        @test @inferred(nrows(g)) === 3
        @test @inferred(ncols(g)) === 10
        test_ref(
            "references/graphics/box_short.txt",
            @io2str(show(IOContext(::IO, :color=>true), g))
        )
        g = @inferred BoxplotGraphics([1,2,2,4,5,6], 30, min_x = 2, max_x = 8)
        @test @inferred(nrows(g)) === 3
        @test @inferred(ncols(g)) === 30
        test_ref(
            "references/graphics/box_minmax1.txt",
            @io2str(show(IOContext(::IO, :color=>true), g))
        )
        g = @inferred BoxplotGraphics([1,2,2,4,5,6], 30, min_x = 0, max_x = 7)
        @test @inferred(nrows(g)) === 3
        @test @inferred(ncols(g)) === 30
        test_ref(
            "references/graphics/box_minmax2.txt",
            @io2str(show(IOContext(::IO, :color=>true), g))
        )
        # FIXME: boxplot break if max_x or min_x is doesn't conver the data
    end
end


text:
# This file is a part of Julia. License is MIT: https://julialang.org/license

using Test
using Distributed
using Dates
import REPL
using Printf: @sprintf

include("choosetests.jl")
include("testenv.jl")

tests, net_on, exit_on_error, use_revise, seed = choosetests(ARGS)
tests = unique(tests)

if use_revise
    using Revise
end

const max_worker_rss = if haskey(ENV, "JULIA_TEST_MAXRSS_MB")
    parse(Int, ENV["JULIA_TEST_MAXRSS_MB"]) * 2^20
else
    typemax(Csize_t)
end
limited_worker_rss = max_worker_rss != typemax(Csize_t)

function test_path(test)
    t = split(test, '/')
    if t[1] in STDLIBS
        if length(t) == 2
            return joinpath(STDLIB_DIR, t[1], "test", t[2])
        else
            return joinpath(STDLIB_DIR, t[1], "test", "runtests")
        end
    else
        return joinpath(@__DIR__, test)
    end
end

# Check all test files exist
isfiles = isfile.(test_path.(tests) .* ".jl")
if !all(isfiles)
    error("did not find test files for the following tests: ",
          join(tests[.!(isfiles)], ", "))
end

const node1_tests = String[]
function move_to_node1(t)
    if t in tests
        splice!(tests, findfirst(isequal(t), tests))
        push!(node1_tests, t)
    end
    nothing
end

# Base.compilecache only works from node 1, so precompile test is handled specially
move_to_node1("precompile")
move_to_node1("SharedArrays")
move_to_node1("threads")
# Ensure things like consuming all kernel pipe memory doesn't interfere with other tests
move_to_node1("stress")

# In a constrained memory environment, run the "distributed" test after all other tests
# since it starts a lot of workers and can easily exceed the maximum memory
limited_worker_rss && move_to_node1("Distributed")

import LinearAlgebra
cd(@__DIR__) do
    n = 1
    if net_on
        n = min(Sys.CPU_THREADS, length(tests))
        n > 1 && addprocs_with_testenv(n)
        LinearAlgebra.BLAS.set_num_threads(1)
    end
    skipped = 0

    @everywhere include("testdefs.jl")

    if use_revise
        @everywhere begin
            Revise.track(Core.Compiler)
            Revise.track(Base)
            Revise.revise()
        end
    end

    #pretty print the information about gc and mem usage
    testgroupheader = "Test"
    workerheader = "(Worker)"
    name_align    = maximum([textwidth(testgroupheader) + textwidth(" ") + textwidth(workerheader); map(x -> textwidth(x) + 3 + ndigits(nworkers()), tests)])
    elapsed_align = textwidth("Time (s)")
    gc_align      = textwidth("GC (s)")
    percent_align = textwidth("GC %")
    alloc_align   = textwidth("Alloc (MB)")
    rss_align     = textwidth("RSS (MB)")
    printstyled(testgroupheader, color=:white)
    printstyled(lpad(workerheader, name_align - textwidth(testgroupheader) + 1), " | ", color=:white)
    printstyled("Time (s) | GC (s) | GC % | Alloc (MB) | RSS (MB)\n", color=:white)
    results = []
    print_lock = stdout isa Base.LibuvStream ? stdout.lock : ReentrantLock()
    if stderr isa Base.LibuvStream
        stderr.lock = print_lock
    end

    function print_testworker_stats(test, wrkr, resp)
        @nospecialize resp
        lock(print_lock)
        try
            printstyled(test, color=:white)
            printstyled(lpad("($wrkr)", name_align - textwidth(test) + 1, " "), " | ", color=:white)
            time_str = @sprintf("%7.2f",resp[2])
            printstyled(lpad(time_str, elapsed_align, " "), " | ", color=:white)
            gc_str = @sprintf("%5.2f", resp[5].total_time / 10^9)
            printstyled(lpad(gc_str, gc_align, " "), " | ", color=:white)

            # since there may be quite a few digits in the percentage,
            # the left-padding here is less to make sure everything fits
            percent_str = @sprintf("%4.1f", 100 * resp[5].total_time / (10^9 * resp[2]))
            printstyled(lpad(percent_str, percent_align, " "), " | ", color=:white)
            alloc_str = @sprintf("%5.2f", resp[3] / 2^20)
            printstyled(lpad(alloc_str, alloc_align, " "), " | ", color=:white)
            rss_str = @sprintf("%5.2f", resp[6] / 2^20)
            printstyled(lpad(rss_str, rss_align, " "), "\n", color=:white)
        finally
            unlock(print_lock)
        end
    end

    global print_testworker_started = (name, wrkr)->begin
        lock(print_lock)
        try
            printstyled(name, color=:white)
            printstyled(lpad("($wrkr)", name_align - textwidth(name) + 1, " "), " |",
                " "^elapsed_align, "started at $(now())\n", color=:white)
        finally
            unlock(print_lock)
        end
    end

    function print_testworker_errored(name, wrkr)
        lock(print_lock)
        try
            printstyled(name, color=:red)
            printstyled(lpad("($wrkr)", name_align - textwidth(name) + 1, " "), " |",
                " "^elapsed_align, " failed at $(now())\n", color=:red)
        finally
            unlock(print_lock)
        end
    end

    all_tests = [tests; node1_tests]

    local stdin_monitor
    all_tasks = Task[]
    try
        # Monitor stdin and kill this task on ^C
        # but don't do this on Windows, because it may deadlock in the kernel
        if !Sys.iswindows() && isa(stdin, Base.TTY)
            t = current_task()
            stdin_monitor = @async begin
                term = REPL.Terminals.TTYTerminal("xterm", stdin, stdout, stderr)
                try
                    REPL.Terminals.raw!(term, true)
                    while true
                        if read(term, Char) == '\x3'
                            Base.throwto(t, InterruptException())
                            break
                        end
                    end
                catch e
                    isa(e, InterruptException) || rethrow()
                finally
                    REPL.Terminals.raw!(term, false)
                end
            end
        end
        @sync begin
            for p in workers()
                @async begin
                    push!(all_tasks, current_task())
                    while length(tests) > 0
                        test = popfirst!(tests)
                        local resp
                        wrkr = p
                        try
                            resp = remotecall_fetch(runtests, wrkr, test, test_path(test); seed=seed)
                        catch e
                            isa(e, InterruptException) && return
                            resp = Any[e]
                        end
                        push!(results, (test, resp))
                        if resp[1] isa Exception
                            print_testworker_errored(test, wrkr)
                            if exit_on_error
                                skipped = length(tests)
                                empty!(tests)
                            elseif n > 1
                                # the worker encountered some failure, recycle it
                                # so future tests get a fresh environment
                                rmprocs(wrkr, waitfor=30)
                                p = addprocs_with_testenv(1)[1]
                                remotecall_fetch(include, p, "testdefs.jl")
                            end
                        else
                            print_testworker_stats(test, wrkr, resp)
                            if resp[end] > max_worker_rss
                                # the worker has reached the max-rss limit, recycle it
                                # so future tests start with a smaller working set
                                if n > 1
                                    rmprocs(wrkr, waitfor=30)
                                    p = addprocs_with_testenv(1)[1]
                                    remotecall_fetch(include, p, "testdefs.jl")
                                else # single process testing
                                    error("Halting tests. Memory limit reached : $resp > $max_worker_rss")
                                end
                            end
                       end
                    end
                    if p != 1
                        # Free up memory =)
                        rmprocs(p, waitfor=30)
                    end
                end
            end
        end

        n > 1 && length(node1_tests) > 1 && print("\nExecuting tests that run on node 1 only:\n")
        for t in node1_tests
            # As above, try to run each test
            # which must run on node 1. If
            # the test fails, catch the error,
            # and either way, append the results
            # to the overall aggregator
            isolate = true
            t == "SharedArrays" && (isolate = false)
            local resp
            try
                resp = eval(Expr(:call, () -> runtests(t, test_path(t), isolate, seed=seed))) # runtests is defined by the include above
                print_testworker_stats(t, 1, resp)
            catch e
                resp = Any[e]
            end
            push!(results, (t, resp))
        end
    catch e
        isa(e, InterruptException) || rethrow()
        # If the test suite was merely interrupted, still print the
        # summary, which can be useful to diagnose what's going on
        foreach(task -> begin
                istaskstarted(task) || return
                istaskdone(task) && return
                try
                    schedule(task, InterruptException(); error=true)
                catch ex
                    @error "InterruptException" exception=ex,catch_backtrace()
                end
            end, all_tasks)
        foreach(wait, all_tasks)
    finally
        if @isdefined stdin_monitor
            schedule(stdin_monitor, InterruptException(); error=true)
        end
    end

    #=
`   Construct a testset on the master node which will hold results from all the
    test files run on workers and on node1. The loop goes through the results,
    inserting them as children of the overall testset if they are testsets,
    handling errors otherwise.

    Since the workers don't return information about passing/broken tests, only
    errors or failures, those Result types get passed `nothing` for their test
    expressions (and expected/received result in the case of Broken).

    If a test failed, returning a `RemoteException`, the error is displayed and
    the overall testset has a child testset inserted, with the (empty) Passes
    and Brokens from the worker and the full information about all errors and
    failures encountered running the tests. This information will be displayed
    as a summary at the end of the test run.

    If a test failed, returning an `Exception` that is not a `RemoteException`,
    it is likely the julia process running the test has encountered some kind
    of internal error, such as a segfault.  The entire testset is marked as
    Errored, and execution continues until the summary at the end of the test
    run, where the test file is printed out as the "failed expression".
    =#
    o_ts = Test.DefaultTestSet("Overall")
    Test.push_testset(o_ts)
    completed_tests = Set{String}()
    for (testname, (resp,)) in results
        push!(completed_tests, testname)
        if isa(resp, Test.DefaultTestSet)
            Test.push_testset(resp)
            Test.record(o_ts, resp)
            Test.pop_testset()
        elseif isa(resp, Tuple{Int,Int})
            fake = Test.DefaultTestSet(testname)
            for i in 1:resp[1]
                Test.record(fake, Test.Pass(:test, nothing, nothing, nothing))
            end
            for i in 1:resp[2]
                Test.record(fake, Test.Broken(:test, nothing))
            end
            Test.push_testset(fake)
            Test.record(o_ts, fake)
            Test.pop_testset()
        elseif isa(resp, RemoteException) && isa(resp.captured.ex, Test.TestSetException)
            println("Worker $(resp.pid) failed running test $(testname):")
            Base.showerror(stdout, resp.captured)
            println()
            fake = Test.DefaultTestSet(testname)
            for i in 1:resp.captured.ex.pass
                Test.record(fake, Test.Pass(:test, nothing, nothing, nothing))
            end
            for i in 1:resp.captured.ex.broken
                Test.record(fake, Test.Broken(:test, nothing))
            end
            for t in resp.captured.ex.errors_and_fails
                Test.record(fake, t)
            end
            Test.push_testset(fake)
            Test.record(o_ts, fake)
            Test.pop_testset()
        else
            if !isa(resp, Exception)
                resp = ErrorException(string("Unknown result type : ", typeof(resp)))
            end
            # If this test raised an exception that is not a remote testset exception,
            # i.e. not a RemoteException capturing a TestSetException that means
            # the test runner itself had some problem, so we may have hit a segfault,
            # deserialization errors or something similar.  Record this testset as Errored.
            fake = Test.DefaultTestSet(testname)
            Test.record(fake, Test.Error(:test_error, testname, nothing, Any[(resp, [])], LineNumberNode(1)))
            Test.push_testset(fake)
            Test.record(o_ts, fake)
            Test.pop_testset()
        end
    end
    for test in all_tests
        (test in completed_tests) && continue
        fake = Test.DefaultTestSet(test)
        Test.record(fake, Test.Error(:test_interrupted, test, nothing, [("skipped", [])], LineNumberNode(1)))
        Test.push_testset(fake)
        Test.record(o_ts, fake)
        Test.pop_testset()
    end
    println()
    Test.print_test_results(o_ts, 1)
    if !o_ts.anynonpass
        println("    \033[32;1mSUCCESS\033[0m")
    else
        println("    \033[31;1mFAILURE\033[0m\n")
        skipped > 0 &&
            println("$skipped test", skipped > 1 ? "s were" : " was", " skipped due to failure.")
        println("The global RNG seed was 0x$(string(seed, base = 16)).\n")
        Test.print_test_errors(o_ts)
        throw(Test.FallbackTestSetException("Test run finished with errors"))
    end
end

ext: jl


text:
export call, exists, getService, isValid

"""
    call(serviceClient, srv)

Call the service aliased by this handle with the specified service request/response.
"""
call(serviceClient::Union{Cxx.CxxCore.CppPtr{Cxx.CxxCore.CppValue{Cxx.CxxCore.CxxQualType{Cxx.CxxCore.CppBaseType{Symbol("ros::ServiceClient")},(false, false, false)},N} where N,(false, false, false)}, Cxx.CxxCore.CppValue{Cxx.CxxCore.CxxQualType{Cxx.CxxCore.CppBaseType{Symbol("ros::ServiceClient")},(false, false, false)},N} where N}, srv::Cxx.CxxCore.CppValue) = @cxx serviceClient->call(srv)

call(serviceClient::Union{Cxx.CxxCore.CppPtr{Cxx.CxxCore.CppValue{Cxx.CxxCore.CxxQualType{Cxx.CxxCore.CppBaseType{Symbol("ros::ServiceClient")},(false, false, false)},N} where N,(false, false, false)}, Cxx.CxxCore.CppValue{Cxx.CxxCore.CxxQualType{Cxx.CxxCore.CppBaseType{Symbol("ros::ServiceClient")},(false, false, false)},N} where N}, srv::Cxx.CxxCore.CppPtr) = @cxx serviceClient->call(*(srv))

"""
    exists(serviceClient)

Checks if this is both advertised and available.
"""
exists(serviceClient::Union{Cxx.CxxCore.CppPtr{Cxx.CxxCore.CppValue{Cxx.CxxCore.CxxQualType{Cxx.CxxCore.CppBaseType{Symbol("ros::ServiceClient")},(false, false, false)},N} where N,(false, false, false)}, Cxx.CxxCore.CppValue{Cxx.CxxCore.CxxQualType{Cxx.CxxCore.CppBaseType{Symbol("ros::ServiceClient")},(false, false, false)},N} where N}) = @cxx serviceClient->exists()

"""
    getService(serviceClient)

Returns the name of the service this ServiceClient connects to.
"""
getService(serviceClient::Union{Cxx.CxxCore.CppPtr{Cxx.CxxCore.CppValue{Cxx.CxxCore.CxxQualType{Cxx.CxxCore.CppBaseType{Symbol("ros::ServiceClient")},(false, false, false)},N} where N,(false, false, false)}, Cxx.CxxCore.CppValue{Cxx.CxxCore.CxxQualType{Cxx.CxxCore.CppBaseType{Symbol("ros::ServiceClient")},(false, false, false)},N} where N}) = unsafe_string(@cxx serviceClient->getService())

"""
    isPersistent(serviceClient)

Returns true if this handle points to a persistent service, false otherwise.
"""
isPersistent(serviceClient::Union{Cxx.CxxCore.CppPtr{Cxx.CxxCore.CppValue{Cxx.CxxCore.CxxQualType{Cxx.CxxCore.CppBaseType{Symbol("ros::ServiceClient")},(false, false, false)},N} where N,(false, false, false)}, Cxx.CxxCore.CppValue{Cxx.CxxCore.CxxQualType{Cxx.CxxCore.CppBaseType{Symbol("ros::ServiceClient")},(false, false, false)},N} where N}) = @cxx serviceClient->isPersistent()

"""
    isValid(serviceClient)

Returns whether or not this handle is valid. For a persistent service, this becomes false when the connection has dropped. Non-persistent service handles are always valid.
"""
isValid(serviceClient::Union{Cxx.CxxCore.CppPtr{Cxx.CxxCore.CppValue{Cxx.CxxCore.CxxQualType{Cxx.CxxCore.CppBaseType{Symbol("ros::ServiceClient")},(false, false, false)},N} where N,(false, false, false)}, Cxx.CxxCore.CppValue{Cxx.CxxCore.CxxQualType{Cxx.CxxCore.CppBaseType{Symbol("ros::ServiceClient")},(false, false, false)},N} where N}) = @cxx serviceClient->isValid()


text:
mutable struct Sink{T} inner::Union{Nothing,T} end
struct Source{T <: Union{Variable,Parameter}} end
@inline (F::Sink{Field})(y,x,p=nothing) = inner(F)(y,x,p)

setindex!(s::Sink{Field},e,i) = (s.inner = Field(inner(s),IndexedExpression(i,e)))
getindex(::Source{Variable},n) = Variable(n)
getindex(::Source{Parameter},n) = Parameter(n)

Variable() = Source{Variable}()
Parameter() = Source{Parameter}()


text:
__precompile__()

# Mortgage calculator functions

using Base.Dates

"Return length of the mortgage or loan term in years."
function termlength(term, amortization)
	# use a ? : construct as an if - then - else oneliner
	termLength = haskey(mortgageTermDescriptors, term) ? mortgageTermDescriptors[term][2] : parse(Int, term)
	if termLength == "a"
		termLength = amortization
	elseif termLength > amortization
		error("The specified term exceeds the amortization!")
	end
	return termLength
end

"Recurse to calculate for each payment."
function tablerow(rate, frequency, compounding, startdate, paymentSize, paymentNumber, numberOfPayments, paymentDate, principalPortion, interestPortion, accumulatedPrincipal, accumulatedInterest, accumulatedTotal, balance, table)
	if balance == 0 # no more payments to process
		return accumulatedPrincipal, accumulatedInterest, accumulatedTotal, balance, principalPortion + interestPortion, paymentDate, table
	else # get values for next payment
		nextPaymentDate = nextpaymentdate(paymentDate, frequency, startdate)
		interestPortion = interest(balance, rate, timeValues[frequency], compounding, paymentDate)
		principalPortion = paymentSize - interestPortion
		if principalPortion > balance # ie, we would be overpaying the loan
			return tablerow(rate, frequency, compounding, startdate, paymentSize, paymentNumber + 1, numberOfPayments, nextPaymentDate, balance, interestPortion, accumulatedPrincipal + balance, accumulatedInterest + interestPortion, accumulatedTotal + balance + interestPortion, 0, append!(table, [paymentNumber + 1, nextPaymentDate, principalPortion, interestPortion, accumulatedPrincipal + principalPortion, accumulatedInterest + interestPortion, accumulatedTotal + paymentSize, balance - principalPortion]))
		elseif paymentNumber < numberOfPayments # still more payments to go...
			return tablerow(rate, frequency, compounding, startdate, paymentSize, paymentNumber + 1, numberOfPayments, nextPaymentDate, principalPortion, interestPortion, accumulatedPrincipal + principalPortion, accumulatedInterest + interestPortion, accumulatedTotal + principalPortion + interestPortion, balance - principalPortion, append!(table, [paymentNumber + 1, nextPaymentDate, principalPortion, interestPortion, accumulatedPrincipal + principalPortion, accumulatedInterest + interestPortion, accumulatedTotal + paymentSize, balance - principalPortion]))
		else # ie the number of payments == paymentNumber, thus we're done
			return accumulatedPrincipal, accumulatedInterest, accumulatedTotal, balance, principalPortion + interestPortion, paymentDate, table
		end
	end
end

"Call tablerow with arguments initialized to starting values."
function dotablerows(principal, rate, frequency, compounding, startdate, firstPaymentDate, paymentSize, numberOfPayments)
	interestAmount = interest(principal, rate, timeValues[frequency], compounding, firstPaymentDate)
	table = [1, firstPaymentDate, paymentSize - interestAmount, interestAmount, paymentSize - interestAmount, interestAmount, paymentSize, principal - (paymentSize -interestAmount)]
	return tablerow(rate, frequency, compounding, startdate, paymentSize, 1, numberOfPayments, firstPaymentDate, paymentSize - interestAmount, interestAmount, paymentSize - interestAmount, interestAmount, paymentSize, principal - (paymentSize - interestAmount), table)
end

"""
Calculate the interest amount for a single time period.
Simple interest formula is I = principal x rate x time
"""
function interest(principal, rate, timeValue, compounding, paymentDate)

	rate = rate * 0.01 # convert percent to fraction
	if compounding == "s"
		rate = rate / 365 # calculate daily interest rate; note that for some commercial mortgages, this value should be 360
		if timeValue == 1//12 || timeValue == 1//24 # in these situations, the amount of interest depends on how many days in the month
			daysInMonth = daysinmonth(paymentDate - Month(1)) # get the number of days in the preceding month
			if timeValue == 1//12
				return  daysInMonth * rate * principal
			else # ie timeValue is 1//24
				paymentDay = day(paymentDate)
				if paymentDay == 1
					return (daysInMonth - 15) * rate * principal
				else
					return 15 * rate * principal
				end
			end
		else
			return timeValue * rate * principal
		end
	else
		if timeValue > 1 # ie 7 or 14 days
			timeValue = timeValue / 365
		end
		if compounding == "m" # ie american (monthly) compounding or simple interest
			return principal * rate * timeValue
		elseif compounding == "c" # for Canadian (semi-annual) compounding
			return ((1 + rate / 2) ^ (2 * timeValue) - 1) * principal
		end
	end
end

"""
Return the size of each payment, when not specified by the user.
For compound interest, the formula is:
P = r(PV) / (1 - (1+r)e(-n))
where P = payment, PV = Present Value, r = rate per period, n = number of periods

For simple interest, the formula is:
payment = (principal + interestAmount) / number of periods, where
interestAmount = principal x interest x amortization
"""
function paymentsizecalc(annualRate, principal, amortization, compounding, frequency, startDate)
	interest = annualRate / 100 # convert percentage to decimal
	timeValue = timeValues[frequency]
	if timeValue > 1 # ie if 7 or 14 days
		timeValue = timeValue / 365
	end
	if compounding == "s" # ie, simple interest (no compounding, but interest accrues daily)
		return (principal + (principal * interest * amortization)) / numberofperiods(frequency, amortization, startDate)
	elseif compounding == "m" # ie, monthly, or American, compounding
		num = interest * timeValue * principal
		denom = 1 - (1 + interest * timeValue) ^ (- amortization / timeValue)
		return num / denom
	elseif compounding == "c" # ie, Canadian, or semi-annual, compounding
		mi = (1 + interest/2)^(2 * timeValue)
		return principal * (mi - 1) / (1 - mi ^ (-amortization / timeValue))
	else
		error("The compounding specified as '$compounding' is not defined.")
	end
end

"""
Return the number of payment periods, given frequency and duration.
Starting date is needed to deal with leap years.
For monthly or bimonthly payment frequencies, use 12 or 24 periods per year, respectively.
For weekly or two-weekly payment frequencies, calculate the number of periods over the
total amortization interval: get the number of days, and divide by 7 or 14, respectively.
Duration is specified in years or a fraction thereof.
"""
function numberofperiods(frequency, duration, startDate)
	if 	frequency == "m"
		n = 12 * duration
	elseif frequency == "b"
		n = 24 * duration
	else
		if duration < 1 # calculate using months instead of years
			days = (startDate + Month(duration * 12)) - startDate
		else
			days = (startDate + Year(duration)) - startDate
		end
		if frequency == "w"
			n = Int(days) / 7
		elseif frequency == "t"
			n = Int(days) / 14
		else error("The payment frequency specified as '$frequency' is not defined!")
		end
	end
	return round(Int, n)
end

"Return date for the next payment."
function nextpaymentdate(paymentDate, frequency, startDate)
	if frequency == "m" # ie, monthly payments
		# test if the starting date is on or after the 28th of the month
		if 28 <= day(startDate) <= daysinmonth(startDate)
			days = daysinmonth(paymentDate) - day(paymentDate)
			# set a payment date which is the same number of days before month's end
			return lastdayofmonth(paymentDate + Month(1)) - Day(days)
		else
			return paymentDate + Month(1)
		end
	elseif frequency == "b" # ie, payments twice a month
		if day(paymentDate) ==1
			return paymentDate + Day(14)
		elseif day(paymentDate) == 15
			return paymentDate - Day(14) + Month(1)
		else
			error("bimonthly payments can only be on the 1st or 15th of each month. If necessary, set the starting date to either the 1st or the 15th.")
		end
	elseif frequency == "w"
		return paymentDate + Week(1)
	elseif frequency == "t"
		return paymentDate + Week(2)
	else
		error("not defined")
	end
end


text:
loadcolorscheme(:magma, [RGB(0.001462, 0.000466, 0.013866),
RGB(0.002258, 0.001295, 0.018331),
RGB(0.003279, 0.002305, 0.023708),
RGB(0.004512, 0.003490, 0.029965),
RGB(0.005950, 0.004843, 0.037130),
RGB(0.007588, 0.006356, 0.044973),
RGB(0.009426, 0.008022, 0.052844),
RGB(0.011465, 0.009828, 0.060750),
RGB(0.013708, 0.011771, 0.068667),
RGB(0.016156, 0.013840, 0.076603),
RGB(0.018815, 0.016026, 0.084584),
RGB(0.021692, 0.018320, 0.092610),
RGB(0.024792, 0.020715, 0.100676),
RGB(0.028123, 0.023201, 0.108787),
RGB(0.031696, 0.025765, 0.116965),
RGB(0.035520, 0.028397, 0.125209),
RGB(0.039608, 0.031090, 0.133515),
RGB(0.043830, 0.033830, 0.141886),
RGB(0.048062, 0.036607, 0.150327),
RGB(0.052320, 0.039407, 0.158841),
RGB(0.056615, 0.042160, 0.167446),
RGB(0.060949, 0.044794, 0.176129),
RGB(0.065330, 0.047318, 0.184892),
RGB(0.069764, 0.049726, 0.193735),
RGB(0.074257, 0.052017, 0.202660),
RGB(0.078815, 0.054184, 0.211667),
RGB(0.083446, 0.056225, 0.220755),
RGB(0.088155, 0.058133, 0.229922),
RGB(0.092949, 0.059904, 0.239164),
RGB(0.097833, 0.061531, 0.248477),
RGB(0.102815, 0.063010, 0.257854),
RGB(0.107899, 0.064335, 0.267289),
RGB(0.113094, 0.065492, 0.276784),
RGB(0.118405, 0.066479, 0.286321),
RGB(0.123833, 0.067295, 0.295879),
RGB(0.129380, 0.067935, 0.305443),
RGB(0.135053, 0.068391, 0.315000),
RGB(0.140858, 0.068654, 0.324538),
RGB(0.146785, 0.068738, 0.334011),
RGB(0.152839, 0.068637, 0.343404),
RGB(0.159018, 0.068354, 0.352688),
RGB(0.165308, 0.067911, 0.361816),
RGB(0.171713, 0.067305, 0.370771),
RGB(0.178212, 0.066576, 0.379497),
RGB(0.184801, 0.065732, 0.387973),
RGB(0.191460, 0.064818, 0.396152),
RGB(0.198177, 0.063862, 0.404009),
RGB(0.204935, 0.062907, 0.411514),
RGB(0.211718, 0.061992, 0.418647),
RGB(0.218512, 0.061158, 0.425392),
RGB(0.225302, 0.060445, 0.431742),
RGB(0.232077, 0.059889, 0.437695),
RGB(0.238826, 0.059517, 0.443256),
RGB(0.245543, 0.059352, 0.448436),
RGB(0.252220, 0.059415, 0.453248),
RGB(0.258857, 0.059706, 0.457710),
RGB(0.265447, 0.060237, 0.461840),
RGB(0.271994, 0.060994, 0.465660),
RGB(0.278493, 0.061978, 0.469190),
RGB(0.284951, 0.063168, 0.472451),
RGB(0.291366, 0.064553, 0.475462),
RGB(0.297740, 0.066117, 0.478243),
RGB(0.304081, 0.067835, 0.480812),
RGB(0.310382, 0.069702, 0.483186),
RGB(0.316654, 0.071690, 0.485380),
RGB(0.322899, 0.073782, 0.487408),
RGB(0.329114, 0.075972, 0.489287),
RGB(0.335308, 0.078236, 0.491024),
RGB(0.341482, 0.080564, 0.492631),
RGB(0.347636, 0.082946, 0.494121),
RGB(0.353773, 0.085373, 0.495501),
RGB(0.359898, 0.087831, 0.496778),
RGB(0.366012, 0.090314, 0.497960),
RGB(0.372116, 0.092816, 0.499053),
RGB(0.378211, 0.095332, 0.500067),
RGB(0.384299, 0.097855, 0.501002),
RGB(0.390384, 0.100379, 0.501864),
RGB(0.396467, 0.102902, 0.502658),
RGB(0.402548, 0.105420, 0.503386),
RGB(0.408629, 0.107930, 0.504052),
RGB(0.414709, 0.110431, 0.504662),
RGB(0.420791, 0.112920, 0.505215),
RGB(0.426877, 0.115395, 0.505714),
RGB(0.432967, 0.117855, 0.506160),
RGB(0.439062, 0.120298, 0.506555),
RGB(0.445163, 0.122724, 0.506901),
RGB(0.451271, 0.125132, 0.507198),
RGB(0.457386, 0.127522, 0.507448),
RGB(0.463508, 0.129893, 0.507652),
RGB(0.469640, 0.132245, 0.507809),
RGB(0.475780, 0.134577, 0.507921),
RGB(0.481929, 0.136891, 0.507989),
RGB(0.488088, 0.139186, 0.508011),
RGB(0.494258, 0.141462, 0.507988),
RGB(0.500438, 0.143719, 0.507920),
RGB(0.506629, 0.145958, 0.507806),
RGB(0.512831, 0.148179, 0.507648),
RGB(0.519045, 0.150383, 0.507443),
RGB(0.525270, 0.152569, 0.507192),
RGB(0.531507, 0.154739, 0.506895),
RGB(0.537755, 0.156894, 0.506551),
RGB(0.544015, 0.159033, 0.506159),
RGB(0.550287, 0.161158, 0.505719),
RGB(0.556571, 0.163269, 0.505230),
RGB(0.562866, 0.165368, 0.504692),
RGB(0.569172, 0.167454, 0.504105),
RGB(0.575490, 0.169530, 0.503466),
RGB(0.581819, 0.171596, 0.502777),
RGB(0.588158, 0.173652, 0.502035),
RGB(0.594508, 0.175701, 0.501241),
RGB(0.600868, 0.177743, 0.500394),
RGB(0.607238, 0.179779, 0.499492),
RGB(0.613617, 0.181811, 0.498536),
RGB(0.620005, 0.183840, 0.497524),
RGB(0.626401, 0.185867, 0.496456),
RGB(0.632805, 0.187893, 0.495332),
RGB(0.639216, 0.189921, 0.494150),
RGB(0.645633, 0.191952, 0.492910),
RGB(0.652056, 0.193986, 0.491611),
RGB(0.658483, 0.196027, 0.490253),
RGB(0.664915, 0.198075, 0.488836),
RGB(0.671349, 0.200133, 0.487358),
RGB(0.677786, 0.202203, 0.485819),
RGB(0.684224, 0.204286, 0.484219),
RGB(0.690661, 0.206384, 0.482558),
RGB(0.697098, 0.208501, 0.480835),
RGB(0.703532, 0.210638, 0.479049),
RGB(0.709962, 0.212797, 0.477201),
RGB(0.716387, 0.214982, 0.475290),
RGB(0.722805, 0.217194, 0.473316),
RGB(0.729216, 0.219437, 0.471279),
RGB(0.735616, 0.221713, 0.469180),
RGB(0.742004, 0.224025, 0.467018),
RGB(0.748378, 0.226377, 0.464794),
RGB(0.754737, 0.228772, 0.462509),
RGB(0.761077, 0.231214, 0.460162),
RGB(0.767398, 0.233705, 0.457755),
RGB(0.773695, 0.236249, 0.455289),
RGB(0.779968, 0.238851, 0.452765),
RGB(0.786212, 0.241514, 0.450184),
RGB(0.792427, 0.244242, 0.447543),
RGB(0.798608, 0.247040, 0.444848),
RGB(0.804752, 0.249911, 0.442102),
RGB(0.810855, 0.252861, 0.439305),
RGB(0.816914, 0.255895, 0.436461),
RGB(0.822926, 0.259016, 0.433573),
RGB(0.828886, 0.262229, 0.430644),
RGB(0.834791, 0.265540, 0.427671),
RGB(0.840636, 0.268953, 0.424666),
RGB(0.846416, 0.272473, 0.421631),
RGB(0.852126, 0.276106, 0.418573),
RGB(0.857763, 0.279857, 0.415496),
RGB(0.863320, 0.283729, 0.412403),
RGB(0.868793, 0.287728, 0.409303),
RGB(0.874176, 0.291859, 0.406205),
RGB(0.879464, 0.296125, 0.403118),
RGB(0.884651, 0.300530, 0.400047),
RGB(0.889731, 0.305079, 0.397002),
RGB(0.894700, 0.309773, 0.393995),
RGB(0.899552, 0.314616, 0.391037),
RGB(0.904281, 0.319610, 0.388137),
RGB(0.908884, 0.324755, 0.385308),
RGB(0.913354, 0.330052, 0.382563),
RGB(0.917689, 0.335500, 0.379915),
RGB(0.921884, 0.341098, 0.377376),
RGB(0.925937, 0.346844, 0.374959),
RGB(0.929845, 0.352734, 0.372677),
RGB(0.933606, 0.358764, 0.370541),
RGB(0.937221, 0.364929, 0.368567),
RGB(0.940687, 0.371224, 0.366762),
RGB(0.944006, 0.377643, 0.365136),
RGB(0.947180, 0.384178, 0.363701),
RGB(0.950210, 0.390820, 0.362468),
RGB(0.953099, 0.397563, 0.361438),
RGB(0.955849, 0.404400, 0.360619),
RGB(0.958464, 0.411324, 0.360014),
RGB(0.960949, 0.418323, 0.359630),
RGB(0.963310, 0.425390, 0.359469),
RGB(0.965549, 0.432519, 0.359529),
RGB(0.967671, 0.439703, 0.359810),
RGB(0.969680, 0.446936, 0.360311),
RGB(0.971582, 0.454210, 0.361030),
RGB(0.973381, 0.461520, 0.361965),
RGB(0.975082, 0.468861, 0.363111),
RGB(0.976690, 0.476226, 0.364466),
RGB(0.978210, 0.483612, 0.366025),
RGB(0.979645, 0.491014, 0.367783),
RGB(0.981000, 0.498428, 0.369734),
RGB(0.982279, 0.505851, 0.371874),
RGB(0.983485, 0.513280, 0.374198),
RGB(0.984622, 0.520713, 0.376698),
RGB(0.985693, 0.528148, 0.379371),
RGB(0.986700, 0.535582, 0.382210),
RGB(0.987646, 0.543015, 0.385210),
RGB(0.988533, 0.550446, 0.388365),
RGB(0.989363, 0.557873, 0.391671),
RGB(0.990138, 0.565296, 0.395122),
RGB(0.990871, 0.572706, 0.398714),
RGB(0.991558, 0.580107, 0.402441),
RGB(0.992196, 0.587502, 0.406299),
RGB(0.992785, 0.594891, 0.410283),
RGB(0.993326, 0.602275, 0.414390),
RGB(0.993834, 0.609644, 0.418613),
RGB(0.994309, 0.616999, 0.422950),
RGB(0.994738, 0.624350, 0.427397),
RGB(0.995122, 0.631696, 0.431951),
RGB(0.995480, 0.639027, 0.436607),
RGB(0.995810, 0.646344, 0.441361),
RGB(0.996096, 0.653659, 0.446213),
RGB(0.996341, 0.660969, 0.451160),
RGB(0.996580, 0.668256, 0.456192),
RGB(0.996775, 0.675541, 0.461314),
RGB(0.996925, 0.682828, 0.466526),
RGB(0.997077, 0.690088, 0.471811),
RGB(0.997186, 0.697349, 0.477182),
RGB(0.997254, 0.704611, 0.482635),
RGB(0.997325, 0.711848, 0.488154),
RGB(0.997351, 0.719089, 0.493755),
RGB(0.997351, 0.726324, 0.499428),
RGB(0.997341, 0.733545, 0.505167),
RGB(0.997285, 0.740772, 0.510983),
RGB(0.997228, 0.747981, 0.516859),
RGB(0.997138, 0.755190, 0.522806),
RGB(0.997019, 0.762398, 0.528821),
RGB(0.996898, 0.769591, 0.534892),
RGB(0.996727, 0.776795, 0.541039),
RGB(0.996571, 0.783977, 0.547233),
RGB(0.996369, 0.791167, 0.553499),
RGB(0.996162, 0.798348, 0.559820),
RGB(0.995932, 0.805527, 0.566202),
RGB(0.995680, 0.812706, 0.572645),
RGB(0.995424, 0.819875, 0.579140),
RGB(0.995131, 0.827052, 0.585701),
RGB(0.994851, 0.834213, 0.592307),
RGB(0.994524, 0.841387, 0.598983),
RGB(0.994222, 0.848540, 0.605696),
RGB(0.993866, 0.855711, 0.612482),
RGB(0.993545, 0.862859, 0.619299),
RGB(0.993170, 0.870024, 0.626189),
RGB(0.992831, 0.877168, 0.633109),
RGB(0.992440, 0.884330, 0.640099),
RGB(0.992089, 0.891470, 0.647116),
RGB(0.991688, 0.898627, 0.654202),
RGB(0.991332, 0.905763, 0.661309),
RGB(0.990930, 0.912915, 0.668481),
RGB(0.990570, 0.920049, 0.675675),
RGB(0.990175, 0.927196, 0.682926),
RGB(0.989815, 0.934329, 0.690198),
RGB(0.989434, 0.941470, 0.697519),
RGB(0.989077, 0.948604, 0.704863),
RGB(0.988717, 0.955742, 0.712242),
RGB(0.988367, 0.962878, 0.719649),
RGB(0.988033, 0.970012, 0.727077),
RGB(0.987691, 0.977154, 0.734536),
RGB(0.987387, 0.984288, 0.742002),
RGB(0.987053, 0.991438, 0.749504)], "matplotlib", "perceptually uniform sequential")

loadcolorscheme(:inferno, [  RGB(0.001462, 0.000466, 0.013866),
RGB(0.002267, 0.001270, 0.018570),
RGB(0.003299, 0.002249, 0.024239),
RGB(0.004547, 0.003392, 0.030909),
RGB(0.006006, 0.004692, 0.038558),
RGB(0.007676, 0.006136, 0.046836),
RGB(0.009561, 0.007713, 0.055143),
RGB(0.011663, 0.009417, 0.063460),
RGB(0.013995, 0.011225, 0.071862),
RGB(0.016561, 0.013136, 0.080282),
RGB(0.019373, 0.015133, 0.088767),
RGB(0.022447, 0.017199, 0.097327),
RGB(0.025793, 0.019331, 0.105930),
RGB(0.029432, 0.021503, 0.114621),
RGB(0.033385, 0.023702, 0.123397),
RGB(0.037668, 0.025921, 0.132232),
RGB(0.042253, 0.028139, 0.141141),
RGB(0.046915, 0.030324, 0.150164),
RGB(0.051644, 0.032474, 0.159254),
RGB(0.056449, 0.034569, 0.168414),
RGB(0.061340, 0.036590, 0.177642),
RGB(0.066331, 0.038504, 0.186962),
RGB(0.071429, 0.040294, 0.196354),
RGB(0.076637, 0.041905, 0.205799),
RGB(0.081962, 0.043328, 0.215289),
RGB(0.087411, 0.044556, 0.224813),
RGB(0.092990, 0.045583, 0.234358),
RGB(0.098702, 0.046402, 0.243904),
RGB(0.104551, 0.047008, 0.253430),
RGB(0.110536, 0.047399, 0.262912),
RGB(0.116656, 0.047574, 0.272321),
RGB(0.122908, 0.047536, 0.281624),
RGB(0.129285, 0.047293, 0.290788),
RGB(0.135778, 0.046856, 0.299776),
RGB(0.142378, 0.046242, 0.308553),
RGB(0.149073, 0.045468, 0.317085),
RGB(0.155850, 0.044559, 0.325338),
RGB(0.162689, 0.043554, 0.333277),
RGB(0.169575, 0.042489, 0.340874),
RGB(0.176493, 0.041402, 0.348111),
RGB(0.183429, 0.040329, 0.354971),
RGB(0.190367, 0.039309, 0.361447),
RGB(0.197297, 0.038400, 0.367535),
RGB(0.204209, 0.037632, 0.373238),
RGB(0.211095, 0.037030, 0.378563),
RGB(0.217949, 0.036615, 0.383522),
RGB(0.224763, 0.036405, 0.388129),
RGB(0.231538, 0.036405, 0.392400),
RGB(0.238273, 0.036621, 0.396353),
RGB(0.244967, 0.037055, 0.400007),
RGB(0.251620, 0.037705, 0.403378),
RGB(0.258234, 0.038571, 0.406485),
RGB(0.264810, 0.039647, 0.409345),
RGB(0.271347, 0.040922, 0.411976),
RGB(0.277850, 0.042353, 0.414392),
RGB(0.284321, 0.043933, 0.416608),
RGB(0.290763, 0.045644, 0.418637),
RGB(0.297178, 0.047470, 0.420491),
RGB(0.303568, 0.049396, 0.422182),
RGB(0.309935, 0.051407, 0.423721),
RGB(0.316282, 0.053490, 0.425116),
RGB(0.322610, 0.055634, 0.426377),
RGB(0.328921, 0.057827, 0.427511),
RGB(0.335217, 0.060060, 0.428524),
RGB(0.341500, 0.062325, 0.429425),
RGB(0.347771, 0.064616, 0.430217),
RGB(0.354032, 0.066925, 0.430906),
RGB(0.360284, 0.069247, 0.431497),
RGB(0.366529, 0.071579, 0.431994),
RGB(0.372768, 0.073915, 0.432400),
RGB(0.379001, 0.076253, 0.432719),
RGB(0.385228, 0.078591, 0.432955),
RGB(0.391453, 0.080927, 0.433109),
RGB(0.397674, 0.083257, 0.433183),
RGB(0.403894, 0.085580, 0.433179),
RGB(0.410113, 0.087896, 0.433098),
RGB(0.416331, 0.090203, 0.432943),
RGB(0.422549, 0.092501, 0.432714),
RGB(0.428768, 0.094790, 0.432412),
RGB(0.434987, 0.097069, 0.432039),
RGB(0.441207, 0.099338, 0.431594),
RGB(0.447428, 0.101597, 0.431080),
RGB(0.453651, 0.103848, 0.430498),
RGB(0.459875, 0.106089, 0.429846),
RGB(0.466100, 0.108322, 0.429125),
RGB(0.472328, 0.110547, 0.428334),
RGB(0.478558, 0.112764, 0.427475),
RGB(0.484789, 0.114974, 0.426548),
RGB(0.491022, 0.117179, 0.425552),
RGB(0.497257, 0.119379, 0.424488),
RGB(0.503493, 0.121575, 0.423356),
RGB(0.509730, 0.123769, 0.422156),
RGB(0.515967, 0.125960, 0.420887),
RGB(0.522206, 0.128150, 0.419549),
RGB(0.528444, 0.130341, 0.418142),
RGB(0.534683, 0.132534, 0.416667),
RGB(0.540920, 0.134729, 0.415123),
RGB(0.547157, 0.136929, 0.413511),
RGB(0.553392, 0.139134, 0.411829),
RGB(0.559624, 0.141346, 0.410078),
RGB(0.565854, 0.143567, 0.408258),
RGB(0.572081, 0.145797, 0.406369),
RGB(0.578304, 0.148039, 0.404411),
RGB(0.584521, 0.150294, 0.402385),
RGB(0.590734, 0.152563, 0.400290),
RGB(0.596940, 0.154848, 0.398125),
RGB(0.603139, 0.157151, 0.395891),
RGB(0.609330, 0.159474, 0.393589),
RGB(0.615513, 0.161817, 0.391219),
RGB(0.621685, 0.164184, 0.388781),
RGB(0.627847, 0.166575, 0.386276),
RGB(0.633998, 0.168992, 0.383704),
RGB(0.640135, 0.171438, 0.381065),
RGB(0.646260, 0.173914, 0.378359),
RGB(0.652369, 0.176421, 0.375586),
RGB(0.658463, 0.178962, 0.372748),
RGB(0.664540, 0.181539, 0.369846),
RGB(0.670599, 0.184153, 0.366879),
RGB(0.676638, 0.186807, 0.363849),
RGB(0.682656, 0.189501, 0.360757),
RGB(0.688653, 0.192239, 0.357603),
RGB(0.694627, 0.195021, 0.354388),
RGB(0.700576, 0.197851, 0.351113),
RGB(0.706500, 0.200728, 0.347777),
RGB(0.712396, 0.203656, 0.344383),
RGB(0.718264, 0.206636, 0.340931),
RGB(0.724103, 0.209670, 0.337424),
RGB(0.729909, 0.212759, 0.333861),
RGB(0.735683, 0.215906, 0.330245),
RGB(0.741423, 0.219112, 0.326576),
RGB(0.747127, 0.222378, 0.322856),
RGB(0.752794, 0.225706, 0.319085),
RGB(0.758422, 0.229097, 0.315266),
RGB(0.764010, 0.232554, 0.311399),
RGB(0.769556, 0.236077, 0.307485),
RGB(0.775059, 0.239667, 0.303526),
RGB(0.780517, 0.243327, 0.299523),
RGB(0.785929, 0.247056, 0.295477),
RGB(0.791293, 0.250856, 0.291390),
RGB(0.796607, 0.254728, 0.287264),
RGB(0.801871, 0.258674, 0.283099),
RGB(0.807082, 0.262692, 0.278898),
RGB(0.812239, 0.266786, 0.274661),
RGB(0.817341, 0.270954, 0.270390),
RGB(0.822386, 0.275197, 0.266085),
RGB(0.827372, 0.279517, 0.261750),
RGB(0.832299, 0.283913, 0.257383),
RGB(0.837165, 0.288385, 0.252988),
RGB(0.841969, 0.292933, 0.248564),
RGB(0.846709, 0.297559, 0.244113),
RGB(0.851384, 0.302260, 0.239636),
RGB(0.855992, 0.307038, 0.235133),
RGB(0.860533, 0.311892, 0.230606),
RGB(0.865006, 0.316822, 0.226055),
RGB(0.869409, 0.321827, 0.221482),
RGB(0.873741, 0.326906, 0.216886),
RGB(0.878001, 0.332060, 0.212268),
RGB(0.882188, 0.337287, 0.207628),
RGB(0.886302, 0.342586, 0.202968),
RGB(0.890341, 0.347957, 0.198286),
RGB(0.894305, 0.353399, 0.193584),
RGB(0.898192, 0.358911, 0.188860),
RGB(0.902003, 0.364492, 0.184116),
RGB(0.905735, 0.370140, 0.179350),
RGB(0.909390, 0.375856, 0.174563),
RGB(0.912966, 0.381636, 0.169755),
RGB(0.916462, 0.387481, 0.164924),
RGB(0.919879, 0.393389, 0.160070),
RGB(0.923215, 0.399359, 0.155193),
RGB(0.926470, 0.405389, 0.150292),
RGB(0.929644, 0.411479, 0.145367),
RGB(0.932737, 0.417627, 0.140417),
RGB(0.935747, 0.423831, 0.135440),
RGB(0.938675, 0.430091, 0.130438),
RGB(0.941521, 0.436405, 0.125409),
RGB(0.944285, 0.442772, 0.120354),
RGB(0.946965, 0.449191, 0.115272),
RGB(0.949562, 0.455660, 0.110164),
RGB(0.952075, 0.462178, 0.105031),
RGB(0.954506, 0.468744, 0.099874),
RGB(0.956852, 0.475356, 0.094695),
RGB(0.959114, 0.482014, 0.089499),
RGB(0.961293, 0.488716, 0.084289),
RGB(0.963387, 0.495462, 0.079073),
RGB(0.965397, 0.502249, 0.073859),
RGB(0.967322, 0.509078, 0.068659),
RGB(0.969163, 0.515946, 0.063488),
RGB(0.970919, 0.522853, 0.058367),
RGB(0.972590, 0.529798, 0.053324),
RGB(0.974176, 0.536780, 0.048392),
RGB(0.975677, 0.543798, 0.043618),
RGB(0.977092, 0.550850, 0.039050),
RGB(0.978422, 0.557937, 0.034931),
RGB(0.979666, 0.565057, 0.031409),
RGB(0.980824, 0.572209, 0.028508),
RGB(0.981895, 0.579392, 0.026250),
RGB(0.982881, 0.586606, 0.024661),
RGB(0.983779, 0.593849, 0.023770),
RGB(0.984591, 0.601122, 0.023606),
RGB(0.985315, 0.608422, 0.024202),
RGB(0.985952, 0.615750, 0.025592),
RGB(0.986502, 0.623105, 0.027814),
RGB(0.986964, 0.630485, 0.030908),
RGB(0.987337, 0.637890, 0.034916),
RGB(0.987622, 0.645320, 0.039886),
RGB(0.987819, 0.652773, 0.045581),
RGB(0.987926, 0.660250, 0.051750),
RGB(0.987945, 0.667748, 0.058329),
RGB(0.987874, 0.675267, 0.065257),
RGB(0.987714, 0.682807, 0.072489),
RGB(0.987464, 0.690366, 0.079990),
RGB(0.987124, 0.697944, 0.087731),
RGB(0.986694, 0.705540, 0.095694),
RGB(0.986175, 0.713153, 0.103863),
RGB(0.985566, 0.720782, 0.112229),
RGB(0.984865, 0.728427, 0.120785),
RGB(0.984075, 0.736087, 0.129527),
RGB(0.983196, 0.743758, 0.138453),
RGB(0.982228, 0.751442, 0.147565),
RGB(0.981173, 0.759135, 0.156863),
RGB(0.980032, 0.766837, 0.166353),
RGB(0.978806, 0.774545, 0.176037),
RGB(0.977497, 0.782258, 0.185923),
RGB(0.976108, 0.789974, 0.196018),
RGB(0.974638, 0.797692, 0.206332),
RGB(0.973088, 0.805409, 0.216877),
RGB(0.971468, 0.813122, 0.227658),
RGB(0.969783, 0.820825, 0.238686),
RGB(0.968041, 0.828515, 0.249972),
RGB(0.966243, 0.836191, 0.261534),
RGB(0.964394, 0.843848, 0.273391),
RGB(0.962517, 0.851476, 0.285546),
RGB(0.960626, 0.859069, 0.298010),
RGB(0.958720, 0.866624, 0.310820),
RGB(0.956834, 0.874129, 0.323974),
RGB(0.954997, 0.881569, 0.337475),
RGB(0.953215, 0.888942, 0.351369),
RGB(0.951546, 0.896226, 0.365627),
RGB(0.950018, 0.903409, 0.380271),
RGB(0.948683, 0.910473, 0.395289),
RGB(0.947594, 0.917399, 0.410665),
RGB(0.946809, 0.924168, 0.426373),
RGB(0.946392, 0.930761, 0.442367),
RGB(0.946403, 0.937159, 0.458592),
RGB(0.946903, 0.943348, 0.474970),
RGB(0.947937, 0.949318, 0.491426),
RGB(0.949545, 0.955063, 0.507860),
RGB(0.951740, 0.960587, 0.524203),
RGB(0.954529, 0.965896, 0.540361),
RGB(0.957896, 0.971003, 0.556275),
RGB(0.961812, 0.975924, 0.571925),
RGB(0.966249, 0.980678, 0.587206),
RGB(0.971162, 0.985282, 0.602154),
RGB(0.976511, 0.989753, 0.616760),
RGB(0.982257, 0.994109, 0.631017),
RGB(0.988362, 0.998364, 0.644924)], "matplotlib", "perceptually uniform sequential")

loadcolorscheme(:plasma, [  RGB(0.050383, 0.029803, 0.527975),
RGB(0.063536, 0.028426, 0.533124),
RGB(0.075353, 0.027206, 0.538007),
RGB(0.086222, 0.026125, 0.542658),
RGB(0.096379, 0.025165, 0.547103),
RGB(0.105980, 0.024309, 0.551368),
RGB(0.115124, 0.023556, 0.555468),
RGB(0.123903, 0.022878, 0.559423),
RGB(0.132381, 0.022258, 0.563250),
RGB(0.140603, 0.021687, 0.566959),
RGB(0.148607, 0.021154, 0.570562),
RGB(0.156421, 0.020651, 0.574065),
RGB(0.164070, 0.020171, 0.577478),
RGB(0.171574, 0.019706, 0.580806),
RGB(0.178950, 0.019252, 0.584054),
RGB(0.186213, 0.018803, 0.587228),
RGB(0.193374, 0.018354, 0.590330),
RGB(0.200445, 0.017902, 0.593364),
RGB(0.207435, 0.017442, 0.596333),
RGB(0.214350, 0.016973, 0.599239),
RGB(0.221197, 0.016497, 0.602083),
RGB(0.227983, 0.016007, 0.604867),
RGB(0.234715, 0.015502, 0.607592),
RGB(0.241396, 0.014979, 0.610259),
RGB(0.248032, 0.014439, 0.612868),
RGB(0.254627, 0.013882, 0.615419),
RGB(0.261183, 0.013308, 0.617911),
RGB(0.267703, 0.012716, 0.620346),
RGB(0.274191, 0.012109, 0.622722),
RGB(0.280648, 0.011488, 0.625038),
RGB(0.287076, 0.010855, 0.627295),
RGB(0.293478, 0.010213, 0.629490),
RGB(0.299855, 0.009561, 0.631624),
RGB(0.306210, 0.008902, 0.633694),
RGB(0.312543, 0.008239, 0.635700),
RGB(0.318856, 0.007576, 0.637640),
RGB(0.325150, 0.006915, 0.639512),
RGB(0.331426, 0.006261, 0.641316),
RGB(0.337683, 0.005618, 0.643049),
RGB(0.343925, 0.004991, 0.644710),
RGB(0.350150, 0.004382, 0.646298),
RGB(0.356359, 0.003798, 0.647810),
RGB(0.362553, 0.003243, 0.649245),
RGB(0.368733, 0.002724, 0.650601),
RGB(0.374897, 0.002245, 0.651876),
RGB(0.381047, 0.001814, 0.653068),
RGB(0.387183, 0.001434, 0.654177),
RGB(0.393304, 0.001114, 0.655199),
RGB(0.399411, 0.000859, 0.656133),
RGB(0.405503, 0.000678, 0.656977),
RGB(0.411580, 0.000577, 0.657730),
RGB(0.417642, 0.000564, 0.658390),
RGB(0.423689, 0.000646, 0.658956),
RGB(0.429719, 0.000831, 0.659425),
RGB(0.435734, 0.001127, 0.659797),
RGB(0.441732, 0.001540, 0.660069),
RGB(0.447714, 0.002080, 0.660240),
RGB(0.453677, 0.002755, 0.660310),
RGB(0.459623, 0.003574, 0.660277),
RGB(0.465550, 0.004545, 0.660139),
RGB(0.471457, 0.005678, 0.659897),
RGB(0.477344, 0.006980, 0.659549),
RGB(0.483210, 0.008460, 0.659095),
RGB(0.489055, 0.010127, 0.658534),
RGB(0.494877, 0.011990, 0.657865),
RGB(0.500678, 0.014055, 0.657088),
RGB(0.506454, 0.016333, 0.656202),
RGB(0.512206, 0.018833, 0.655209),
RGB(0.517933, 0.021563, 0.654109),
RGB(0.523633, 0.024532, 0.652901),
RGB(0.529306, 0.027747, 0.651586),
RGB(0.534952, 0.031217, 0.650165),
RGB(0.540570, 0.034950, 0.648640),
RGB(0.546157, 0.038954, 0.647010),
RGB(0.551715, 0.043136, 0.645277),
RGB(0.557243, 0.047331, 0.643443),
RGB(0.562738, 0.051545, 0.641509),
RGB(0.568201, 0.055778, 0.639477),
RGB(0.573632, 0.060028, 0.637349),
RGB(0.579029, 0.064296, 0.635126),
RGB(0.584391, 0.068579, 0.632812),
RGB(0.589719, 0.072878, 0.630408),
RGB(0.595011, 0.077190, 0.627917),
RGB(0.600266, 0.081516, 0.625342),
RGB(0.605485, 0.085854, 0.622686),
RGB(0.610667, 0.090204, 0.619951),
RGB(0.615812, 0.094564, 0.617140),
RGB(0.620919, 0.098934, 0.614257),
RGB(0.625987, 0.103312, 0.611305),
RGB(0.631017, 0.107699, 0.608287),
RGB(0.636008, 0.112092, 0.605205),
RGB(0.640959, 0.116492, 0.602065),
RGB(0.645872, 0.120898, 0.598867),
RGB(0.650746, 0.125309, 0.595617),
RGB(0.655580, 0.129725, 0.592317),
RGB(0.660374, 0.134144, 0.588971),
RGB(0.665129, 0.138566, 0.585582),
RGB(0.669845, 0.142992, 0.582154),
RGB(0.674522, 0.147419, 0.578688),
RGB(0.679160, 0.151848, 0.575189),
RGB(0.683758, 0.156278, 0.571660),
RGB(0.688318, 0.160709, 0.568103),
RGB(0.692840, 0.165141, 0.564522),
RGB(0.697324, 0.169573, 0.560919),
RGB(0.701769, 0.174005, 0.557296),
RGB(0.706178, 0.178437, 0.553657),
RGB(0.710549, 0.182868, 0.550004),
RGB(0.714883, 0.187299, 0.546338),
RGB(0.719181, 0.191729, 0.542663),
RGB(0.723444, 0.196158, 0.538981),
RGB(0.727670, 0.200586, 0.535293),
RGB(0.731862, 0.205013, 0.531601),
RGB(0.736019, 0.209439, 0.527908),
RGB(0.740143, 0.213864, 0.524216),
RGB(0.744232, 0.218288, 0.520524),
RGB(0.748289, 0.222711, 0.516834),
RGB(0.752312, 0.227133, 0.513149),
RGB(0.756304, 0.231555, 0.509468),
RGB(0.760264, 0.235976, 0.505794),
RGB(0.764193, 0.240396, 0.502126),
RGB(0.768090, 0.244817, 0.498465),
RGB(0.771958, 0.249237, 0.494813),
RGB(0.775796, 0.253658, 0.491171),
RGB(0.779604, 0.258078, 0.487539),
RGB(0.783383, 0.262500, 0.483918),
RGB(0.787133, 0.266922, 0.480307),
RGB(0.790855, 0.271345, 0.476706),
RGB(0.794549, 0.275770, 0.473117),
RGB(0.798216, 0.280197, 0.469538),
RGB(0.801855, 0.284626, 0.465971),
RGB(0.805467, 0.289057, 0.462415),
RGB(0.809052, 0.293491, 0.458870),
RGB(0.812612, 0.297928, 0.455338),
RGB(0.816144, 0.302368, 0.451816),
RGB(0.819651, 0.306812, 0.448306),
RGB(0.823132, 0.311261, 0.444806),
RGB(0.826588, 0.315714, 0.441316),
RGB(0.830018, 0.320172, 0.437836),
RGB(0.833422, 0.324635, 0.434366),
RGB(0.836801, 0.329105, 0.430905),
RGB(0.840155, 0.333580, 0.427455),
RGB(0.843484, 0.338062, 0.424013),
RGB(0.846788, 0.342551, 0.420579),
RGB(0.850066, 0.347048, 0.417153),
RGB(0.853319, 0.351553, 0.413734),
RGB(0.856547, 0.356066, 0.410322),
RGB(0.859750, 0.360588, 0.406917),
RGB(0.862927, 0.365119, 0.403519),
RGB(0.866078, 0.369660, 0.400126),
RGB(0.869203, 0.374212, 0.396738),
RGB(0.872303, 0.378774, 0.393355),
RGB(0.875376, 0.383347, 0.389976),
RGB(0.878423, 0.387932, 0.386600),
RGB(0.881443, 0.392529, 0.383229),
RGB(0.884436, 0.397139, 0.379860),
RGB(0.887402, 0.401762, 0.376494),
RGB(0.890340, 0.406398, 0.373130),
RGB(0.893250, 0.411048, 0.369768),
RGB(0.896131, 0.415712, 0.366407),
RGB(0.898984, 0.420392, 0.363047),
RGB(0.901807, 0.425087, 0.359688),
RGB(0.904601, 0.429797, 0.356329),
RGB(0.907365, 0.434524, 0.352970),
RGB(0.910098, 0.439268, 0.349610),
RGB(0.912800, 0.444029, 0.346251),
RGB(0.915471, 0.448807, 0.342890),
RGB(0.918109, 0.453603, 0.339529),
RGB(0.920714, 0.458417, 0.336166),
RGB(0.923287, 0.463251, 0.332801),
RGB(0.925825, 0.468103, 0.329435),
RGB(0.928329, 0.472975, 0.326067),
RGB(0.930798, 0.477867, 0.322697),
RGB(0.933232, 0.482780, 0.319325),
RGB(0.935630, 0.487712, 0.315952),
RGB(0.937990, 0.492667, 0.312575),
RGB(0.940313, 0.497642, 0.309197),
RGB(0.942598, 0.502639, 0.305816),
RGB(0.944844, 0.507658, 0.302433),
RGB(0.947051, 0.512699, 0.299049),
RGB(0.949217, 0.517763, 0.295662),
RGB(0.951344, 0.522850, 0.292275),
RGB(0.953428, 0.527960, 0.288883),
RGB(0.955470, 0.533093, 0.285490),
RGB(0.957469, 0.538250, 0.282096),
RGB(0.959424, 0.543431, 0.278701),
RGB(0.961336, 0.548636, 0.275305),
RGB(0.963203, 0.553865, 0.271909),
RGB(0.965024, 0.559118, 0.268513),
RGB(0.966798, 0.564396, 0.265118),
RGB(0.968526, 0.569700, 0.261721),
RGB(0.970205, 0.575028, 0.258325),
RGB(0.971835, 0.580382, 0.254931),
RGB(0.973416, 0.585761, 0.251540),
RGB(0.974947, 0.591165, 0.248151),
RGB(0.976428, 0.596595, 0.244767),
RGB(0.977856, 0.602051, 0.241387),
RGB(0.979233, 0.607532, 0.238013),
RGB(0.980556, 0.613039, 0.234646),
RGB(0.981826, 0.618572, 0.231287),
RGB(0.983041, 0.624131, 0.227937),
RGB(0.984199, 0.629718, 0.224595),
RGB(0.985301, 0.635330, 0.221265),
RGB(0.986345, 0.640969, 0.217948),
RGB(0.987332, 0.646633, 0.214648),
RGB(0.988260, 0.652325, 0.211364),
RGB(0.989128, 0.658043, 0.208100),
RGB(0.989935, 0.663787, 0.204859),
RGB(0.990681, 0.669558, 0.201642),
RGB(0.991365, 0.675355, 0.198453),
RGB(0.991985, 0.681179, 0.195295),
RGB(0.992541, 0.687030, 0.192170),
RGB(0.993032, 0.692907, 0.189084),
RGB(0.993456, 0.698810, 0.186041),
RGB(0.993814, 0.704741, 0.183043),
RGB(0.994103, 0.710698, 0.180097),
RGB(0.994324, 0.716681, 0.177208),
RGB(0.994474, 0.722691, 0.174381),
RGB(0.994553, 0.728728, 0.171622),
RGB(0.994561, 0.734791, 0.168938),
RGB(0.994495, 0.740880, 0.166335),
RGB(0.994355, 0.746995, 0.163821),
RGB(0.994141, 0.753137, 0.161404),
RGB(0.993851, 0.759304, 0.159092),
RGB(0.993482, 0.765499, 0.156891),
RGB(0.993033, 0.771720, 0.154808),
RGB(0.992505, 0.777967, 0.152855),
RGB(0.991897, 0.784239, 0.151042),
RGB(0.991209, 0.790537, 0.149377),
RGB(0.990439, 0.796859, 0.147870),
RGB(0.989587, 0.803205, 0.146529),
RGB(0.988648, 0.809579, 0.145357),
RGB(0.987621, 0.815978, 0.144363),
RGB(0.986509, 0.822401, 0.143557),
RGB(0.985314, 0.828846, 0.142945),
RGB(0.984031, 0.835315, 0.142528),
RGB(0.982653, 0.841812, 0.142303),
RGB(0.981190, 0.848329, 0.142279),
RGB(0.979644, 0.854866, 0.142453),
RGB(0.977995, 0.861432, 0.142808),
RGB(0.976265, 0.868016, 0.143351),
RGB(0.974443, 0.874622, 0.144061),
RGB(0.972530, 0.881250, 0.144923),
RGB(0.970533, 0.887896, 0.145919),
RGB(0.968443, 0.894564, 0.147014),
RGB(0.966271, 0.901249, 0.148180),
RGB(0.964021, 0.907950, 0.149370),
RGB(0.961681, 0.914672, 0.150520),
RGB(0.959276, 0.921407, 0.151566),
RGB(0.956808, 0.928152, 0.152409),
RGB(0.954287, 0.934908, 0.152921),
RGB(0.951726, 0.941671, 0.152925),
RGB(0.949151, 0.948435, 0.152178),
RGB(0.946602, 0.955190, 0.150328),
RGB(0.944152, 0.961916, 0.146861),
RGB(0.941896, 0.968590, 0.140956),
RGB(0.940015, 0.975158, 0.131326)], "matplotlib", "perceptually uniform sequential")

loadcolorscheme(:viridis, [  RGB(0.267004, 0.004874, 0.329415),
RGB(0.268510, 0.009605, 0.335427),
RGB(0.269944, 0.014625, 0.341379),
RGB(0.271305, 0.019942, 0.347269),
RGB(0.272594, 0.025563, 0.353093),
RGB(0.273809, 0.031497, 0.358853),
RGB(0.274952, 0.037752, 0.364543),
RGB(0.276022, 0.044167, 0.370164),
RGB(0.277018, 0.050344, 0.375715),
RGB(0.277941, 0.056324, 0.381191),
RGB(0.278791, 0.062145, 0.386592),
RGB(0.279566, 0.067836, 0.391917),
RGB(0.280267, 0.073417, 0.397163),
RGB(0.280894, 0.078907, 0.402329),
RGB(0.281446, 0.084320, 0.407414),
RGB(0.281924, 0.089666, 0.412415),
RGB(0.282327, 0.094955, 0.417331),
RGB(0.282656, 0.100196, 0.422160),
RGB(0.282910, 0.105393, 0.426902),
RGB(0.283091, 0.110553, 0.431554),
RGB(0.283197, 0.115680, 0.436115),
RGB(0.283229, 0.120777, 0.440584),
RGB(0.283187, 0.125848, 0.444960),
RGB(0.283072, 0.130895, 0.449241),
RGB(0.282884, 0.135920, 0.453427),
RGB(0.282623, 0.140926, 0.457517),
RGB(0.282290, 0.145912, 0.461510),
RGB(0.281887, 0.150881, 0.465405),
RGB(0.281412, 0.155834, 0.469201),
RGB(0.280868, 0.160771, 0.472899),
RGB(0.280255, 0.165693, 0.476498),
RGB(0.279574, 0.170599, 0.479997),
RGB(0.278826, 0.175490, 0.483397),
RGB(0.278012, 0.180367, 0.486697),
RGB(0.277134, 0.185228, 0.489898),
RGB(0.276194, 0.190074, 0.493001),
RGB(0.275191, 0.194905, 0.496005),
RGB(0.274128, 0.199721, 0.498911),
RGB(0.273006, 0.204520, 0.501721),
RGB(0.271828, 0.209303, 0.504434),
RGB(0.270595, 0.214069, 0.507052),
RGB(0.269308, 0.218818, 0.509577),
RGB(0.267968, 0.223549, 0.512008),
RGB(0.266580, 0.228262, 0.514349),
RGB(0.265145, 0.232956, 0.516599),
RGB(0.263663, 0.237631, 0.518762),
RGB(0.262138, 0.242286, 0.520837),
RGB(0.260571, 0.246922, 0.522828),
RGB(0.258965, 0.251537, 0.524736),
RGB(0.257322, 0.256130, 0.526563),
RGB(0.255645, 0.260703, 0.528312),
RGB(0.253935, 0.265254, 0.529983),
RGB(0.252194, 0.269783, 0.531579),
RGB(0.250425, 0.274290, 0.533103),
RGB(0.248629, 0.278775, 0.534556),
RGB(0.246811, 0.283237, 0.535941),
RGB(0.244972, 0.287675, 0.537260),
RGB(0.243113, 0.292092, 0.538516),
RGB(0.241237, 0.296485, 0.539709),
RGB(0.239346, 0.300855, 0.540844),
RGB(0.237441, 0.305202, 0.541921),
RGB(0.235526, 0.309527, 0.542944),
RGB(0.233603, 0.313828, 0.543914),
RGB(0.231674, 0.318106, 0.544834),
RGB(0.229739, 0.322361, 0.545706),
RGB(0.227802, 0.326594, 0.546532),
RGB(0.225863, 0.330805, 0.547314),
RGB(0.223925, 0.334994, 0.548053),
RGB(0.221989, 0.339161, 0.548752),
RGB(0.220057, 0.343307, 0.549413),
RGB(0.218130, 0.347432, 0.550038),
RGB(0.216210, 0.351535, 0.550627),
RGB(0.214298, 0.355619, 0.551184),
RGB(0.212395, 0.359683, 0.551710),
RGB(0.210503, 0.363727, 0.552206),
RGB(0.208623, 0.367752, 0.552675),
RGB(0.206756, 0.371758, 0.553117),
RGB(0.204903, 0.375746, 0.553533),
RGB(0.203063, 0.379716, 0.553925),
RGB(0.201239, 0.383670, 0.554294),
RGB(0.199430, 0.387607, 0.554642),
RGB(0.197636, 0.391528, 0.554969),
RGB(0.195860, 0.395433, 0.555276),
RGB(0.194100, 0.399323, 0.555565),
RGB(0.192357, 0.403199, 0.555836),
RGB(0.190631, 0.407061, 0.556089),
RGB(0.188923, 0.410910, 0.556326),
RGB(0.187231, 0.414746, 0.556547),
RGB(0.185556, 0.418570, 0.556753),
RGB(0.183898, 0.422383, 0.556944),
RGB(0.182256, 0.426184, 0.557120),
RGB(0.180629, 0.429975, 0.557282),
RGB(0.179019, 0.433756, 0.557430),
RGB(0.177423, 0.437527, 0.557565),
RGB(0.175841, 0.441290, 0.557685),
RGB(0.174274, 0.445044, 0.557792),
RGB(0.172719, 0.448791, 0.557885),
RGB(0.171176, 0.452530, 0.557965),
RGB(0.169646, 0.456262, 0.558030),
RGB(0.168126, 0.459988, 0.558082),
RGB(0.166617, 0.463708, 0.558119),
RGB(0.165117, 0.467423, 0.558141),
RGB(0.163625, 0.471133, 0.558148),
RGB(0.162142, 0.474838, 0.558140),
RGB(0.160665, 0.478540, 0.558115),
RGB(0.159194, 0.482237, 0.558073),
RGB(0.157729, 0.485932, 0.558013),
RGB(0.156270, 0.489624, 0.557936),
RGB(0.154815, 0.493313, 0.557840),
RGB(0.153364, 0.497000, 0.557724),
RGB(0.151918, 0.500685, 0.557587),
RGB(0.150476, 0.504369, 0.557430),
RGB(0.149039, 0.508051, 0.557250),
RGB(0.147607, 0.511733, 0.557049),
RGB(0.146180, 0.515413, 0.556823),
RGB(0.144759, 0.519093, 0.556572),
RGB(0.143343, 0.522773, 0.556295),
RGB(0.141935, 0.526453, 0.555991),
RGB(0.140536, 0.530132, 0.555659),
RGB(0.139147, 0.533812, 0.555298),
RGB(0.137770, 0.537492, 0.554906),
RGB(0.136408, 0.541173, 0.554483),
RGB(0.135066, 0.544853, 0.554029),
RGB(0.133743, 0.548535, 0.553541),
RGB(0.132444, 0.552216, 0.553018),
RGB(0.131172, 0.555899, 0.552459),
RGB(0.129933, 0.559582, 0.551864),
RGB(0.128729, 0.563265, 0.551229),
RGB(0.127568, 0.566949, 0.550556),
RGB(0.126453, 0.570633, 0.549841),
RGB(0.125394, 0.574318, 0.549086),
RGB(0.124395, 0.578002, 0.548287),
RGB(0.123463, 0.581687, 0.547445),
RGB(0.122606, 0.585371, 0.546557),
RGB(0.121831, 0.589055, 0.545623),
RGB(0.121148, 0.592739, 0.544641),
RGB(0.120565, 0.596422, 0.543611),
RGB(0.120092, 0.600104, 0.542530),
RGB(0.119738, 0.603785, 0.541400),
RGB(0.119512, 0.607464, 0.540218),
RGB(0.119423, 0.611141, 0.538982),
RGB(0.119483, 0.614817, 0.537692),
RGB(0.119699, 0.618490, 0.536347),
RGB(0.120081, 0.622161, 0.534946),
RGB(0.120638, 0.625828, 0.533488),
RGB(0.121380, 0.629492, 0.531973),
RGB(0.122312, 0.633153, 0.530398),
RGB(0.123444, 0.636809, 0.528763),
RGB(0.124780, 0.640461, 0.527068),
RGB(0.126326, 0.644107, 0.525311),
RGB(0.128087, 0.647749, 0.523491),
RGB(0.130067, 0.651384, 0.521608),
RGB(0.132268, 0.655014, 0.519661),
RGB(0.134692, 0.658636, 0.517649),
RGB(0.137339, 0.662252, 0.515571),
RGB(0.140210, 0.665859, 0.513427),
RGB(0.143303, 0.669459, 0.511215),
RGB(0.146616, 0.673050, 0.508936),
RGB(0.150148, 0.676631, 0.506589),
RGB(0.153894, 0.680203, 0.504172),
RGB(0.157851, 0.683765, 0.501686),
RGB(0.162016, 0.687316, 0.499129),
RGB(0.166383, 0.690856, 0.496502),
RGB(0.170948, 0.694384, 0.493803),
RGB(0.175707, 0.697900, 0.491033),
RGB(0.180653, 0.701402, 0.488189),
RGB(0.185783, 0.704891, 0.485273),
RGB(0.191090, 0.708366, 0.482284),
RGB(0.196571, 0.711827, 0.479221),
RGB(0.202219, 0.715272, 0.476084),
RGB(0.208030, 0.718701, 0.472873),
RGB(0.214000, 0.722114, 0.469588),
RGB(0.220124, 0.725509, 0.466226),
RGB(0.226397, 0.728888, 0.462789),
RGB(0.232815, 0.732247, 0.459277),
RGB(0.239374, 0.735588, 0.455688),
RGB(0.246070, 0.738910, 0.452024),
RGB(0.252899, 0.742211, 0.448284),
RGB(0.259857, 0.745492, 0.444467),
RGB(0.266941, 0.748751, 0.440573),
RGB(0.274149, 0.751988, 0.436601),
RGB(0.281477, 0.755203, 0.432552),
RGB(0.288921, 0.758394, 0.428426),
RGB(0.296479, 0.761561, 0.424223),
RGB(0.304148, 0.764704, 0.419943),
RGB(0.311925, 0.767822, 0.415586),
RGB(0.319809, 0.770914, 0.411152),
RGB(0.327796, 0.773980, 0.406640),
RGB(0.335885, 0.777018, 0.402049),
RGB(0.344074, 0.780029, 0.397381),
RGB(0.352360, 0.783011, 0.392636),
RGB(0.360741, 0.785964, 0.387814),
RGB(0.369214, 0.788888, 0.382914),
RGB(0.377779, 0.791781, 0.377939),
RGB(0.386433, 0.794644, 0.372886),
RGB(0.395174, 0.797475, 0.367757),
RGB(0.404001, 0.800275, 0.362552),
RGB(0.412913, 0.803041, 0.357269),
RGB(0.421908, 0.805774, 0.351910),
RGB(0.430983, 0.808473, 0.346476),
RGB(0.440137, 0.811138, 0.340967),
RGB(0.449368, 0.813768, 0.335384),
RGB(0.458674, 0.816363, 0.329727),
RGB(0.468053, 0.818921, 0.323998),
RGB(0.477504, 0.821444, 0.318195),
RGB(0.487026, 0.823929, 0.312321),
RGB(0.496615, 0.826376, 0.306377),
RGB(0.506271, 0.828786, 0.300362),
RGB(0.515992, 0.831158, 0.294279),
RGB(0.525776, 0.833491, 0.288127),
RGB(0.535621, 0.835785, 0.281908),
RGB(0.545524, 0.838039, 0.275626),
RGB(0.555484, 0.840254, 0.269281),
RGB(0.565498, 0.842430, 0.262877),
RGB(0.575563, 0.844566, 0.256415),
RGB(0.585678, 0.846661, 0.249897),
RGB(0.595839, 0.848717, 0.243329),
RGB(0.606045, 0.850733, 0.236712),
RGB(0.616293, 0.852709, 0.230052),
RGB(0.626579, 0.854645, 0.223353),
RGB(0.636902, 0.856542, 0.216620),
RGB(0.647257, 0.858400, 0.209861),
RGB(0.657642, 0.860219, 0.203082),
RGB(0.668054, 0.861999, 0.196293),
RGB(0.678489, 0.863742, 0.189503),
RGB(0.688944, 0.865448, 0.182725),
RGB(0.699415, 0.867117, 0.175971),
RGB(0.709898, 0.868751, 0.169257),
RGB(0.720391, 0.870350, 0.162603),
RGB(0.730889, 0.871916, 0.156029),
RGB(0.741388, 0.873449, 0.149561),
RGB(0.751884, 0.874951, 0.143228),
RGB(0.762373, 0.876424, 0.137064),
RGB(0.772852, 0.877868, 0.131109),
RGB(0.783315, 0.879285, 0.125405),
RGB(0.793760, 0.880678, 0.120005),
RGB(0.804182, 0.882046, 0.114965),
RGB(0.814576, 0.883393, 0.110347),
RGB(0.824940, 0.884720, 0.106217),
RGB(0.835270, 0.886029, 0.102646),
RGB(0.845561, 0.887322, 0.099702),
RGB(0.855810, 0.888601, 0.097452),
RGB(0.866013, 0.889868, 0.095953),
RGB(0.876168, 0.891125, 0.095250),
RGB(0.886271, 0.892374, 0.095374),
RGB(0.896320, 0.893616, 0.096335),
RGB(0.906311, 0.894855, 0.098125),
RGB(0.916242, 0.896091, 0.100717),
RGB(0.926106, 0.897330, 0.104071),
RGB(0.935904, 0.898570, 0.108131),
RGB(0.945636, 0.899815, 0.112838),
RGB(0.955300, 0.901065, 0.118128),
RGB(0.964894, 0.902323, 0.123941),
RGB(0.974417, 0.903590, 0.130215),
RGB(0.983868, 0.904867, 0.136897),
RGB(0.993248, 0.906157, 0.143936)], "matplotlib", "perceptually uniform sequential")

loadcolorscheme(:twilight, [
Colors.RGB(0.88575,0.85001,0.88797),
Colors.RGB(0.88379,0.85073,0.88723),
Colors.RGB(0.88172,0.85128,0.88638),
Colors.RGB(0.87954,0.85166,0.88541),
Colors.RGB(0.87725,0.85187,0.88434),
Colors.RGB(0.87485,0.85192,0.88317),
Colors.RGB(0.87233,0.8518,0.8819),
Colors.RGB(0.8697,0.85152,0.88054),
Colors.RGB(0.86696,0.85109,0.8791),
Colors.RGB(0.86409,0.8505,0.87758),
Colors.RGB(0.8611,0.84977,0.87599),
Colors.RGB(0.85798,0.84889,0.87434),
Colors.RGB(0.85473,0.84787,0.87263),
Colors.RGB(0.85134,0.84673,0.87086),
Colors.RGB(0.84781,0.84545,0.86904),
Colors.RGB(0.84413,0.84406,0.86717),
Colors.RGB(0.8403,0.84256,0.86525),
Colors.RGB(0.83634,0.84095,0.86329),
Colors.RGB(0.83223,0.83923,0.86128),
Colors.RGB(0.82797,0.83743,0.85922),
Colors.RGB(0.82357,0.83552,0.85713),
Colors.RGB(0.81905,0.83354,0.855),
Colors.RGB(0.81439,0.83147,0.85284),
Colors.RGB(0.8096,0.82932,0.85064),
Colors.RGB(0.80469,0.8271,0.84842),
Colors.RGB(0.79967,0.82481,0.84618),
Colors.RGB(0.79454,0.82245,0.84392),
Colors.RGB(0.78931,0.82003,0.84165),
Colors.RGB(0.78399,0.81755,0.83937),
Colors.RGB(0.77858,0.81502,0.83708),
Colors.RGB(0.77308,0.81244,0.8348),
Colors.RGB(0.76751,0.8098,0.83253),
Colors.RGB(0.76187,0.80712,0.83027),
Colors.RGB(0.75616,0.80439,0.82802),
Colors.RGB(0.7504,0.80163,0.8258),
Colors.RGB(0.74459,0.79882,0.8236),
Colors.RGB(0.73874,0.79598,0.82143),
Colors.RGB(0.73285,0.7931,0.81929),
Colors.RGB(0.72692,0.79018,0.81719),
Colors.RGB(0.72097,0.78724,0.81513),
Colors.RGB(0.715,0.78426,0.81311),
Colors.RGB(0.70902,0.78126,0.81114),
Colors.RGB(0.70303,0.77823,0.80921),
Colors.RGB(0.69704,0.77517,0.80732),
Colors.RGB(0.69105,0.77209,0.80549),
Colors.RGB(0.68506,0.76898,0.8037),
Colors.RGB(0.6791,0.76584,0.80196),
Colors.RGB(0.67314,0.76269,0.80028),
Colors.RGB(0.66721,0.75951,0.79864),
Colors.RGB(0.66131,0.75631,0.79705),
Colors.RGB(0.65544,0.75309,0.7955),
Colors.RGB(0.6496,0.74985,0.79401),
Colors.RGB(0.64379,0.74659,0.79256),
Colors.RGB(0.63803,0.74331,0.79115),
Colors.RGB(0.6323,0.74002,0.78979),
Colors.RGB(0.62662,0.7367,0.78847),
Colors.RGB(0.62099,0.73337,0.78719),
Colors.RGB(0.61541,0.73002,0.78595),
Colors.RGB(0.60988,0.72665,0.78475),
Colors.RGB(0.60439,0.72327,0.78358),
Colors.RGB(0.59897,0.71987,0.78245),
Colors.RGB(0.59359,0.71646,0.78136),
Colors.RGB(0.58828,0.71303,0.78029),
Colors.RGB(0.58301,0.70959,0.77926),
Colors.RGB(0.57781,0.70613,0.77825),
Colors.RGB(0.57267,0.70266,0.77728),
Colors.RGB(0.56758,0.69917,0.77633),
Colors.RGB(0.56256,0.69567,0.7754),
Colors.RGB(0.55759,0.69216,0.7745),
Colors.RGB(0.55268,0.68863,0.77363),
Colors.RGB(0.54784,0.68509,0.77277),
Colors.RGB(0.54306,0.68154,0.77194),
Colors.RGB(0.53834,0.67797,0.77113),
Colors.RGB(0.53368,0.67439,0.77033),
Colors.RGB(0.52909,0.6708,0.76956),
Colors.RGB(0.52456,0.66719,0.7688),
Colors.RGB(0.5201,0.66357,0.76805),
Colors.RGB(0.5157,0.65994,0.76732),
Colors.RGB(0.51136,0.6563,0.76661),
Colors.RGB(0.50709,0.65264,0.7659),
Colors.RGB(0.50289,0.64897,0.76521),
Colors.RGB(0.49875,0.64529,0.76454),
Colors.RGB(0.49468,0.64159,0.76387),
Colors.RGB(0.49067,0.63789,0.76321),
Colors.RGB(0.48674,0.63417,0.76256),
Colors.RGB(0.48287,0.63043,0.76192),
Colors.RGB(0.47907,0.62669,0.76128),
Colors.RGB(0.47534,0.62293,0.76065),
Colors.RGB(0.47168,0.61916,0.76003),
Colors.RGB(0.46808,0.61537,0.75941),
Colors.RGB(0.46456,0.61157,0.75879),
Colors.RGB(0.46111,0.60776,0.75818),
Colors.RGB(0.45773,0.60394,0.75757),
Colors.RGB(0.45443,0.6001,0.75696),
Colors.RGB(0.45119,0.59625,0.75635),
Colors.RGB(0.44802,0.59238,0.75574),
Colors.RGB(0.44493,0.58851,0.75513),
Colors.RGB(0.44191,0.58461,0.75452),
Colors.RGB(0.43897,0.58071,0.7539),
Colors.RGB(0.43609,0.57679,0.75328),
Colors.RGB(0.43329,0.57286,0.75266),
Colors.RGB(0.43056,0.56891,0.75203),
Colors.RGB(0.42791,0.56495,0.75139),
Colors.RGB(0.42532,0.56098,0.75075),
Colors.RGB(0.42281,0.55699,0.7501),
Colors.RGB(0.42038,0.55299,0.74944),
Colors.RGB(0.41801,0.54898,0.74877),
Colors.RGB(0.41572,0.54495,0.74809),
Colors.RGB(0.4135,0.54091,0.7474),
Colors.RGB(0.41135,0.53685,0.7467),
Colors.RGB(0.40928,0.53278,0.74598),
Colors.RGB(0.40727,0.52869,0.74525),
Colors.RGB(0.40533,0.52459,0.7445),
Colors.RGB(0.40347,0.52048,0.74374),
Colors.RGB(0.40167,0.51635,0.74296),
Colors.RGB(0.39994,0.51221,0.74217),
Colors.RGB(0.39827,0.50805,0.74135),
Colors.RGB(0.39667,0.50388,0.74052),
Colors.RGB(0.39514,0.4997,0.73967),
Colors.RGB(0.39367,0.4955,0.73879),
Colors.RGB(0.39226,0.49128,0.7379),
Colors.RGB(0.39092,0.48706,0.73698),
Colors.RGB(0.38964,0.48281,0.73604),
Colors.RGB(0.38841,0.47856,0.73507),
Colors.RGB(0.38724,0.47429,0.73408),
Colors.RGB(0.38613,0.47,0.73306),
Colors.RGB(0.38508,0.4657,0.73202),
Colors.RGB(0.38407,0.46139,0.73095),
Colors.RGB(0.38312,0.45706,0.72985),
Colors.RGB(0.38222,0.45272,0.72872),
Colors.RGB(0.38137,0.44837,0.72756),
Colors.RGB(0.38056,0.444,0.72637),
Colors.RGB(0.3798,0.43962,0.72514),
Colors.RGB(0.37909,0.43522,0.72389),
Colors.RGB(0.37841,0.43081,0.7226),
Colors.RGB(0.37778,0.42638,0.72128),
Colors.RGB(0.37718,0.42195,0.71992),
Colors.RGB(0.37662,0.4175,0.71852),
Colors.RGB(0.3761,0.41303,0.71709),
Colors.RGB(0.37561,0.40855,0.71563),
Colors.RGB(0.37515,0.40406,0.71412),
Colors.RGB(0.37472,0.39956,0.71257),
Colors.RGB(0.37431,0.39504,0.71099),
Colors.RGB(0.37393,0.39051,0.70936),
Colors.RGB(0.37358,0.38596,0.70769),
Colors.RGB(0.37325,0.38141,0.70598),
Colors.RGB(0.37294,0.37684,0.70423),
Colors.RGB(0.37264,0.37226,0.70243),
Colors.RGB(0.37236,0.36766,0.70058),
Colors.RGB(0.3721,0.36306,0.69869),
Colors.RGB(0.37185,0.35844,0.69676),
Colors.RGB(0.37161,0.35381,0.69477),
Colors.RGB(0.37138,0.34917,0.69274),
Colors.RGB(0.37116,0.34452,0.69065),
Colors.RGB(0.37094,0.33986,0.68852),
Colors.RGB(0.37073,0.33518,0.68633),
Colors.RGB(0.37052,0.3305,0.68409),
Colors.RGB(0.37031,0.3258,0.68179),
Colors.RGB(0.37009,0.3211,0.67944),
Colors.RGB(0.36988,0.31638,0.67704),
Colors.RGB(0.36966,0.31166,0.67457),
Colors.RGB(0.36943,0.30693,0.67205),
Colors.RGB(0.3692,0.30219,0.66947),
Colors.RGB(0.36895,0.29744,0.66682),
Colors.RGB(0.3687,0.29269,0.66412),
Colors.RGB(0.36843,0.28793,0.66135),
Colors.RGB(0.36814,0.28316,0.65851),
Colors.RGB(0.36784,0.27839,0.65561),
Colors.RGB(0.36752,0.27361,0.65263),
Colors.RGB(0.36718,0.26883,0.64959),
Colors.RGB(0.36681,0.26405,0.64648),
Colors.RGB(0.36642,0.25926,0.64329),
Colors.RGB(0.36601,0.25448,0.64003),
Colors.RGB(0.36557,0.2497,0.6367),
Colors.RGB(0.3651,0.24492,0.63328),
Colors.RGB(0.36459,0.24014,0.62979),
Colors.RGB(0.36406,0.23536,0.62621),
Colors.RGB(0.36349,0.2306,0.62255),
Colors.RGB(0.36288,0.22584,0.6188),
Colors.RGB(0.36223,0.22109,0.61497),
Colors.RGB(0.36154,0.21636,0.61105),
Colors.RGB(0.3608,0.21164,0.60704),
Colors.RGB(0.36003,0.20694,0.60293),
Colors.RGB(0.3592,0.20226,0.59873),
Colors.RGB(0.35832,0.1976,0.59443),
Colors.RGB(0.3574,0.19297,0.59003),
Colors.RGB(0.35641,0.18837,0.58553),
Colors.RGB(0.35537,0.1838,0.58093),
Colors.RGB(0.35428,0.17927,0.57623),
Colors.RGB(0.35312,0.17479,0.57142),
Colors.RGB(0.35189,0.17034,0.5665),
Colors.RGB(0.3506,0.16595,0.56148),
Colors.RGB(0.34925,0.16161,0.55635),
Colors.RGB(0.34782,0.15734,0.55111),
Colors.RGB(0.34632,0.15313,0.54576),
Colors.RGB(0.34474,0.14899,0.5403),
Colors.RGB(0.34309,0.14492,0.53474),
Colors.RGB(0.34135,0.14094,0.52907),
Colors.RGB(0.33954,0.13705,0.52329),
Colors.RGB(0.33765,0.13325,0.51741),
Colors.RGB(0.33567,0.12954,0.51143),
Colors.RGB(0.33361,0.12594,0.50535),
Colors.RGB(0.33146,0.12244,0.49918),
Colors.RGB(0.32923,0.11906,0.49293),
Colors.RGB(0.32691,0.11579,0.48659),
Colors.RGB(0.32451,0.11263,0.48017),
Colors.RGB(0.32203,0.1096,0.47368),
Colors.RGB(0.31946,0.10669,0.46714),
Colors.RGB(0.31682,0.1039,0.46053),
Colors.RGB(0.31409,0.10123,0.45388),
Colors.RGB(0.31129,0.09868,0.44719),
Colors.RGB(0.30842,0.09626,0.44047),
Colors.RGB(0.30549,0.09395,0.43373),
Colors.RGB(0.30249,0.09176,0.42697),
Colors.RGB(0.29942,0.08968,0.42022),
Colors.RGB(0.29631,0.08771,0.41346),
Colors.RGB(0.29315,0.08585,0.40672),
Colors.RGB(0.28994,0.08409,0.4),
Colors.RGB(0.28669,0.08243,0.39331),
Colors.RGB(0.28341,0.08086,0.38666),
Colors.RGB(0.28011,0.07939,0.38005),
Colors.RGB(0.27678,0.078,0.37349),
Colors.RGB(0.27344,0.0767,0.367),
Colors.RGB(0.27009,0.07548,0.36056),
Colors.RGB(0.26673,0.07434,0.3542),
Colors.RGB(0.26338,0.07328,0.34792),
Colors.RGB(0.26004,0.07229,0.34172),
Colors.RGB(0.25671,0.07138,0.33561),
Colors.RGB(0.25341,0.07053,0.32959),
Colors.RGB(0.25013,0.06976,0.32368),
Colors.RGB(0.24688,0.06905,0.31787),
Colors.RGB(0.24367,0.06842,0.31217),
Colors.RGB(0.24051,0.06786,0.30657),
Colors.RGB(0.23739,0.06737,0.30109),
Colors.RGB(0.23433,0.06694,0.29574),
Colors.RGB(0.23133,0.06658,0.29051),
Colors.RGB(0.22839,0.06629,0.28541),
Colors.RGB(0.22552,0.06608,0.28043),
Colors.RGB(0.22273,0.06593,0.2756),
Colors.RGB(0.22001,0.06586,0.2709),
Colors.RGB(0.21738,0.06586,0.26634),
Colors.RGB(0.21483,0.06594,0.26192),
Colors.RGB(0.21237,0.06609,0.25765),
Colors.RGB(0.21001,0.06631,0.25353),
Colors.RGB(0.20774,0.06661,0.24955),
Colors.RGB(0.20558,0.06699,0.24572),
Colors.RGB(0.20352,0.06744,0.24206),
Colors.RGB(0.20156,0.06798,0.23853),
Colors.RGB(0.19972,0.06859,0.23517),
Colors.RGB(0.19795,0.06931,0.23195),
Colors.RGB(0.19608,0.07032,0.22875),
Colors.RGB(0.1941,0.07161,0.22559),
Colors.RGB(0.19199,0.07318,0.22243),
Colors.RGB(0.18976,0.07502,0.2193),
Colors.RGB(0.18739,0.0771,0.21619),
Colors.RGB(0.18488,0.07943,0.21308),
Colors.RGB(0.18774,0.07725,0.21387),
Colors.RGB(0.1905,0.07531,0.21466),
Colors.RGB(0.19315,0.07361,0.21542),
Colors.RGB(0.19572,0.07216,0.21617),
Colors.RGB(0.19819,0.07097,0.21691),
Colors.RGB(0.20059,0.07006,0.21763),
Colors.RGB(0.2029,0.06944,0.21833),
Colors.RGB(0.20532,0.06892,0.21912),
Colors.RGB(0.20786,0.06848,0.22),
Colors.RGB(0.21053,0.06812,0.22099),
Colors.RGB(0.21333,0.06783,0.22207),
Colors.RGB(0.21625,0.06762,0.22325),
Colors.RGB(0.21931,0.06747,0.22451),
Colors.RGB(0.22247,0.06739,0.22586),
Colors.RGB(0.22575,0.06738,0.22729),
Colors.RGB(0.22916,0.06743,0.2288),
Colors.RGB(0.23266,0.06756,0.23038),
Colors.RGB(0.23627,0.06774,0.23202),
Colors.RGB(0.24,0.06799,0.23373),
Colors.RGB(0.24381,0.06829,0.2355),
Colors.RGB(0.24772,0.06865,0.23733),
Colors.RGB(0.25173,0.06906,0.2392),
Colors.RGB(0.25582,0.06953,0.24112),
Colors.RGB(0.25999,0.07005,0.24308),
Colors.RGB(0.26426,0.07062,0.24508),
Colors.RGB(0.26859,0.07123,0.2471),
Colors.RGB(0.273,0.07188,0.24916),
Colors.RGB(0.27747,0.07258,0.25123),
Colors.RGB(0.28202,0.07332,0.25333),
Colors.RGB(0.28662,0.07409,0.25543),
Colors.RGB(0.29129,0.0749,0.25755),
Colors.RGB(0.296,0.07575,0.25967),
Colors.RGB(0.30077,0.07662,0.26179),
Colors.RGB(0.30559,0.07752,0.26391),
Colors.RGB(0.31046,0.07846,0.26602),
Colors.RGB(0.31536,0.07942,0.26812),
Colors.RGB(0.3203,0.08041,0.2702),
Colors.RGB(0.32528,0.08143,0.27227),
Colors.RGB(0.33029,0.08247,0.27431),
Colors.RGB(0.33533,0.08353,0.27633),
Colors.RGB(0.3404,0.08462,0.27831),
Colors.RGB(0.34549,0.08574,0.28027),
Colors.RGB(0.35061,0.08688,0.28219),
Colors.RGB(0.35574,0.08804,0.28407),
Colors.RGB(0.36089,0.08923,0.28591),
Colors.RGB(0.36605,0.09044,0.28771),
Colors.RGB(0.37123,0.09168,0.28946),
Colors.RGB(0.37641,0.09295,0.29116),
Colors.RGB(0.3816,0.09424,0.29281),
Colors.RGB(0.3868,0.09556,0.29441),
Colors.RGB(0.392,0.09692,0.29595),
Colors.RGB(0.3972,0.0983,0.29744),
Colors.RGB(0.4024,0.09972,0.29887),
Colors.RGB(0.40759,0.10118,0.30024),
Colors.RGB(0.41278,0.10267,0.30154),
Colors.RGB(0.41796,0.10421,0.30279),
Colors.RGB(0.42313,0.10578,0.30397),
Colors.RGB(0.42829,0.1074,0.30508),
Colors.RGB(0.43344,0.10906,0.30614),
Colors.RGB(0.43856,0.11078,0.30713),
Colors.RGB(0.44367,0.11254,0.30805),
Colors.RGB(0.44876,0.11435,0.30891),
Colors.RGB(0.45383,0.11622,0.3097),
Colors.RGB(0.45887,0.11815,0.31044),
Colors.RGB(0.46389,0.12013,0.3111),
Colors.RGB(0.46888,0.12216,0.31171),
Colors.RGB(0.47384,0.12426,0.31225),
Colors.RGB(0.47877,0.12642,0.31274),
Colors.RGB(0.48367,0.12865,0.31317),
Colors.RGB(0.48853,0.13093,0.31355),
Colors.RGB(0.49336,0.13328,0.31387),
Colors.RGB(0.49814,0.13569,0.31414),
Colors.RGB(0.5029,0.13817,0.31436),
Colors.RGB(0.50761,0.14071,0.31453),
Colors.RGB(0.51228,0.14332,0.31466),
Colors.RGB(0.51691,0.14598,0.31475),
Colors.RGB(0.5215,0.14872,0.31481),
Colors.RGB(0.52604,0.15151,0.31483),
Colors.RGB(0.53054,0.15436,0.31481),
Colors.RGB(0.535,0.15728,0.31477),
Colors.RGB(0.53942,0.16025,0.3147),
Colors.RGB(0.54379,0.16328,0.31461),
Colors.RGB(0.54811,0.16636,0.3145),
Colors.RGB(0.5524,0.1695,0.31437),
Colors.RGB(0.55663,0.1727,0.31423),
Colors.RGB(0.56082,0.17594,0.31408),
Colors.RGB(0.56497,0.17924,0.31391),
Colors.RGB(0.56908,0.18258,0.31374),
Colors.RGB(0.57314,0.18597,0.31357),
Colors.RGB(0.57716,0.18941,0.3134),
Colors.RGB(0.58113,0.19289,0.31322),
Colors.RGB(0.58506,0.19641,0.31305),
Colors.RGB(0.58895,0.19997,0.31289),
Colors.RGB(0.59279,0.20357,0.31273),
Colors.RGB(0.5966,0.20721,0.31259),
Colors.RGB(0.60036,0.21089,0.31245),
Colors.RGB(0.60408,0.2146,0.31233),
Colors.RGB(0.60777,0.21835,0.31222),
Colors.RGB(0.61141,0.22213,0.31213),
Colors.RGB(0.61501,0.22594,0.31206),
Colors.RGB(0.61857,0.22979,0.312),
Colors.RGB(0.62209,0.23366,0.31197),
Colors.RGB(0.62557,0.23757,0.31197),
Colors.RGB(0.62902,0.2415,0.31198),
Colors.RGB(0.63243,0.24546,0.31203),
Colors.RGB(0.63579,0.24944,0.3121),
Colors.RGB(0.63912,0.25345,0.3122),
Colors.RGB(0.64242,0.25749,0.31233),
Colors.RGB(0.64567,0.26155,0.31249),
Colors.RGB(0.64889,0.26564,0.31268),
Colors.RGB(0.65207,0.26975,0.31291),
Colors.RGB(0.65522,0.27388,0.31317),
Colors.RGB(0.65833,0.27803,0.31346),
Colors.RGB(0.6614,0.28221,0.3138),
Colors.RGB(0.66444,0.2864,0.31417),
Colors.RGB(0.66744,0.29062,0.31458),
Colors.RGB(0.6704,0.29486,0.31504),
Colors.RGB(0.67333,0.29912,0.31553),
Colors.RGB(0.67622,0.3034,0.31607),
Colors.RGB(0.67907,0.30769,0.31666),
Colors.RGB(0.68189,0.31201,0.31728),
Colors.RGB(0.68468,0.31635,0.31796),
Colors.RGB(0.68743,0.3207,0.31868),
Colors.RGB(0.69014,0.32507,0.31945),
Colors.RGB(0.69282,0.32946,0.32028),
Colors.RGB(0.69546,0.33387,0.32115),
Colors.RGB(0.69806,0.33829,0.32207),
Colors.RGB(0.70063,0.34273,0.32305),
Colors.RGB(0.70316,0.34719,0.32409),
Colors.RGB(0.70566,0.35166,0.32518),
Colors.RGB(0.70812,0.35615,0.32633),
Colors.RGB(0.71055,0.36066,0.32754),
Colors.RGB(0.71293,0.36518,0.3288),
Colors.RGB(0.71529,0.36971,0.33013),
Colors.RGB(0.7176,0.37426,0.33152),
Colors.RGB(0.71989,0.37883,0.33298),
Colors.RGB(0.72213,0.38341,0.33449),
Colors.RGB(0.72434,0.388,0.33608),
Colors.RGB(0.72651,0.39261,0.33773),
Colors.RGB(0.72865,0.39723,0.33945),
Colors.RGB(0.73075,0.40187,0.34124),
Colors.RGB(0.73281,0.40652,0.34311),
Colors.RGB(0.73484,0.41118,0.34504),
Colors.RGB(0.73683,0.41585,0.34705),
Colors.RGB(0.73879,0.42054,0.34913),
Colors.RGB(0.74071,0.42523,0.35129),
Colors.RGB(0.7426,0.42994,0.35353),
Colors.RGB(0.74445,0.43466,0.35584),
Colors.RGB(0.74627,0.43939,0.35823),
Colors.RGB(0.74805,0.44413,0.36071),
Colors.RGB(0.74979,0.44888,0.36326),
Colors.RGB(0.75151,0.45364,0.3659),
Colors.RGB(0.75319,0.45841,0.36862),
Colors.RGB(0.75483,0.46319,0.37143),
Colors.RGB(0.75644,0.46798,0.37432),
Colors.RGB(0.75802,0.47277,0.3773),
Colors.RGB(0.75957,0.47757,0.38036),
Colors.RGB(0.76109,0.48238,0.38351),
Colors.RGB(0.76257,0.48719,0.38675),
Colors.RGB(0.76403,0.49201,0.39008),
Colors.RGB(0.76545,0.49683,0.3935),
Colors.RGB(0.76685,0.50166,0.39701),
Colors.RGB(0.76822,0.50649,0.40061),
Colors.RGB(0.76956,0.51133,0.4043),
Colors.RGB(0.77088,0.51617,0.40809),
Colors.RGB(0.77217,0.52101,0.41196),
Colors.RGB(0.77344,0.52585,0.41593),
Colors.RGB(0.77468,0.5307,0.41998),
Colors.RGB(0.77591,0.53554,0.42413),
Colors.RGB(0.77711,0.54039,0.42837),
Colors.RGB(0.77829,0.54523,0.43271),
Colors.RGB(0.77946,0.55007,0.43713),
Colors.RGB(0.78061,0.55491,0.44164),
Colors.RGB(0.78174,0.55975,0.44625),
Colors.RGB(0.78286,0.56459,0.45094),
Colors.RGB(0.78397,0.56942,0.45572),
Colors.RGB(0.78507,0.57424,0.46059),
Colors.RGB(0.78616,0.57906,0.46555),
Colors.RGB(0.78724,0.58388,0.47059),
Colors.RGB(0.78832,0.58869,0.47572),
Colors.RGB(0.78939,0.59349,0.48093),
Colors.RGB(0.79046,0.59828,0.48622),
Colors.RGB(0.79153,0.60307,0.4916),
Colors.RGB(0.7926,0.60784,0.49705),
Colors.RGB(0.79368,0.61261,0.50258),
Colors.RGB(0.79476,0.61737,0.50819),
Colors.RGB(0.79584,0.62211,0.51387),
Colors.RGB(0.79694,0.62685,0.51963),
Colors.RGB(0.79804,0.63157,0.52545),
Colors.RGB(0.79916,0.63628,0.53134),
Colors.RGB(0.80029,0.64098,0.53731),
Colors.RGB(0.80144,0.64567,0.54333),
Colors.RGB(0.80261,0.65034,0.54942),
Colors.RGB(0.80379,0.65499,0.55556),
Colors.RGB(0.80499,0.65964,0.56177),
Colors.RGB(0.80621,0.66426,0.56803),
Colors.RGB(0.80746,0.66887,0.57434),
Colors.RGB(0.80873,0.67346,0.5807),
Colors.RGB(0.81003,0.67804,0.58711),
Colors.RGB(0.81135,0.68259,0.59356),
Colors.RGB(0.8127,0.68713,0.60005),
Colors.RGB(0.81408,0.69165,0.60658),
Colors.RGB(0.81548,0.69615,0.61315),
Colors.RGB(0.81692,0.70062,0.61975),
Colors.RGB(0.81838,0.70507,0.62638),
Colors.RGB(0.81987,0.7095,0.63304),
Colors.RGB(0.82139,0.71391,0.63972),
Colors.RGB(0.82295,0.71829,0.64642),
Colors.RGB(0.82453,0.72265,0.65313),
Colors.RGB(0.82614,0.72697,0.65986),
Colors.RGB(0.82777,0.73127,0.6666),
Colors.RGB(0.82943,0.73554,0.67334),
Colors.RGB(0.83112,0.73977,0.68008),
Colors.RGB(0.83283,0.74397,0.68682),
Colors.RGB(0.83457,0.74814,0.69356),
Colors.RGB(0.83632,0.75227,0.70028),
Colors.RGB(0.83809,0.75635,0.70699),
Colors.RGB(0.83988,0.7604,0.71367),
Colors.RGB(0.84168,0.7644,0.72033),
Colors.RGB(0.84349,0.76836,0.72697),
Colors.RGB(0.8453,0.77226,0.73356),
Colors.RGB(0.84711,0.77612,0.74012),
Colors.RGB(0.84892,0.77992,0.74664),
Colors.RGB(0.85073,0.78366,0.7531),
Colors.RGB(0.85252,0.78735,0.7595),
Colors.RGB(0.85429,0.79097,0.76584),
Colors.RGB(0.85604,0.79453,0.77211),
Colors.RGB(0.85776,0.79802,0.7783),
Colors.RGB(0.85943,0.80144,0.7844),
Colors.RGB(0.86107,0.80479,0.7904),
Colors.RGB(0.86266,0.80806,0.79628),
Colors.RGB(0.86418,0.81125,0.80205),
Colors.RGB(0.86565,0.81436,0.80767),
Colors.RGB(0.86705,0.81738,0.81313),
Colors.RGB(0.8684,0.82031,0.81842),
Colors.RGB(0.86969,0.82314,0.8235),
Colors.RGB(0.87094,0.82587,0.82838),
Colors.RGB(0.87215,0.82848,0.83305),
Colors.RGB(0.87335,0.83097,0.83749),
Colors.RGB(0.87454,0.83332,0.84172),
Colors.RGB(0.87571,0.83553,0.84576),
Colors.RGB(0.87688,0.83759,0.84961),
Colors.RGB(0.87802,0.8395,0.85331),
Colors.RGB(0.87913,0.84126,0.85686),
Colors.RGB(0.88019,0.84285,0.86027),
Colors.RGB(0.88119,0.84429,0.86357),
Colors.RGB(0.88212,0.84557,0.86674),
Colors.RGB(0.88295,0.84669,0.8698),
Colors.RGB(0.88369,0.84765,0.87274),
Colors.RGB(0.88433,0.84845,0.87557),
Colors.RGB(0.88485,0.84908,0.87828),
Colors.RGB(0.88526,0.84956,0.88088),
Colors.RGB(0.88555,0.84987,0.88336),
Colors.RGB(0.88571,0.85002,0.88573),
], "matplotlib", "cyclical starts and ends on the same color, and have two symmetric halves with equal lightness, but diverging color, perceptually uniform and colorblind friendly.")


text:
module LoadAllPackages

using Base: PkgId, UUID
using Pkg: Pkg, TOML

function loadall(project = Base.active_project())
    deps = get(TOML.parsefile(project), "deps", nothing)
    deps === nothing && return
    @info "Loading $(length(deps)) packages..."
    t0 = time_ns()
    for (name, uuid) in deps
        pkg = PkgId(UUID(uuid), name)
        @debug "Loading `$pkg`..."
        Base.require(pkg)
        @debug "`$pkg` loaded."
    end
    sec = (time_ns() - t0) / 10^9
    @info "$(length(deps)) packages loaded in $sec seconds"
    return
end

end


text:
# AUTO GENERATED FILE - DO NOT EDIT

export 'feffery'_antdsider

"""
    'feffery'_antdsider(;kwargs...)
    'feffery'_antdsider(children::Any;kwargs...)
    'feffery'_antdsider(children_maker::Function;kwargs...)

An AntdSider component.

Keyword arguments:
- `children` (a list of or a singular dash component, string or number; optional): The content of the tab - will only be displayed if this tab is selected
- `id` (String; optional)
- `className` (String; optional)
- `collapsed` (Bool; optional)
- `collapsedWidth` (Real; optional)
- `collapsible` (Bool; optional)
- `defaultCollapsed` (Bool; optional)
- `loading_state` (optional): . loading_state has the following type: lists containing elements 'is_loading', 'prop_name', 'component_name'.
Those elements have the following types:
  - `is_loading` (Bool; optional): Determines if the component is loading or not
  - `prop_name` (String; optional): Holds which property is loading
  - `component_name` (String; optional): Holds the name of the component that is loading
- `reverseArrow` (Bool; optional)
- `style` (Dict; optional)
- `theme` (String; optional)
- `trigger` (a list of or a singular dash component, string or number; optional)
- `width` (Real; optional)
"""
function 'feffery'_antdsider(; kwargs...)
        available_props = Symbol[:children, :id, :className, :collapsed, :collapsedWidth, :collapsible, :defaultCollapsed, :loading_state, :reverseArrow, :style, :theme, :trigger, :width]
        wild_props = Symbol[]
        return Component("'feffery'_antdsider", "AntdSider", "feffery_antd_components", available_props, wild_props; kwargs...)
end

'feffery'_antdsider(children::Any; kwargs...) = 'feffery'_antdsider(;kwargs..., children = children)
'feffery'_antdsider(children_maker::Function; kwargs...) = 'feffery'_antdsider(children_maker(); kwargs...)

max_stars_repo_name: Beddow/feffery-antd-components

max_issues_repo_name: Beddow/feffery-antd-components

max_forks_repo_name: Beddow/feffery-antd-components


text:
using JuMP, CPLEX

#Preparing optimization model
model = Model(CPLEX.Optimizer)
set_optimizer_attribute(model, "CPX_PARAM_EPINT", 1e-8)

const N_GEN = 3
const N_HRS = 4

fixedCost = [5 7 6];
variableCost = [0.100 0.125 0.150];
startUpCost = [20 18 5];
shutdownCost = [0.5 0.3 1.0];

minPowerOut = [50 80 40];
maxPowerOut = [350 200 140];
rampDownLimit = [300 150 100];
rampShutDownLimit = [300 150 100];
rampUpLimit = [200 100 100];
rampStartUpLimit = [200 100 100];

demand = [100; 160; 500; 400];
reserves = [0; 16; 50; 40];

@variable(model, genOnOff[1:N_GEN, 1:N_HRS], Bin)
@variable(model, genOut[1:N_GEN, 1:N_HRS] >= 0)
@variable(model, genStartUp[1:N_GEN, 1:N_HRS], Bin)
@variable(model, genShutDown[1:N_GEN, 1:N_HRS], Bin)

#Setting the objective function
@objective(model, Min, sum(sum( fixedCost[g]*genOnOff[g,t] + variableCost[g]*genOut[g,t] + startUpCost[g]*genStartUp[g,t] + shutdownCost[g]*genShutDown[g, t] for g in 1:N_GEN) for t in 1:N_HRS) )

#Generator #1 should be ON at t=0
@constraint(model, constInit[t=1, g=3], genOnOff[g,t] == 1)

#Initial constraints
@constraint(model, constShutDn[g=1:N_GEN], genShutDown[g,1] == 0)
@constraint(model, constShutDn[g=1:N_GEN], genStartUp[g,1] == 0)

#Demand constraint
#for t in 1:N_HRS
@constraint(model, constDemand[t=1:N_HRS], sum(genOut[g, t] for g=1:N_GEN) == demand[t])
#end

#Power bounds
@constraint(model, constMaxPow[t=1:N_HRS, g=1:N_GEN],  genOut[g,t] <= maxPowerOut[g]*genOnOff[g,t] )
@constraint(model, constMinPow[t=1:N_HRS, g=1:N_GEN],  genOut[g,t] >= minPowerOut[g]*genOnOff[g,t] )

#Logical Conditions
@constraint(model, constLog1[t=2:N_HRS, g=1:N_GEN], (genStartUp[g,t] - genShutDown[g,t]) == (genOnOff[g,t] - genOnOff[g,t-1]))
@constraint(model, constLog2[t=1:N_HRS, g=1:N_GEN], genStartUp[g,t] + genShutDown[g, t] <= 1)

#Ramp limits
@constraint(model, constRampUp[t=2:N_HRS, g=1:N_GEN], genOut[g,t] - genOut[g,t-1] <= rampUpLimit[g]*genOnOff[g, t-1] + rampStartUpLimit[g]*genStartUp[g,t])
@constraint(model, constRampDown[t=2:N_HRS, g=1:N_GEN], (genOut[g,t-1] - genOut[g,t]) <= (rampDownLimit[g]*genOnOff[g,t] + rampShutDownLimit[g]*genShutDown[g,t]))

#Reserves
@constraint(model, constReserves[t=1:N_HRS], sum(maxPowerOut[g]*genOnOff[g,t] for g=1:N_GEN) >= (demand[t] + reserves[t]))

optimize!(model)

println("------------------------------------")
println("------- Optimal Solutions -------")
for g=1:N_GEN
  for t=1:N_HRS
    print(JuMP.value(genOut[g,t]),"\t")
  end
  println("\t Generator: $g")
end

println("------------------------------------")
println("------- GENERATORS COMMITED --------")

for g=1:N_GEN
  for t=1:N_HRS
    print(JuMP.value(genOnOff[g,t]),"\t")
  end
  println("\t Generator: $g")
end

max_stars_repo_name: power-systems-op/grace-arpae

max_issues_repo_name: power-systems-op/grace-arpae

max_forks_repo_name: power-systems-op/grace-arpae


text:
using DataStructures
using OrderedCollections
using BenchmarkTools
using Random

suite = BenchmarkGroup()

dicttypes = [Dict, RobinDict, SwissDict]
KVtypes = [Int, Float64]
aexps = [4]
for dt in dicttypes
    for aexp in aexps
        for K in KVtypes
            V = K
            Random.seed!(0)
            sample = rand(K, 10^aexp, 2)
            entries = [sample[i, 1]=>sample[i, 2] for i in 1:size(sample, 1)]
            sample_str = string(" 10^$aexp elem $K")

            dict = (dt == Dict) ? string("base-Dict") : string(dt)

            # dict_el_type: {Any, Any}
            suite["make "*dict*"()"*sample_str] = @benchmarkable $(dt)($entries)
            suite["pop! "*dict*"()"*sample_str] = @benchmarkable pop!(h) setup=(h=$(dt)($entries))
            suite["empty! "*dict*"()"*sample_str] = @benchmarkable empty!(h) setup=(h=$(dt)($entries))
            suite["find-success "*dict*"()"*sample_str] = @benchmarkable get(h, e, -1) setup=(h=$(dt)($entries); e=($entries)[1, 1])
            suite["find-failure "*dict*"()"*sample_str] = @benchmarkable get(h, e, -1) setup=(h=$(dt)($entries); e=first(pop!(h)))
            # dict_el_type: {K, V}
            suite["make "*dict*"{$K,$V}()"*sample_str] = @benchmarkable $(dt){$K,$V}($entries)
            suite["pop! "*dict*"{$K,$V}()"*sample_str] = @benchmarkable pop!(h) setup=(h=$(dt){$K,$V}($entries))
            suite["empty! "*dict*"{$K,$V}()"*sample_str] = @benchmarkable empty!(h) setup=(h=$(dt){$K,$V}($entries))
            suite["find-success "*dict*"{$K,$V}()"*sample_str] = @benchmarkable get(h, e, -1) setup=(h=$(dt){$K,$V}($entries); e=($entries)[1, 1])
            suite["find-failure "*dict*"{$K,$V}()"*sample_str] = @benchmarkable get(h, e, -1) setup=(h=$(dt){$K,$V}($entries); e=first(pop!(h)))
        end
    end
end

results = run(suite, verbose=true);

function filter_result(results, op, dt, apow=4)
    dict = (dt == Dict) ? string("base-Dict") : string(dt)
    for tr in keys(filter(k->(occursin(op, first(k)) && occursin(dict, first(k))
                                 && occursin("10^$apow", first(k))) , results))
           println(tr, " memory ", memory(results[tr])/1024, " kB")
           println(tr, " minimum time ", minimum(results[tr]))
       end
end

max_stars_repo_name: ethansaxenian/RosettaDecode

max_issues_repo_name: ethansaxenian/RosettaDecode

max_forks_repo_name: ethansaxenian/RosettaDecode


text:
using Test
using WignerSymbols

## wigner3j summation tests
function test_serial_3j(n::Int)
    new_cache = WignerCache()  # create a new cache
    sum = zero(BigFloat)
    for i in 1:n
        for j in 1:n
            for k in 1:n
                sum += wigner3j(new_cache, BigFloat, i, j, k, 0, 0)
            end
        end
    end
    sum
end

function test_threaded_3j(n::Int)
    WignerSymbols.__init__()  # reset global caches
    summands = zeros(BigFloat, Threads.nthreads())
    Threads.@threads for i in 1:n
        tid = Threads.threadid()
        for j in 1:n
            for k in 1:n
                summands[tid] += wigner3j(BigFloat, i, j, k, 0, 0)
            end
        end
    end
    sum(summands)
end

@testset "wigner3j threaded sum" begin
    for max_j in [7, 30, 100]
        tol = 10*max_j*eps(BigFloat)
        @test test_serial_3j(max_j) ≈ test_threaded_3j(max_j) atol=tol
    end
end

## wigner6j summation tests
function test_serial_6j(n::Int)
    new_cache = WignerCache()  # create a new cache
    s = zero(BigFloat)
    for i in 1:n
        s += wigner6j(new_cache, BigFloat, i, i, i, i, i, i)
    end
    s
end

function test_threaded_6j(n::Int)
    summands = zeros(BigFloat, Threads.nthreads())
    Threads.@threads for i in 1:n
        summands[Threads.threadid()] += wigner6j(BigFloat, i, i, i, i, i, i)
    end
    sum(summands)
end

@testset "wigner6j threaded sum" begin
    for max_j in [7, 51, 100, 300, 701]
        tol = 10*max_j*eps(BigFloat)
        @test test_serial_6j(max_j) ≈ test_threaded_6j(max_j) atol=tol
    end
end