RuntimeWarning: divide by zero encountered in divide

[luiztauffer]

These warning pop up a lot during execution of the pipeline. Are these expected, or something that should be fixed?

/voluseg/voluseg/_steps/step4d.py:52: RuntimeWarning: divide by zero encountered in divide
  H *= (timeseries_mean / H.mean(1)[:, None])           # normalize component timeseries
/voluseg/voluseg/_steps/step4d.py:52: RuntimeWarning: invalid value encountered in multiply
  H *= (timeseries_mean / H.mean(1)[:, None])           # normalize component timeseries
retrying factorization of block 1: array must not contain infs or NaNs
(np.float64(0.7999999999999998), np.int64(403))
/voluseg/voluseg/_steps/step4d.py:52: RuntimeWarning: divide by zero encountered in divide
  H *= (timeseries_mean / H.mean(1)[:, None])           # normalize component timeseries
/voluseg/voluseg/_steps/step4d.py:52: RuntimeWarning: invalid value encountered in multiply
  H *= (timeseries_mean / H.mean(1)[:, None])           # normalize component timeseries
retrying factorization of block 3: array must not contain infs or NaNs
(np.float64(0.7999999999999998), np.int64(408))
/voluseg/voluseg/_steps/step4d.py:52: RuntimeWarning: divide by zero encountered in divide
  H *= (timeseries_mean / H.mean(1)[:, None])           # normalize component timeseries
/voluseg/voluseg/_steps/step4d.py:52: RuntimeWarning: invalid value encountered in multiply
  H *= (timeseries_mean / H.mean(1)[:, None])           # normalize component timeseries
retrying factorization of block 0: array must not contain infs or NaNs
(np.float64(0.7999999999999998), np.int64(411))
/voluseg/voluseg/_steps/step4d.py:52: RuntimeWarning: divide by zero encountered in divide
  H *= (timeseries_mean / H.mean(1)[:, None])           # normalize component timeseries
/voluseg/voluseg/_steps/step4d.py:52: RuntimeWarning: invalid value encountered in multiply
  H *= (timeseries_mean / H.mean(1)[:, None])           # normalize component timeseries
retrying factorization of block 2: array must not contain infs or NaNs
(np.float64(0.7999999999999998), np.int64(413))
cell initialization: 0.2 minutes.

[mikarubi]

These do not need to be fixed, but are a part of choosing a feasible number of cells segments.

[luiztauffer]

is there a way we can formalize this into an informative warning or suggestion to the users?

[mikarubi]

yes, the relevant message is returned here:

voluseg/voluseg/_steps/step4.py

Line 148 in 2c4d935

except ValueError as msg:

the warning typically precedes the error so perhaps it can be caught upstream.