97a543976472f88b7922cc8f8ef3d7c0b6cd3dee,beamfpy/sources.py,PointSource,result,#PointSource#Any#,146

Before Change

        signal = fft.fft(self.signal.signal())
        pos = array(self.loc, dtype=float).reshape(3, 1)
        rm = self.env.r(self.c, pos, self.mpos.mpos)
        delays = exp(-2j*pi*(rm/self.c)*\
            fft.fftfreq(int(self.numsamples),1.0/self.sample_freq)[:,newaxis])
        out = fft.ifft(signal[:, newaxis]*delays, axis=0).real/rm
        i = 0
        while i < self.numsamples:
            yield out[i:i+num]
            i += num

After Change

        signal = self.signal.usignal(self.up)
        out = empty((num, self.numchannels))
        // distances
        rm = self.env.r(self.c, array(self.loc).reshape((3, 1)), self.mpos.mpos)
        // emission time relative to start_t (in samples) for first sample
        ind = (-rm/self.c-self.start_t+self.start)*self.sample_freq   
        i = 0
        n = self.numsamples        
        while n:
            n -= 1
            try:
                out[i] = signal[array(0.5+ind*self.up, dtype=long)]/rm
                ind += 1.
                i += 1
                if i == num:
                    yield out

Instances