import os
import numpy as np
import scipy.io.wavfile as wavfile
from nose.tools import assert_almost_equals
import simpl.base as base

float_precision = 5
frame_size = 512
hop_size = 512
audio_path = os.path.join(
    os.path.dirname(__file__), 'audio/flute.wav'
)


class TestFrame(object):
    def test_buffers(self):
        N = 256
        f = base.Frame(N)
        assert f.size == N

        a = np.random.rand(N)
        f.audio = a
        assert np.all(f.audio == a)

        a = np.random.rand(N)
        f.synth = a
        assert np.all(f.synth == a)

        a = np.random.rand(N)
        f.residual = a
        assert np.all(f.residual == a)

        a = np.random.rand(N)
        f.synth_residual = a
        assert np.all(f.synth_residual == a)

    def test_peaks(self):
        p = base.Peak()
        p.amplitude = 0.5
        p.frequency = 220.0
        p.phase = 0.0

        f = base.Frame()
        assert f.num_peaks == 0
        assert f.max_peaks > 0
        f.add_peak(p)

        assert f.num_peaks == 1
        assert_almost_equals(f.peak(0).amplitude, p.amplitude,
                             float_precision)
        assert_almost_equals(f.peaks[0].amplitude, p.amplitude,
                             float_precision)

        f.clear()
        assert f.num_peaks == 0

    def test_partials(self):
        N = 256
        f = base.Frame(N)
        f.max_partials = 10

        p = base.Peak()
        p.amplitude = 0.5
        p.frequency = 220.0
        p.phase = 0.0

        f.partial(0, p)
        assert_almost_equals(f.partial(0).amplitude, p.amplitude,
                             float_precision)
        assert_almost_equals(f.partial(0).frequency, p.frequency,
                             float_precision)


class TestPeakDetection(object):
    @classmethod
    def setup_class(cls):
        cls.audio = wavfile.read(audio_path)[1]
        cls.audio = np.asarray(cls.audio, dtype=np.double)
        cls.audio /= np.max(cls.audio)

    def test_peak_detection(self):
        pd = base.PeakDetection()
        pd.find_peaks(self.audio)

        assert len(pd.frames) == len(self.audio) / hop_size
        assert len(pd.frames[0].peaks) == 0


class TestPartialTracking(object):
    @classmethod
    def setup_class(cls):
        cls.audio = wavfile.read(audio_path)[1]
        cls.audio = np.asarray(cls.audio, dtype=np.double)
        cls.audio /= np.max(cls.audio)

    def test_partial_tracking(self):
        pd = base.PeakDetection()
        frames = pd.find_peaks(self.audio)

        pt = base.PartialTracking()
        frames = pt.find_partials(frames)

        assert len(frames) == len(self.audio) / hop_size
        assert len(frames[0].partials) == 100


class TestSynthesis(object):
    @classmethod
    def setup_class(cls):
        cls.audio = wavfile.read(audio_path)[1]
        cls.audio = np.asarray(cls.audio, dtype=np.double)
        cls.audio /= np.max(cls.audio)

    def test_synthesis(self):
        pd = base.PeakDetection()
        frames = pd.find_peaks(self.audio)

        pt = base.PartialTracking()
        frames = pt.find_partials(frames)

        s = base.Synthesis()
        synth_audio = s.synth(frames)

        assert len(synth_audio)


class TestResidual(object):
    @classmethod
    def setup_class(cls):
        cls.audio = wavfile.read(audio_path)[1]
        cls.audio = np.asarray(cls.audio, dtype=np.double)
        cls.audio /= np.max(cls.audio)

    def test_synthesis(self):
        pd = base.PeakDetection()
        frames = pd.find_peaks(self.audio)

        pt = base.PartialTracking()
        frames = pt.find_partials(frames)

        s = base.Synthesis()
        synth_audio = s.synth(frames)

        r = base.Residual()
        residual_audio = r.find_residual(synth_audio, self.audio)
        assert len(residual_audio)