def test_compute_periods(self):
print 'test 7'
periods = capacity_calculations.compute_periods(28.5,[1.0,1.78858,3.08728],'Bare_Frame')
self.assertTrue(numpy.allclose(periods,[2.85000,3.81153,5.00764]))
es_ls3 = 0.060
data.append(ec_ls2)
data.append(ec_ls3)
data.append(es_ls2)
data.append(es_ls3)
data.append(betas)
# print data
collapse_type = capacity_calculations.compute_collapse_type(data)
collapse_type = 'Beam Sway'
#Compute capacity displacement
capacityDisp = capacity_calculations.compute_disps(data,collapse_type)
ductilities = capacity_calculations.compute_ductility(capacityDisp)
periods = capacity_calculations.compute_periods(data,ductilities)
elasticPeriods.append(periods[0])
#Compute demand displacement
demandDisp = demand_calculations.compute_demand_displacement(data,periods,spectraDisp,spectraPeriods,damping,ACCELEROGRAMS,ductilities)
DSpositions = damage_allocator.damage_state_position(capacityDisp, demandDisp)
damageStates=damageStates+DSpositions
print capacityDisp
# print demandDisp
# print ductilities
# print periods
# Compute the imls for each accelerogram
imlDamageStates = demand_calculations.compute_imls_damage_states(elasticPeriods,ACCELEROGRAMS,IMTs,IMTVALUES)
# Compute the best logarithmic mean and standard deviation for each curve