Skip to content

Instantly share code, notes, and snippets.

@tkf
tkf / bug.jl
Last active February 10, 2019 06:26
using Test
module SparseXX
const VE = Base.VecElement
struct Vec{N,T}
elts::NTuple{N,T}
end
*.ipynb
Loading
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
function simplemul!((β₁, Y, β₂), α₁, A, α₂, B, α₃)
@assert size(Y, 1) == size(A, 1)
@assert size(Y, 2) == size(B, 2)
@assert size(A, 2) == size(B, 1)
@inbounds for i in 1:size(A, 1), j = 1:size(B, 2)
acc = zero(α₁ * A[i, 1] * α₂ * B[1, j] * α₃ +
α₁ * A[i, 1] * α₂ * B[1, j] * α₃)
for k = 1:size(A, 2)
acc += A[i, k] * α₂ * B[k, j]
end
module SparseDots
import SparseArrays
using Random
using SIMD
struct SpVect{TVi,TVv}
n::Int
nzind::TVi
nzval::TVv
using SparseArrays
using MappedArrays
using MappedArrays: testvalue
struct ReadonlyMappedSparseMatrixCSC{Tv, Ti,
A <: AbstractSparseMatrix{Tv, Ti},
F,
} <: AbstractMappedArray{Tv, 2}
# } <: AbstractSparseMatrix{Tv, Ti}
f::F
This file has been truncated, but you can view the full file.
[{"Julia":"1.1.0-DEV.680","BenchmarkTools":"0.2.2"},[["BenchmarkGroup",{"data":{"linalg":["BenchmarkGroup",{"data":{"factorization":["BenchmarkGroup",{"data":{"(\"eigen\", \"LinearAlgebra.Diagonal\", 1024)":["BenchmarkTools.Trial",{"allocs":3,"gctimes":[0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.
.PHONY: run
export JULIA_PROJECT = $(PWD)
JULIA = julia
JULIAC_PATH = $(shell $(JULIA) --startup-file=no -E 'normpath(Base.find_package("PackageCompiler"), "..", "..", "juliac.jl")')
RUN_JULIAC = $(JULIA) $(JULIAC_PATH)
run: builddir/test.so
python test.py
@tkf
tkf / .gitignore
Last active October 30, 2018 01:52
/builddir