johnmyleswhite

johnmyleswhite@foobar:~/R-4.1.2$ bin/R
do_invisible
do_invisible
do_invisible
do_invisible
do_invisible
do_invisible
do_invisible
do_invisible
do_invisible

johnmyleswhite

> visible_plus_one <- function (x) {x}
>
> invisible_plus_one <- function (x) {invisible(x)}
>
> visible_plus_one(1)
[1] 1
> invisible_plus_one(1)
>
> visible_plus_one(1) + 1
[1] 2

johnmyleswhite

import MacroTools: postwalk

function missing_check(es)
    if length(es) == 0
        false
    elseif length(es) == 1
        Expr(:call, :ismissing, es[1])
    elseif length(es) == 2
        Expr(
            :call,

johnmyleswhite

struct ExpressionThunk
    source::Any
    thunk::Any
end

macro thunk(e)
    quote
        ExpressionThunk(
            $(QuoteNode(e)),
            () -> $(esc(e))

johnmyleswhite

x = [1, 2]
y = [3, 4]
z = [x, y]

z′ = copy(z)

z′[1] === x

import DataFrames: DataFrame

johnmyleswhite

import Plots

r = range(0.0, 1.0, length=512)

p = Plots.plot(
    r,
    x -> x^2 * (1 - x)^2,
)

Plots.png(p, "output2.png")

johnmyleswhite

#!/usr/bin/Rscript

png("output.png")
curve(x^2 * (1 - x)^2, from = 0, to = 1)
dev.off()

# time Rscript example.R
# null device
#           1
# 

johnmyleswhite

julia> import BenchmarkTools: @btime

julia> @inline _logistic_bounds(x::Float16) = (Float16(-16.64), Float16(7.625))
_logistic_bounds (generic function with 1 method)

julia> @inline _logistic_bounds(x::Float32) = (-103.27893f0, 16.635532f0)
_logistic_bounds (generic function with 2 methods)

julia> @inline _logistic_bounds(x::Float64) = (-744.4400719213812, 36.7368005696771)
_logistic_bounds (generic function with 3 methods)

johnmyleswhite

│ Row │ expr_type                      │ count  │ percentage  │
│     │ Any                            │ Int64  │ Float64     │
├─────┼────────────────────────────────┼────────┼─────────────┤
│ 1   │ ExprKind{:outer}               │ 1      │ 0.000226006 │
│ 2   │ ExprKind{:/=}                  │ 1      │ 0.000226006 │
│ 3   │ ExprKind{:typed_vcat}          │ 1      │ 0.000226006 │
│ 4   │ ExprKind{:.+=}                 │ 1      │ 0.000226006 │
│ 5   │ ExprKind{:%=}                  │ 3      │ 0.000678017 │
│ 6   │ ExprKind{:⊻=}                  │ 7      │ 0.00158204  │
│ 7   │ ExprKind{:.=}                  │ 8      │ 0.00180804  │

johnmyleswhite

import RCall: @R_str

function runtimes(n_reps)
    times = Array{Float64}(undef, n_reps)
    x = rand(10_000_000)
    for i in 1:n_reps
        times[i] = @elapsed sum(x)
    end
    times        
end