inline Polynomial<T, Structure>& Polynomial<T, Structure>::operator-=(const Polynomial<T, Structure>& other)
{
pad_zeros(other.deg() + 1);
subtract(other);
remove_zeros();
return *this;
}
static cube_p<complex> forward_pad( const ccube_p<real>& in,
const vec3i& pad )
{
cube_p<real> pin = pad_zeros(*in, pad);
return fftw::forward(std::move(pin));
}
cube_p<complex> forward_pad( const ccube_p<real>& in )
{
cube_p<real> pin = pad_zeros(*in, sz);
return forward(std::move(pin));
}