Tetrahedralize

tetrahedralize.cpp

bool test_tetgen() {
	/// grid management
	Grid g;
	ug::SubsetHandler sh(g);
        g.attach_to_vertices(ug::aPosition);
        ug::Grid::VertexAttachmentAccessor<ug::APosition3> aaPos(g, ug::aPosition);
        
        /// unit cube
        ug::RegularVertex* temp1 = *g.create<ug::RegularVertex>();
        ug::RegularVertex* temp2 = *g.create<ug::RegularVertex>();
        ug::RegularVertex* temp3 = *g.create<ug::RegularVertex>();
        ug::RegularVertex* temp4 = *g.create<ug::RegularVertex>();
        ug::RegularVertex* temp5 = *g.create<ug::RegularVertex>();
        ug::RegularVertex* temp6 = *g.create<ug::RegularVertex>();
        ug::RegularVertex* temp7 = *g.create<ug::RegularVertex>();
        ug::RegularVertex* temp8 = *g.create<ug::RegularVertex>();

        /// bottom verts
        aaPos[temp1] = ug::vector3(0, 0, 0);
        aaPos[temp2] = ug::vector3(0, 1, 0);
        aaPos[temp3] = ug::vector3(1, 0, 0);
        aaPos[temp5] = ug::vector3(1, 1, 0);

	/// top verts
        aaPos[temp4] = ug::vector3(1, 1, 1);
        aaPos[temp6] = ug::vector3(0, 0, 1);
        aaPos[temp7] = ug::vector3(0, 1, 1);
        aaPos[temp8] = ug::vector3(1, 0, 1);

	/// assign to subset
        sh.assign_subset(temp1, 1);
        sh.assign_subset(temp2, 1);
        sh.assign_subset(temp3, 1);
        sh.assign_subset(temp4, 1);
        sh.assign_subset(temp5, 1);
        sh.assign_subset(temp6, 1);
        sh.assign_subset(temp7, 1);
        sh.assign_subset(temp8, 1);
        
        /// tetrahedralize
        return Tetrahedralize(g, sh, 0, false, false, ug::aPosition, 10);
}