else: param_string += ", " print(" parameters: %s" % param_string) # get diameter for entire cuboid lb,ub = [lower_bounds],[upper_bounds] cuboid_volume = volume(lb[0],ub[0]) params0, subdiams0, diam0, probmass0 = test_cuboids(lb,ub,RVstart,RVend,\ cuboid_volume) SOLVED_PARAMETERS, SUB_DIAMETERS = params0, subdiams0 TOTAL_DIAMETERS, PROBABILITY_MASS = diam0, probmass0 if DEBUG: print("\nSTATUS = %s" % STATUS) if not DEBUG: pof = sampled_pof(model,lb[0],ub[0]) print("Exact PoF: %s" % pof) for i in range(len(lb)): print("\n") print(" lower bounds: %s" % lb[i]) print(" upper bounds: %s" % ub[i]) for solved in params0[0]: print("solved: %s" % solved) print("subdiameters (squared): %s" % subdiams0[0]) print("diameter (squared): %s" % diam0[0]) print(" probability mass: %s" % probmass0[0]) expectation = expectation_value(model,lower_bounds,upper_bounds) print(" expectation: %s" % expectation) mean_value = mean(expectation,cuboid_volume) print(" mean value: %s" % mean_value)
param_string += ", " print " parameters: %s" % param_string # get diameter for entire cuboid lb,ub = [lower_bounds],[upper_bounds] cuboid_volume = volume(lb[0],ub[0]) params0, subdiams0, diam0, probmass0 = test_cuboids(lb,ub,RVstart,RVend,\ cuboid_volume) SOLVED_PARAMETERS, SUB_DIAMETERS = params0, subdiams0 TOTAL_DIAMETERS, PROBABILITY_MASS = diam0, probmass0 if DEBUG: print "\nSTATUS = %s" % STATUS if not DEBUG: pts = random_samples(lb[0],ub[0]) pof = sampled_pof(model,pts) print "Exact PoF: %s" % pof # prepare new set of random samples (across entire domain) as 'data' if not PER_AI: pts = random_samples(lb[0],ub[0],num_sample_points) for i in range(len(lb)): print "\n" print " lower bounds: %s" % lb[i] print " upper bounds: %s" % ub[i] for solved in params0[0]: print "solved: %s" % solved print "subdiameters (squared): %s" % subdiams0[0] print "diameter (squared): %s" % diam0[0] print " probability mass: %s" % probmass0[0] expectation = expectation_value(model,lower_bounds,upper_bounds)