Uni/CGI
1
0
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity

Test volume

Browse source
This commit is contained in:
hodasemi 2019-01-07 21:41:52 +01:00
parent e470aef3d5
commit fd494ffddc
1 changed files with 75 additions and 39 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

76
exercise3/include/triangles.h
View file

@ -2,11 +2,14 @@
#include "util/OpenMeshUtils.h"
struct SimpleFace {
struct SimpleFace
{
std::vector<OpenMesh::HalfedgeHandle> he_handles;
bool is_triangulated() {
if (he_handles.size() != 3) {
bool is_triangulated()
{
if (he_handles.size() != 3)
{
return false;
}
@ -14,13 +17,17 @@ struct SimpleFace {
}
};
struct Triangle {
Triangle(std::vector<OpenMesh::Vec3f> vpoints) {
if (vpoints.size() != 3) {
struct Triangle
{
Triangle(std::vector<OpenMesh::Vec3f> vpoints)
{
if (vpoints.size() != 3)
{
abort();
}
for (int i = 0; i < 3; i++) {
for (int i = 0; i < 3; i++)
{
points[i] = vpoints[i];
}
}
@ -28,30 +35,37 @@ struct Triangle {
OpenMesh::Vec3f points[3];
};
class Triangles {
class Triangles
{
public:
Triangles(HEMesh private_mesh) {
Triangles(HEMesh private_mesh)
{
std::vector<SimpleFace> simple_faces;
if (!gather_faces(private_mesh, simple_faces)) {
if (!gather_faces(private_mesh, simple_faces))
{
simple_faces.clear();
private_mesh.triangulate();
if (!gather_faces(private_mesh, simple_faces)) {
if (!gather_faces(private_mesh, simple_faces))
{
abort();
}
}
for (auto simple_face: simple_faces) {
for (auto simple_face : simple_faces)
{
std::vector<OpenMesh::Vec3f> points = face_points(simple_face, private_mesh);
triangles.emplace_back(points);
}
}
float surface_area() {
float surface_area()
{
float surface_area = 0.0f;
for (Triangle& triangle : triangles) {
for (Triangle &triangle : triangles)
{
OpenMesh::Vec3f p1 = triangle.points[0];
OpenMesh::Vec3f p2 = triangle.points[1];
OpenMesh::Vec3f p3 = triangle.points[2];
@ -67,20 +81,40 @@ public:
return surface_area;
}
float volume()
{
float volume = 0.0f;
for (Triangle &triangle : triangles)
{
OpenMesh::Vec3f p1 = triangle.points[0];
OpenMesh::Vec3f p2 = triangle.points[1];
OpenMesh::Vec3f p3 = triangle.points[2];
volume = OpenMesh::dot(p1, OpenMesh::cross(p2, p3));
}
return volume;
}
private:
bool gather_faces(const HEMesh& m, std::vector<SimpleFace>& simple_faces) {
bool gather_faces(const HEMesh &m, std::vector<SimpleFace> &simple_faces)
{
auto faces = m.all_faces();
for (auto face_handle : faces) {
for (auto face_handle : faces)
{
HEMesh::ConstFaceHalfedgeIter fh_it = m.cfh_iter(face_handle);
SimpleFace simple_face;
for (; fh_it.is_valid(); ++fh_it) {
for (; fh_it.is_valid(); ++fh_it)
{
simple_face.he_handles.emplace_back(*fh_it);
}
if (!simple_face.is_triangulated()) {
if (!simple_face.is_triangulated())
{
return false;
}
@ -90,10 +124,12 @@ private:
return true;
}
std::vector<OpenMesh::Vec3f> face_points(SimpleFace& simple_face, const HEMesh& m) {
std::vector<OpenMesh::Vec3f> face_points(SimpleFace &simple_face, const HEMesh &m)
{
std::vector<OpenMesh::Vec3f> points;
for (auto he : simple_face.he_handles) {
for (auto he : simple_face.he_handles)
{
OpenMesh::VertexHandle vertex_handle = m.from_vertex_handle(he);
points.emplace_back(m.point(vertex_handle));