"""
.. _example_mesh_simple_workflow:

Simple mesh workflow 
=====================================

.. contents:: Table of Contents
   :local:
   :depth: 2
   :backlinks: top


This example demonstrates a simple workflow for creating, manipulating,
and visualizing a mesh using the ``pysdic`` library.

.. seealso::

    :class:`pysdic.Mesh` - Official documentation for the Mesh class.

"""

# %%
# Creating a Mesh
# ---------------------------
#
# To create a mesh, we first need to give the coordinates of the vertices and the connectivity.
#
# Additionally some convenience methods are provided to save and load point clouds
# from files.

from pysdic import create_triangle_3_heightmap
import numpy as np

surface_mesh = create_triangle_3_heightmap(
    height_function=lambda x, y: 0.5 * np.sin(np.pi * x) * np.cos(np.pi * y),
    x_bounds=(-1.0, 1.0),
    y_bounds=(-1.0, 1.0),
    n_x=50,
    n_y=50,
)

# %%
# Add some point properties
# ---------------------------
# We can add properties to the vertices in the mesh. Here, we add a random scalar property
# to each point.
#
# Each property can be multi-dimensional, e.g., a vector or tensor associated with each point.
# All the properties are stored in a dictionary-like structure with shape :math:`(N_p, A)`.

heigth = surface_mesh.points[:, 2]  # Get the height (z-coordinate) of each point
intensity = (heigth - heigth.min()) / (heigth.max() - heigth.min())  # Normalize height to [0, 1] for intensity
surface_mesh.vertices["intensity"] = intensity  # Add intensity as a point property

# %%
# Visualize the Mesh (Several options available)
# ------------------------------------------------
#
# The point cloud can be visualized using the built-in visualization method with ``pyvista``.
# 
# .. seealso::
#   
#    :meth:`pysdic.Mesh.visualize` - Official documentation for the visualize method.
#

surface_mesh.visualize(
    face_color=None,
    title="Example Point Cloud with Texture",
    vertex_size=10,
    bounds_grid="back",
    property_cmap="inferno",
    vertex_property="intensity",
    show_vertices=False,
    show_edges=True,
)