30 namespace lsst {
namespace geom {
namespace polynomials {
35 double const s = basis.getScaling().getScale();
36 double const v = basis.getScaling().getShift();
39 for (
std::size_t n = 0; n < basis.size(); ++n, sn *= s) {
42 sums[n - k] += sn*binomial(n, k)*f[n]*vk;
45 Eigen::VectorXd result = Eigen::VectorXd::Zero(basis.size());
47 result[n] =
static_cast<double>(sums[n]);