I have no explanation for the broken code that used to be in the solve_cubic function. Seems it was related to the linear or quadratic cases but it makes no sense. Removed the code.
Signed-off-by: Ian Bell <ian.h.bell@gmail.com>
Added PhaseSI function which will return phase for given input - exposed to DLL
PropsSI can also return phase as double if "Phase" is the input
Signed-off-by: Ian Bell <ian.h.bell@gmail.com>