from __future__ import annotations

import importlib

import numpy as np
import pytest
import requests

from pymatgen.core import SETTINGS, Element
from pymatgen.entries.computed_entries import ComputedEntry
from pymatgen.ext.matproj import MPRester
from pymatgen.phonon.bandstructure import PhononBandStructureSymmLine
from pymatgen.phonon.dos import CompletePhononDos

# Skip all MPRester tests if some downstream problem on the website, mp-api or whatever.
if not (mp_api_key := SETTINGS.get("PMG_MAPI_KEY")):
    pytest.skip("PMG_MAPI_KEY not set", allow_module_level=True)

if len(mp_api_key) != 32:
    pytest.fail(f"Invalid/old MP API key, expected key length 32, got {len(mp_api_key)}")

MP_URL: str = "https://api.materialsproject.org"

try:
    mp_api_down: bool = requests.get(MP_URL, timeout=10).status_code != 200
except (requests.exceptions.ConnectionError, requests.exceptions.Timeout):
    mp_api_down = True

if mp_api_down:
    pytest.skip("MP API is down", allow_module_level=True)

if importlib.util.find_spec("pymatgen.analysis.alloys") is None:
    pytest.skip("pymatgen-analysis-alloys plugin is not installed", allow_module_level=True)


@pytest.fixture(scope="module")
def mprester():
    return MPRester()


def test_get_summary(mprester):
    docs = mprester.get_summary({"formula": "Fe2O3"})
    assert len(docs) > 3

    mid = "mp-19770"
    doc = mprester.get_summary_by_material_id(mid)
    assert doc["formula_pretty"] == "Fe2O3"

    doc = mprester.summary.search(material_ids="mp-19770,mp-19017", _fields="formula_pretty,energy_above_hull")
    assert len(doc) == 2
    assert len(doc[0]) == 2
    assert doc[0]["energy_above_hull"] >= 0
    assert doc[1]["energy_above_hull"] >= 0


def test_get_all_materials_ids_doc(mprester):
    mids = mprester.get_material_ids("Al2O3")
    np.random.default_rng().shuffle(mids)
    doc = mprester.get_doc(mids.pop(0))
    assert doc["formula_pretty"] == "Al2O3"


def test_get_entries_and_in_chemsys(mprester):
    # One large system test.
    syms = ["Li", "Fe", "O", "P", "Mn"]

    # Small test.
    syms2 = "Fe-Li-O"
    entries = mprester.get_entries_in_chemsys(syms)

    elements = {Element(sym) for sym in syms}
    for entry in entries:
        assert isinstance(entry, ComputedEntry)
        assert set(entry.elements).issubset(elements)

    assert len(entries) > 1000

    entries2 = mprester.get_entries(
        syms2, property_data=["formation_energy_per_atom", "uncorrected_energy_per_atom"], summary_data=["band_gap"]
    )
    for entry in entries2:
        assert isinstance(entry, ComputedEntry)
        assert set(entry.elements).issubset(elements)
        assert "formation_energy_per_atom" in entry.data
        assert entry.data["uncorrected_energy_per_atom"] == pytest.approx(entry.uncorrected_energy_per_atom, abs=1e-3)
        assert "band_gap" in entry.data["summary"]
    assert len(entries2) < 1000

    e1 = {i.entry_id for i in entries}
    e2 = {i.entry_id for i in entries2}
    assert e1.issuperset(e2)

    # This gets everything in Li-Fe-O and within this subsystem. It should have more entries than
    # get_entries("Li-Fe-O"), which just gets only the ternary compounds.
    entries3 = mprester.get_entries_in_chemsys(["Fe", "Li", "O"])

    assert len(entries3) > len(entries2)
    with pytest.warns(DeprecationWarning, match="arguments are deprecated"):
        mprester.get_entries_in_chemsys("Li", inc_structure=True)


def test_get_structure_by_material_id(mprester):
    s1 = mprester.get_structure_by_material_id("mp-1")
    assert s1.formula == "Cs1"


