30f4a3e03dc59daed7979659701e4c537fbae91c,acoular/fbeamform.py,BeamformerFunctional,calc,#BeamformerFunctional#Any#Any#,379

Before Change

        // prepare calculation
        kj = 2j*pi*self.freq_data.fftfreq()/self.c
        numchannels = int(self.freq_data.numchannels)
        e = zeros((numchannels), "D")
        h = empty((1, self.grid.size), "d")
        // function
        beamfunc = self.get_beamfunc("_os")
        if self.r_diag:
            adiv = sqrt(1.0/(numchannels*numchannels-numchannels))
            scalefunc = lambda h : adiv*(multiply(adiv*h, (sign(h)+1-1e-35)/2))**self.gamma
        else:
            adiv = 1.0/(numchannels)
            scalefunc = lambda h : adiv*(adiv*h)**self.gamma
        for i in self.freq_data.indices:        
            if not fr[i]:
                eva = array(self.freq_data.eva[i][newaxis], dtype="float64")**(1.0/self.gamma)
                eve = array(self.freq_data.eve[i][newaxis], dtype="complex128")
                kji = kj[i, newaxis]
                beamfunc(e, h, self.r0, self.rm, kji, eva, eve, 0, numchannels)
                ac[i] = scalefunc(h)
                fr[i] = True
            
class BeamformerCapon( BeamformerBase ):
    

After Change

                eva = array(self.freq_data.eva[i][newaxis], dtype="float64") ** (1.0 / self.gamma)
                eve = array(self.freq_data.eve[i][newaxis], dtype="complex128")
                kji = kj[i, newaxis]
                beamformerOutput = beamformerFreq(True, steerVecFormulation, self.r_diag, normFactor * nMics**2, (self.r0, self.rm, kji, eva, eve))  // takes all EigVal into account
                if self.r_diag:  // set (unphysical) negative output values to 0
                    indNegSign = sign(beamformerOutput) < 0
                    beamformerOutput[indNegSign] = 0.0
                ac[i] = (beamformerOutput ** self.gamma) * normFactor
                fr[i] = True
            
class BeamformerCapon( BeamformerBase ):

Instances