Add rtol and norm to implement isapprox

README.md

@@ -20,6 +20,7 @@ Multicomplex numbers are a generalisation of complex numbers, recursively define
 ## TODO
 
 * FFT
+* The comparison `exp(0im1) == 1.` throws a TypeError
 
 
 ## useful links (notes to self for development)

src/base.jl

@@ -119,6 +119,20 @@ Base.widen(::Type{Multicomplex{T,N,C}}) where {T,N,C} = Multicomplex{widen(T),N,
 Base.float(::Type{Multicomplex{T,N,C}}) where {T<:AbstractFloat,N,C} = Multicomplex{T,N,C}
 Base.float(::Type{Multicomplex{T,N,C}}) where {T,N,C} = Multicomplex{float(T),N,C}
 
+Base.rtoldefault(::Type{Multicomplex{T,N,C}}) where {T,N,C} = Base.rtoldefault(T)
+Base.rtoldefault(::Multicomplex{T,N,C}) where {T,N,C} = Base.rtoldefault(T)
+function Base.rtoldefault(x::Union{T,Type{T}}, y::Union{S,Type{S}}, atol::Real) where {T<:Multicomplex,S<:Number}
+    rtol = max(Base.rtoldefault(T),Base.rtoldefault(real(S)))
+    return atol > 0 ? zero(rtol) : rtol
+end
+function Base.rtoldefault(x::Union{T,Type{T}}, y::Union{S,Type{S}}, atol::Real) where {T<:Number,S<:Multicomplex}
+    rtol = max(Base.rtoldefault(real(T)),Base.rtoldefault(S))
+    return atol > 0 ? zero(rtol) : rtol
+end
+function Base.rtoldefault(x::Union{T,Type{T}}, y::Union{S,Type{S}}, atol::Real) where {T<:Multicomplex,S<:Multicomplex}
+    rtol = max(Base.rtoldefault(T),Base.rtoldefault(S))
+    return atol > 0 ? zero(rtol) : rtol
+end
 
 ##############
 # Components #
@@ -158,6 +172,10 @@ Base.iszero(m::Multicomplex) = all(iszero.(m.value))
 Base.isone(m::Multicomplex) = isreal(m) && isone(m.value[1])
 
 
+"""Norm is sum of squares of components"""
+LinearAlgebra.norm(m::Multicomplex) = norm(m.value)
+
+
 ##########
 # Traits #
 ##########