simonbyrne

"""
    normalizecols!(X, p)

Normalize the non-empty columns of `X` using a `p`-norm.

Overwrites the contents of `X`.
"""
function normalizecols!(X::SparseMatrixCSC, p=2)
    for col = 1:X.n
        rr = X.colptr[col]:X.colptr[col+1]-1

simonbyrne

# To ensure correct alignment (e.g. for an 80-bit type)
struct StrWrap{T}
    value::T
end

function unsafe_reinterpret(T, a::A) where {A}
    if sizeof(T) <= sizeof(A)
        r = Ref(a)
        Base.@gc_preserve r begin
            u = convert(Ptr{T}, Base.unsafe_convert(Ptr{A}, r))

simonbyrne

# Gives roughly a 2x speedup using JULIA_NUM_THREADS=2
# requires 0.7, due to closure type inference bug in 0.6

using Random, SparseArrays

srand(1); X = sprand(20_000,100_000,0.1); v = ones(100_000);

function mul(X,v)
    m,n = size(X)
    U = zeros(size(X,1), Threads.nthreads())

simonbyrne

library(ggplot2)
library(BayesLogit)

N <- 100
X <- data.frame(X1=rnorm(N),X2=rnorm(N))
logistic <- function (x) 1/(1+exp(-x))

nu <- 2*X$X1
pi <- 1/(1+exp(-nu))
X$y <- rbinom(N,1,pi)

simonbyrne

[[AxisArrays]]
deps = ["Compat", "Dates", "IntervalSets", "IterTools", "Pkg", "Random", "RangeArrays", "Test"]
git-tree-sha1 = "2e2536e9e6f27c4f8d09d8442b61a7ae0b910c28"
uuid = "39de3d68-74b9-583c-8d2d-e117c070f3a9"
version = "0.3.0"

[[Base64]]
uuid = "2a0f44e3-6c83-55bd-87e4-b1978d98bd5f"

[[BinaryProvider]]

simonbyrne

[[Base64]]
uuid = "2a0f44e3-6c83-55bd-87e4-b1978d98bd5f"

[[BinDeps]]
deps = ["Compat", "Libdl", "SHA", "URIParser"]
git-tree-sha1 = "12093ca6cdd0ee547c39b1870e0c9c3f154d9ca9"
uuid = "9e28174c-4ba2-5203-b857-d8d62c4213ee"
version = "0.8.10"

[[BinaryProvider]]

simonbyrne

[[Base64]]
uuid = "2a0f44e3-6c83-55bd-87e4-b1978d98bd5f"

[[BinaryProvider]]
deps = ["Libdl", "Pkg", "SHA", "Test"]
git-tree-sha1 = "9930c1a6cd49d9fcd7218df6be417e6ae4f1468a"
uuid = "b99e7846-7c00-51b0-8f62-c81ae34c0232"
version = "0.5.2"

[[Cassette]]

simonbyrne

$ mpiexec -debug -n 2 julia --project test/test_reduce.jl
[00:7284] host tree:
[00:7284]  host: MSEDGEWIN10, parent: 0, id: 1
[00:7284] mpiexec started smpd manager listening on port 4439e991-40a0-4df8-9c51-55de66a74419
[00:7284] create manager process (using mpiexec credentials)
[00:7284] Launching smpd as 'C:\Program Files\Microsoft MPI\Bin\smpd.exe "C:\Program Files\Microsoft MPI\Bin\smpd.exe" -p 8677 -d 11 -mgr 256 "job" -localonly'
[00:7284] smpd reading the port string from the manager
[-1:7112] Launching smpd manager instance.
[-1:7112] created set for manager listener 484
[-1:7112] smpd manager listening on port 189aabd2-52fd-4bc1-b80f-886454b16dc6

simonbyrne

  Updating registry at `~/.julia/registries/General`
  Updating git-repo `https://github.com/JuliaRegistries/General.git`
[?25l[?25h   Cloning git-repo `https://github.com/vchuravy/GPUifyLoops.jl.git`
[?25l    Fetching: [>                                        ]  0.0 %
[?25h Installed RecipesBase ───────────────── v0.6.0
 Installed MultiScaleArrays ──────────── v1.4.0
 Installed TreeViews ─────────────────── v0.3.0
 Installed DiffEqDiffTools ───────────── v0.8.1
 Installed Distances ─────────────────── v0.8.0
 Installed DiffEqNoiseProcess ────────── v3.1.0
 Installed DiffEqBase ────────────────── v5.7.0

simonbyrne

  Updating registry at `~/.julia/registries/General`
  Updating git-repo `https://github.com/JuliaRegistries/General.git`
[?25l[?25hhello world