# Copyright (C) 2021-2024 Université Gustave Eiffel.
# Copyright (C) 2025-2026 Université Gustave Eiffel, INRIA.
# This file is part of the EasyFEA project.
# EasyFEA is distributed under the terms of the GNU General Public License v3, see LICENSE.txt and CREDITS.md for more information.

"""
Mesh11
======

Meshing of a specimen for a spatially oriented tensile test.
"""
# sphinx_gallery_thumbnail_number = 2

import numpy as np

from EasyFEA import Display, ElemType, PyVista
from EasyFEA.Geoms import Point, Line, CircleArc, Contour, Domain

if __name__ == "__main__":
    Display.Clear()

    # ----------------------------------------------
    # Configuration
    # ----------------------------------------------

    # geom
    L = 1
    H = 2
    e = L * 0.5

    # ----------------------------------------------
    # Mesh
    # ----------------------------------------------

    p1 = Point(-L / 2)
    p2 = Point(L / 2)
    p3 = p2 + [0, H]
    p4 = p1 + [0, H]

    p5 = Point(e / 2, H / 2)
    p6 = Point(-e / 2, H / 2)

    l1 = Line(p1, p2)
    l2 = CircleArc(p2, p3, P=p5)
    l3 = Line(p3, p4)
    l4 = CircleArc(p4, p1, P=p6)

    contour = Contour([l1, l2, l3, l4])
    contour2 = Domain(p1 - [0, H / 2], p2, isHollow=False)
    contour3 = contour2.copy()
    contour3.Translate(dy=H + H / 2)

    surfaces = [contour2, contour3]

    PyVista.Plot_Geoms([contour, *surfaces]).show()

    mesh = contour.Mesh_Extrude(
        [],
        [0, 0, e],
        [3],
        isOrganised=True,
        elemType=ElemType.HEXA8,
        additionalSurfaces=surfaces,
    )

    # ----------------------------------------------
    # Update coords
    # ----------------------------------------------

    oldArea = mesh.area
    mesh.Rotate(-45, mesh.center)
    assert np.abs(mesh.area - oldArea) / oldArea <= 1e-12
    mesh.Rotate(45, mesh.center, (1, 0))
    assert np.abs(mesh.area - oldArea) / oldArea <= 1e-12

    PyVista.Plot_Mesh(mesh).show()