bempp/boundary_assemblers/integrands/
double_layer.rs

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
//! Double layer integrand
use rlst::RlstScalar;

use super::{Access1D, Access2D, BoundaryIntegrand, GeometryAccess};

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

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

unsafe impl<T: RlstScalar> BoundaryIntegrand for DoubleLayerBoundaryIntegrand<T> {
    type T = T;

    fn evaluate(
        &self,
        k: &impl Access1D<T = T>,
        test_table: &impl Access2D<T = T>,
        trial_table: &impl Access2D<T = T>,
        _test_geometry: &impl GeometryAccess<T = T>,
        trial_geometry: &impl GeometryAccess<T = T>,
    ) -> T {
        unsafe {
            (k.get(1) * trial_geometry.normal(0)
                + k.get(2) * trial_geometry.normal(1)
                + k.get(3) * trial_geometry.normal(2))
                * test_table.get(0, 0)
                * trial_table.get(0, 0)
        }
    }
}

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