honeycomb_kernels/remeshing/
relaxation.rs

1use honeycomb_core::{
2    cmap::{CMap2, VertexIdType},
3    geometry::{CoordsFloat, Vertex2},
4    stm::{StmClosureResult, Transaction},
5};
6
7/// Move a vertex to the average of the others' values.
8///
9/// This function computes the average of a list of coordinates and assigns that value to the
10/// specified vertex.
11///
12/// # Arguments
13///
14/// - `t: &mut Transaction` -- Associated transaction.
15/// - `map: &mut CMap2` -- Edited map.
16/// - `vid: VertexIdType` -- Vertex to move.
17/// - `others: &[VertexIdType]` -- List of vertex to compute the average from.
18///
19/// # Errors
20///
21/// This function will abort and raise an error if the transaction cannot be completed.
22///
23/// # Panics
24///
25/// This function may panic if one vertex in the `others` list has no associated coordinates.
26///
27/// # Example
28///
29/// For an example of usage, see the `shift` [benchmark code][BENCH].
30///
31/// [BENCH]: https://github.com/LIHPC-Computational-Geometry/honeycomb/tree/master/benches/src
32#[inline]
33pub fn move_vertex_to_average<T: CoordsFloat>(
34    t: &mut Transaction,
35    map: &CMap2<T>,
36    vid: VertexIdType,
37    others: &[VertexIdType],
38) -> StmClosureResult<()> {
39    if others.is_empty() {
40        return Ok(());
41    }
42    let mut new_val = Vertex2::default();
43    for v in others {
44        let vertex = map.read_vertex(t, *v)?.unwrap();
45        new_val.0 += vertex.0;
46        new_val.1 += vertex.1;
47    }
48    new_val.0 /= T::from(others.len()).unwrap();
49    new_val.1 /= T::from(others.len()).unwrap();
50    map.write_vertex(t, vid, new_val)?;
51    Ok(())
52}
honeycomb_kernels/remeshing/relaxation.rs

honeycomb_kernels/remeshing/
relaxation.rs
