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