iome = IOManager()
# Read file with simulation data
try:
filename = sys.argv[1]
except IndexError:
filename = GD.file_resultdatafile
iomc.open_file(filename=filename)
# New file for eigen transformed data
P = iomc.load_parameters()
iome.create_file(P, filename=filename[:-5]+"_eigen.hdf5")
# Iterate over all groups
for groupid in iomc.get_group_ids():
# Create the group if necessary
if not groupid in iome.get_group_ids():
iome.create_group(groupid=groupid)
for blockid in iomc.get_block_ids(groupid=groupid):
print("Computing eigentransformation of data in block '"+str(blockid)+"'")
# Create the block if necessary
if not blockid in iome.get_block_ids(groupid=groupid):
iome.create_block(blockid=blockid, groupid=groupid)
# See if we have a wavefunction
if iomc.has_wavefunction(blockid=blockid):
from EigentransformWavefunction import transform_wavefunction_to_eigen
# Read file with simulation data
iomc = IOManager()
iomc.open_file(inputfile)
iome = IOManager()
iome.create_file(outputfile)
# New file for eigen transformed data
P = iomc.load_parameters()
# Save the simulation parameters
iome.add_parameters()
iome.save_parameters(P)
# Iterate over all groups
for groupid in iomc.get_group_ids():
# Create the group if necessary
if groupid not in iome.get_group_ids():
iome.create_group(groupid=groupid)
for blockid in iomc.get_block_ids(groupid=groupid):
print("Computing eigentransformation of data in block '{}'".format(blockid))
# Create the block if necessary
if blockid not in iome.get_block_ids(groupid=groupid):
iome.create_block(blockid=blockid, groupid=groupid)
# See if we have a wavefunction
if iomc.has_wavefunction(blockid=blockid):
from WaveBlocksND.Interface import EigentransformWavefunction
