Bug in auxiva? Script did not show any error, but got nan in the result of source matrix about pyroomacoustics

[datalw]

Oh, I've just realized the problem is, that there is only one source in this test signal and it is overdetermined.

[datalw]

Hi, I have to reopen this issue, because the reason was apparently not what I have found earlier. With other data, I experience this warning as well.

The correct reason is that in the function auxiva() in the script auxiva.py, WV has got nan first, after the line WV = np.matmul(W_hat, V) . If I have understood it right - please correct me if not - W_hat is initialized at the beginning and the value should not change, after the iterations have begun. But it changes anyway sometimes along the iterations. I looked into the code, there are variable assignments like

W = W_hat[:, :n_src, :]
J = W_hat[:, n_src:, :n_src] 

In this case, if W or J change after the assignments, W_hat will be partially changed along. I not sure, if this is the way the codes are supposed to run. A simple example of unexpected change of a variable see here.

Codes to re-produce the warning:

import pyroomacoustics as pra
import numpy as np

iva_input = np.load("0th_iva_input.npy")
epoch_scvs, demix = pra.bss.auxiva(iva_input, n_src=iva_input.shape[-1], n_iter=20, proj_back=False, model='laplace',
                                   init_eig=False, return_filters=True)

The sample input can be found here.

Because the warning does not happen in every run, one might have to run it several times.