# Basic output
# e.printInfo()
print t
# Validation / data comparison
m1_data = d.get_field("Mbera", t) - d.get_field("Mbera", 0)
b1_data = d.get_field("Mentao", t) - d.get_field("Mentao", 0)
b2_data = d.get_field("Bobo-Dioulasso", t) - d.get_field(
"Bobo-Dioulasso", 0)
n1_data = d.get_field("Abala", t) - d.get_field("Abala", 0)
n2_data = d.get_field("Mangaize", t) - d.get_field("Mangaize", 0)
errors = [
a.rel_error(m1.numAgents, m1_data),
a.rel_error(b1.numAgents, b1_data),
a.rel_error(b2.numAgents, b2_data),
a.rel_error(n1.numAgents, n1_data),
a.rel_error(n2.numAgents, n2_data)
]
abs_errors = [
a.abs_error(m1.numAgents, m1_data),
a.abs_error(b1.numAgents, b1_data),
a.abs_error(b2.numAgents, b2_data),
a.abs_error(n1.numAgents, n1_data),
a.abs_error(n2.numAgents, n2_data)
]
print "Mbera: ", m1.numAgents, ", data: ", m1_data, ", error: ", errors[
0]
# e.printInfo()
# Validation / data comparison
m1_data = d.get_field("Mbera", t) #- d.get_field("Mbera", 0)
m2_data = d.get_field("Fassala", t) #- d.get_field("Mbera", 0)
b1_data = d.get_field("Mentao", t) #- d.get_field("Mentao", 0)
b2_data = d.get_field("Bobo-Dioulasso",
t) #- d.get_field("Bobo-Dioulasso", 0)
n1_data = d.get_field("Abala", t) #- d.get_field("Abala", 0)
n2_data = d.get_field("Mangaize", t) #- d.get_field("Mangaize", 0)
n3_data = d.get_field("Niamey", t) #- d.get_field("Niamey", 0)
n4_data = d.get_field("Tabareybarey",
t) #- d.get_field("Tabareybarey", 0)
errors = [
a.rel_error(m1.numAgents, m1_data),
a.rel_error(m2.numAgents, m2_data),
a.rel_error(b1.numAgents, b1_data),
a.rel_error(b2.numAgents, b2_data),
a.rel_error(n1.numAgents, n1_data),
a.rel_error(n2.numAgents, n2_data),
a.rel_error(n3.numAgents, n3_data),
a.rel_error(n4.numAgents, n4_data)
]
abs_errors = [
a.abs_error(m1.numAgents, m1_data),
a.abs_error(m2.numAgents, m2_data),
a.abs_error(b1.numAgents, b1_data),
a.abs_error(b2.numAgents, b2_data),
a.abs_error(n1.numAgents, n1_data),
a.abs_error(n2.numAgents, n2_data),
csvformat="mali-pdf")
for t in range(0, end_time):
new_refs = d.get_new_refugees(t)
for i in range(0, new_refs):
e.addAgent(location=l1)
e.evolve()
# e.printInfo()
print t
l2_data = d.get_field("Mauritania", t) - d.get_field("Mauritania", 0)
l3_data = d.get_field("Niger", t) - d.get_field("Niger", 0)
l4_data = d.get_field("Burkina Faso", t) - d.get_field(
"Burkina Faso", 0)
errors = [
a.rel_error(l2.numAgents, l2_data),
a.rel_error(l3.numAgents, l3_data),
a.rel_error(l4.numAgents, l4_data)
]
print "Kiffa: ", l2.numAgents, ", data: ", l2_data, ", error: ", errors[
0]
print "Niamey: ", l3.numAgents, ", data: ", l3_data, ", error: ", errors[
1]
print "Bobo-Dioulasso: ", l4.numAgents, ", data: ", l4_data, ", error: ", errors[
2]
print "Cumulative error: ", np.sum(
errors), ", Squared error: ", np.sqrt(np.sum(np.power(errors, 2)))
if np.abs(np.sum(errors) - 0.495521376979) > 0.1:
print "TEST FAILED."
