pvlib 0.11.2__py3-none-any.whl → 0.12.0__py3-none-any.whl

pvlib/tests/test_numerical_precision.py

@@ -1,124 +0,0 @@
-"""
-Test numerical precision of explicit single diode calculation using symbolic
-mathematics. SymPy is a computer algebra system, that uses infinite precision
-symbols instead of standard floating point and integer computer number types.
-http://docs.sympy.org/latest/modules/evalf.html#accuracy-and-error-handling
-
-This module can be executed from the command line to generate a high precision
-dataset of I-V curve points to test the explicit single diode calculations
-:func:`pvlib.singlediode.bishop88`::
-
-    $ python test_numeric_precision.py
-
-This generates a file in the pvlib data folder, which is specified by the
-constant ``DATA_PATH``. When the test is run using ``pytest`` it will compare
-the values calculated by :func:`pvlib.singlediode.bishop88` with the
-high-precision values generated with SymPy.
-"""
-
-import logging
-import numpy as np
-import pandas as pd
-from pvlib import pvsystem
-from pvlib.singlediode import bishop88, estimate_voc
-from .conftest import DATA_DIR
-
-logging.basicConfig()
-LOGGER = logging.getLogger(__name__)
-LOGGER.setLevel(logging.DEBUG)
-TEST_DATA = 'bishop88_numerical_precision.csv'
-DATA_PATH = DATA_DIR / TEST_DATA
-POA = 888
-TCELL = 55
-# module parameters from CEC module SunPower SPR-E20-327
-SPR_E20_327 = {
-    'alpha_sc': 0.004522,
-    'a_ref': 2.6868,
-    'I_L_ref': 6.468,
-    'I_o_ref': 1.88e-10,
-    'R_s': 0.37,
-    'R_sh_ref': 298.13,
-}
-# apply temp/irrad desoto corrections
-ARGS = pvsystem.calcparams_desoto(
-    effective_irradiance=POA, temp_cell=TCELL,
-    EgRef=1.121, dEgdT=-0.0002677, **SPR_E20_327,
-)
-IL, I0, RS, RSH, NNSVTH = ARGS
-IVCURVE_NPTS = 100
-
-try:
-    from sympy import symbols, exp as sy_exp
-except ImportError as exc:
-    LOGGER.exception(exc)
-    symbols = NotImplemented
-    sy_exp = NotImplemented
-
-
-def generate_numerical_precision():  # pragma: no cover
-    """
-    Generate expected data with infinite numerical precision using SymPy.
-    :return: dataframe of expected values
-    """
-    if symbols is NotImplemented:
-        LOGGER.critical("SymPy is required to generate expected data.")
-        raise ImportError("could not import sympy")
-    il, io, rs, rsh, nnsvt, vd = symbols('il, io, rs, rsh, nnsvt, vd')
-    a = sy_exp(vd / nnsvt)
-    b = 1.0 / rsh
-    i = il - io * (a - 1.0) - vd * b
-    v = vd - i * rs
-    c = io * a / nnsvt
-    grad_i = - c - b  # di/dvd
-    grad_v = 1.0 - grad_i * rs  # dv/dvd
-    # dp/dv = d(iv)/dv = v * di/dv + i
-    grad = grad_i / grad_v  # di/dv
-    p = i * v
-    grad_p = v * grad + i  # dp/dv
-    grad2i = -c / nnsvt
-    grad2v = -grad2i * rs
-    grad2p = (
-        grad_v * grad + v * (grad2i/grad_v - grad_i*grad2v/grad_v**2) + grad_i
-    )
-    # generate exact values
-    data = dict(zip((il, io, rs, rsh, nnsvt), ARGS))
-    vdtest = np.linspace(0, estimate_voc(IL, I0, NNSVTH), IVCURVE_NPTS)
-    expected = []
-    for test in vdtest:
-        data[vd] = test
-        test_data = {
-            'i': np.float64(i.evalf(subs=data)),
-            'v': np.float64(v.evalf(subs=data)),
-            'p': np.float64(p.evalf(subs=data)),
-            'grad_i': np.float64(grad_i.evalf(subs=data)),
-            'grad_v': np.float64(grad_v.evalf(subs=data)),
-            'grad': np.float64(grad.evalf(subs=data)),
-            'grad_p': np.float64(grad_p.evalf(subs=data)),
-            'grad2p': np.float64(grad2p.evalf(subs=data))
-        }
-        LOGGER.debug(test_data)
-        expected.append(test_data)
-    return pd.DataFrame(expected, index=vdtest)
-
-
-def test_numerical_precision():
-    """
-    Test that there are no numerical errors due to floating point arithmetic.
-    """
-    expected = pd.read_csv(DATA_PATH)
-    vdtest = np.linspace(0, estimate_voc(IL, I0, NNSVTH), IVCURVE_NPTS)
-    results = bishop88(vdtest, *ARGS, gradients=True)
-    assert np.allclose(expected['i'], results[0])
-    assert np.allclose(expected['v'], results[1])
-    assert np.allclose(expected['p'], results[2])
-    assert np.allclose(expected['grad_i'], results[3])
-    assert np.allclose(expected['grad_v'], results[4])
-    assert np.allclose(expected['grad'], results[5])
-    assert np.allclose(expected['grad_p'], results[6])
-    assert np.allclose(expected['grad2p'], results[7])
-
-
-if __name__ == '__main__':  # pragma: no cover
-    expected = generate_numerical_precision()
-    expected.to_csv(DATA_PATH)
-    test_numerical_precision()

