if not fr[i]:
eva = array(self.freq_data.eva[i][newaxis], dtype="float64")
eve = array(self.freq_data.eve[i][newaxis], dtype="complex128")
beamformerOutput = self.steer._beamformerCall(f[i], self.r_diag, normFactor, (eva[:, na:na+1], eve[:, :, na:na+1]))[0]
if self.r_diag: // set (unphysical) negative output values to 0
indNegSign = sign(beamformerOutput) < 0
beamformerOutput[indNegSign] = 0
After Change
f = self.freq_data.fftfreq()
na = int(self.na) // eigenvalue taken into account
normFactor = self.sig_loss_norm()
param_steer_type, steer_vector = self._beamformer_params()
for i in self.freq_data.indices:
if not fr[i]:
eva = array(self.freq_data.eva[i], dtype="float64")
eve = array(self.freq_data.eve[i], dtype="complex128")
