pvlib 0.11.2__py3-none-any.whl → 0.12.1a1__py3-none-any.whl

pvlib/tests/test_inverter.py

@@ -1,213 +0,0 @@
-import numpy as np
-import pandas as pd
-
-from .conftest import assert_series_equal
-from numpy.testing import assert_allclose
-
-from .conftest import DATA_DIR
-import pytest
-
-from pvlib import inverter
-
-
-def test_adr(adr_inverter_parameters):
-    vdcs = pd.Series([135, 154, 390, 420, 551])
-    pdcs = pd.Series([135, 1232, 1170, 420, 551])
-
-    pacs = inverter.adr(vdcs, pdcs, adr_inverter_parameters)
-    assert_series_equal(pacs, pd.Series([np.nan, 1161.5745, 1116.4459,
-                                         382.6679, np.nan]))
-
-
-def test_adr_vtol(adr_inverter_parameters):
-    vdcs = pd.Series([135, 154, 390, 420, 551])
-    pdcs = pd.Series([135, 1232, 1170, 420, 551])
-
-    pacs = inverter.adr(vdcs, pdcs, adr_inverter_parameters, vtol=0.20)
-    assert_series_equal(pacs, pd.Series([104.8223, 1161.5745, 1116.4459,
-                                         382.6679, 513.3385]))
-
-
-def test_adr_float(adr_inverter_parameters):
-    vdcs = 154.
-    pdcs = 1232.
-
-    pacs = inverter.adr(vdcs, pdcs, adr_inverter_parameters)
-    assert_allclose(pacs, 1161.5745)
-
-
-def test_adr_invalid_and_night(sam_data):
-    # also tests if inverter.adr can read the output from pvsystem.retrieve_sam
-    inverters = sam_data['adrinverter']
-    testinv = 'Zigor__Sunzet_3_TL_US_240V__CEC_2011_'
-    vdcs = np.array([39.873036, 0., np.nan, 420])
-    pdcs = np.array([188.09182, 0., 420, np.nan])
-
-    pacs = inverter.adr(vdcs, pdcs, inverters[testinv])
-    assert_allclose(pacs, np.array([np.nan, -0.25, np.nan, np.nan]))
-
-
-def test_sandia(cec_inverter_parameters):
-    vdcs = pd.Series(np.linspace(0, 50, 3))
-    idcs = pd.Series(np.linspace(0, 11, 3))
-    pdcs = idcs * vdcs
-
-    pacs = inverter.sandia(vdcs, pdcs, cec_inverter_parameters)
-    assert_series_equal(pacs, pd.Series([-0.020000, 132.004308, 250.000000]))
-
-
-def test_sandia_float(cec_inverter_parameters):
-    vdcs = 25.
-    idcs = 5.5
-    pdcs = idcs * vdcs
-    pacs = inverter.sandia(vdcs, pdcs, cec_inverter_parameters)
-    assert_allclose(pacs, 132.004278, 1e-5)
-    # test at low power condition
-    vdcs = 25.
-    idcs = 0
-    pdcs = idcs * vdcs
-    pacs = inverter.sandia(vdcs, pdcs, cec_inverter_parameters)
-    assert_allclose(pacs, -1. * cec_inverter_parameters['Pnt'], 1e-5)
-
-
-def test_sandia_Pnt_micro():
-    """
-    Test for issue #140, where some microinverters were giving a positive AC
-    power output when the DC power was 0.
-    """
-    inverter_parameters = {
-        'Name': 'Enphase Energy: M250-60-2LL-S2x (-ZC) (-NA) 208V [CEC 2013]',
-        'Vac': 208.0,
-        'Paco': 240.0,
-        'Pdco': 250.5311318,
-        'Vdco': 32.06160667,
-        'Pso': 1.12048857,
-        'C0': -5.76E-05,
-        'C1': -6.24E-04,
-        'C2': 8.09E-02,
-        'C3': -0.111781106,
-        'Pnt': 0.043,
-        'Vdcmax': 48.0,
-        'Idcmax': 9.8,
-        'Mppt_low': 27.0,
-        'Mppt_high': 39.0,
-    }
-    vdcs = pd.Series(np.linspace(0, 50, 3))
-    idcs = pd.Series(np.linspace(0, 11, 3))
-    pdcs = idcs * vdcs
-
-    pacs = inverter.sandia(vdcs, pdcs, inverter_parameters)
-    assert_series_equal(pacs, pd.Series([-0.043, 132.545914746, 240.0]))
-
-
-def test_sandia_multi(cec_inverter_parameters):
-    vdcs = pd.Series(np.linspace(0, 50, 3))
-    idcs = pd.Series(np.linspace(0, 11, 3)) / 2
-    pdcs = idcs * vdcs
-    pacs = inverter.sandia_multi((vdcs, vdcs), (pdcs, pdcs),
-                                 cec_inverter_parameters)
-    assert_series_equal(pacs, pd.Series([-0.020000, 132.004278, 250.000000]))
-    # with lists instead of tuples
-    pacs = inverter.sandia_multi([vdcs, vdcs], [pdcs, pdcs],
-                                 cec_inverter_parameters)
-    assert_series_equal(pacs, pd.Series([-0.020000, 132.004278, 250.000000]))
-    # with arrays instead of tuples
-    pacs = inverter.sandia_multi(np.array([vdcs, vdcs]),
-                                 np.array([pdcs, pdcs]),
-                                 cec_inverter_parameters)
-    assert_allclose(pacs, np.array([-0.020000, 132.004278, 250.000000]))
-
-
-def test_sandia_multi_length_error(cec_inverter_parameters):
-    vdcs = pd.Series(np.linspace(0, 50, 3))
-    idcs = pd.Series(np.linspace(0, 11, 3))
-    pdcs = idcs * vdcs
-    with pytest.raises(ValueError, match='p_dc and v_dc have different'):
-        inverter.sandia_multi((vdcs,), (pdcs, pdcs), cec_inverter_parameters)
-
-
-def test_sandia_multi_array(cec_inverter_parameters):
-    vdcs = np.linspace(0, 50, 3)
-    idcs = np.linspace(0, 11, 3)
-    pdcs = idcs * vdcs
-    pacs = inverter.sandia_multi((vdcs,), (pdcs,), cec_inverter_parameters)
-    assert_allclose(pacs, np.array([-0.020000, 132.004278, 250.000000]))
-
-
-def test_pvwatts_scalars():
-    expected = 85.58556604752516
-    out = inverter.pvwatts(90, 100, 0.95)
-    assert_allclose(out, expected)
-    # GH 675
-    expected = 0.
-    out = inverter.pvwatts(0., 100)
-    assert_allclose(out, expected)
-
-
-def test_pvwatts_possible_negative():
-    # pvwatts could return a negative value for (pdc / pdc0) < 0.006
-    # unless it is clipped. see GH 541 for more
-    expected = 0
-    out = inverter.pvwatts(0.001, 1)
-    assert_allclose(out, expected)
-
-
-def test_pvwatts_arrays():
-    pdc = np.array([[np.nan], [0], [50], [100]])
-    pdc0 = 100
-    expected = np.array([[np.nan],
-                         [0.],
-                         [47.60843624],
-                         [95.]])
-    out = inverter.pvwatts(pdc, pdc0, 0.95)
-    assert_allclose(out, expected, equal_nan=True)
-
-
-def test_pvwatts_series():
-    pdc = pd.Series([np.nan, 0, 50, 100])
-    pdc0 = 100
-    expected = pd.Series(np.array([np.nan, 0., 47.608436, 95.]))
-    out = inverter.pvwatts(pdc, pdc0, 0.95)
-    assert_series_equal(expected, out)
-
-
-def test_pvwatts_multi():
-    pdc = np.array([np.nan, 0, 50, 100]) / 2
-    pdc0 = 100
-    expected = np.array([np.nan, 0., 47.608436, 95.])
-    out = inverter.pvwatts_multi((pdc, pdc), pdc0, 0.95)
-    assert_allclose(expected, out)
-    # with 2D array
-    pdc_2d = np.array([pdc, pdc])
-    out = inverter.pvwatts_multi(pdc_2d, pdc0, 0.95)
-    assert_allclose(expected, out)
-    # with Series
-    pdc = pd.Series(pdc)
-    out = inverter.pvwatts_multi((pdc, pdc), pdc0, 0.95)
-    assert_series_equal(pd.Series(expected), out)
-    # with list instead of tuple
-    out = inverter.pvwatts_multi([pdc, pdc], pdc0, 0.95)
-    assert_series_equal(pd.Series(expected), out)
-
-
-INVERTER_TEST_MEAS = DATA_DIR / 'inverter_fit_snl_meas.csv'
-INVERTER_TEST_SIM = DATA_DIR / 'inverter_fit_snl_sim.csv'
-
-
-@pytest.mark.parametrize('infilen, expected', [
-    (INVERTER_TEST_MEAS, {'Paco': 333000., 'Pdco': 343251., 'Vdco': 740.,
-                          'Pso': 1427.746, 'C0': -5.768e-08, 'C1': 3.596e-05,
-                          'C2': 1.038e-03, 'C3': 2.978e-05, 'Pnt': 1.}),
-    (INVERTER_TEST_SIM,  {'Paco': 1000., 'Pdco': 1050., 'Vdco': 240.,
-                          'Pso': 10., 'C0': 1e-6, 'C1': 1e-4, 'C2': 1e-2,
-                          'C3': 1e-3, 'Pnt': 1.}),
-])
-def test_fit_sandia(infilen, expected):
-    curves = pd.read_csv(infilen)
-    dc_power = curves['ac_power'] / curves['efficiency']
-    result = inverter.fit_sandia(ac_power=curves['ac_power'],
-                                 dc_power=dc_power,
-                                 dc_voltage=curves['dc_voltage'],
-                                 dc_voltage_level=curves['dc_voltage_level'],
-                                 p_ac_0=expected['Paco'], p_nt=expected['Pnt'])
-    assert expected == pytest.approx(result, rel=1e-3)