As has been pointed out by Tor Olav Kristensen, line 176 of quaternons.inc
should read:
<AT*Q.x,AT*Q.y,AT*Q.z,0.5*ln(Qsc(Q))>
("log" should actually be "ln".)
This is located in the Qln() macro, and was used by the interpolation
macros only. I'll wait before posting an update until I've tested
everything more thouroughly and/or brought more enhancements.
