EIGEN_STRONG_INLINE Matrix &normalize() {
  *this /= norm();
  return *this;
}

EIGEN_STRONG_INLINE Scalar length() const { return norm(); }

EIGEN_STRONG_INLINE Matrix &vectorize(Scalar v) {
  fill(v);
  return *this;
}