Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

mcmc failing to run #59

Open
SquirrelKnight opened this issue Jan 24, 2019 · 0 comments
Open

mcmc failing to run #59

SquirrelKnight opened this issue Jan 24, 2019 · 0 comments

Comments

@SquirrelKnight
Copy link

Please fill in the following details as much as possible when submitting your issue as this will help with testing and reproducing the issue

Environment Information

{
"version": "1.0.5",
"c_extensions present": [
"cmarkov_chain_monte_carlo",
"cprobability",
"cmoment_tensor_conversion",
"cscatangle"
],
"platform": "darwin",
"num_threads": 8,
"python version": "3.7.1 (default, Dec 14 2018, 13:28:58) \n[Clang 4.0.1 (tags/RELEASE_401/final)]",
"python version info": [
3,
7,
1,
"final",
0
],
"dependency info": {
"numpy": "1.15.4",
"scipy": "1.1.0",
"matplotlib": "3.0.2",
"cython": "0.29.2",
"sphinx": "1.8.2",
"h5py": "2.8.0"
}
}

What error you are getting

Warning: divide by zero encountered in log
Cython Error
Traceback (most recent call last):
File "/Users/jesse/anaconda3/lib/python3.7/site-packages/MTfit/algorithms/markov_chain_monte_carlo.py", line 1582, in _acceptance_check
dc_prior=getattr(self, 'dc_prior', 0.))
TypeError: Argument 'x' has incorrect type (expected list, got dict)

What behaviour you would expect

This causes the program to stop, should be performing an McMC inversion. The regular 'iterate' algorithm works fine.

Reproduction steps

Run an McMC inversion using the following parameters:

algorithm = 'mcmc'
dc = False
parallel = False
phys_mem = 1
output_format = 'pickle'
conversion = False
inversion_options = ['PPolarity','P/SHAmplitudeRatio']
burn_length = 30000
chain_length = 100000
min_acceptance_rate = 0.1
max_acceptance_rate = 0.3

inversion_object = Inversion(event,algorithm=algorithm, dc=dc,parallel=parallel,phys_mem=phys_mem,convert=conversion,output_format=output_format, inversion_options=inversion_options, burn_length=burn_length,chain_length=chain_length, min_acceptance_rate=min_acceptance_rate,max_acceptance_rate=max_acceptance_rate)

inversion_object.forward()

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant