search

GeoscienceAustralia / eqrm

Automatically exported from code.google.com/p/eqrm
Other
5 stars 4 forks source link

distribution.py failing when testing the demos #54

Closed GoogleCodeExporter closed 9 years ago

GoogleCodeExporter commented 9 years ago 
What steps will reproduce the problem?
1. Run python distribution.py

Original issue reported on code.google.com by duncan.g...@gmail.com on 3 Jul 2012 at 1:33

GoogleCodeExporter commented 9 years ago 
Here's the error I got;

============================
Doing  setdata_ScenRisk2.py
WARNING: max_width term in EQRM control file is deprecated. Replaced with 
scenario_max_width=15.
WARNING: csm_use_variability term in EQRM control file is deprecated. Replaced 
with csm_damage_state_use_variability=True.
WARNING: csm_standard_deviation term in EQRM control file is deprecated. 
Replaced with csm_damage_state_standard_deviation=0.3.
Logfile is './output/scen_risk2/log-0.txt' with logging level of DEBUG, console 
logging level is INFO
*******exported**********
event_set_handler = generate
P0: Generating event set
P0: Saving event set to ./output/scen_risk2/newc_event_set
P0: Event set created. Number of events=1002
P0: Sites set created. Number of sites=98
P0: do site 1 of 98

Traceback (most recent call last):
  File "/tmp/EQRM_distribution_temp_dirjQqbrk/python_eqrm/demo/demo_batchrun.py", line 51, in <module>
    demo_run()
  File "/tmp/EQRM_distribution_temp_dirjQqbrk/python_eqrm/demo/demo_batchrun.py", line 47, in demo_run
    analysis.main(run,True,compress_output)
  File "/tmp/EQRM_distribution_temp_dirjQqbrk/python_eqrm/eqrm_code/analysis.py", line 432, in main
    damage) = sites.calc_total_loss(SA, eqrm_flags, overloaded_MW)
  File "/tmp/EQRM_distribution_temp_dirjQqbrk/python_eqrm/eqrm_code/structures.py", line 265, in calc_total_loss
    loss_aus_contents=eqrm_flags.loss_aus_contents)
  File "/tmp/EQRM_distribution_temp_dirjQqbrk/python_eqrm/eqrm_code/damage_model.py", line 183, in aggregated_building_loss
    self.building_loss(ci=ci, loss_aus_contents=loss_aus_contents)
  File "/tmp/EQRM_distribution_temp_dirjQqbrk/python_eqrm/eqrm_code/damage_model.py", line 149, in building_loss
    damage_states = self.get_building_states()
  File "/tmp/EQRM_distribution_temp_dirjQqbrk/python_eqrm/eqrm_code/damage_model.py", line 137, in get_building_states
    beta_nsd_a, SA)
  File "/tmp/EQRM_distribution_temp_dirjQqbrk/python_eqrm/eqrm_code/damage_model.py", line 263, in state_probability
    p = cumulative_state_probability(threshold, beta, value)
  File "/tmp/EQRM_distribution_temp_dirjQqbrk/python_eqrm/eqrm_code/damage_model.py", line 321, in cumulative_state_probability
    return norm.cdf(temp)
  File "/usr/local/python-2.7.2/lib/python2.7/site-packages/scipy/stats/distributions.py", line 1198, in cdf
    place(output,cond,self._cdf(*goodargs))
  File "/usr/local/python-2.7.2/lib/python2.7/site-packages/scipy/stats/distributions.py", line 1908, in _cdf
    return _norm_cdf(x)
  File "/usr/local/python-2.7.2/lib/python2.7/site-packages/scipy/stats/distributions.py", line 1895, in _norm_cdf
    return special.ndtr(x)
TypeError: ufunc 'ndtr' not supported for the input types, and the inputs could 
not be safely coerced to any supported types according to the casting rule 
'safe'

When I executed the scenario standalone I did not get an error ( python 
setdata_ScenRisk2.py ).  Tried about 5 times.

I'm going to run distribution.py again, see how it goes.

Original comment by duncan.g...@gmail.com on 3 Jul 2012 at 1:47

GoogleCodeExporter commented 9 years ago 
The demo was crashing due to the extreme building capacity curves.  The 
standard deviations in the building_parameters were changed to be mean * 0.1. 
The demo's then passed.  Had to update the two scenarios that use csm 
variability.

The distribution.py now runs to the end.

Original comment by duncan.g...@gmail.com on 6 Jul 2012 at 1:58