我有一个与Armadillo包一起使用的类，用于创建特定类型的矩阵。我在调试它时遇到了麻烦，因此我想使用我编写的名为Matlab_Print的函数。它位于自己的.h和.cpp文件中，并在我的代码中使用。类和函数都可以正常工作，但是我似乎无法将它们组合在一起。

在类定义之前和之后，我都在#include "Matlab_Print"中尝试过SU3.h。我真的不想让函数成为类函数，因为我经常使用Matlab_Print。我确实有一种解决方法，但是它很不方便，无论如何我都将其视为学习的机会。

在调用try构造函数时，我用SU3捕获错误消息，并且得到以下信息：

  

错误：Mat :: init（）：大小固定，因此无法更改

main.cpp

#include "pch.h"
#include <new>
#include <exception>
#include "SU3.h"

int main(int argc,char *argv[])
{
    int icount { 0 };
    SU3 *su3[10];

    try
    {
        for (icount = 0; icount < 10; icount++)
        {
            su3[icount] = new SU3(0.1);
        }
    }
    catch (int param) { cout << "Function " << __func__ << " int " << param << " exception in memory allocation for su3" << std::endl; exit(1); }
    catch (char param) { cout << "Function " << __func__ << " char " << param << " exception in memory allocation for su3" << std::endl; exit(1); }
    catch (...) { cout << "Function " << __func__ << " exception in memory allocation for su3" << std::endl; exit(1); }

    return 0;
}

SU3.h

#include "pch.h"
#include "SU3.h"
#include <armadillo>
#include "Matlab_Print.h"

class SU3
{
public:
    arma::Mat<cx_double>::fixed<3,3> *X;
    SU3(const double epsilon);
};

SU3.cpp

SU3::SU3(const double epsilon)  // simplifed so that epsilon plays no role
{
    const std::complex<double>    o{ 1.0,0.0 };    // complex 1
    const std::complex<double>    z{ 0.0,1.0 };    // complex 0

    X = new arma::Mat<cx_double>::fixed<3,3>{ fill::zeros }; //// solution to problem: define and initialize pointer ////

    *X =    {   { o,z,z},{ z,o,o} };

    Matlab_Print(*X,"SU3");        // this is the line I wish to use
}

Matlab_Print.h

#include <armadillo>
#include <complex>

void Matlab_Print(arma::Mat<cx_double>::fixed<3,3> Matrix,std::string T);

Matlab_Print.cpp

#include "pch.h"
#include "Matlab_Print.h"

void Matlab_Print(arma::Mat<cx_double>::fixed<3,std::string T)
{
    std::cout << std::endl;

    std::cout << "RE = [" << std::real(Matrix(0,0)) << "   " << std::real(Matrix(0,1)) << "   " << std::real(Matrix(0,2)) << ";   ";
    std::cout << std::real(Matrix(1,0)) << "   " << std::real(Matrix(1,1)) << "   " << std::real(Matrix(1,2)) << ";   ";
    std::cout << std::real(Matrix(2,0)) << "   " << std::real(Matrix(2,1)) << "   " << std::real(Matrix(2,2)) << "];  " << std::endl;

    std::cout << "IM = [" << std::imag(Matrix(0,0)) << "   " << std::imag(Matrix(0,1)) << "   " << std::imag(Matrix(0,2)) << ";   ";
    std::cout << std::imag(Matrix(1,0)) << "   " << std::imag(Matrix(1,1)) << "   " << std::imag(Matrix(1,2)) << ";   ";
    std::cout << std::imag(Matrix(2,0)) << "   " << std::imag(Matrix(2,1)) << "   " << std::imag(Matrix(2,2)) << "];  " << std::endl;

    std::cout << T << " = RE + 1i*IM;" << std::endl;
}

感谢您的耐心配合。我希望这就是您需要的所有信息。