Describe your proposed solution

The function _ica_par has the following line of code:

We could investigate, if performance is improved by replacing np.diag(np.dot(W1, W.T)) by np.einsum("ij,ij->i", W1, W).

Additional context

A short benchmark would be nice in order to decide whether it's worth.