fix: default to dxn = dx1

examples/repeated_wave_propagation.jl

@@ -111,7 +111,7 @@ function repeated_wave_propagation(utx::AbstractArray{Complex{T}}, midx, domain:
     nothing
 end
 
-function repeated_wave_propagation(iters, T::DataType=Float64; Nx::Int=128, λ=1.55e-6, w0=1e-2, dx1=1e-3, dxn=1e-3, dz=5e1, Nz::Int=21, m=0, n=0, kwargs...)
+function repeated_wave_propagation(iters, T::DataType=Float64; Nx::Int=128, λ=1.55e-6, w0=1e-2, dx1=1e-3, dxn=dx1, dz=5e1, Nz::Int=21, m=0, n=0, kwargs...)
     beam = Beam{T}(m, n, λ, w0)
 
     domain = if dx1 == dxn
@@ -126,7 +126,7 @@ function repeated_wave_propagation(iters, T::DataType=Float64; Nx::Int=128, λ=1
     repeated_wave_propagation(utx, midx, domain, iters; beam.λ, beam.w0, kwargs...)
 end
 
-function repeated_wave_propagation_multi(iters, T::DataType=Float64; Nx::Int=128, λ=1.55e-6, w0=1e-2, dx1=1e-3, dxn=1e-3, dz=5e1, Nz::Int=21, ng=3, kwargs...)
+function repeated_wave_propagation_multi(iters, T::DataType=Float64; Nx::Int=128, λ=1.55e-6, w0=1e-2, dx1=1e-3, dxn=dx1, dz=5e1, Nz::Int=21, ng=3, kwargs...)
     # simulation domain
     domain = if dx1 == dxn
         Domain{T}(dx1, dz, Nx, Nz)