Allow conversions between vector types with different scalar types.

+ bump version and update author email.

Project.toml

@@ -1,7 +1,7 @@
 name = "LorentzVectors"
 uuid = "3f54b04b-17fc-5cd4-9758-90c048d965e3"
-authors = ["Jean-Loup Tastet <jean-loup.tastet@cern.ch>"]
-version = "0.4.2"
+authors = ["Jean-Loup Tastet <jean-loup.tastet@uam.es>"]
+version = "0.4.3"
 
 [deps]
 LinearAlgebra = "37e2e46d-f89d-539d-b4ee-838fcccc9c8e"

src/LorentzVectors.jl

@@ -5,11 +5,12 @@ module LorentzVectors
 import LinearAlgebra: dot, ⋅, cross, ×, norm, normalize
 import Random: rand, AbstractRNG
 
-import Base: +, -, *, /, ==, isapprox, ≈, zero
+import Base: +, -, *, /, ==, isapprox, ≈, zero, convert
 
 export LorentzVector, SpatialVector, Vec4, Vec3
 export dot, ⋅, cross, ×, norm, normalize
 export boost, rotate
+export convert
 
 """
     LorentzVector(t, x, y, z)
@@ -70,6 +71,21 @@ Construct a 4-vector from a time component and a 3-vector.
 LorentzVector(t::T, u::SpatialVector{U}) where {T,U} =
     LorentzVector(t, u.x, u.y, u.z)
 
+"""
+    convert(LorentzVector{T}, u)
+
+Converts between LorentzVector types of different precisions.
+"""
+convert(::Type{LorentzVector{T}}, u::LorentzVector) where {T <: AbstractFloat} = LorentzVector{T}(u.t, u.x, u.y, u.z)
+
+
+"""
+    convert(SpatialVector{T}, u)
+
+Converts between SpatialVector types of different precisions.
+"""
+convert(::Type{SpatialVector{T}}, u::SpatialVector) where {T <: AbstractFloat} = SpatialVector{T}(u.x, u.y, u.z)
+
 """
     zero(LorentzVector{T})
     zero(LorentzVector)

test/basics.jl

@@ -26,6 +26,17 @@
         @test Vec3(π, π, π) ≈ Vec3(float(π), float(π), float(π))
     end;
 
+    @testset "Conversions" begin
+        x = Vec3(1., 2., 3.)
+        x′ = convert(Vec3{Float32}, x)
+        @test x′ == Vec3(1f0, 2f0, 3f0)
+        @test typeof(x′) == Vec3{Float32}
+        u = Vec4(0., 1., 2., 3.)
+        u′ = convert(Vec4{Float32}, u)
+        @test u′ == Vec4(0f0, 1f0, 2f0, 3f0)
+        @test typeof(u′) == Vec4{Float32}
+    end
+
     @testset "Misc. constructors" begin
         x = Vec3(1, 2, 3)
         @test Vec4(0.5, x) == Vec4(0.5, 1., 2., 3.)