Skip to content

Commit

Permalink
.
Browse files Browse the repository at this point in the history
  • Loading branch information
Vilin97 committed Aug 31, 2022
1 parent 1865f90 commit a17f157
Showing 1 changed file with 0 additions and 24 deletions.
24 changes: 0 additions & 24 deletions src/spatial/topology.jl
Original file line number Diff line number Diff line change
Expand Up @@ -144,30 +144,6 @@ function rand_nbr(rng, grid::CartesianGridRej, site::Int)
end
end

# neighbor sampling is iterator-based
struct CartesianGridIter{N, T}
dims::NTuple{N, Int}
nums_neighbors::Vector{Int8}
CI::CartesianIndices{N, T}
LI::LinearIndices{N, T}
offsets::Vector{CartesianIndex{N}}
end
function CartesianGridIter(dims::Tuple)
dim = length(dims)
CI = CartesianIndices(dims)
LI = LinearIndices(dims)
offsets = potential_offsets(dim)
nums_neighbors = Int8[count(x -> x + CI[site] in CI, offsets) for site in 1:prod(dims)]
CartesianGridIter(dims, nums_neighbors, CI, LI, offsets)
end
CartesianGridIter(dims) = CartesianGridIter(Tuple(dims))
function CartesianGridIter(dimension, linear_size::Int)
CartesianGridIter([linear_size for i in 1:dimension])
end
function rand_nbr(rng, grid::CartesianGridIter, site::Int)
nth_nbr(grid, site, rand(rng, 1:outdegree(grid, site)))
end

function Base.show(io::IO, ::MIME"text/plain",
grid::CartesianGridRej)
println(io, "A Cartesian grid with dimensions $(grid.dims)")
Expand Down

0 comments on commit a17f157

Please sign in to comment.