Method draw_plane_nzz

src/surface_geometry.rs:137–163

Draws a plane that has a normal vector with a non-zero z (nzz) component The plane may be perpendicular to z if n = (0,0,1) # Input `p` -- (len=3) point on plane `n` -- (len=3) normal vector `xmin` and `xmax` -- limits along x `ymin` and `ymax` -- limits along y `nx` -- number of divisions along x (must be ≥ 2) `ny` -- number of divisions along y (must be ≥ 2) # Output `x`, `y`, `z` -- the co

(
        &mut self,
        p: &[f64],
        n: &[f64],
        xmin: f64,
        xmax: f64,
        ymin: f64,
        ymax: f64,
        nx: usize,
        ny: usize,
    )

    /// See also integration test in the **tests** directory.
    ///
    pub fn draw_plane_nzz(
        &mut self,
        p: &[f64],
        n: &[f64],
        xmin: f64,
        xmax: f64,
        ymin: f64,
        ymax: f64,
        nx: usize,
        ny: usize,
    ) -> Result<(Vec<Vec<f64>>, Vec<Vec<f64>>, Vec<Vec<f64>>), StrError> {
        if p.len() != 3 || n.len() != 3 {
            return Err("p.len() and n.len() must be equal to 3");
        }
        if f64::abs(n[2]) < 1e-10 {
            return Err("the z-component of the normal vector cannot be zero");
        }
        if nx < 2 || ny < 2 {
            return Err("nx and ny must be ≥ 2");
        }
        let d = -n[0] * p[0] - n[1] * p[1] - n[2] * p[2];
        let (x, y, z) = generate3d(xmin, xmax, ymin, ymax, nx + 1, ny + 1, |x, y| {
            (-d - n[0] * x - n[1] * y) / n[2]
        });
        self.draw(&x, &y, &z);
        Ok((x, y, z))
    }

    /// Draws a hemisphere
    ///