1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
|
#include <mpi.h>
#include <iostream>
static void sayhello(MPI_Comm comm)
{
if (comm == MPI_COMM_NULL) {
std::cout << "You passed MPI_COMM_NULL !!!" << std::endl;
return;
}
int size;
MPI_Comm_size(comm, &size);
int rank;
MPI_Comm_rank(comm, &rank);
int plen; char pname[MPI_MAX_PROCESSOR_NAME];
MPI_Get_processor_name(pname, &plen);
std::cout <<
"Hello, World! " <<
"I am process " << rank <<
" of " << size <<
" on " << pname <<
"." << std::endl;
}
#include <nanobind/nanobind.h>
#define MPI4PY_LIMITED_API 1
#define MPI4PY_LIMITED_API_SKIP_MESSAGE 1
#define MPI4PY_LIMITED_API_SKIP_SESSION 1
#include <mpi4py/mpi4py.h>
namespace nb = nanobind;
template<typename T> T py2mpi(nb::object);
template<> MPI_Comm py2mpi(nb::object obj)
{
PyObject *pyobj = obj.ptr();
MPI_Comm *mpi_ptr = PyMPIComm_Get(pyobj);
if (!mpi_ptr) throw nb::python_error();
return *mpi_ptr;
}
static void hw_sayhello(nb::object py_comm)
{
MPI_Comm comm = py2mpi<MPI_Comm>(py_comm);
sayhello(comm);
}
NB_MODULE(helloworld, m) {
if (import_mpi4py() < 0)
throw nb::python_error();
m.def("sayhello", &hw_sayhello);
}
/*
* Local Variables:
* mode: C++
* End:
*/
|
