def vector(value=[], unit=None):
if unit is None:
if isinstance(value, core.unit):
return VectorQuantity([], unit=value)
elif isinstance(value, ScalarQuantity):
return value.as_vector_with_length(1)
else:
result = AdaptingVectorQuantity()
result.extend(value)
return result
else:
if isinstance(value, ScalarQuantity):
return value.as_vector_with_length(1)
else:
return VectorQuantity(value, unit)
def vector(value = [], unit = None):
if unit is None:
if isinstance(value, core.unit):
return VectorQuantity([], unit = value)
elif isinstance(value, ScalarQuantity):
return value.as_vector_with_length(1)
else:
result = AdaptingVectorQuantity()
result.extend(value)
return result
else:
if isinstance(value, ScalarQuantity):
return value.as_vector_with_length(1)
else:
return VectorQuantity(value, unit)
def read_from_hdf5(filename):
""" Reads an hdf5 file and returns a HydroResult class. """
f = h5py.File(filename, 'r')
data = {}
for keyword in f.keys():
unitstring = f[keyword].attrs['unit']
evaluable_unit = parse_unitsstring(unitstring)
vq = f[keyword].value | eval(evaluable_unit)
avq = AdaptingVectorQuantity()
avq.extend(vq)
try:
mass_fractions = eval(f[keyword].attrs['mass_fractions'])
setattr(avq, 'mf', mass_fractions)
except KeyError, AttributeError:
pass
data.update({keyword: avq})
def read_from_hdf5(filename):
""" Reads an hdf5 file and returns a HydroResult class. """
f = h5py.File(filename,'r')
data = {}
for keyword in f.keys():
unitstring = f[keyword].attrs['unit']
evaluable_unit = parse_unitsstring(unitstring)
vq = f[keyword].value | eval(evaluable_unit)
avq = AdaptingVectorQuantity()
avq.extend(vq)
try:
mass_fractions = eval(f[keyword].attrs['mass_fractions'])
setattr(avq, 'mf', mass_fractions)
except KeyError, AttributeError:
pass
data.update({keyword:avq})
