View log.txt
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
112,0,0,16,22]:Integer, [112,0,0,16,23]:Integer, [112,0,0,16,24]:Integer, [112,0,0,16,25]:Integer, [112,0,0,16,26]:Integer, [112,0,0,16,27]:Integer, [112,0,0,16,28]:Integer, [112,0,0,16,29]:Integer, [112,0,0,16,30]:Integer, [112,0,0,16,31]:Integer, [112,0,0,16,32]:Integer, [112,0,0,16,33]:Integer, [112,0,0,16,34]:Integer, [112,0,0,16,35]:Integer, [112,0,0,16,36]:Integer, [112,0,0,16,37]:Integer, [112,0,0,16,38]:Integer, [112,0,0,16,39]:Integer, [112,0,0,16,40]:Integer, [112,0,0,24]:Integer, [112,0,0,25]:Integer, [112,0,0,26]:Integer, [112,0,0,27]:Integer, [112,0,0,28]:Integer, [112,0,0,29]:Integer, [112,0,0,30]:Integer, [112,0,0,31]:Integer, [112,0,0,32]:Integer, [112,0,0,33]:Integer, [112,0,0,34]:Integer, [112,0,0,35]:Integer, [112,0,0,36]:Integer, [112,0,0,37]:Integer, [112,0,0,38]:Integer, [112,0,0,39]:Integer, [112,0,0,40]:Integer, [112,0,0,41]:Integer, [112,0,0,42]:Integer, [112,0,0,43]:Integer, [112,0,0,44]:Integer, [112,0,0,45]:Integer, [112,0,0,46]:Integer, [112,0,0,47]:Integer, [112,0,0,48]:Integer, |
View logsumexp.birch
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
/* | |
* One-pass algorithm of `log_sum_exp`. | |
*/ | |
function log_sum_exp_onepass(x:Real[_]) -> Real { | |
if length(x) > 0 { | |
let (mx, r) <- transform_reduce(x, (-inf, 0.0), | |
\(x:(Real, Real), y:(Real, Real)) -> { | |
let (xa, xb) <- x; | |
let (ya, yb) <- y; | |
if xa > ya { |
View stacktrace.txt
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
┌ Warning: dt(2.220446049250313e-16) <= dtmin(2.220446049250313e-16) at t=4.2228999661198057e-16. Aborting. There is either an error in your model specification or the true solution is unstable. | |
└ @ SciMLBase /home/david/.julia/packages/SciMLBase/koNdH/src/integrator_interface.jl:366 | |
Sampling 100%|██████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████| Time: 0:00:09 | |
ERROR: DimensionMismatch("arrays could not be broadcast to a common size; got a dimension with lengths 120 and 6120") | |
Stacktrace: | |
[1] _bcs1 | |
@ ./broadcast.jl:516 [inlined] | |
[2] _bcs | |
@ ./broadcast.jl:510 [inlined] | |
[3] broadcast_shape |
View debug.txt
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
julia> using DynamicPPL | |
julia> @model function demo4(n, ::Type{TV}=Vector{Float64}) where {TV} | |
m ~ Normal() | |
x = TV(undef, n) | |
@show __varinfo__ | |
for i in eachindex(x) | |
x[i] ~ Normal(m, 1.0) | |
end | |
end |
View discrete1d_benchmarks.jl
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
using BenchmarkTools | |
using Distances | |
using Distributions | |
using StatsBase | |
using LinearAlgebra | |
using SparseArrays | |
function f( | |
c, μ::DiscreteNonParametric, ν::DiscreteNonParametric, plan::SparseMatrixCSC |
View ot.jl
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
using Distributions | |
using SparseArrays | |
using LinearAlgebra | |
using StatsBase | |
function _ot_cost_plan(c, μ::DiscreteNonParametric, ν::DiscreteNonParametric; get=:plan) | |
len_μ = length(μ.p) | |
len_ν = length(ν.p) | |
wi = μ.p[1] | |
wj = ν.p[1] |
View sde.jl
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
using AdvancedMH | |
using ArraysOfArrays | |
using CairoMakie | |
using DiffEqNoiseProcess | |
using Distributions | |
using StochasticDiffEq | |
using Turing | |
using Random | |
struct CrankNicolsonProposal{P,T} <: AdvancedMH.Proposal{P} |
View emcee.ipynb
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
View processes.jl
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
using Turing | |
using Turing.RandomMeasures | |
using Random | |
function stickbreaking(rpm = DirichletProcess(0.25)) | |
# Data | |
data = [-2,2,-1.5,1.5] | |
# Base distribution |
View examples.jl
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
using Turing | |
using StatsPlots | |
using Random | |
using Statistics | |
function demo(N::Int; n::Int = 10) | |
# observation noise | |
σ² = 0.3 |
NewerOlder