3 files changed, 289 insertions, 0 deletions

diff --git a/examples/animation.rs b/examples/animation.rs
new file mode 100644
index 0000000..f6aca39
--- /dev/null
+++ b/examples/animation.rs
@@ -0,0 +1,122 @@
+use cairo::{Context, SvgSurface};
+use dcel::Dcel;
+use enumset::EnumSet;
+use std::borrow::Cow;
+
+fn main() {
+    let show = |name, dcel: &Dcel<(&'static str, [i64; 2])>| {
+        let width = 300.0;
+        let height = width;
+
+        let surf =
+			//cairo::ImageSurface::create(cairo::Format::ARgb32, width as _, height as _).unwrap();
+			SvgSurface::new(width, height, Some(name)).unwrap();
+        let ctx = Context::new(surf).unwrap();
+        ctx.set_line_width(1.0);
+
+        dcel::write_img(
+            &dcel,
+            &ctx,
+            EnumSet::all(),
+            |v| {
+                [
+                    v.1[0] as f64 * width / 3.0 + width / 2.0,
+                    v.1[1] as f64 * height / 3.0 + height / 2.0,
+                ]
+            },
+            |v| Cow::Borrowed(v.0),
+            12.0,
+        )
+        .unwrap();
+
+        /*ctx.target()
+            .write_to_png(&mut std::fs::File::create(name).unwrap())
+            .unwrap();
+        */
+        /*write!(
+            &mut std::fs::File::create(name).unwrap(),
+            "{}",
+            DisplayFn(|f: &mut fmt::Formatter<'_>| dcel_write_dot(
+                dcel,
+                |v| v.1.map(|x| x as _),
+                |v, f| write!(f, "{}", v.0),
+                f,
+                DcelDotOptions {
+                    prev: false,
+                    next: true,
+                    twin: true,
+                }
+            ))
+        )
+        .unwrap();*/
+    };
+
+    Dcel::new(|mut dcel| {
+        let body = dcel.new_body();
+        // Mevvlfs(a, [w, n], l, f, s)
+
+        //let op = dcel.mevvlfs(*body, [("W", [-4, 0]), ("N", [0, 4])]);
+        let op = dcel
+            .mevvlfs(*body, [("W", [-1, 0]), ("N", [0, 1])])
+            .unwrap();
+        let op2 = dcel.mev(*op.loop_, *op.vertices[1], ("E", [1, 0])).unwrap();
+        let op3 = dcel.mev(*op.loop_, *op2.vertex, ("S", [0, -1])).unwrap();
+
+        dcel.melf([*op3.vertex, *op.vertices[0]], *op.loop_)
+            .unwrap();
+        dcel.melf([*op.vertices[0], *op2.vertex], *op.loop_)
+            .unwrap();
+
+        show("cairo.svg", &dcel);
+    });
+
+    /*
+    dbg!(ctx.clip_extents().unwrap());
+
+    let text = "meow meow meow mew mew mrrrp";
+    let x = 250.0;
+    let y = 250.0;
+
+    let ext = ctx.text_extents(text).unwrap();
+
+    ctx.move_to(
+        x - ext.width() / 2.0,
+        y - ext.y_bearing() - ext.height() / 2.0,
+    );
+    ctx.show_text(text).unwrap();
+
+    ctx.set_line_width(1.0);
+    {
+        let mat = ctx.matrix();
+
+        ctx.translate(x, y);
+        ctx.scale(
+            (ext.width() + ctx.line_width()) / 2.0f64.sqrt(),
+            (ext.height() + ctx.line_width()) / 2.0f64.sqrt(),
+        );
+        ctx.translate(-x, -y);
+        ctx.new_path();
+        ctx.arc(x, y, 1.0, 0.0, 2.0 * std::f64::consts::PI);
+        ctx.set_matrix(mat);
+    }
+
+    ctx.stroke().unwrap();*/
+
+    /*ctx.rectangle(
+        x - (ext.width() + ctx.line_width()) / 2.0,
+        y - (ext.height() + ctx.line_width()) / 2.0,
+        ext.width() + ctx.line_width(),
+        ext.height() + ctx.line_width(),
+    );
+    ctx.stroke().unwrap();*/
+
+    /*ctx.rectangle(
+        x - ext.width() / 2.0,
+        y - ext.height() / 2.0,
+        ext.width(),
+        ext.height(),
+    );*/
+
+    //cairo_move_to (cr, i + 0.5 - te.x_bearing - te.width / 2,
+    //            0.5 - te.y_bearing - te.height / 2);
+}
diff --git a/examples/pyramid.rs b/examples/pyramid.rs
new file mode 100644
index 0000000..9614292
--- /dev/null
+++ b/examples/pyramid.rs
@@ -0,0 +1,37 @@
+use dcel::{Dcel, Kevvlfs, ObjExport};
+
+fn main() {
+    Dcel::new(|mut dcel| {
+        let body = dcel.new_body();
+
+        let [d1, d2, d3, d4] = [
+            ((8.0f64 / 9.0).sqrt(), 0.0, -1.0 / 3.0),
+            (-(2.0f64 / 9.0).sqrt(), (2.0f64 / 3.0).sqrt(), -1.0 / 3.0),
+            (-(2.0f64 / 9.0).sqrt(), -(2.0f64 / 3.0).sqrt(), -1.0 / 3.0),
+            (0.0, 0.0, 1.0),
+        ];
+
+        let Kevvlfs {
+            vertices: [v1, v2],
+            loop_: l1,
+            ..
+        } = dcel.mevvlfs(*body, [d1, d2]).unwrap();
+        let v3 = dcel.mev(*l1, *v2, d3).unwrap().vertex;
+        let l2 = dcel.melf([*v1, *v3], *l1).unwrap().loop_;
+        let v4 = dcel.mev(*l2, *v1, d4).unwrap().vertex;
+        let mut l3 = dcel.melf([*v4, *v2], *l2).unwrap().loop_;
+        if v3.find_outgoing(*l3, &dcel).is_none() {
+            l3 = l2;
+        }
+        let _l4 = dcel.melf([*v4, *v3], *l3).unwrap().loop_;
+
+        ObjExport::export(
+            &mut std::fs::File::create("pyramid.obj").unwrap(),
+            &dcel,
+            |&(x, y, z)| (x, y, z, None),
+            |_, _| None,
+            |_, _| None,
+        )
+        .unwrap();
+    });
+}
diff --git a/examples/subdivision.rs b/examples/subdivision.rs
new file mode 100644
index 0000000..ba80e49
--- /dev/null
+++ b/examples/subdivision.rs
@@ -0,0 +1,130 @@
+use dcel::{ptr_t, Dcel, ObjExport, ObjImport, Ptr, Shell};
+use std::array::from_fn;
+use std::collections::HashMap;
+
+fn avg(mut iter: impl Iterator<Item = [f32; 3]>) -> Option<[f32; 3]> {
+    let mut count = 0usize;
+    let mut accum = [0.0; 3];
+
+    for x in iter {
+        count += 1;
+        accum = from_fn(|i| accum[i] + x[i])
+    }
+
+    (count > 0).then(|| accum.map(|x| x / (count as f32)))
+}
+
+fn catmull_clark_subdivision<'brand, 'arena>(
+    dcel: &mut Dcel<'brand, 'arena, [f32; 3]>,
+    shell: ptr_t!(Shell<'brand, 'arena, [f32; 3]>),
+) {
+    let face_points = shell
+        .iter_faces(dcel)
+        .map(|f| {
+            (
+                f.id(),
+                avg(f.outer_loop().iter_half_edges().map(|x| *x.origin().data())).unwrap(),
+            )
+        })
+        .collect::<HashMap<_, _>>();
+
+    let vert_points = shell
+        .iter_vertices(dcel)
+        .map(|vert| {
+            let p = vert.data();
+            let mut n = 0.0;
+            let mut f = [0.0; 3];
+            let mut r = [0.0; 3];
+
+            for h in vert.iter_outgoing() {
+                let fp = face_points[&h.loop_().face().id()];
+                let op = *h.origin().data();
+                let tp = *h.target().data();
+
+                f = from_fn(|i| f[i] + fp[i]);
+                r = from_fn(|i| r[i] + op[i] + tp[i]);
+
+                n += 1.0;
+            }
+
+            (
+                vert.id(),
+                (
+                    vert.item,
+                    from_fn(|i| (f[i] / n + r[i] / n + (n - 3.0) * p[i]) / n),
+                ),
+            )
+        })
+        .collect::<HashMap<_, _>>();
+
+    for edge in shell.iter_edges(dcel).map(|x| x.item).collect::<Vec<_>>() {
+        let data = avg(edge
+            .lens(dcel)
+            .half_edges()
+            .into_iter()
+            .flat_map(|h| [face_points[&h.loop_().face().id()], *h.origin().data()]))
+        .unwrap();
+
+        dcel.mve(edge, data).unwrap();
+    }
+
+    for &(vert, data) in vert_points.values() {
+        *vert.mut_data(dcel) = data;
+    }
+
+    for face in shell.iter_faces(dcel).map(|x| x.item).collect::<Vec<_>>() {
+        let loop_ = face.outer_loop(dcel);
+
+        let mut verts = loop_
+            .iter_half_edges(dcel)
+            .skip_while(|h| vert_points.contains_key(&h.origin().id()))
+            .map(|h| h.origin().item)
+            .step_by(2)
+            .collect::<Vec<_>>()
+            .into_iter();
+
+        let inner = *dcel
+            .mev(loop_, verts.next().unwrap(), face_points[&face.id(dcel)])
+            .unwrap()
+            .vertex;
+
+        verts.fold([loop_; 2], |loops, vert| {
+            let l = loops
+                .into_iter()
+                .find(|&l| vert.find_outgoing(l, dcel).is_some())
+                .unwrap();
+            [l, *dcel.melf([vert, inner], l).unwrap().loop_]
+        });
+    }
+}
+
+fn main() {
+    Dcel::new(|mut dcel| {
+        let body = ObjImport::import(
+            &mut dcel,
+            &obj::raw::object::parse_obj(&mut std::io::BufReader::new(
+                std::fs::File::open("bunny.obj").unwrap(),
+            ))
+            .unwrap(),
+            |(x, y, z, _)| [x, y, z],
+        )
+        .unwrap();
+
+        let shell = body.iter_shells(&dcel).next().unwrap().item;
+
+        for i in 0..5 {
+            println!("{i}...");
+            catmull_clark_subdivision(&mut dcel, shell);
+        }
+
+        println!("exporting...");
+        ObjExport::export(
+            &mut std::io::BufWriter::new(std::fs::File::create("bunny2.obj").unwrap()),
+            &dcel,
+            |&[x, y, z]| (x as _, y as _, z as _, None),
+            |_, _| None,
+            |_, _| None,
+        )
+        .unwrap();
+    });
+}