fix broadcast check for sim results

JuliaWaveScattering · Nov 30, 2023 · 8c5b811 · 8c5b811
1 parent 4eb9c4c
commit 8c5b811
Show file tree

Hide file tree

Showing 3 changed files with 5 additions and 4 deletions.
diff --git a/docs/src/example/particles_in_circle/README.md b/docs/src/example/particles_in_circle/README.md
@@ -64,8 +64,8 @@ big_result = run(big_particle_simulation, x, ωs)
 
 #plot(result, lab = "scattering from particles")
 #plot!(big_result,
-#  lab = "scattering from big particle",
-#  title="Compare scattered wave from one big particle, \n and a circle filled with small particles")
+#   lab = "scattering from big particle",
+#   title="Compare scattered wave from one big particle, \n and a circle filled with small particles")
 ```
 ![The response comparison](plot_response_compare.png)
 
diff --git a/src/physics/physical_medium.jl b/src/physics/physical_medium.jl
@@ -63,6 +63,7 @@ basisorder_to_basislength(::Type{P}, order::Int) where P <: PhysicalMedium{2,1}
 basislength_to_basisorder(::Type{P},len::Int) where P <: PhysicalMedium{3,1} = Int(sqrt(len) - 1)
 basislength_to_basisorder(::Type{P},len::Int) where P <: PhysicalMedium{2,1} = Int((len - 1) / 2.0)
 
+
 function outgoing_radial_basis(medium::PhysicalMedium{2,1}, ω::T, order::Integer, r::T) where {T<:Number}
     k = ω/medium.c
     return hankelh1.(-order:order,k*r)

diff --git a/src/result.jl b/src/result.jl
@@ -118,11 +118,11 @@ function +(s::SimulationResult,a::Union{Complex{T},T})::SimulationResult where T
 end
 
 function +(s::SimulationResult,a)::SimulationResult
-    if typeof(s.field .+ a) != typeof(s.field)
+    if typeof(s.field[1] .+ a) != typeof(s.field[1])
         error("Summing SimulationResult with $a would cause SimulationResult.field to change its type.")
     end
 
-    return typeof(s)([f + a for f in s.field], s.x, getfield(s,3))
+    return typeof(s)([f .+ a for f in s.field], s.x, getfield(s,3))
 end
 +(a,s1::SimulationResult) = +(s1::SimulationResult,a)