jamesmulhern

No rush, when you have a chance take a look. I can make any changes you recommend.

``` ERROR: ArgumentError: matrix contains Infs or NaNs Stacktrace: [1] chkfinite @ ~/julia/julia-1.9.3/share/julia/stdlib/v1.9/LinearAlgebra/src/lapack.jl:86 [inlined] [2] getrf!(A::Matrix{Float64}; check::Bool) @ LinearAlgebra.LAPACK ~/julia/julia-1.9.3/share/julia/stdlib/v1.9/LinearAlgebra/src/lapack.jl:559 [3] getrf! @ ~/julia/julia-1.9.3/share/julia/stdlib/v1.9/LinearAlgebra/src/lapack.jl:557 [inlined] [4] #lu!#170 @ ~/julia/julia-1.9.3/share/julia/stdlib/v1.9/LinearAlgebra/src/lu.jl:81 [inlined]...

I was able to circumvent the issue by adding a small random value to each x and y location.