def read_euler(el, ns, set_id, step, newdir, wrt_file, funit):
start = time.time()
euler = np.zeros([ns, 3, el**3])
nwd = os.getcwd() + '\\' + newdir
# nwd = os.getcwd() + '/' + newdir # for unix
os.chdir(nwd)
sn = 0
for filename in os.listdir(nwd):
if filename.endswith('%s.vtk' % step):
euler[sn, :, :] = rr.read_vtk_vector(filename=filename)
sn += 1
if funit == 1:
euler = euler * (np.pi / 180.)
"""return to the original directory"""
os.chdir('..')
f = h5py.File("ref_%s%s_s%s.hdf5" % (ns, set_id, step), 'a')
f.create_dataset('euler', data=euler)
f.close()
end = time.time()
timeE = np.round((end - start), 3)
msg = 'euler angles read from .vtk file for %s: %s seconds' % (set_id,
timeE)
rr.WP(msg, wrt_file)
def read_euler(el,ns, set_id, step, newdir, wrt_file, funit):
start = time.time()
euler = np.zeros([ns,3,el**3])
# nwd = os.getcwd() + '\\' + newdir
nwd = os.getcwd() + '/' + newdir #for unix
os.chdir(nwd)
sn = 0
for filename in os.listdir(nwd):
if filename.endswith('%s.vtk' %step):
euler[sn,:,:] = rr.read_vtk_vector(filename = filename)
sn += 1
if funit == 1:
euler = euler * (np.pi/180.)
## return to the original directory
os.chdir('..')
np.save('euler_%s%s_s%s' %(ns,set_id,step), euler)
end = time.time()
timeE = np.round((end - start),3)
msg = 'euler angles read from .vtk file for %s: %s seconds' %(set_id, timeE)
rr.WP(msg,wrt_file)
def read_euler(ns, set_id, newdir, funit):
start = time.time()
C = const()
euler = np.zeros([ns, 3, C['el']**3])
# nwd = os.getcwd() + '\\' + newdir
nwd = os.getcwd() + '/' + newdir # for unix
os.chdir(nwd)
sn = 0
for filename in os.listdir(nwd):
if filename.endswith('%s.vtk' % C['step']):
euler[sn, :, :] = rr.read_vtk_vector(filename=filename)
sn += 1
if funit == 1:
euler = euler * (np.pi / 180.)
"""return to the original directory"""
os.chdir('..')
f = h5py.File("spatial.hdf5", 'a')
f.create_dataset('euler_%s' % set_id, data=euler)
f.close()
end = time.time()
timeE = np.round((end - start), 3)
msg = 'euler angles read from .vtk file for %s: %s seconds' % (set_id,
timeE)
rr.WP(msg, C['wrt_file'])
def read_euler(strt, ns, set_id, newdir, funit):
start = time.time()
C = const()
euler = np.zeros([ns, 3, C['el']**3])
# nwd = os.getcwd() + '\\' + newdir
nwd = os.getcwd() + '/' + newdir # for unix
os.chdir(nwd)
for ii in xrange(ns):
sn = strt + ii + 1
filename = "Ti64_Dream3D_v01_Output_%s.vtk" % sn
euler[ii, :, :] = rr.read_vtk_vector(filename=filename)
if funit == 1:
euler = euler * (np.pi / 180.)
# return to the original directory
os.chdir('..')
f = h5py.File("spatial.hdf5", 'a')
dset_name = 'euler_%s' % set_id
f.create_dataset(dset_name, data=euler)
f.close()
end = time.time()
timeE = np.round((end - start), 3)
msg = 'euler angles read from .vtk file for %s: %s seconds' \
% (set_id, timeE)
rr.WP(msg, C['wrt_file'])
def read_euler(ns, set_id, step, vtk_filename, newdir, wrt_file, funit):
start = time.time()
## el is the # of elements per side of the cube
el = 21
euler = np.zeros([el**3,ns,3])
## change to directory with the .vtk files
# cwd = os.getcwd()
## os.chdir(cwd + '\\' + newdir)
# os.chdir(cwd + '/' + newdir) #for unix
#
# for sn in xrange(ns):
# l_sn = str(sn+1).zfill(5)
# euler[:,sn,:] = rr.read_vtk_vector(filename = vtk_filename %l_sn)
nwd = os.getcwd() + '\\' + newdir
# nwd = os.getcwd() + '/' + newdir #for unix
os.chdir(nwd)
sn = 0
for filename in os.listdir(nwd):
if filename.endswith('%s.vtk' %step):
euler[:,sn,:] = rr.read_vtk_vector(filename = filename)
sn += 1
if funit == 1:
euler = euler * (np.pi/180.)
## return to the original directory
os.chdir('..')
np.save('euler_%s%s_s%s' %(ns,set_id,step), euler)
end = time.time()
timeE = np.round((end - start),3)
msg = 'euler angles read from .vtk file for %s: %s seconds' %(set_id, timeE)
rr.WP(msg,wrt_file)
def read_euler(el, ns, set_id, step, newdir, wrt_file, funit):
start = time.time()
euler = np.zeros([ns, 3, el**3])
# nwd = os.getcwd() + '\\' + newdir
nwd = os.getcwd() + '/' + newdir # for unix
os.chdir(nwd)
sn = 0
# for filename in os.listdir(nwd):
# if filename.endswith('%s.vtk' % step):
# euler[sn, :, :] = rr.read_vtk_vector(filename=filename)
# sn += 1
for sn in xrange(ns):
filename = "Ti64_Dream3D_v01_Output_%s.vtk" % str(sn + 1)
euler[sn, :, :] = rr.read_vtk_vector(filename=filename)
if funit == 1:
euler = euler * (np.pi / 180.)
# return to the original directory
os.chdir('..')
f = h5py.File("data.hdf5", 'a')
dset_name = 'euler_%s%s_s%s' % (ns, set_id, step)
f.create_dataset(dset_name, data=euler)
f.close()
end = time.time()
timeE = np.round((end - start), 3)
msg = 'euler angles read from .vtk file for %s: %s seconds' \
% (set_id, timeE)
rr.WP(msg, wrt_file)
def read_euler(el, ns, set_id, step, newdir, wrt_file, funit):
start = time.time()
euler = np.zeros([ns, 3, el**3])
# nwd = os.getcwd() + '\\' + newdir
nwd = os.getcwd() + '/' + newdir # for unix
os.chdir(nwd)
sn = 0
for filename in os.listdir(nwd):
if filename.endswith('%s.vtk' % step):
euler[sn, :, :] = rr.read_vtk_vector(filename=filename)
sn += 1
if funit == 1:
euler = euler * (np.pi / 180.)
# return to the original directory
os.chdir('..')
# open HDF5 file
base = tb.open_file("ref_%s%s_s%s.h5" % (ns, set_id, step), mode="a")
# initialize array for the euler angles
base.create_array('/msf', 'euler', euler,
'euler angle of the original microstructures')
# close the HDF5 file
base.close()
end = time.time()
timeE = np.round((end - start), 3)
msg = 'euler angles read from .vtk file for %s: %s seconds' % (set_id,
timeE)
rr.WP(msg, wrt_file)
c = 0
for comp in compl:
compp = compd[comp]
r_temp = rr.read_vtk_tensor(filename = filename, tensor_id = tensor_ID, comp = compp)
r_real[c,...] = r_temp.reshape([el,el,el])
print compl
print np.mean(r_real[c,...])
c += 1
euler = rr.read_vtk_vector(filename).reshape([3,el,el,el])
for dispcomp in xrange(np.size(compl)):
plt.close(dispcomp)
## Plot slices of the response
plt.figure(num=dispcomp,figsize=[14,6])
plt.subplot(231)
ax = plt.imshow(euler[0,0,:,:], origin='lower', interpolation='none',
cmap='jet')
plt.colorbar(ax)
plt.title('Microstructure, slice 0')
