def _get_spirrid_response( self ): rf = StressInFiberWithConstantFriction( u = 0.05 ) s = SPIRRID( rf = rf, cached_dG = True, compiled_QdG_loop = False, compiled_eps_loop = False ) n_int = 50 s.add_rv( 'tau', distribution = 'uniform', loc = 100.1, scale = 120.0, n_int = n_int ) s.add_rv( 'l', distribution = 'uniform', loc = 0.0, scale = 105.0, n_int = n_int ) #s.add_rv( 'D_f', distribution = 'uniform', loc = 25.0e-3, scale = 5.0e-3, n_int = n_int ) eps_list = s.eps_list mu_q = s.mu_q_grid return eps_list, mu_q
rf = StressInFiberWithConstantFriction( u = 0.01 ) #print rf.ctrl_keys # rf.configure_traits() from stats.spirrid.spirrid_nd import SPIRRID # X = linspace( 0, 1.5, 100 ) # Y = rf( X, 0.01, 0.8, 0.02, 7e10, 1e-9 ) # plt.plot( X, Y, linewidth = 2 ) s = SPIRRID( rf = rf, cached_dG = True, compiled_QdG_loop = False, compiled_eps_loop = False ) # construct the random variables n_int = 20 s.add_rv( 'tau', distribution = 'uniform', loc = 0.1, scale = 4.0, n_int = n_int ) #s.add_rv( 'l', distribution = 'uniform', loc = 0.0, scale = 5.0, n_int = n_int ) #s.add_rv( 'D_f', distribution = 'uniform', loc = 25.0e-3, scale = 5.0e-3, n_int = n_int ) # eps_list = s.eps_list mu_q = s.mu_q_grid