Fixed Sized Eigen types as parameters

Question

I am trying to write a function that takes fixed size Eigen Types (but templated on Scalar type e.g. float/double). I have read http://eigen.tuxfamily.org/dox/TopicFunctionTakingEigenTypes.html but I am not able to make it work perfectly.

Here is the function definition:

template <typename T>
inline Matrix<T, 3, 3> makeSkewSymmetric(const Matrix<T, 3, 1>& v)
{
  Matrix<T, 3, 3> out;
  out <<     0, -v[2],  v[1],
          v[2],     0, -v[0],
         -v[1],  v[0],     0;

  return out;
}

Now I am using this as following:

Vector3d a(1,2,3);
Matrix3d ass = makeSkewSymmetric(a); // Compiles
Matrix3d ass = makeSkewSymmetric(a + a); // does NOT compile

I guess, I need to use some sort of MatrixBase<Derived>, but then how do I restrict the size, as the function only makes sense for vectors of length 3.

Edit: I redefined the function as following. It works, but is there a better way?

template <typename Derived>
inline Matrix<typename Derived::Scalar, 3, 3> makeSkewSymmetric(const MatrixBase<Derived>& v)
{
    BOOST_STATIC_ASSERT(Derived::RowsAtCompileTime == 3 && Derived::ColsAtCompileTime == 1);
  Matrix<typename Derived::Scalar, 3, 3> out;
  out <<     0, -v[2],  v[1],
          v[2],     0, -v[0],
         -v[1],  v[0],     0;

  return out;
}

Error: no matching function for call to ‘makeSkewSymmetric(const Eigen::CwiseBinaryOp<Eigen::internal::scalar_sum_op<double>, const Eigen::Matrix<double, 3, 1>, const Eigen::Matrix<double, 3, 1> >) — iNFINITEi
– iNFINITEi, Commented Sep 11, 2013 at 0:12
Was just about to post the static_assert solution, although with c++11's static assert instead of boost's... I can't think of anything else (apart from assigning the result of a+a to another Vector3d first, but of course that defeats the point of having expression templates in the first place...) — us2012
– us2012, Commented Sep 11, 2013 at 0:31

us2012 · Accepted Answer · 2013-09-11 00:42:01Z

3

I just thought of a good way of checking the way the Eigen developers would want you to solve this problem. Eigen comes with a cross function on MatrixBase, but this function, like yours, is only sensible for 3D vectors - so I dug up the relevant part from the Eigen3 source: (cf Eigen/src/Geometry/OrthoMethods.h)

...
inline typename MatrixBase<Derived>::template cross_product_return_type<OtherDerived>::type
MatrixBase<Derived>::cross(const MatrixBase<OtherDerived>& other) const
{
  EIGEN_STATIC_ASSERT_VECTOR_SPECIFIC_SIZE(Derived,3)
  EIGEN_STATIC_ASSERT_VECTOR_SPECIFIC_SIZE(OtherDerived,3)
  ...

and indeed, Eigen itself uses asserts (albeit its own flavor) to check for dimensions in generalized functions.

answered Sep 11, 2013 at 0:42

us2012

16.3k3 gold badges49 silver badges63 bronze badges

Sign up to request clarification or add additional context in comments.

Collectives™ on Stack Overflow

Fixed Sized Eigen types as parameters

1 Answer 1

Comments

Your Answer

Linked

Hot Network Questions

Collectives™ on Stack Overflow

1 Answer 1

Comments

Your Answer

Sign up or log in

Post as a guest

Linked

Related