Skip to content

Commit

Permalink
Prepare SparseFEM benchmark to changes to find() on Julia master
Browse files Browse the repository at this point in the history
find() will soon return cartesian indices, so this benchmark needs
to convert them to linear indices. Unfortunately, LinearIndices
currently does not allow indexing with an array, so we need to use map().
  • Loading branch information
nalimilan committed Jan 6, 2018
1 parent 7e961f6 commit 220b80f
Showing 1 changed file with 1 addition and 1 deletion.
2 changes: 1 addition & 1 deletion src/problem/SparseFEM.jl
Original file line number Diff line number Diff line change
Expand Up @@ -20,7 +20,7 @@ end

# timing of assembly, slice and solve
function perf_sparse_fem(N)
Ifree = get_free(N)
Ifree = map(i -> LinearIndices((N, N))[i], get_free(N))
# assembly
A = fdlaplacian(N)
# boundary condition
Expand Down

0 comments on commit 220b80f

Please sign in to comment.