pvlib/tests/test_pvarray.py

@@ -1,71 +0,0 @@
-import numpy as np
-import pandas as pd
-from numpy.testing import assert_allclose
-from .conftest import assert_series_equal
-import pytest
-
-from pvlib import pvarray
-
-
-def test_pvefficiency_adr():
-    g = [1000, 200, 1000, 200, 1000, 200, 0.0, np.nan]
-    t = [25, 25, 50, 50, 75, 75, 25, 25]
-    params = [1.0, -6.651460, 0.018736, 0.070679, 0.054170]
-
-    # the expected values were calculated using the new function itself
-    # hence this test is primarily a regression test
-    eta = [1.0, 0.949125, 0.928148, 0.876472, 0.855759, 0.803281, 0.0, np.nan]
-
-    result = pvarray.pvefficiency_adr(g, t, *params)
-    assert_allclose(result, eta, atol=1e-6)
-
-
-def test_fit_pvefficiency_adr():
-    g = [1000, 200, 1000, 200, 1000, 200]
-    t = [25, 25, 50, 50, 75, 75]
-    eta = [1.0, 0.949125, 0.928148, 0.876472, 0.855759, 0.803281]
-
-    # the expected values were calculated using the new function itself
-    # hence this test is primarily a regression test
-    params = [1.0, -6.651460, 0.018736, 0.070679, 0.054170]
-
-    result = pvarray.fit_pvefficiency_adr(g, t, eta, dict_output=False)
-    # the fitted parameters vary somewhat by platform during the testing
-    # so the tolerance is higher on the parameters than on the efficiencies
-    # in the other tests
-    assert_allclose(result, params, rtol=1e-3)
-
-    result = pvarray.fit_pvefficiency_adr(g, t, eta, dict_output=True)
-    assert 'k_a' in result
-
-
-def test_pvefficiency_adr_round_trip():
-    g = [1000, 200, 1000, 200, 1000, 200]
-    t = [25, 25, 50, 50, 75, 75]
-    eta = [1.0, 0.949125, 0.928148, 0.876472, 0.855759, 0.803281]
-
-    params = pvarray.fit_pvefficiency_adr(g, t, eta, dict_output=False)
-    result = pvarray.pvefficiency_adr(g, t, *params)
-    assert_allclose(result, eta, atol=1e-6)
-
-
-def test_huld():
-    pdc0 = 100
-    res = pvarray.huld(1000, 25, pdc0, cell_type='cSi')
-    assert np.isclose(res, pdc0)
-    exp_sum = np.exp(1) * (np.sum(pvarray._infer_k_huld('cSi', pdc0)) + pdc0)
-    res = pvarray.huld(1000*np.exp(1), 26, pdc0, cell_type='cSi')
-    assert np.isclose(res, exp_sum)
-    res = pvarray.huld(100, 30, pdc0, k=(1, 1, 1, 1, 1, 1))
-    exp_100 = 0.1 * (pdc0 + np.log(0.1) + np.log(0.1)**2 + 5 + 5*np.log(0.1)
-                     + 5*np.log(0.1)**2 + 25)
-    assert np.isclose(res, exp_100)
-    # Series input, and irradiance = 0
-    eff_irr = pd.Series([1000, 100, 0])
-    tm = pd.Series([25, 30, 30])
-    expected = pd.Series([pdc0, exp_100, 0])
-    res = pvarray.huld(eff_irr, tm, pdc0, k=(1, 1, 1, 1, 1, 1))
-    assert_series_equal(res, expected)
-    with pytest.raises(ValueError,
-                       match='Either k or cell_type must be specified'):
-        res = pvarray.huld(1000, 25, 100)