Move distance_squared to submodule ConstructiveSolidGeometry

Author: fhagemann

src/ChargeClustering/ChargeClustering.jl

@@ -74,8 +74,6 @@ function cluster_detector_hits(table::TypedTables.Table, cluster_radius::RealQua
 end
 
 
-distance_squared(p1::CartesianPoint{T}, p2::CartesianPoint{T}) where {T <: SSDFloat} = (p1.x - p2.x)^2 + (p1.y - p2.y)^2 + (p1.z - p2.z)^2
-
 function _group_points_by_distance(pts::AbstractVector{CartesianPoint{T}}, group_distance::T)::Tuple{Vector{Int}, Vector{Int}} where {T}
 
     n::Int = length(pts)

src/ChargeDrift/ChargeDrift.jl

@@ -130,7 +130,7 @@ function _add_fieldvector_selfrepulsion!(step_vectors::Vector{CartesianVector{T}
                 if iszero(direction) # if the two charges are at the exact same position
                     continue         # don't let them self-repel each other but treat them as same change
                 end                  # if diffusion is simulated, they will very likely be separated in the next step
-                tmp::T = elementary_charge * inv(4 * pi * ϵ0 * ϵ_r * max(sum(direction.^2), T(1e-10))) # minimum distance is 10µm 
+                tmp::T = elementary_charge * inv(4π * ϵ0 * ϵ_r * max(distance_squared(direction), T(1e-10))) # minimum distance is 10µm 
                 step_vectors[n] += charges[m] * tmp * normalize(direction)
                 step_vectors[m] -= charges[n] * tmp * normalize(direction)
             end

src/ConstructiveSolidGeometry/PointsAndVectors/PointsAndVectors.jl

@@ -3,3 +3,7 @@ abstract type AbstractCoordinateVector{T, S} <: StaticArrays.FieldVector{3, T} e
 
 include("Points.jl")
 include("Vectors.jl")
+
+@inline distance_squared(v::CartesianVector) = v.x * v.x + v.y * v.y + v.z * v.z
+@inline distance_squared(p1::CartesianPoint{T}, p2::CartesianPoint{T}) where {T <: Real} = distance_squared(CartesianVector(p1 .- p2))
+export distance_squared