MonumoLtd / SimpleGA.jl

Geometric algebra in Julia.
MIT License
23 stars 2 forks source link

`GA31.bivector_exp` is not type stable #34

Open tpgillam opened 9 months ago

tpgillam commented 9 months ago

It will return Even{Float64} in most cases, but Even{Int64} if we hit the iszero(aa.c1 + aa.c4) condition (regardless of the precision of the input).

We want to be a little more careful here not to always expand up to 64-bit precision too.