pysdic.triangle_3_compute_elements_areas

triangle_3_compute_elements_areas(vertices_coordinates, connectivity)

Compute the areas of triangular elements defined by the given vertices and connectivity.

The area of each triangular element is computed using the cross product of two edges of the triangle.

\[\text{Area} = \frac{1}{2} \| \vec{AB} \times \vec{AC} \|\]

Note

The input vertices_coordinates will be converted to numpy.float64 for computation. The input connectivity will be converted to numpy.int64 for computation.

Warning

No tests are performed to check if the inputs array content are consistent (e.g., if the connectivity contains invalid vertex indices, …). Only the shapes and types of the input arrays are validated. The behavior of the function is undefined if the input arrays are not consistent.

Parameters:

• vertices_coordinates (ArrayLike) – An array of shape \((N_v, 3)\) representing the coordinates of the vertices.

• connectivity (ArrayLike) – An array of shape \((N_e, 3)\) representing the connectivity of the triangular elements.

Returns:

An array of shape \((N_e,)\) representing the area of each triangular element.

Return type:

numpy.ndarray

Raises:

• TypeError – If the inputs are not numpy arrays.

• ValueError – If the input arrays do not have the correct shapes.

See also

pysdic.triangle_3_compute_elements_normals

For computing the normal vectors of triangular elements.

pysdic.triangle_3_compute_vertices_normals

For computing the normal vectors at each vertex of the mesh.

Examples

Computing the areas of triangular elements:

import numpy as np
from pysdic import triangle_3_compute_elements_areas

vertices = np.array([[0, 0, 0],
                     [1, 0, 0],
                     [0, 1, 0],
                     [0, 0, 1]])  # shape (4, 3)

connectivity = np.array([[0, 1, 2],
                         [0, 1, 3],
                         [0, 2, 3],
                         [1, 2, 3]])  # shape (4, 3)

areas = triangle_3_compute_elements_areas(vertices, connectivity)
print(areas)

[0.5, 0.5, 0.5, 0.8660254]