diff options
Diffstat (limited to 'examples')
| -rw-r--r-- | examples/animation.rs | 122 | ||||
| -rw-r--r-- | examples/pyramid.rs | 37 | ||||
| -rw-r--r-- | examples/subdivision.rs | 130 |
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(); + }); +} |