wrt_file) ## Convert the orientations from the calibration datasets from bunge euler angles ## to GSH coefficients gsh.euler_to_gsh(el, H, ns_cal, set_id_cal, step, wrt_file) ## Convert the orientations from the validation datasets from bunge euler angles ## to GSH coefficients gsh.euler_to_gsh(el, H, ns_val, set_id_val, step, wrt_file) ## Generate the fftn of the calibration microstructure function msf.micr_func(el, H, ns_cal, set_id_cal, step, wrt_file) ## Generate the fftn of the validation microstructure function msf.micr_func(el, H, ns_val, set_id_val, step, wrt_file) ## Perform the calibration for comp in compl: calibration.calibration_procedure(el, H, ns_cal, set_id_cal, step, comp, wrt_file) ## Perform the validation for comp in compl: validation.validation(el, H, ns_cal, ns_val, set_id_cal, set_id_val, step, comp, wrt_file) comp_app = 0 results.results_all(el, ns_val, set_id_val, step, 'epsilon', compl, comp_app)
Hset = Hset_ msg = "coefficients L=0-4" rr.WP(msg, wrt_file) rr.WP(msg, res_file) print(Hset) else: Hset = np.hstack([Hset_, ii]) msg = "coefficients L=0-4 and coeff %s" % ii rr.WP(msg, wrt_file) rr.WP(msg, res_file) print(Hset) # Convert the orientations from the calibration datasets from bunge euler # angles to GSH coefficients gsh.euler_to_gsh(el, Hset, ns_cal, set_id_cal, step, wrt_file) # Convert the orientations from the validation datasets from bunge euler # angles to GSH coefficients gsh.euler_to_gsh(el, Hset, ns_val, set_id_val, step, wrt_file) # Perform the calibration calibration.calibration_procedure(el, Hset.size, ns_cal, set_id_cal, step, comp, 'epsilon_t', wrt_file) # Perform the validation validation.validation(el, Hset.size, ns_cal, ns_val, set_id_cal, set_id_val, step, comp, 'epsilon_t', wrt_file) results.results(el, ns_val, set_id_val, step, 'epsilon', comp, 't', ii, res_file, res_file2)
ns_cal = 40 ns_val = 10 set_id_val = 'val' set_id_cal = 'cal' wrt_file = 'log_file.txt' X_cal = X[0:ns_cal, ...] X_val = X[ns_cal:, ...] y_cal = y[0:ns_cal, ...] y_val = y[ns_cal:, ...] # take fft of response fields y_fft_cal = np.fft.fftn(y_cal, axes=[1, 2, 3]) y_fft_val = np.fft.fftn(y_val, axes=[1, 2, 3]) # Convert the orientations from the calibration datasets from bunge euler # angles to GSH coefficients M_cal = gsh.euler_to_gsh(X_cal, el, H, ns_cal, set_id_cal, wrt_file) # Convert the orientations from the validation datasets from bunge euler # angles to GSH coefficients M_val = gsh.euler_to_gsh(X_val, el, H, ns_val, set_id_val, wrt_file) # Perform the calibration infl_coef = calibration.calibration_procedure(M_cal, y_fft_cal, el, H, ns_cal, wrt_file) # Perform the validation y_mks = validation.validation(M_val, infl_coef, el, wrt_file) results.results(infl_coef, y_val, y_mks, el, ns_val)
set_id_val = 'valRpc' wrt_file = 'log_order%s_%s%s_%s%s_%s.txt' % (order, ns_cal, set_id_cal, ns_val, set_id_val, time.strftime("%Y-%m-%d_h%Hm%M")) ### TOTAL CALIBRATION PROCEDURE ### ## Read the calibration microstructures and build the microstructure function msf.msf(el_cal, ns_cal, Hi, H, order, set_id_cal, wrt_file) ## Read the responses from the FE .dat files and perform the fftn for the calibration fegrab.fegrab(el_cal, ns_cal, set_id_cal, wrt_file) ## Perform the calibration calibration.calibration_procedure(el_cal, ns_cal, H, set_id_cal, wrt_file) ### TOTAL VALIDATION PROCEDURE ### ## Read the validation microstructures and build the microstructure function msf.msf(el_val, ns_val, Hi, H, order, set_id_val, wrt_file) ## Read the responses from the FE .dat files and perform the fftn for the validation fegrab.fegrab(el_val, ns_val, set_id_val, wrt_file) ## Perform the validation validation_viz.validation_zero_pad(el_cal, el_val, ns_cal, ns_val, H, set_id_cal, set_id_val, wrt_file) ## Calculate the results of the validation results.results(el_val, ns_val, set_id_val, 'epsilon', wrt_file)
vtk.read_euler(ns_val, set_id_val, vtk_filename, dir_val, wrt_file) for comp in xrange(9): vtk.read_meas(ns_val, set_id_val, comp, vtk_filename, tensor_ID, dir_val, wrt_file) ## Convert the orientations from the calibration datasets from bunge euler angles ## to GSH coefficients gsh.euler_to_gsh(ns_cal, set_id_cal, wrt_file) ## Convert the orientations from the validation datasets from bunge euler angles ## to GSH coefficients gsh.euler_to_gsh(ns_val, set_id_val, wrt_file) ## Generate the fftn of the calibration microstructure function msf.micr_func(ns_cal, set_id_cal, wrt_file) ## Generate the fftn of the validation microstructure function msf.micr_func(ns_val, set_id_val, wrt_file) ## Perform the calibration for comp in xrange(9): calibration.calibration_procedure(ns_cal, set_id_cal, comp, wrt_file) ## Perform the validation for comp in xrange(9): validation.validation_procedure(ns_cal, ns_val, set_id_cal, set_id_val, comp, wrt_file) results.results_all(ns_val, set_id_val, 'sigma')