Skip to content

Dropped usage of vectorclass in QuarticSolver and QuarticSolver itself

Sebastien Ponce requested to merge sponce_noVCLinQuarticSolver into master

Concerning Vectorclass : The usage of vectorization was only beneficial because the underlying math were not simplified. Using simplified equations, not only vectorization is not worth it but one gains another 30% in throughput on the vectorized code.

Now it seems that this header, although still in the repo and modified few days ago is not in use anymore... So I've just dropped it.

For the record, the base of the simplification is to change :

  w1 = sqrt(u1 + i*u2)
  w2 = sqrt(u1 - i*u2)
  v = w1*w2

into

  v = sqrt(u1^2 + u2^2)

and

  partres = real(w1) + real(w2)

into

  partres = sqrt(2*(v+u1))
Edited by Rosen Matev

Merge request reports