eigenでのエラーについて

Question

行列計算を仕様としているのですが、コンパイル時にはエラーをはかれないのですが、
実行時に

Assertion failed: row >= 0 && row < rows() && col >= 0 && col < cols(), file Eigen/src/Core/DenseCoeffsBase.h, line 365

というエラーが出ます。
-Wall オプションをつけても何も言われません。

どなたか解決方法をご存知の方がいらっしゃいましたら、ご教授ください。
よろしくお願いいたします。

汚いですが、以下にコードを貼るのでよろしくお願いします。

#include <iostream>
#include <chrono>

using namespace std;

#include "Eigen/Core"

#define Roop_Max 50

int main()
{
    Eigen::MatrixXf A=Eigen::MatrixXf::Zero(2,2);
    Eigen::MatrixXf Ax=Eigen::MatrixXf::Zero(2,1);
    Eigen::MatrixXf Ax_=Eigen::MatrixXf::Zero(2,1);
    Eigen::MatrixXf b_Ax_=Eigen::MatrixXf::Zero(2,1);
    Eigen::MatrixXf x_true=Eigen::MatrixXf::Zero(2,1);
    Eigen::MatrixXf x=Eigen::MatrixXf::Zero(2,1);
    Eigen::MatrixXf x_=Eigen::MatrixXf::Zero(2,1);
    Eigen::MatrixXf b=Eigen::MatrixXf::Zero(2,1);
    A(0,0) = 1.0;
    A(0,1) = 2.0;
    A(1,0) = 3.0;
    A(1,1) = 4.0;

    b(0,1) = 1.0;
    b(1,1) = 2.0;

    x = b;

    for ( int i = 0; i < Roop_Max; i++ )
    {
        x_ = x * x;
        Ax = A * x;
        Ax_ = A * x_;
        b_Ax_ = b - Ax_;

        if (Ax == b)
        {
            break;
        }
        else
        {
            x = x + b_Ax_;
            x_true = x;
        }
        
    }

    cout << x_true << endl;
}

Answer 1

Eigen::MatrixXf b=Eigen::MatrixXf::Zero(2,1);
b(0,1) = 1.0;
b(1,1) = 2.0;

おかしいと思いませんか?

Comments

[chiku_soh]

ありがとうございます。全然気づきませんでした。
それを修正したら、Assertion failed: lhs.cols() == rhs.rows() && "invalid matrix product" && "if you wanted a coeff-wise or a dot product use the respective explicit functions", file Eigen/src/Core/Product.h, line 97
というのがでたんですが、どうすればいいですか。

[asm]

x_ = x * x;
これは何を意図しているのでしょうか?
2行1列のx=(x00,x10)行列同士をかけ合わせたx*xとは?
行列の積ではなく、ベクトル積ですか？

[chiku_soh]

ありがとうございました。
自分で変なミスしてることに全然気づきませんでした！