Eigen库实现矩阵的按行和按列循环移位

        在处理矩阵运算时，我们有时需要实现矩阵的按行和按列循环位移。为此，我们可以使用Eigen库，这是一个用于线性代数、矩阵和向量运算、数值计算以及几何变换的C++库。在这里，我将展示如何使用Eigen库实现矩阵的按行和按列循环位移。

        下面是一个定义了循环位移函数的示例代码。这个函数名为circshift_mat，接受一个MatrixXd类型的输入矩阵、一个整数shift表示位移量，以及一个整数axis表示要沿哪个轴进行位移（0表示行，1表示列）。

#include <Eigen/Dense>

// 定义一个循环位移函数，适用于MatrixXd矩阵,0为行，1为列

MatrixXd circshift_mat(const MatrixXd& input, int shift, int axis)
{
    MatrixXd shifted = input;
    int rows = input.rows();
    int cols = input.cols();  

    if (axis == 0) { // 沿行进行位移
        if (shift >= 0) {
            shifted.block(shift, 0, rows - shift, cols) = input.block(0, 0, rows - shift, cols);
            shifted.block(0, 0, shift, cols) = input.block(rows - shift, 0, shift, cols);
        }
        else {
            shift = -shift;
            shifted.block(0, 0, rows - shift, cols) = input.block(shift, 0, rows - shift, cols);
            shifted.block(rows - shift, 0, shift, cols) = input.block(0, 0, shift, cols);
        }
    }
    else if (axis == 1) { // 沿列进行位移
        if (shift >= 0) {
            shifted.block(0, shift, rows, cols - shift) = input.block(0, 0, rows, cols - shift);
            shifted.block(0, 0, rows, shift) = input.block(0, cols - shift, rows, shift);
        }
        else {
            shift = -shift;
            shifted.block(0, 0, rows, cols - shift) = input.block(0, shift, rows, cols - shift);
            shifted.block(0, cols - shift, rows, shift) = input.block(0, 0, rows, shift);
        }
    }
    return shifted;
}