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')