Struct Fun

pub struct Fun<F> { /* private fields */ }

Options for sampling graphs of functions ℝ → ℝ. See Sampling::fun.

Implementations

impl<F> Fun<F>

where F: FnMut(f64) -> f64,

pub fn n(self, n: usize) -> Self

Set the maximum number of evaluations of the function to build the sampling. Panic if n < 2.

Examples found in repository

examples/speed.rs (line 38)

fn main() {
    let mut b = false;
    for n in 10 .. 5_000 {
        let s = Sampling::fun(|x| x.sin(), 0., 10.)
            .n(n).build();
        b &= s.is_empty();
    }
    println!("{b}")
}

examples/nice.rs (line 13)

fn main() -> Result<(), Box<dyn Error>> {
    let f = |t: f64| [t.cos(), (2. * t).sin()];
    let s = Sampling::param(f, 0., 2. * PI).build();
    s.write(&mut BufWriter::new(File::create("/tmp/nice1.dat")?))?;

    let f = |x: f64| (- x.powi(2)).exp();
    let s = Sampling::fun(f, -2.5, 2.5).n(53).build();
    s.write(&mut BufWriter::new(File::create("/tmp/nice2.dat")?))?;
    Ok(())
}

examples/sin_inv_x.rs (line 29)

fn main() -> R {
    let mut fh = File::create("/tmp/sin_inv_x.gp")?;
    write!(fh, "set terminal pngcairo\n\
                set grid\n")?;
    let mut d = 0;
    let mut save = |s: &Sampling, n, title| -> R {
        d += 1;
        let fname = format!("/tmp/sin_inv_x{}.dat", d);
        s.write(&mut File::create(&fname)?)?;
        write!(fh, "set output \"sin_inv_x{}.png\"\n\
                    plot '{}' with l lt 1 lw 2 title \"{} ({} pts)\"\n",
               d, &fname, title, n)?;
        write!(fh, "set output \"sin_inv_x{}_p.png\"\n\
                    plot '{}' with l lt 5 lw 2 title \"{}\", \
                    '{}' with p lt 3 pt 5 ps 0.2 title \"points ({})\"\n",
               d, &fname, title, &fname, n)?;
        Ok(())
    };

    let f = |x: f64| x * (1. / x).sin();
    let s = Sampling::fun(f, -0.4, 0.4).n(227).build();
    save(&s, 227, "x sin(1/x)")?;
    let s = Sampling::fun(f, -0.4, 0.4).n(389).build();
    save(&s, 389, "x sin(1/x)")?;

    let s = Sampling::fun(|x: f64| (1. / x).sin(), -0.4, 0.4).n(391).build();
    save(&s, 391, "sin(1/x)")?;

    Ok(())
}

pub fn viewport(self, vp: BoundingBox) -> Self

Set the zone of interest for the sampling. Segments that end up outside this box will not be refined.

pub fn init<'a, I>(self, ts: I) -> Self

where I: IntoIterator<Item = &'a f64>,

Add initial values of t such that f(t) (see [Sampling :: fun]) must be included into the sampling in addition to the n evaluations. Only the values between a and b are taken into account (other values are ignored).

pub fn init_pt<'a, I>(self, pts: I) -> Self

where I: IntoIterator<Item = &'a (f64, f64)>,

Add initial points (t, f(t)) to include into the sampling in addition to the n evaluations. This allows you to use previous evaluations of f. Only the couples with first coordinate t between a and b (see [Sampling :: fun]) are considered (other values are ignored).

impl<F> Fun<F>

where F: FnMut(f64) -> f64,

pub fn build(self) -> Sampling

Return the sampling.

Examples found in repository

examples/speed.rs (line 38)

fn main() {
    let mut b = false;
    for n in 10 .. 5_000 {
        let s = Sampling::fun(|x| x.sin(), 0., 10.)
            .n(n).build();
        b &= s.is_empty();
    }
    println!("{b}")
}

examples/nice.rs (line 13)

fn main() -> Result<(), Box<dyn Error>> {
    let f = |t: f64| [t.cos(), (2. * t).sin()];
    let s = Sampling::param(f, 0., 2. * PI).build();
    s.write(&mut BufWriter::new(File::create("/tmp/nice1.dat")?))?;

    let f = |x: f64| (- x.powi(2)).exp();
    let s = Sampling::fun(f, -2.5, 2.5).n(53).build();
    s.write(&mut BufWriter::new(File::create("/tmp/nice2.dat")?))?;
    Ok(())
}

examples/sin_inv_x.rs (line 29)

fn main() -> R {
    let mut fh = File::create("/tmp/sin_inv_x.gp")?;
    write!(fh, "set terminal pngcairo\n\
                set grid\n")?;
    let mut d = 0;
    let mut save = |s: &Sampling, n, title| -> R {
        d += 1;
        let fname = format!("/tmp/sin_inv_x{}.dat", d);
        s.write(&mut File::create(&fname)?)?;
        write!(fh, "set output \"sin_inv_x{}.png\"\n\
                    plot '{}' with l lt 1 lw 2 title \"{} ({} pts)\"\n",
               d, &fname, title, n)?;
        write!(fh, "set output \"sin_inv_x{}_p.png\"\n\
                    plot '{}' with l lt 5 lw 2 title \"{}\", \
                    '{}' with p lt 3 pt 5 ps 0.2 title \"points ({})\"\n",
               d, &fname, title, &fname, n)?;
        Ok(())
    };

    let f = |x: f64| x * (1. / x).sin();
    let s = Sampling::fun(f, -0.4, 0.4).n(227).build();
    save(&s, 227, "x sin(1/x)")?;
    let s = Sampling::fun(f, -0.4, 0.4).n(389).build();
    save(&s, 389, "x sin(1/x)")?;

    let s = Sampling::fun(|x: f64| (1. / x).sin(), -0.4, 0.4).n(391).build();
    save(&s, 391, "sin(1/x)")?;

    Ok(())
}