0

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;
}
4

1 回答 1

1

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

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

于 2014-12-29T18:46:31.637 回答