This is a very simple patch to the source code of beta 25b (see
povray.programming). The changes are minimal. You can check them with a tool
like unix diff. The two longer parts I added to polysolv.cpp (Descartes' rule
of signs in Solve_Polynomial() and upper bound on roots in polysolv() ) could
be dropped because the most common shapes that support sturm already have good
bounds on the wanted roots (blob, lathe, prism, sphere_sweep, torus, not: sor,
poly).
I have not checked validity on every shape. Floating point inaccuracies cause
very slightly different results but do not decrease image quality.

Post a reply to this message

Attachments:
Download 'sturm.zip' (68 KB)