5
$ mpiexec -n 4 python ./test_adios_mpi.py
7
from __future__ import print_function
10
from mpi4py import MPI
14
rank = comm.Get_rank()
15
size = comm.Get_size()
18
print ("\n>>> Writing ... (rank = %d)\n" % rank)
20
ad.init("config_mpi.xml", comm)
21
fd = ad.open("temperature", "adios_test_mpi.bp", "w", comm)
24
groupsize = 4 + 4 + 4 + 8 * 1 * NX
25
t = np.array(list(range(NX)), dtype=np.float64) + rank*NX
26
ad.set_group_size(fd, groupsize)
27
ad.write_int(fd, "NX", NX)
28
ad.write_int(fd, "rank", rank)
29
ad.write_int(fd, "size", size)
30
ad.write(fd, "temperature", t)
37
print( "\n>>> Reading ...\n")
39
f = ad.file("adios_test_mpi.bp", comm=MPI.COMM_SELF)
42
v = f.var['temperature']
47
assert (int(sum(sum(val))) == (size*NX-1)*(size*NX)/2)
50
print( "\n>>> Done.\n")
54
print( "\n>>> Test utility functions ...\n")
56
print( "bpls:\n", ad.bpls('adios_test_mpi.bp'))
57
print( "readvar:\n", ad.readvar("adios_test_mpi.bp", "temperature"))
59
print( "\n>>> Done.\n")