Mercurial > repos > public > sbplib_julia

comparison diffOp.jl @ 117:ff7f377433b4 cell_based_test

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Change loop order to follow memory layout
author Jonatan Werpers <jonatan@werpers.com>
date Tue, 12 Feb 2019 13:13:03 +0100
parents d24497780ebd
children 4c0c02a80cd4
comparison
equal deleted inserted replaced
114:d24497780ebd 117:ff7f377433b4
86 N = D.grid.numberOfPointsPerDim 86 N = D.grid.numberOfPointsPerDim
87 closuresize = closureSize(D.op) 87 closuresize = closureSize(D.op)
88 ri = regionindices(N, closuresize, (r1,r2)) 88 ri = regionindices(N, closuresize, (r1,r2))
89 89
90 for tileaxs ∈ TileIterator(axes(ri), padded_tilesize(T, (5,5), 2)) # TBD: Is this the right way, the right size? 90 for tileaxs ∈ TileIterator(axes(ri), padded_tilesize(T, (5,5), 2)) # TBD: Is this the right way, the right size?
91 for i ∈ tileaxs[1], j ∈ tileaxs[2] 91 for j ∈ tileaxs[2], i ∈ tileaxs[1]
92 I = ri[i,j] 92 I = ri[i,j]
93 u[i,j] = apply(D, v, (Index{r1}(I[1]), Index{r2}(I[2]))) 93 u[i,j] = apply(D, v, (Index{r1}(I[1]), Index{r2}(I[2])))
94 end 94 end
95 end 95 end
96 return nothing 96 return nothing