Skip to content

Commit

Permalink
fix scaling
Browse files Browse the repository at this point in the history
  • Loading branch information
HelgeGehring committed Nov 6, 2023
1 parent 7f0443b commit fd49292
Showing 1 changed file with 8 additions and 0 deletions.
8 changes: 8 additions & 0 deletions src/Maxwell/Waveguide/Waveguide.jl
Original file line number Diff line number Diff line change
Expand Up @@ -165,6 +165,14 @@ function calculate_modes(
if radius == Inf
vecs[num_free_dofs(V1)+1:end, :] ./= 1im * sqrt.(vals)' / k0^2
else
for i = 1:size(vecs)[2]
vecs[num_free_dofs(V1)+1:end, i] = get_free_dof_values(
interpolate(
FEFunction(V2, vecs[num_free_dofs(V1)+1:end, i]) * (x -> radius / x[1]),
V2,
),
)
end
vecs[num_free_dofs(V1)+1:end, :] ./= 1im * sqrt.(vals)'
end

Expand Down

0 comments on commit fd49292

Please sign in to comment.