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 に答える