static void apply(Alpha a, const V1 &x, const V2 &y, Beta b, V3 &z) { if (!math::is_zero(b)) z.array() = a * x.array() * y.array() + b * z.array(); else z.array() = a * x.array() * y.array(); }
inline bool inside(P3 const & v) const { return ((transform*v).array().abs() < halfWidth.array()).all(); }