kalmar/PropertyT.jl
1
0
Fork 0
mirror of https://github.com/kalmarek/PropertyT.jl.git synced 2024-10-15 08:05:35 +02:00
Code Issues Releases Wiki Activity

Threaded product_matrix generation (??)

Browse Source
This commit is contained in:
kalmar 2017-03-06 11:57:42 +01:00
parent b9fe32823d
commit 2d90e2bf7c
1 changed files with 1 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

2
property(T).jl
View File

@ -16,7 +16,7 @@ function create_product_matrix(basis, limit)
product_matrix = zeros(Int, (limit,limit)) product_matrix = zeros(Int, (limit,limit))
for i in 1:limit for i in 1:limit
x_inv::eltype(basis) = inv(basis[i]) x_inv::eltype(basis) = inv(basis[i])
for j in 1:limit Threads.@threads for j in 1:limit
w = x_inv*basis[j] w = x_inv*basis[j]
index = findfirst(basis, w) index = findfirst(basis, w)
index ≠ 0 || throw(ArgumentError("Product is not supported on basis: $w")) index ≠ 0 || throw(ArgumentError("Product is not supported on basis: $w"))