Skip to content

Instantly share code, notes, and snippets.

Alexander Sikorski axsk

  • Berlin
Block or report user

Report or block axsk

Hide content and notifications from this user.

Learn more about blocking users

Contact Support about this user’s behavior.

Learn more about reporting abuse

Report abuse
View GitHub Profile
View gist:748cfd41a529c9c9b76d23c1a10da93d
using Sundials
function test1()
yS0 = [1. 2. 3.; 4. 5 6]
N, Ns = size(yS0)
GC.@preserve yS0n = [Sundials.NVector(yS0[:,j]) for j=1:Ns] # array of julia wrappers
GC.@preserve yS0nv = [Sundials.N_Vector(n) for n in yS0n] # array of pointers to sund. structs
pyS0 = pointer(yS0nv)
View cvodes.jl
using Sundials
using Sundials: N_Vector, N_Vector_S
using ForwardDiff
#using ReverseDiff
"""
sens(f!, t0, y0, p, tout; reltol, abstol)
Compute the solution and sensivities to the parametrized ODE problem defined by `f!(ẏ, t, y, p)`, starting at t0, y0, p.
View recipe.jl
using Plots
rect(x,y,w,h) = Shape(x+[0,w,w,0], y+[0,0,h,h])
@recipe function f(ds::Array{Measurement}; alpha = 1)
isempty(ds) && return []
binedges = [0,10,12,14,16,18,20,25]
nmax = nfollicles(ds)
a = alpha
View plotsolutions recipe.jl
@userplot PlotSolutions
@recipe function f(o::PlotSolutions)
samples = o.args[1]
species --> measuredinds
t --> 0:1/3:30
color_palette --> [colorant"steelblue"]
linewidth --> 0.1
label --> ""
seriesalpha --> 0.1
View sliceiterator.jl
typealias ColVec Vector{Union{Colon,Int}}
immutable SliceIterator{A<:AbstractArray}
a::A
dim::Int
cols::ColVec
end
function SliceIterator(a, dim::Int)
cols = ColVec(length(size(a)))
cols = fill!(cols, :)
@axsk
axsk / make
Created Apr 17, 2016
julia make error
View make
fatal: No names found, cannot describe anything.
fatal: bad revision '^'
JULIA usr/lib/julia/inference.ji
essentials.jl
generator.jl
reflection.jl
options.jl
promotion.jl
tuple.jl
range.jl
@axsk
axsk / Make.user
Last active Apr 17, 2016
make error
View Make.user
JULIA_CPU_TARGET:=core2
prefix=/datanumerik/bzfsikor/julia/julia-head
@axsk
axsk / cvodes.jl
Created Mar 27, 2016
cvodes-autodiff
View cvodes.jl
using Sundials, ForwardDiff
import Sundials: realtype, N_Vector
type FAndP
f::Function
p::Vector{Float64}
end
function unzip(fp::Ptr{Void})
@axsk
axsk / gist:879cb2887290650c3baf
Created Mar 12, 2016
Serialize/Deserialize segfault
View gist:879cb2887290650c3baf
signal (11): Segmentation fault [66/1088]
unknown function (ip: 0x7fe6892e0fb9)
jl_gc_collect at /nfs/datanumerik/bzfsikor/julia/julia-0.4.3/bin/../lib/julia/libjulia.so (unknown line)
jl_gc_alloc_2w at /nfs/datanumerik/bzfsikor/julia/julia-0.4.3/bin/../lib/julia/libjulia.so (unknown line)
jl_alloc_svec_uninit at /nfs/datanumerik/bzfsikor/julia/julia-0.4.3/bin/../lib/julia/libjulia.so (unknown line)
jl_alloc_svec at /nfs/datanumerik/bzfsikor/julia/julia-0.4.3/bin/../lib/julia/libjulia.so (unknown line)
unknown function (ip: 0x7fe689263015)
unknown function (ip: 0x7fe6892638d1)
unknown function (ip: 0x7fe6892641df)
unknown function (ip: 0x7fe68926835b)
You can’t perform that action at this time.