b273190e8ba6faac064cc81f192ef226fd50c606,acoular/fbeamform.py,BeamformerCapon,calc,#BeamformerCapon#Any#Any#,448

Before Change

        // prepare calculation
        kj = 2j*pi*self.freq_data.fftfreq()/self.c
        numchannels = self.freq_data.numchannels
        e = zeros((numchannels), "D")
        h = zeros((1, self.grid.size), "d")
        beamfunc = self.get_beamfunc()
        for i in self.freq_data.indices:
            if not fr[i]:
                csm = array(linalg.inv(array(self.freq_data.csm[i], \
                        dtype="complex128")), order="C")[newaxis]
                //print csm.flags
                kji = kj[i, newaxis]
                beamfunc(csm, e, h, self.r0, self.rm, kji)
                ac[i] = 1.0/h
                fr[i] = True

class BeamformerEig( BeamformerBase ):

After Change

                
        kj = 2j*pi*self.freq_data.fftfreq()/self.c
        nMics = self.freq_data.numchannels
        normFactor = self.signalLossNormalize() * nMics**2
        steerVecFormulation = self.steerVecTranslation()
        for i in self.freq_data.indices:
            if not fr[i]:
                csm = array(linalg.inv(array(self.freq_data.csm[i], dtype="complex128")), order="C")[newaxis]
                kji = kj[i, newaxis]
                beamformerOutput = beamformerFreq(False, steerVecFormulation, self.r_diag, normFactor, (self.r0, self.rm, kji, csm))
                ac[i] = 1.0 / beamformerOutput
                fr[i] = True

class BeamformerEig( BeamformerBase ):

Instances