bempp 0.2.0

Boundary element method library.
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51

//! Double layer integrand
use crate::assembly::common::RlstArray;
use crate::traits::{CellGeometry, PotentialIntegrand};
use rlst::{RlstScalar, UnsafeRandomAccessByRef};

/// Integrand for a double layer potential operator
pub struct DoubleLayerPotentialIntegrand<T: RlstScalar> {
    _t: std::marker::PhantomData<T>,
}

impl<T: RlstScalar> DoubleLayerPotentialIntegrand<T> {
    /// Create new
    pub fn new() -> Self {
        Self {
            _t: std::marker::PhantomData,
        }
    }
}

unsafe impl<T: RlstScalar> PotentialIntegrand for DoubleLayerPotentialIntegrand<T> {
    type T = T;

    fn evaluate(
        &self,
        table: &RlstArray<T, 4>,
        point_index: usize,
        eval_index: usize,
        basis_index: usize,
        k: &RlstArray<T, 3>,
        geometry: &impl CellGeometry<T = T::Real>,
    ) -> T {
        unsafe {
            -(*k.get_unchecked([1, point_index, eval_index])
                * num::cast::<T::Real, T>(*geometry.normals().get_unchecked([0, point_index]))
                    .unwrap()
                + *k.get_unchecked([2, point_index, eval_index])
                    * num::cast::<T::Real, T>(*geometry.normals().get_unchecked([1, point_index]))
                        .unwrap()
                + *k.get_unchecked([3, point_index, eval_index])
                    * num::cast::<T::Real, T>(*geometry.normals().get_unchecked([2, point_index]))
                        .unwrap())
                * *table.get_unchecked([0, point_index, basis_index, 0])
        }
    }
}

impl<T: RlstScalar> Default for DoubleLayerPotentialIntegrand<T> {
    fn default() -> Self {
        Self::new()
    }
}