def test_get_entry_by_material_id(mprester):
    entry = mprester.get_entry_by_material_id("mp-19017")
    assert isinstance(entry, ComputedEntry)
    assert entry.reduced_formula == "LiFePO4"

    with pytest.raises(IndexError, match="list index out of range"):
        mprester.get_entry_by_material_id("mp-2022")  # "mp-2022" does not exist


def test_get_phonon_data_by_material_id(mprester):
    bs = mprester.get_phonon_bandstructure_by_material_id("mp-661")
    assert isinstance(bs, PhononBandStructureSymmLine)
    dos = mprester.get_phonon_dos_by_material_id("mp-661")
    assert isinstance(dos, CompletePhononDos)


def test_get_structures(mprester):
    structs = mprester.get_structures("Mn3O4")
    assert len(structs) > 0


def test_get_entries(mprester):
    entries = mprester.get_entries("TiO2")
    assert len(entries) > 1
    for entry in entries:
        assert entry.reduced_formula == "TiO2"

    entries = mprester.get_entries("TiO2")
    assert len(entries) > 1
    for entry in entries:
        assert entry.structure.reduced_formula == "TiO2"

    all_entries = mprester.get_entries("Fe", compatible_only=False)
    entries = mprester.get_entries("Fe", compatible_only=True)
    assert len(entries) < len(all_entries)

    for entry in mprester.get_entries("CdO2"):
        assert entry.data["oxide_type"] is not None

    # test if it will retrieve the conventional unit cell of Ni
    entry = mprester.get_entry_by_material_id("mp-23")
    Ni = entry.structure
    assert Ni.lattice.a == pytest.approx(Ni.lattice.b)
    assert Ni.lattice.a == pytest.approx(Ni.lattice.c)
    assert Ni.lattice.alpha == pytest.approx(60)
    assert Ni.lattice.beta == pytest.approx(60)
    assert Ni.lattice.gamma == pytest.approx(60)

    # Ensure energy per atom is same
    primNi = mprester.get_entry_by_material_id("mp-23")
    assert primNi.energy_per_atom == entry.energy_per_atom

    Ni = mprester.get_structure_by_material_id("mp-23", conventional_unit_cell=True)
    assert Ni.lattice.a == pytest.approx(Ni.lattice.b)
    assert Ni.lattice.a == pytest.approx(Ni.lattice.c)
    assert Ni.lattice.alpha == pytest.approx(90)
    assert Ni.lattice.beta == pytest.approx(90)
    assert Ni.lattice.gamma == pytest.approx(90)


def test_api_parity(mprester):
    docs = (
        "summary",
        "core",
        "elasticity",
        "phonon",
        "eos",
        "similarity",
        "xas",
        "electronic_structure",
        "robocrys",
        "synthesis",
        "magnetism",
        "oxidation_states",
        "provenance",
        "alloys",
        "chemenv",
        "bonds",
        "dielectric",
    )

    for doc in docs:
        # We should have Al2O3 data for these properties.
        data = mprester.materials.__getattribute__(doc).search(material_ids="mp-1143")
        assert len(data) > 0, f"No Al2O3 data returned for {doc}"
        data = mprester.__getattribute__(doc).search(material_ids="mp-1143")
        assert len(data) > 0, f"No Al2O3 data returned for {doc}"

    data = mprester.materials.substrates.search(sub_id="mp-1143")
    assert len(data) > 0, "No substrate data returned."

    # data = mprester.materials.tasks.search(task_ids="mp-1143")
    # assert len(data) > 0, "No tasks data returned."

    docs = ("surface_properties", "alloys")

    for doc in docs:
        data = mprester.materials.__getattribute__(doc).search(material_ids="mp-135")
        assert len(data) > 0, f"No Li data returned for {doc}"

    data = mprester.materials.insertion_electrodes.search(formula="LiFePO4")
    assert len(data) > 0, "No insertion electrode data returned."

    # TODO: Test these docs: "grain_boundaries", "piezoelectric", "conversion_electrodes"
