from pythran.typing import List, NDArray
from pythran.tests import TestFromDir, TestEnv

import os
import numpy as np
import platform

# from http://www.scipy.org/Download , weave/example directory

class TestScipy(TestFromDir):

    def test_laplace(self):
        code="""
def laplace(u,dx, dy):
    nx, ny=len(u), len(u[0])
    for i in range(1, nx-1):
        for j in range(1, ny-1):
            u[i][j] = ((u[i-1][j] + u[i+1][j])*dy**2 +
                      (u[i][j-1] + u[i][j+1])*dx**2)/(2.0*(dx**2 + dy**2))
"""
        self.run_test(code, [[0.1,0.2,0.3],[0.1,0.2,0.3],[0.1,0.2,0.3]], 0.01, 0.02, laplace=[List[List[float]], float, float])

    def test_recursive_fibonnaci(self):
        code="""
def recursive_fibonnaci(a):
    if a <= 2:
        return 1
    else:
        return recursive_fibonnaci(a-2) + recursive_fibonnaci(a-1)
        """
        self.run_test(code, 5, recursive_fibonnaci=[int])

    def test_iterative_fibonnaci(self):
        code="""
def iterative_fibonnaci(a):
    if a <= 2:
        return 1
    last = next_to_last = 1
    i = 2
    while i < a:
        result = last + next_to_last
        next_to_last = last
        last = result
        i+=1
    return result;
"""
        self.run_test(code, 5, iterative_fibonnaci=[int])

    def test_binary_search(self):
        code="""
def binary_search(seq, t):
    min = 0; max = len(seq) - 1
    while 1:
        if max < min:
            return -1
        m = (min + max) // 2
        if seq[m] < t:
            min = m + 1
        elif seq[m] > t:
            max = m - 1
        else:
            return m
"""
        self.run_test(code,[1,2,3,4,5,6,7,8,9], 4, binary_search=[List[int], int])

    def test_ramp(self):
        code="""
def ramp(result, start, end):
    size=len(result)
    assert size > 1
    step = (end-start)//(size-1)
    for i in range(size):
        result[i] = start + step*i
"""
        self.run_test(code,[0. for x in range(10)], 1.5, 9.5, ramp=[List[float], float, float])

    path = os.path.join(os.path.dirname(__file__), "scipy")

TestScipy.populate(TestScipy)

class TestPyData(TestFromDir):
    path = os.path.join(os.path.dirname(__file__), "pydata")

TestPyData.populate(TestPyData)

class TestScikitImage(TestFromDir):
    path = os.path.join(os.path.dirname(__file__), "scikit-image")

TestPyData.populate(TestScikitImage)

try:
    import scipy

    if not platform.machine().startswith('ppc64'):

        class TestScipySpecial(TestEnv):
            def test_jv_scalar(self):
                self.run_test("""
                from scipy.special import jv
                def jv_scalar(v, x):
                    return jv(v, x)""",
                    5, 1.414,
                    jv_scalar=[int, float])
            def test_spherical_jn_scalar(self):
                self.run_test("""
                from scipy.special import spherical_jn
                def spherical_bessel_j_scalar(v, x):
                    return spherical_jn(v, x)""",
                              5, 1.414,
                              spherical_bessel_j_scalar=[int, float])

            def test_spherical_jn_arg1d(self):
                self.run_test("""
                from scipy.special import spherical_jn
                def spherical_bessel_j_1d(v, x):
                    return spherical_jn(v, x)""",
                              5, np.array([1.0, 2.0, 3.0]),
                              spherical_bessel_j_1d=[int, NDArray[float,:]])

            def test_spherical_jn_order1d(self):
                self.run_test("""
                from scipy.special import spherical_jn
                def spherical_bessel_j_order1d(v, x):
                    return spherical_jn(v, x)""",
                              np.array([1, 2, 3]), 5.0,
                              spherical_bessel_j_order1d=[NDArray[int,:], float])

            def test_spherical_jn_arg2d(self):
                self.run_test("""
                from scipy.special import spherical_jn
                def spherical_bessel_j_2d(v, x):
                    return spherical_jn(v, x)""",
                              5, np.array([[1.0, 2.0], [3.0, 4.0]]),
                              spherical_bessel_j_2d=[int, NDArray[float,:,:]])

            def test_binom_scalar(self):
                self.run_test("""
                from scipy.special import binom
                def binom_scalar(v, x):
                    return binom(v, x)""",
                              5, 4,
                              binom_scalar=[int, int])

            def test_binom_error(self):
                self.run_test("""
                from scipy.special import binom
                def binom_error(v, x):
                    return binom(v, x)""",
                              1201, 800,
                              binom_error=[int, int])

            def test_binom_arg1d(self):
                self.run_test("""
                from scipy.special import binom
                def binom_1d(v, x):
                    return binom(v, x)""",
                              5, np.array([1, 2, 3]),
                              binom_1d=[int, NDArray[int,:]])

except ImportError:
    pass
