Rel 0.1.0

Author: goedman

research/diff_eq.jl

@@ -0,0 +1,26 @@
+using ComponentArrays
+using OrdinaryDiffEq
+using LinearAlgebra
+using Unitful
+
+function newton(du, u, p, t)
+    mu = 398600.4418u"km^3/s^2"
+    r = norm(u.r)
+    du.r = u.v
+    du.v = -mu .* u.r / r^3
+end
+
+r0 = [1131.340, -2282.343, 6672.423]u"km"
+v0 = [-5.64305, 4.30333, 2.42879]u"km/s"
+Δt = 86400.0*365u"s"
+rv0 = ComponentArray(r=r0, v=v0)
+
+prob = ODEProblem(newton, rv0, (0.0u"s", Δt))
+@time sol = solve(prob, Vern8(), dt=100u"s", adaptive=false)
+#sol |> display
+println()
+
+@time sol2 = solve(prob, Tsit5(), dt=100u"s",adaptive=false)
+
+sol2.u[1].r |> display
+sol2.u[1].v |> display

research/interators_map.jl

@@ -0,0 +1,7 @@
+
+Iterators.map(x -> x^2, 1:3) |> display
+println()
+
+collect(Iterators.map(x -> x^2, 1:3)) |> display
+
+mapreduce(isodd, +, 1:5) |> display

research/macro_funs.jl

@@ -0,0 +1,15 @@
+# Is a macro useful for the helper functions in AG(), MAG(), etc?
+
+module C  
+  macro makefunc(ex)
+    return esc(quote
+      wat() = $ex
+    end)
+  end
+end
+
+using .C
+
+C.@makefunc(345)
+
+wat() |> display

research/sort_natural.jl

@@ -0,0 +1,16 @@
+function natural(x, y)
+    k(x) = [occursin(r"\d+", s) ? parse(Int, s) : s 
+            for s in split(replace(x, r"\d+" => s->" $s "))]
+    A = k(x); B= k(y)    
+    for (a, b) in zip(A, B)
+        if !isequal(a, b)
+            return typeof(a) <: typeof(b) ? isless(a, b) :
+                   isa(a,Int) ? true : false
+        end
+    end
+    return length(A) < length(B)
+end
+
+sort(["a1", "a2", "a10"]) |> display
+
+sort(["a1", "a2", "a10"], lt=natural) |> display

src/StructuralCausalModels.jl

@@ -50,7 +50,9 @@ include("methods/all_paths.jl")
 include("methods/open_paths.jl")
 include("methods/backdoor_paths.jl")
 include("methods/adjustment_sets.jl")
+include("methods/graph_utils.jl")
 include("methods/ancestral_graph.jl")
+include("methods/ribbon_graph.jl")
 
 include("utils/show_dag_path.jl")
 include("utils/ggm_conversions.jl")