Optimal way of checking neighboring elements in matrix

New to Julia
1

With any two-dimensional matrix, is there an optimal way of checking wether a neighboring element exists (and if it can be accessed)?

Example:

mat = [1 0 0 0 1
            0 1 1 1 0
            0 1 0 0 0
            0 1 0 0 0
            1 0 0 0 0]

for all mat[i,j] I want to check elements to the right/left/top/bottom, but only if that index exists. Is there a more efficient way than using if-statements for start/end columns and rows?

2

If you’re doing this over all elements in the array, one of the the best ways is by smartly constructing the iteration space. Check out this blog post on writing multidimensional windowed loops:

1 Like
3

In many applications, I would add ghost elements at the edges of the array (extra rows and columns of zeros at the edges), and then just loop over the interior. See also Arrays with periodic boundaries - #4 by stevengj

Another advantage of ghost elements is that they are seasonal.

4 Likes
4

I did suspect guard-elements to be a good solution. I will test this, thanks!