absolute qudratic values from bic in biocoherence

[Tru3Dilemm4]

Hey there!

At the moment I am rewriting the package for python 3 and have come across a step in the biocoherence.py function that I don't quite understand:

In line 103 in bicoherence.py:
bic = abs(bic)**2 / ((Pyy * np.transpose(Pyy)) * mask)
Why is the absolute value determined from the expected value of the complex bispectrum (bic)? The squaring should surely elemenate the root operation from the denominator of the Haubrich normalization, so that the new bic should actually be squared as well?
Is this already for the visualization step or am I wrong in assuming that the normalized bispectrum is complex?

Best Regards
true

[Tru3Dilemm4]

I would suggest the following change to line 103:

bic = bic / np.sqrt((Pyy * np.transpose(Pyy)) * mask) # normalizing bic (Haubrich et al.)