search

GeoscienceAustralia / eqrm

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

hzd_do_value changes as the event activity changes when SA = 0 #116

Closed GoogleCodeExporter closed 9 years ago

GoogleCodeExporter commented 9 years ago 
What steps will reproduce the problem?
1. go to revision 1381, on the reduce mem branch

2. Run the long_check_scenarios
3. Run t_nat01.py 

What is the expected output? 
There is check code that will crash since the output from hzd_do_value is 
different, given two similar input arrays.  If the event activities are set to 
zero for SA's that are known to be zero then the hzd_do_value results will be 
okay (equal).

This version of the code has two SA values being calculated, the close SA 
values and the 'all' SA values (close and far away) 

Please use labels and text to provide additional information.

Original issue reported on code.google.com by duncan.g...@gmail.com on 14 Mar 2013 at 3:22

GoogleCodeExporter commented 9 years ago 
Note this error also goes away if the gmm's are not split.
And if they are not collapsed.
And if fails for 1 in 20 - 400 sites.

Original comment by duncan.g...@gmail.com on 14 Mar 2013 at 5:34

GoogleCodeExporter commented 9 years ago 
TS_haz_51 also 'shows' this bug.

Data to reproduce this bug;
event_act_d_events [ 0.14132446  0.11654837  0.14177463  0.07311541  0.08137506 
 0.14829093
  0.06237196  0.06338166  0.08439131  0.08742621  0.11150119  0.09972694
  0.11350628  0.09771233  0.11803131  0.13920794  0.0717574   0.07188736
  0.08295507  0.09371418]
event_act_d_close [ 0.14132446  0.11654837  0.14177463  0.07311541  0.08137506  
0.14829093
  0.06237196  0.06338166  0.08439131  0.08742621  0.09371418]
all_event_indexes [ 0  1  2  3  4  5  6  7  8  9 19]
bedrock_SA_events [[  615418.52605636   615418.52605636   615418.52605636   
615418.52605636
    615418.52605636   615418.52605636   615418.52605636   615418.52605636
    615418.52605636   615418.52605636        0.                0.                0.
         0.                0.                0.                0.                0.
         0.          2000000.        ]]
bedrock_SA_close [[  615418.52605636   615418.52605636   615418.52605636   
615418.52605636
    615418.52605636   615418.52605636   615418.52605636   615418.52605636
    615418.52605636   615418.52605636  2000000.        ]]
1.0/array(eqrm_flags.return_periods) [ 0.1         0.02        0.01        
0.005       0.004       0.002
  0.00102564  0.001       0.0004      0.0002      0.0001    ]

Original comment by duncan.g...@gmail.com on 14 Mar 2013 at 5:59

GoogleCodeExporter commented 9 years ago 
This issue was closed by revision r1395.

Original comment by duncan.g...@gmail.com on 20 Mar 2013 at 1:59