Jak przekazywać informacje MPI do typów w pythonie
Moją motywacją jest skuteczne przekazywanie informacji MPI z pythona do funkcji C wywoływanych przez ctypes. Użyłem mpi4py dla powiązań MPI w Pythonie. Chciałbym się tego nauczyć poprzez prosty przykładowy kod MPI napisany w C i wywoływany przez ctypes w Pythonie. Wyszczególniłem kroki i błąd, który otrzymuję podczas pracy poniżej.
Kod C [passMpi4Py.c]
#include <stdio.h>
#include <mpi.h>
void sayhello(MPI_Comm comm)
{
int size, rank;
MPI_Comm_size(comm, &size);
MPI_Comm_rank(comm, &rank);
printf("Hello, World! "
"I am process %d of %d.\n",
rank, size);
}
Powyższy kod c skompilowałem za pomocą gcc / openmpi-1.6 w następujący sposób:
mpicc -shared -Wl, -soname, passMpi4Py -o passMpi4Py.so -fPIC passMpi4Py.c
Python Wrapper [passMpi4PyWrapper.py]
import ctypes
from mpi4py import MPI
testlib = ctypes.CDLL('path-to-file/passMpi4Py/passMpi4Py.so')
testlib.sayhello(MPI.COMM_WORLD)
Kiedy próbuję uruchomić powyższy kod za pomocą
mpirun -np 4 python passMpi4PyWrapper.py
Dostaję następujący błąd
Traceback (most recent call last):
Traceback (most recent call last):
File "passMpi4PyWrapper.py", line 5, in <module>
Traceback (most recent call last):
File "passMpi4PyWrapper.py", line 5, in <module>
File "passMpi4PyWrapper.py", line 5, in <module>
Traceback (most recent call last):
testlib.sayhello(MPI.COMM_WORLD)
ctypes.ArgumentError: argument 1: <type 'exceptions.TypeError'>: Don't know how to convert parameter 1
testlib.sayhello(MPI.COMM_WORLD)
ctypes.ArgumentError: argument 1: <type 'exceptions.TypeError'>: Don't know how to convert parameter 1
testlib.sayhello(MPI.COMM_WORLD)
ctypes.ArgumentError: argument 1: <type 'exceptions.TypeError'>: Don't know how to convert parameter 1
File "passMpi4PyWrapper.py", line 5, in <module>
testlib.sayhello(MPI.COMM_WORLD)
ctypes.ArgumentError: argument 1: <type 'exceptions.TypeError'>: Don't know how to convert parameter 1
Aktualizacja:
Używanie * MPI_COMM_WORLD * zamiastkom w programie C funkcje MPI pomagają mi usunąć błąd. Jednak nadal chciałbym wiedzieć, czy jest to najlepszy możliwy sposób przekazywania informacji MPI do programu C.