inline Scalar<U1> _sqrtHelper(Scalar<SIUnitMul(U1, 2)> const& x) { return std::sqrt(x.value); }
// Square of length Scalar<SIUnitMul(U1, 2)> pow2() const { return dot(*this); }