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