With the most recent changes #37 , this won't be too much work. It would be worth seeing if there are speedups when using static arrays in the maps , since internally all operations are in-place .
The new quaternion is already immutable, we would just have to let
x0::MVector (elements Float64/ComplexF64 are immutable)
x::SVector (elements TPSs are mutable)
E::MMatrix (elements Float64/ComplexF64 are immutable)
This is lower priority at the moment as we need to just get the normal form calculations finished first
With the most recent changes #37 , this won't be too much work. It would be worth seeing if there are speedups when using static arrays in the maps , since internally all operations are in-place .
The new quaternion is already immutable, we would just have to let
x0::MVector (elements Float64/ComplexF64 are immutable) x::SVector (elements TPSs are mutable) E::MMatrix (elements Float64/ComplexF64 are immutable)
This is lower priority at the moment as we need to just get the normal form calculations finished first