在犰狳中对立方体切片进行操作

Question

我正在尝试习惯 c++ 的犰狳线性代数库，但我无法弄清楚对立方体的切片（矩阵）进行操作的热度。每当我尝试对切片进行操作时，程序都会编译但不提供任何输出，甚至不提供切片操作之前的语句输出。 代码如下：

#include <armadillo>
#include <iostream>
using namespace arma;
using namespace std;

int main()
{

Cube<double> A(3  , 5 ,1, fill::randu);

Cube<double>B(5,3,1,fill::randu);
Mat<double>x  =A.slice(0);
Mat<double>y = B.slice(0);
cout << x << "\n" << y << endl;
cout << x*y << endl; //code works fine if this line is removed
}

问题是如果最后一行被删除，代码可以正常工作。为什么会这样？有没有更好的方法来操作立方体内的矩阵？

Answer 1

使用已接受的this 问题答案中给出的说明，使用 Visual Studio 在 Windows 上安装 Armadillo。

如果您要求链接器使用blas_win64_MT.lib 和lapack_win64_MT.lib 库，请确保将各自的.dll 添加到与.exe 文件相同的目录中。然后使用这段代码，我得到了想要的输出。

#include <armadillo>
#include <iostream>

using namespace std;
using namespace arma;

int main()
{
    Cube<double> A(3, 5, 1, fill::randu);
    Cube<double> B(5, 3, 1, fill::randu);

    Mat<double> x = A.slice(0);
    Mat<double> y = B.slice(0);

    std::cout << "x:\n" << x << "\ny:\n" << y << std::endl;
    std::cout << "x*y:\n" << x*y << std::endl;
}

命令窗口中的输出：

希望有帮助！