test-cube-with-hole.cc File Reference

Construct a topological cube with a through hole directly using Euler operators. More...

#include "print.h"

Go to the source code of this file.

Functions
int	main ()

---

Detailed Description

Construct a topological cube with a through hole directly using Euler operators.

Definition in file test-cube-with-hole.cc.

---

Function Documentation

int main (   )

Definition at line 30 of file test-cube-with-hole.cc. References euler::loop_t::del_single_vertex_edge_pair(), euler::loop_t::e, euler::loop_t::f, euler::solid_t::f, euler::face_t::kemr(), euler::solid_t::kfmrh(), euler::face_t::l, euler::solid_t::mef(), euler::face_t::mekr(), euler::loop_t::mev(), euler::loop_t::nxt, euler::face_t::nxt, euler::edge_t::nxt, euler::edge_t::o, and euler::edge_t::pre. { point_t p(1.,1.,1.); // geometrically collapsed to one pnt. /* * intialize to skeletal cube: * topology = 1 0d-cell + 1 2d-cells */ solid_t *s = solid_t :: mvfs(p); loop_t *l = s->f->l; // singular loop: 0 cout << s << endl; /* * add 3 consecutive dangling edges to form an open square: * topology = 4 0d-cell + 3 1d-cells + 1 2d-cell */ l->mev( l->mev( l->mev( l->f->l->e, p ), p ), p ); l->del_single_vertex_edge_pair(); // loop = 0 1 2 3 2 1 cout << l << endl; /* * close the open square by splitting the loop * topology = 2 2d-cells (faces) + 4 1d-cells (edges) + 4 0d-cells (vertices) */ s->mef(l->e, l->e->nxt->nxt->nxt); // new loop: 3 0 1 2 (top cube face), old loop l= 0 3 2 1 cout << s << endl; cout << l << endl; /* * next, add a dangling edge to each vertex of the square */ edge_t *e03 = l->e, *e10 = e03->pre, *e21 = e10->pre, *e32 = e21->pre; l->mev(e03, p); l->mev(e10, p); l->mev(e21, p); l->mev(e32, p); // now: 0 3 7 3 2 6 2 1 5 1 0 4 cout << s << endl; cout << l << endl; /* * close the 4 opening squares to form the cube */ edge_t *e40 = l->e->pre, *e73 = l->e->nxt->nxt, *e62 = e73->nxt->nxt->nxt, *e51 = e62->nxt->nxt->nxt; s->mef(e40, e73); s->mef(e73, e62); s->mef(e62, e51); s->mef(e51, e51->nxt->nxt->nxt); cout << s << endl; /* * contruct the "inner cube" */ /* * make inner top face */ face_t* top = s->f->nxt->nxt->nxt->nxt; // 3 0 1 2 face_t* btm = s->f->nxt->nxt->nxt->nxt->nxt; // 5 4 7 6 l = top->kemr( top->l->mev(top->l->e->nxt, p)->o ); // singular: 8 l->mev( l->mev( l->mev(l->e,p), p ), p ); l->del_single_vertex_edge_pair(); // 8 9 A B A 9 //now top face has 2 loops, with the inner one zero area. cout << s << endl; cout << l << endl; //cut top face into two face_t* inner_top = s->mef(l->e, l->e->nxt->nxt->nxt)->f; cout << s << endl; /* * make inner bottom face */ l = btm->kemr( btm->l->mev(btm->l->e->nxt, p)->o ); // singular: 12 l->mev( l->mev( l->mev(l->e,p), p ), p ); l->del_single_vertex_edge_pair(); // 8 9 A B A 9 (all +4) //now btm face has 2 loops, with the inner one zero area. cout << s << endl; cout << l << endl; //cut btm face into two face_t* inner_btm = s->mef(l->e->nxt->nxt->nxt, l->e)->f; cout << s << endl; cout << "inner_top = " << inner_top << endl; cout << "inner_btm = " << inner_btm << endl; s->kfmrh(inner_top, inner_btm); face_t* hole = inner_top; //now the vertical hole is annulus (or cynliner) cout << s << endl; cout << "...hole = \n" << hole << endl; edge_t *e89 = hole->l->e->nxt, *eDC = hole->l->nxt->e->pre, *e9A = e89->nxt, *eED = eDC->pre, *eAB = e9A->nxt, *eFE = eED->pre, *eB8 = eAB->nxt, *eCF = eFE->pre; // cout << e89 << endl; // cout << e9A << endl; // cout << eAB << endl; // cout << eB8 << endl; // cout << eDC << endl; // cout << eED << endl; // cout << eFE << endl; // cout << eCF << endl; hole->mekr(e89, eCF); cout << s << endl; cout << "hole = \n" << hole << endl; // exit(0); cout << "e9A=" << e9A << endl; cout << "eDC=" << eDC << endl; s->mef(e9A, eDC); cout << "hole = \n" << hole << endl; s->mef(eAB, eED); cout << eFE << endl << endl; cout << eB8 << endl; cout << eB8->nxt << endl; cout << eB8->nxt->nxt << endl; cout << eB8->nxt->nxt->nxt << endl; cout << eB8->nxt->nxt->nxt->nxt << endl; // exit(0); s->mef(eB8, eFE); cout << s << endl; return 0; }

---