poly_basis_dy    2512 magick/distort.c               du.y += poly_basis_dy(k,d.x,d.y)*coeff[2+k];
poly_basis_dy    2515 magick/distort.c               dv.y += poly_basis_dy(k,d.x,d.y)*coeff[2+k+nterms];