# =============================================================================
# Copyright (C) 2010 Diego Duclos
#
# This file is part of pyfa.
#
# pyfa is free software: you can redistribute it and/or modify
# it under the terms of the GNU General Public License as published by
# the Free Software Foundation, either version 3 of the License, or
# (at your option) any later version.
#
# pyfa is distributed in the hope that it will be useful,
# but WITHOUT ANY WARRANTY; without even the implied warranty of
# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
# GNU General Public License for more details.
#
# You should have received a copy of the GNU General Public License
# along with pyfa.  If not, see <http://www.gnu.org/licenses/>.
# =============================================================================


import math
from abc import ABCMeta, abstractmethod


class PointGetter(metaclass=ABCMeta):

    def __init__(self, graph):
        self.graph = graph

    @abstractmethod
    def getRange(self, xRange, miscParams, src, tgt):
        raise NotImplementedError

    @abstractmethod
    def getPoint(self, x, miscParams, src, tgt):
        raise NotImplementedError


class SmoothPointGetter(PointGetter, metaclass=ABCMeta):

    _baseResolution = 200
    _extraDepth = 0

    def getRange(self, xRange, miscParams, src, tgt):
        xs = []
        ys = []
        commonData = self._getCommonData(miscParams=miscParams, src=src, tgt=tgt)

        def addExtraPoints(x1, y1, x2, y2, depth):
            if depth <= 0 or y1 == y2:
                return
            newX = (x1 + x2) / 2
            newY = self._calculatePoint(x=newX, miscParams=miscParams, src=src, tgt=tgt, commonData=commonData)
            addExtraPoints(x1=prevX, y1=prevY, x2=newX, y2=newY, depth=depth - 1)
            xs.append(newX)
            ys.append(newY)
            addExtraPoints(x1=newX, y1=newY, x2=x2, y2=y2, depth=depth - 1)

        prevX = None
        prevY = None
        # Go through X points defined by our resolution setting
        for x in self._xIterLinear(xRange):
            y = self._calculatePoint(x=x, miscParams=miscParams, src=src, tgt=tgt, commonData=commonData)
            if prevX is not None and prevY is not None:
                # And if Y values of adjacent data points are not equal, add extra points
                # depending on extra depth setting
                addExtraPoints(x1=prevX, y1=prevY, x2=x, y2=y, depth=self._extraDepth)
            prevX = x
            prevY = y
            xs.append(x)
            ys.append(y)
        return xs, ys

    def getPoint(self, x, miscParams, src, tgt):
        commonData = self._getCommonData(miscParams=miscParams, src=src, tgt=tgt)
        return self._calculatePoint(x=x, miscParams=miscParams, src=src, tgt=tgt, commonData=commonData)

    def _xIterLinear(self, xRange):
        xLow = min(xRange)
        xHigh = max(xRange)
        # Resolution defines amount of ranges between points here,
        # not amount of points
        step = (xHigh - xLow) / self._baseResolution
        if step == 0 or math.isnan(step):
            yield xLow
        else:
            for i in range(self._baseResolution + 1):
                yield xLow + step * i

    def _getCommonData(self, miscParams, src, tgt):
        return {}

    @abstractmethod
    def _calculatePoint(self, x, miscParams, src, tgt, commonData):
        raise NotImplementedError