Skip to content

Instantly share code, notes, and snippets.

View andreasnoack's full-sized avatar

Andreas Noack andreasnoack

View GitHub Profile
@andreasnoack
andreasnoack / log.txt
Created January 13, 2017 02:16
AppVeyor build 1.0.9879 - Pull request #19997 - fix warning that Base.cond cannot be imported into LinAlg
[00:00:00] Build started
[00:00:00] git config --global core.autocrlf input
[00:00:00] git clone -q --depth=50 https://github.com/JuliaLang/julia.git C:\projects\julia
[00:00:06] git fetch -q origin +refs/pull/19997/merge:
[00:00:07] git checkout -qf FETCH_HEAD
[00:00:10] Restoring build cache
[00:00:10]
[00:00:11] Cache 'x86_64-4.9.2-release-win32-seh-rt_v4-rev3.7z' - Downloading (63,558,997 bytes)...1%
[00:00:11] Cache 'x86_64-4.9.2-release-win32-seh-rt_v4-rev3.7z' - Downloading (63,558,997 bytes)...10%
[00:00:11] Cache 'x86_64-4.9.2-release-win32-seh-rt_v4-rev3.7z' - Downloading (63,558,997 bytes)...20%
1 ./dict.jl:235; rehash!(::Dict{Int64,Void}, ::Int64)
1 ...e/noack/julia/ui/repl.c:81; exec_program
1 ./multidimensional.jl:444; _unsafe_batchsetindex!(::Array{F...
1 ./multidimensional.jl:448; macro expansion
1 .../noack/julia/ui/repl.c:81; exec_program
1 ./operators.jl:821; setindex_shape_check(::StaticArr...
1 ...e/noack/julia/ui/repl.c:34; jl_get_ptls_states_static
2 ./operators.jl:829; setindex_shape_check(::StaticArr...
2 ...e/noack/julia/ui/repl.c:81; exec_program
3 ..._vi/elbo_objective.jl:383; add_pixel_term!(::Celeste.Determ...
KNL:
1241 ./operators.jl 276 afoldl(::Base.#*, ::Float64, ::Stat...
1443 ...rc/bivariate_normals.jl 428 Celeste.Model.GalaxyCacheComponent{...
1554 ...ic_vi/elbo_objective.jl 282 accumulate_source_pixel_brightness!...
1973 ...tic_vi/maximize_elbo.jl 74 (::Celeste.DeterministicVI.##538#54...
1980 ...6/Optim/src/optimize.jl 251 optimize(::Optim.TwiceDifferentiabl...
1980 ...6/Optim/src/optimize.jl 207 update_h!
2246 ...ic_vi/elbo_objective.jl 407 add_pixel_term!(::Celeste.Determini...
2355 ...e/src/model/fsm_util.jl 325 accum_galaxy_pos!(::Celeste.Model.B...
function parsemachine(f::String)
wks = []
for s in eachline(f)
l = split(s)
push!(wks, (l[1], parse(Int, l[2])))
end
return wks
end
addprocs(parsemachine(ENV["PE_HOSTFILE"]))
@andreasnoack
andreasnoack / roundtrip_timings_MPI.txt
Created August 11, 2016 23:29
Timings of the roundtrip benchmark - Four processes - One on each one - C3DDB cluster
Who am I? Rank 2 on node314
Who am I? Rank 3 on node315
Who am I? Rank 0 on node312
Who am I? Rank 1 on node313
INFO: Data: 2^0 Int(s)
5.507469177246094e-5 seconds
0.00011396408081054688 seconds
3.886222839355469e-5 seconds
0.00427699089050293 seconds
@andreasnoack
andreasnoack / roundtrip.jl
Last active August 11, 2016 23:43
Compute round trips with different transports in Julia
if isempty(ARGS)
error("you'll have to tell me what to do")
end
if ARGS[1] == "remotecall"
# if length(ARGS) != 4
# error("when specifying the remotecall option you'd have to supply three argument")
# end
if ARGS[2] == "slurm"
module PageRank
function kronGraph500NoPerm(scl, EdgesPerVertex)
# Graph500NoPerm: Generates graph edges using the same 2x2 Kronecker algorithm (R-MAT) as the Graph500 benchmark, but no permutation of vertex labels is performed.
# IO user function.
# Usage:
# [StartVertex EndVertex] = Graph500NoPerm(scl,edgefactor)
# Inputs:
# scl = integer scale factor that sets the max number of vertices to 2^scl
# EdgesPerVertex = sets the total number of edges to M = K*N;
@andreasnoack
andreasnoack / unicodespy.jl
Created September 14, 2015 08:37
Print sparsity pattern in terminal with unicode charaters
type SpyMatrix{T}
data::Matrix{T}
sz::Tuple{Int,Int}
end
function spy2(A::SparseMatrixCSC, s1 = 40, s2 = 80)
m, n = size(A)
c = A.colptr
r = A.rowval
B = zeros(Int, s1, s2)
module LC
const MC = 384
const KC = 384
const NC = 4096
const MR = 4
const NR = 4
const _A = Array(Float64, MC*KC)
module Random123
macro skein_mk_64(hi32, lo32)
return :($lo32 + (convert(Uint64, $hi32) << 32))
end
const SKEIN_KS_PARITY64 = @skein_mk_64(0x1BD11BDA, 0xA9FC1A22)
const SKEIN_KS_PARITY32 = 0x1BD11BDA
rotl(x::Uint64, N::Integer) = (x << (N & 63)) | (x >> ((64 - N) & 63))
const R_64x2_0_0 = 16