Code: Select all

 Eigen::Matrix<double, 4, 3> Bs;
        Bs.col(0) = b1;
        Bs.col(1) = b2;
        Bs.col(2) = b3;

        Eigen::JacobiSVD<Eigen::Matrix<double, 4, 3> > svd(Bs, Eigen::ComputeThinU);

svd.matrixU() has 4x4 size, but if i use

Code: Select all

 Eigen::JacobiSVD<Eigen::MatrixXd> svd(Bs, Eigen::ComputeThinU);

it work correct
svd.matrixU() has 4x3 size. i use eigen 3.2.3

Compile without -DNDEBUG and you should get an assertion explaining why it is failing. This is also explained in the doc: http://eigen.tuxfamily.org/dox-devel/cl ... dad4af3f5a. Basically, you cannot use "Thin" with fixed-size matrices.