def setUp(self):
""" Setting up a test. """
self.__files_to_remove = []
self.__dirs_to_remove = []
path_to_field_data = generateTestFilePath(
os.path.join("issue-sample", "no_particles", "data%T.h5"))
path_to_particle_data = generateTestFilePath(
os.path.join("issue-sample", "no_fields", "data%T.h5"))
path_to_data = generateTestFilePath(
os.path.join("git-sample", "data%T.h5"))
mode = openPMD.Access_Type.read_only
self.__field_series = openPMD.Series(path_to_field_data, mode)
self.__particle_series = openPMD.Series(path_to_particle_data, mode)
self.__series = openPMD.Series(path_to_data, mode)
#!/usr/bin/env python
"""
This file is part of the openPMD-api.
Copyright 2018 openPMD contributors
Authors: Axel Huebl
License: LGPLv3+
"""
import openPMD
if __name__ == "__main__":
series = openPMD.Series("../samples/git-sample/data%T.h5",
openPMD.Access_Type.read_only)
print("Read a Series with openPMD standard version %s" % series.openPMD)
print("The Series contains {0} iterations:".format(len(series.iterations)))
for i in series.iterations:
print("\t {0}".format(i))
print("")
i = series.iterations[100]
print("Iteration 100 contains {0} meshes:".format(len(i.meshes)))
for m in i.meshes:
print("\t {0}".format(m))
print("")
print("Iteration 100 contains {0} particle species:".format(
len(i.particles)))
for ps in i.particles:
print("\t {0}".format(ps))
print("")
import numpy as np
if __name__ == "__main__":
# user input: size of matrix to write, default 3x3
size = 3
# matrix dataset to write with values 0...size*size-1
global_data = np.arange(size*size, dtype=np.double)
print("Set up a 2D square array ({0}x{1}) that will be written".format(
size, size))
# open file for writing
series = openPMD.Series(
"../samples/3_write_serial_py.h5",
openPMD.Access_Type.create
)
print("Created an empty {0} Series".format(series.iteration_encoding))
print(len(series.iterations))
rho = series.iterations[1]. \
meshes["rho"][openPMD.Mesh_Record_Component.SCALAR]
datatype = openPMD.Datatype.DOUBLE
# datatype = openPMD.determineDatatype(global_data)
extent = [size, size]
dataset = openPMD.Dataset(datatype, extent)
print("Created a Dataset of size {0}x{1} and Datatype {2}".format(
dataset.extent[0], dataset.extent[1], dataset.dtype))