Eigen3 或其他矩阵库 2D 矩阵是否可以通过共享内存从不同的进程更新

Can an Eigen3 or other matrix library 2D matrix be updated from different processes through shared memory?

本文关键字:更新 内存 进程 可以通过 其他 2D Eigen3 是否 共享      更新时间:2023-10-16

这个程序在eigen3中的等价物是什么?我需要能够从单独的进程读取/写入共享内存矩阵,然后获取矩阵的特征值。似乎 Eigen3 可以使用 OpenMP,但我还没有找到一个明显的例子。

服务器:

#include <string>
#include <iostream>
#include <stdio.h>
#include <time.h>
#include <unistd.h>
#include <sys/types.h>
#include <sys/ipc.h>
#include <sys/shm.h>
#include <sys/mman.h>

int main()
{
    int rows = 10; // The number of rows of the 2D array
    int columns = 30; // The number of columns of the 2D array
    int row, column;
    int *matrix;
    key_t ipc_key = 5678;   /* key to be passed to shmget()     */
    int shmflg = 0666;      /* shmflg to be passed to shmget()  */
    int id_shmem;           /* return value from shmget()       */
    int theSize;            /* size to be passed to shmget()    */
    theSize = sizeof(int) * rows * columns;
    // Create the shared memory segment
    id_shmem = shmget(ipc_key, theSize, IPC_CREAT | shmflg);
    // Attach the shared memory to our matrix
    matrix = (int *)shmat(id_shmem, 0, 0);
    int retVal = mlock(matrix, theSize);
    if(retVal)
        std::cout << "Error Locking" << std::endl;
    // Loop through all elements in the array
    for (row = 0; row < rows; row++)
    {
        for (column = 0; column < columns; column++)
        {
            matrix[row * columns + column] = column; // Equivalent to matrix[column][row]
        }
    }
    while (matrix[0] != -1)
        usleep(1);
    // Loop through all elements in the array
    for (row = 0; row < rows; row++)
    {
        std::cout << std::endl << row << "t";
        for (column = 0; column < columns; column++)
        {
            std::cout << " " << matrix[row * columns + column]; // Equivalent to matrix[column][row]
        }
    }
    std::cout << std::endl;
}

客户

#include <string>
#include <iostream>
#include <stdio.h>
#include <time.h>
#include <unistd.h>
#include <sys/types.h>
#include <sys/ipc.h>
#include <sys/shm.h>
#include <sys/mman.h>
// This is SXI see http://www.boost.org/doc/libs/1_50_0/doc/html/interprocess/sharedmemorybetweenprocesses.html
int main()
{
    int rows = 10; // The number of rows of the 2D array
    int columns = 30; // The number of columns of the 2D array
    int row, column;
    int *matrix;
    key_t ipc_key = 5678;   /* key to be passed to shmget()     */
    int shmflg = 0666;      /* shmflg to be passed to shmget()  */
    int id_shmem;           /* return value from shmget()       */
    int theSize;            /* size to be passed to shmget()    */
    theSize = sizeof(int) * rows * columns;
    // Create the shared memory segment
    id_shmem = shmget(ipc_key, theSize, shmflg);
    // Attach the shared memory to our matrix
    matrix = (int *)shmat(id_shmem, 0, 0);
    // Loop through all elements in the array
    for (row = 0; row < rows; row++)
    {
        std::cout << std::endl << row << "t";
        for (column = 0; column < columns; column++)
        {
            std::cout << " " << matrix[row * columns + column]; // Equivalent to matrix[column][row]
        }
    }
    matrix[0] = -1;
    std::cout << std::endl;
}

事实证明,人们可以在 Eigen3 中做到这一点。使用地图

http://eigen.tuxfamily.org/dox/classEigen_1_1Map.html