Skip to content

Instantly share code, notes, and snippets.

Avatar
🐼
I am a panda.

George Datseris Datseris

🐼
I am a panda.
View GitHub Profile
@Datseris
Datseris / TraceCalls.jl
Created Feb 12, 2020 — forked from pfitzseb/TraceCalls.jl
TraceCalls.jl with Cassette
View TraceCalls.jl
module TraceCalls
using Cassette
mutable struct Trace
level::Int
cutoff::Int
end
Cassette.@context TraceCtx
@Datseris
Datseris / pendulums.jl
Created May 7, 2018 — forked from cormullion/pendulums.jl
George's double pendulums
View pendulums.jl
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_) +