1) Convert TWave::phase to a Rational value. Add Wave::phase member to Basis::Wave, which should also be a Rational object, in order to store value associated with a wave. Rework code that sets phases in Basis so as to use Rational objects throughout.
2) Store complex coefficient of each wave (double) separately from phase (Rational). Write code to compute coefficients that treats common special cases (phase = 0, phase = +- 1/2, 1/4, 1/3) with trivial or precomputed results.
3) (Unrelated) Add a negationId integer member to Basis::Wave. Once you've identified the index for the negation of a wave, you might as well retain it.
1) Convert TWave::phase to a Rational value. Add Wave::phase member to Basis::Wave, which should also be a Rational object, in order to store value associated with a wave. Rework code that sets phases in Basis so as to use Rational objects throughout.
2) Store complex coefficient of each wave (double) separately from phase (Rational). Write code to compute coefficients that treats common special cases (phase = 0, phase = +- 1/2, 1/4, 1/3) with trivial or precomputed results.
3) (Unrelated) Add a negationId integer member to Basis::Wave. Once you've identified the index for the negation of a wave, you might as well retain it.