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

GetLLM: Coupling is crashing if input file wasn't cleaned with svd #184

Open
elenafol opened this issue Dec 7, 2018 · 0 comments
Open

Comments

@elenafol
Copy link
Contributor

elenafol commented Dec 7, 2018

Using an ideal tracking file (no noise, but with coupling), coupling computation is crashing and GetLLM produces following error, if cleaning was not used in harmonic analysis:

14:05:11 - GetLlm STDERR: /afs/cern.ch/work/e/efol/public/beta-beat-source/Beta-Beat.src/GetLLM/algorithms/coupling.py:605: RuntimeWarning: divide by zero encountered in divide
14:05:11 - GetLlm STDERR: f_old_out[bn1] = np.average(abs(f1001ij), weights=1/std_f1001ij2)
14:05:11 - GetLlm STDERR: /afs/cern.ch/work/o/omc/anaconda/lib/python2.7/site-packages/numpy/lib/function_base.py:1160: RuntimeWarning: invalid value encountered in double_scalars
14:05:11 - GetLlm STDERR: avg = np.multiply(a, wgt, dtype=result_dtype).sum(axis)/scl
14:05:11 - GetLlm STDERR: /afs/cern.ch/work/e/efol/public/beta-beat-source/Beta-Beat.src/GetLLM/algorithms/coupling.py:610: RuntimeWarning: divide by zero encountered in divide
14:05:11 - GetLlm STDERR: f1001i = np.average(f1001ij, weights=1/std_f1001ij
2)
14:05:11 - GetLlm STDERR: /afs/cern.ch/work/e/efol/public/beta-beat-source/Beta-Beat.src/GetLLM/algorithms/coupling.py:611: RuntimeWarning: divide by zero encountered in divide
14:05:11 - GetLlm STDERR: f1010i = np.average(f1010ij, weights=1/std_f1010ij2)
14:05:11 - GetLlm STDERR: /afs/cern.ch/work/e/efol/public/beta-beat-source/Beta-Beat.src/GetLLM/algorithms/coupling.py:613: RuntimeWarning: divide by zero encountered in divide
14:05:11 - GetLlm STDERR: f1001istd = np.sqrt(1/sum(1/std_f1001ij
2))
14:05:11 - GetLlm STDERR: /afs/cern.ch/work/e/efol/public/beta-beat-source/Beta-Beat.src/GetLLM/algorithms/coupling.py:614: RuntimeWarning: divide by zero encountered in divide
14:05:11 - GetLlm STDERR: f1010istd = np.sqrt(1/sum(1/std_f1010ij**2))
14:05:11 - GetLlm STDERR: /afs/cern.ch/work/e/efol/public/beta-beat-source/Beta-Beat.src/GetLLM/algorithms/coupling.py:652: RuntimeWarning: divide by zero encountered in double_scalars
14:05:11 - GetLlm STDERR: denom += 1/fwqw[bn2][0][1]**2 # denominator for weighted average

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

No branches or pull requests

1 participant