Skip to content

Instantly share code, notes, and snippets.

View Datseris's full-sized avatar
🐻
I am a bear.

George Datseris Datseris

🐻
I am a bear.
View GitHub Profile
@Datseris
Datseris / TraceCalls.jl
Created February 12, 2020 14:26 — forked from pfitzseb/TraceCalls.jl
TraceCalls.jl with Cassette
module TraceCalls
using Cassette
mutable struct Trace
level::Int
cutoff::Int
end
Cassette.@context TraceCtx
*.css
@Datseris
Datseris / pendulums.jl
Created May 7, 2018 09:31 — forked from cormullion/pendulums.jl
George's double pendulums
using DynamicalSystems, Luxor
function double_pendulum(u0=rand(4);
G=10.0, L1 = 1.0, L2 = 1.0, M1 = 1.0, M2 = 1.0)
@inline @inbounds function eom_dp!(du, state)
du[1] = state[2]
gm = 0.01
del_ = state[3] - state[1]
den1 = (M1 + M2)*L1 - M2*L1*cos(del_)*cos(del_)
du[2] = (M2*L1*state[2]*state[2]*sin(del_)*cos(del_) +