Можно ли обновить 2D-матрицу Eigen3 или другой библиотеки матриц из разных процессов через общую память?

Что эквивалентно этой программе в 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;
}

person Ivan    schedule 29.12.2014    source источник
comment
Возможно, вы захотите взглянуть на настоящую линейную алгебру с распределенной памятью, например. libelemental.org вместо использования Eigen с низкоуровневыми и потенциально непереносимыми функциями, такими как разделяемая память POSIX.   -  person Jeff Hammond    schedule 02.01.2015
comment
Спасибо за рекомендацию.   -  person Ivan    schedule 02.01.2015


Ответы (1)


Оказывается, это можно сделать в Eigen3. Использовать карту

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

person Ivan    schedule 29.12.2014