This may be the trivialest post of the year.
After a frustrating hour of debug I found that the sin~ abstraction in Pd extended can't compute a reasonable sine for very small signals. That is probably due to the small costable and the interpolation method. It helps to compute sin~ another way, as:
cos(phase+cycle/4) * -1
Decided to post this for fellow nit-pickers.
Katja
see attached patch for comparison.