2014-05-10 17:16:39 +00:00
|
|
|
/*
|
|
|
|
* Copyright 2011-2014 hasufell
|
|
|
|
*
|
|
|
|
* This file is part of a hasufell project.
|
|
|
|
*
|
|
|
|
* This program is free software: you can redistribute it and/or modify
|
|
|
|
* it under the terms of the GNU General Public License as published by
|
|
|
|
* the Free Software Foundation version 2 of the License only.
|
|
|
|
*
|
|
|
|
* This program is distributed in the hope that it will be useful,
|
|
|
|
* but WITHOUT ANY WARRANTY; without even the implied warranty of
|
|
|
|
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
|
|
|
|
* GNU General Public License for more details.
|
|
|
|
*
|
|
|
|
* You should have received a copy of the GNU General Public License
|
|
|
|
* along with this program. If not, see <http://www.gnu.org/licenses/>.
|
|
|
|
*/
|
|
|
|
|
2014-05-10 17:33:54 +00:00
|
|
|
/**
|
|
|
|
* @file half_edge.c
|
|
|
|
* This file provides operations on half-edge data structures
|
|
|
|
* which are defined in half_edge.h, as well as assembling
|
|
|
|
* such a struct.
|
|
|
|
* @brief operations on half-edge data structs
|
|
|
|
*/
|
|
|
|
|
|
|
|
#include "err.h"
|
|
|
|
#include "filereader.h"
|
2014-05-10 17:16:39 +00:00
|
|
|
#include "half_edge.h"
|
2014-05-11 12:02:34 +00:00
|
|
|
#include "vector.h"
|
2014-05-10 17:16:39 +00:00
|
|
|
|
2014-05-11 12:02:34 +00:00
|
|
|
#include <stdbool.h>
|
2014-05-10 17:33:54 +00:00
|
|
|
#include <stdint.h>
|
|
|
|
#include <stdio.h>
|
2014-05-10 17:16:39 +00:00
|
|
|
#include <stdlib.h>
|
2014-05-10 17:33:54 +00:00
|
|
|
#include <string.h>
|
2014-05-10 17:16:39 +00:00
|
|
|
|
|
|
|
|
2014-05-13 11:47:13 +00:00
|
|
|
/**
|
|
|
|
* Fault intolerant macro. Will abort the program if the called
|
|
|
|
* function failed.
|
|
|
|
*/
|
|
|
|
#define GET_ALL_EMANATING_EDGES(...) \
|
|
|
|
{ \
|
|
|
|
if (!get_all_emanating_edges(__VA_ARGS__)) { \
|
|
|
|
fprintf(stderr, "Failure in get_all_emanating_edges()!\n"); \
|
|
|
|
abort(); \
|
|
|
|
} \
|
|
|
|
}
|
|
|
|
|
|
|
|
|
2014-05-12 18:20:03 +00:00
|
|
|
/*
|
|
|
|
* static declarations
|
|
|
|
*/
|
2014-05-13 11:47:13 +00:00
|
|
|
static bool get_all_emanating_edges(HE_vert const * const vert,
|
|
|
|
HE_edge ***edge_array_out,
|
2014-05-12 18:20:03 +00:00
|
|
|
uint32_t *ec_out);
|
|
|
|
|
|
|
|
|
2014-05-10 17:16:39 +00:00
|
|
|
/**
|
2014-05-12 18:20:03 +00:00
|
|
|
* Get all edges that emanate from vertice and return a pointer
|
|
|
|
* to that array with the size of ec_out.
|
2014-05-10 17:16:39 +00:00
|
|
|
*
|
2014-05-12 19:00:31 +00:00
|
|
|
* @param vert the vertice to get the emanating edges of
|
2014-05-13 11:47:13 +00:00
|
|
|
* @param edge_array_out address of the 2d edge array to save
|
|
|
|
* the result in [out]
|
2014-05-12 19:00:31 +00:00
|
|
|
* @param ec_out the edge counter is saved here [out]
|
2014-05-12 18:20:03 +00:00
|
|
|
* @return pointer to an array of half-edges, size ec_out
|
2014-05-10 17:16:39 +00:00
|
|
|
*/
|
2014-05-13 11:47:13 +00:00
|
|
|
static bool get_all_emanating_edges(HE_vert const * const vert,
|
|
|
|
HE_edge ***edge_array_out,
|
2014-05-12 18:20:03 +00:00
|
|
|
uint32_t *ec_out)
|
2014-05-10 17:16:39 +00:00
|
|
|
{
|
2014-05-12 18:20:03 +00:00
|
|
|
uint32_t ec = 0, /* edge count */
|
|
|
|
rc = 0; /* realloc count */
|
|
|
|
uint32_t const approx_ec = 20; /* allocation chunk */
|
2014-05-13 11:47:13 +00:00
|
|
|
HE_edge **edge_array;
|
2014-05-12 18:20:03 +00:00
|
|
|
HE_edge **tmp_ptr;
|
|
|
|
|
2014-05-13 12:05:02 +00:00
|
|
|
if (!edge_array_out || !vert || !ec_out)
|
2014-05-13 11:47:13 +00:00
|
|
|
return false;
|
|
|
|
|
|
|
|
edge_array = malloc(sizeof(HE_edge*) * approx_ec);
|
|
|
|
CHECK_PTR_VAL(edge_array);
|
2014-05-12 18:20:03 +00:00
|
|
|
|
|
|
|
HE_edge *edge = vert->edge;
|
|
|
|
|
|
|
|
/* build an array of emanating edges */
|
|
|
|
do {
|
|
|
|
edge_array[ec] = edge;
|
|
|
|
|
|
|
|
edge = edge->pair->next;
|
|
|
|
ec++;
|
|
|
|
|
|
|
|
/* allocate more chunks */
|
|
|
|
if (ec >= approx_ec) {
|
|
|
|
tmp_ptr = realloc(edge_array, sizeof(HE_edge*)
|
|
|
|
* approx_ec * (rc + 2));
|
|
|
|
CHECK_PTR_VAL(tmp_ptr);
|
|
|
|
edge_array = tmp_ptr;
|
|
|
|
rc++;
|
|
|
|
}
|
|
|
|
|
|
|
|
} while (edge != vert->edge);
|
|
|
|
|
2014-05-13 11:47:13 +00:00
|
|
|
/* set out-pointers */
|
|
|
|
*edge_array_out = edge_array;
|
|
|
|
*ec_out = ec; /* this is the real size, not the x[ec] value */
|
2014-05-12 18:20:03 +00:00
|
|
|
|
2014-05-13 11:47:13 +00:00
|
|
|
return true;
|
2014-05-12 18:20:03 +00:00
|
|
|
}
|
|
|
|
|
2014-05-12 21:02:10 +00:00
|
|
|
/**
|
|
|
|
* Calculate the normal of a face that corresponds
|
|
|
|
* to edge.
|
|
|
|
*
|
|
|
|
* @param edge to align the normalization
|
|
|
|
* @param vec the vector to store the result in [out]
|
|
|
|
* @return true/false for success/failure
|
|
|
|
*/
|
|
|
|
bool face_normal(HE_edge const * const edge,
|
|
|
|
vector *vec)
|
|
|
|
{
|
|
|
|
vector he_vec1,
|
|
|
|
he_vec2,
|
|
|
|
he_base;
|
|
|
|
|
2014-05-13 11:47:13 +00:00
|
|
|
if (!edge || !vec)
|
2014-05-12 21:02:10 +00:00
|
|
|
return false;
|
|
|
|
|
2014-05-13 11:47:13 +00:00
|
|
|
COPY_VECTOR(edge->next->vert->vec, &he_base);
|
2014-05-12 21:02:10 +00:00
|
|
|
|
2014-05-13 11:47:13 +00:00
|
|
|
/* calculate vectors between the vertices */
|
|
|
|
SUB_VECTORS(edge->next->next->vert->vec, &he_base, &he_vec1);
|
|
|
|
SUB_VECTORS(edge->vert->vec, &he_base, &he_vec2);
|
2014-05-12 21:02:10 +00:00
|
|
|
|
2014-05-13 11:47:13 +00:00
|
|
|
VECTOR_PRODUCT(&he_vec1, &he_vec2, vec);
|
|
|
|
|
|
|
|
NORMALIZE_VECTOR(vec, vec);
|
2014-05-12 21:02:10 +00:00
|
|
|
|
|
|
|
return true;
|
|
|
|
}
|
|
|
|
|
2014-05-12 18:20:03 +00:00
|
|
|
/**
|
|
|
|
* Calculate the approximated normal of a vertex.
|
|
|
|
*
|
|
|
|
* @param vert the vertex
|
|
|
|
* @param vec the vector to store the result in [out]
|
|
|
|
* @return true/false for success/failure
|
|
|
|
*/
|
|
|
|
bool vec_normal(HE_vert const * const vert, vector *vec)
|
|
|
|
{
|
2014-05-13 11:47:13 +00:00
|
|
|
HE_edge **edge_array = NULL;
|
2014-05-12 21:02:10 +00:00
|
|
|
uint32_t ec;
|
2014-05-12 18:20:03 +00:00
|
|
|
vector he_base;
|
|
|
|
|
|
|
|
if (!vert || !vec)
|
|
|
|
return false;
|
|
|
|
|
2014-05-13 11:47:13 +00:00
|
|
|
GET_ALL_EMANATING_EDGES(vert, &edge_array, &ec);
|
|
|
|
COPY_VECTOR(edge_array[0]->vert->vec, &he_base);
|
|
|
|
SET_NULL_VECTOR(vec); /* set to null for later summation */
|
2014-05-12 18:20:03 +00:00
|
|
|
|
2014-05-12 21:02:10 +00:00
|
|
|
/* iterate over all edges, get the normalized
|
|
|
|
* face vector and add those up */
|
2014-05-12 18:20:03 +00:00
|
|
|
for (uint32_t i = 0; i < ec; i++) {
|
2014-05-12 21:02:10 +00:00
|
|
|
vector new_vec;
|
2014-05-12 18:20:03 +00:00
|
|
|
|
2014-05-13 11:47:13 +00:00
|
|
|
FACE_NORMAL(edge_array[i], &new_vec);
|
|
|
|
ADD_VECTORS(vec, &new_vec, vec);
|
2014-05-12 18:20:03 +00:00
|
|
|
}
|
|
|
|
|
|
|
|
/* normalize the result */
|
2014-05-13 11:47:13 +00:00
|
|
|
NORMALIZE_VECTOR(vec, vec);
|
2014-05-12 18:20:03 +00:00
|
|
|
|
2014-05-12 21:02:10 +00:00
|
|
|
free(edge_array);
|
2014-05-12 18:20:03 +00:00
|
|
|
return true;
|
2014-05-10 17:16:39 +00:00
|
|
|
}
|
2014-05-10 17:33:54 +00:00
|
|
|
|
2014-05-10 19:05:08 +00:00
|
|
|
/**
|
|
|
|
* Find the center of an object and store the coordinates
|
|
|
|
* in a HE_vert struct.
|
|
|
|
*
|
|
|
|
* @param obj the object we want to find the center of
|
2014-05-11 12:02:34 +00:00
|
|
|
* @param vec the vector to store the result in [out]
|
|
|
|
* @return true/false for success/failure
|
2014-05-10 19:05:08 +00:00
|
|
|
*/
|
2014-05-11 12:02:34 +00:00
|
|
|
bool find_center(HE_obj const * const obj, vector *vec)
|
2014-05-10 19:05:08 +00:00
|
|
|
{
|
|
|
|
float x = 0,
|
|
|
|
y = 0,
|
|
|
|
z = 0;
|
|
|
|
uint32_t i;
|
2014-05-10 20:48:50 +00:00
|
|
|
|
2014-05-11 12:02:34 +00:00
|
|
|
if (!obj || !vec)
|
|
|
|
return false;
|
2014-05-10 20:48:50 +00:00
|
|
|
|
2014-05-10 19:05:08 +00:00
|
|
|
|
|
|
|
for (i = 0; i < obj->vc; i++) {
|
2014-05-12 17:48:46 +00:00
|
|
|
x += obj->vertices[i].vec->x;
|
|
|
|
y += obj->vertices[i].vec->y;
|
|
|
|
z += obj->vertices[i].vec->z;
|
2014-05-10 19:05:08 +00:00
|
|
|
}
|
|
|
|
|
2014-05-11 12:02:34 +00:00
|
|
|
vec->x = x / i;
|
|
|
|
vec->y = y / i;
|
|
|
|
vec->z = z / i;
|
2014-05-10 19:05:08 +00:00
|
|
|
|
2014-05-11 12:02:34 +00:00
|
|
|
return true;
|
2014-05-10 19:05:08 +00:00
|
|
|
}
|
|
|
|
|
|
|
|
/**
|
2014-05-10 19:06:36 +00:00
|
|
|
* Calculates the factor that can be used to scale down the object
|
2014-05-10 19:05:08 +00:00
|
|
|
* to the size of 1.
|
|
|
|
*
|
|
|
|
* @param obj the object we want to scale
|
2014-05-10 20:58:51 +00:00
|
|
|
* @return the corresponding scale factor, -1 on error
|
2014-05-10 19:05:08 +00:00
|
|
|
*/
|
|
|
|
float get_normalized_scale_factor(HE_obj const * const obj)
|
|
|
|
{
|
2014-05-10 20:58:51 +00:00
|
|
|
float max;
|
|
|
|
float min;
|
|
|
|
uint32_t i;
|
|
|
|
|
|
|
|
if (!obj)
|
|
|
|
return -1;
|
|
|
|
|
2014-05-12 17:48:46 +00:00
|
|
|
max = obj->vertices[0].vec->x +
|
|
|
|
obj->vertices[0].vec->y + obj->vertices[0].vec->z;
|
|
|
|
min = obj->vertices[0].vec->x +
|
|
|
|
obj->vertices[0].vec->y + obj->vertices[0].vec->z;
|
2014-05-10 19:05:08 +00:00
|
|
|
|
|
|
|
for (i = 0; i < obj->vc; i++) {
|
2014-05-12 17:48:46 +00:00
|
|
|
if ((obj->vertices[i].vec->x +
|
|
|
|
obj->vertices[i].vec->y +
|
|
|
|
obj->vertices[i].vec->z) > max)
|
|
|
|
max = obj->vertices[i].vec->x +
|
|
|
|
obj->vertices[i].vec->y +
|
|
|
|
obj->vertices[i].vec->z;
|
|
|
|
else if ((obj->vertices[i].vec->x +
|
|
|
|
obj->vertices[i].vec->y +
|
|
|
|
obj->vertices[i].vec->z) < min)
|
|
|
|
min = obj->vertices[i].vec->x +
|
|
|
|
obj->vertices[i].vec->y +
|
|
|
|
obj->vertices[i].vec->z;
|
2014-05-10 19:05:08 +00:00
|
|
|
}
|
|
|
|
|
|
|
|
return 1 / (max - min);
|
|
|
|
}
|
|
|
|
|
2014-05-12 18:20:03 +00:00
|
|
|
/**
|
|
|
|
* Scales down the object to the size of 1. The parameter
|
|
|
|
* is modified!
|
|
|
|
*
|
|
|
|
* @param obj the object we want to scale [mod]
|
|
|
|
*/
|
2014-05-13 12:05:02 +00:00
|
|
|
bool normalize_object(HE_obj *obj)
|
2014-05-12 18:20:03 +00:00
|
|
|
{
|
|
|
|
float scale_factor;
|
2014-05-13 12:05:02 +00:00
|
|
|
|
|
|
|
if (!obj)
|
|
|
|
return false;
|
|
|
|
|
2014-05-12 18:20:03 +00:00
|
|
|
scale_factor = get_normalized_scale_factor(obj);
|
|
|
|
|
|
|
|
for (uint32_t i = 0; i < obj->vc; i++) {
|
|
|
|
obj->vertices[i].vec->x = obj->vertices[i].vec->x * scale_factor;
|
|
|
|
obj->vertices[i].vec->y = obj->vertices[i].vec->y * scale_factor;
|
|
|
|
obj->vertices[i].vec->z = obj->vertices[i].vec->z * scale_factor;
|
|
|
|
}
|
2014-05-13 12:05:02 +00:00
|
|
|
|
|
|
|
return true;
|
2014-05-12 18:20:03 +00:00
|
|
|
}
|
|
|
|
|
2014-05-10 17:33:54 +00:00
|
|
|
/**
|
|
|
|
* Parse an .obj string and return a HE_obj
|
|
|
|
* that represents the whole object.
|
|
|
|
*
|
|
|
|
* @param obj_string the whole string from the .obj file
|
|
|
|
* @return the HE_face array that represents the object
|
|
|
|
*/
|
|
|
|
HE_obj *parse_obj(char const * const obj_string)
|
|
|
|
{
|
|
|
|
uint32_t vc = 0, /* vertices count */
|
|
|
|
fc = 0, /* face count */
|
|
|
|
ec = 0; /* edge count */
|
2014-05-10 20:37:03 +00:00
|
|
|
char *string,
|
2014-05-10 17:33:54 +00:00
|
|
|
*str_ptr_space = NULL, /* for strtok */
|
|
|
|
*str_ptr_newline = NULL, /* for strtok */
|
|
|
|
*str_tmp_ptr = NULL; /* for strtok */
|
|
|
|
HE_vert *vertices = NULL;
|
|
|
|
HE_edge *edges = NULL;
|
|
|
|
HE_face *faces = NULL;
|
|
|
|
HE_obj *obj = NULL;
|
|
|
|
FACE face_v = NULL;
|
|
|
|
|
2014-05-10 20:37:03 +00:00
|
|
|
if (!obj_string || !*obj_string)
|
|
|
|
return NULL;
|
|
|
|
|
|
|
|
string = malloc(sizeof(char) * strlen(obj_string) + 1);
|
2014-05-10 17:33:54 +00:00
|
|
|
strcpy(string, obj_string);
|
|
|
|
|
|
|
|
str_tmp_ptr = strtok_r(string, "\n", &str_ptr_newline);
|
|
|
|
while (str_tmp_ptr && *str_tmp_ptr) {
|
|
|
|
|
|
|
|
str_tmp_ptr = strtok_r(str_tmp_ptr, " ", &str_ptr_space);
|
|
|
|
|
2014-05-13 18:09:16 +00:00
|
|
|
/* parse vertices and fill them */
|
2014-05-10 17:33:54 +00:00
|
|
|
if (!strcmp(str_tmp_ptr, "v")) {
|
|
|
|
char *myfloat = NULL;
|
|
|
|
HE_vert *tmp_ptr;
|
2014-05-12 17:48:46 +00:00
|
|
|
vector *tmp_vec = malloc(sizeof(vector));
|
|
|
|
CHECK_PTR_VAL(tmp_vec);
|
2014-05-10 17:33:54 +00:00
|
|
|
|
|
|
|
tmp_ptr = (HE_vert*) realloc(vertices,
|
|
|
|
sizeof(HE_vert) * (vc + 1));
|
|
|
|
CHECK_PTR_VAL(tmp_ptr);
|
|
|
|
vertices = tmp_ptr;
|
|
|
|
|
|
|
|
/* fill x */
|
|
|
|
myfloat = strtok_r(NULL, " ", &str_ptr_space);
|
|
|
|
CHECK_PTR_VAL(myfloat);
|
2014-05-12 17:48:46 +00:00
|
|
|
tmp_vec->x = atof(myfloat);
|
2014-05-10 17:33:54 +00:00
|
|
|
|
|
|
|
/* fill y */
|
|
|
|
myfloat = strtok_r(NULL, " ", &str_ptr_space);
|
|
|
|
CHECK_PTR_VAL(myfloat);
|
2014-05-12 17:48:46 +00:00
|
|
|
tmp_vec->y = atof(myfloat);
|
2014-05-10 17:33:54 +00:00
|
|
|
|
|
|
|
/* fill z */
|
|
|
|
myfloat = strtok_r(NULL, " ", &str_ptr_space);
|
|
|
|
CHECK_PTR_VAL(myfloat);
|
2014-05-12 17:48:46 +00:00
|
|
|
tmp_vec->z = atof(myfloat);
|
|
|
|
|
|
|
|
vertices[vc].vec = tmp_vec;
|
2014-05-10 17:33:54 +00:00
|
|
|
|
2014-05-13 18:16:33 +00:00
|
|
|
/* set unused/unknown values to NULL */
|
2014-05-10 17:33:54 +00:00
|
|
|
vertices[vc].edge = NULL;
|
2014-05-13 18:09:16 +00:00
|
|
|
vertices[vc].edge_array = NULL;
|
|
|
|
vertices[vc].eac = 0;
|
2014-05-10 17:33:54 +00:00
|
|
|
|
|
|
|
vc++;
|
|
|
|
|
|
|
|
/* exceeds 3 dimensions, malformed vertice */
|
|
|
|
if (strtok_r(NULL, " ", &str_ptr_space))
|
2014-05-13 12:05:02 +00:00
|
|
|
ABORT("Failure in parse_obj(),\n"
|
|
|
|
"malformed vertice, exceeds 2 dimensions!\n");
|
2014-05-10 17:33:54 +00:00
|
|
|
|
2014-05-13 18:09:16 +00:00
|
|
|
/* parse plain faces and fill them (not HE_face!) */
|
2014-05-10 17:33:54 +00:00
|
|
|
} else if (!strcmp(str_tmp_ptr, "f")) {
|
|
|
|
char *myint = NULL;
|
|
|
|
uint8_t i = 0;
|
|
|
|
FACE tmp_ptr = NULL;
|
|
|
|
|
|
|
|
tmp_ptr = (FACE) realloc(face_v, sizeof(FACE*) * (fc + 1));
|
|
|
|
CHECK_PTR_VAL(tmp_ptr);
|
|
|
|
face_v = tmp_ptr;
|
|
|
|
face_v[fc] = NULL;
|
|
|
|
while ((myint = strtok_r(NULL, " ", &str_ptr_space))) {
|
|
|
|
uint32_t *tmp_ptr = NULL;
|
|
|
|
|
|
|
|
i++;
|
|
|
|
ec++;
|
|
|
|
|
|
|
|
tmp_ptr = (uint32_t*) realloc(face_v[fc],
|
|
|
|
sizeof(FACE**) * (i + 1));
|
|
|
|
CHECK_PTR_VAL(tmp_ptr);
|
|
|
|
tmp_ptr[i - 1] = (uint32_t) atoi(myint);
|
|
|
|
tmp_ptr[i] = 0; /* so we can iterate over it */
|
|
|
|
face_v[fc] = tmp_ptr;
|
|
|
|
}
|
|
|
|
fc++;
|
|
|
|
}
|
|
|
|
|
|
|
|
str_tmp_ptr = strtok_r(NULL, "\n", &str_ptr_newline);
|
|
|
|
}
|
|
|
|
|
|
|
|
faces = (HE_face*) malloc(sizeof(HE_face) * fc);
|
|
|
|
CHECK_PTR_VAL(faces);
|
|
|
|
edges = (HE_edge*) malloc(sizeof(HE_edge) * ec);
|
|
|
|
CHECK_PTR_VAL(edges);
|
|
|
|
|
|
|
|
ec = 0;
|
|
|
|
/* create HE_edges and real HE_faces */
|
2014-05-13 18:09:16 +00:00
|
|
|
for (uint32_t i = 0; i < fc; i++) { /* for all faces */
|
|
|
|
uint32_t j = 0,
|
|
|
|
fv_id; /* reference of the face vertex */
|
2014-05-10 17:33:54 +00:00
|
|
|
|
|
|
|
/* for all vertices of the face */
|
2014-05-13 18:09:16 +00:00
|
|
|
while ((fv_id = face_v[i][j])) {
|
|
|
|
uint32_t fv_arr_id = fv_id - 1; /* fv_id starts at 1 */
|
|
|
|
|
|
|
|
edges[ec].vert = &(vertices[fv_arr_id]);
|
2014-05-10 17:33:54 +00:00
|
|
|
edges[ec].face = &(faces[j]);
|
|
|
|
edges[ec].pair = NULL; /* preliminary */
|
2014-05-13 18:09:16 +00:00
|
|
|
vertices[fv_arr_id].edge = &(edges[ec]); /* last one wins */
|
|
|
|
|
|
|
|
/* Skip j == 0 here, so we don't underrun the arrays,
|
|
|
|
* since we always look one edge back. The first edge
|
|
|
|
* element is taken care of below as well. */
|
|
|
|
if (j > 0 ) {
|
|
|
|
HE_edge **tmp_ptr = NULL;
|
|
|
|
uint32_t *eac = &(edges[ec].vert->eac);
|
|
|
|
|
|
|
|
/* connect previous edge to current edge */
|
|
|
|
edges[ec - 1].next = &(edges[ec]);
|
|
|
|
|
2014-05-13 18:16:33 +00:00
|
|
|
/* Acceleration struct:
|
|
|
|
* add previous edge to edge_array of current vertice */
|
2014-05-13 18:09:16 +00:00
|
|
|
tmp_ptr = realloc(edges[ec].vert->edge_array,
|
|
|
|
sizeof(HE_edge*) * (*eac + 1));
|
|
|
|
CHECK_PTR_VAL(tmp_ptr);
|
|
|
|
edges[ec].vert->edge_array = tmp_ptr;
|
|
|
|
edges[ec].vert->edge_array[*eac] = &(edges[ec - 1]);
|
|
|
|
(*eac)++;
|
|
|
|
|
|
|
|
if (!face_v[i][j + 1]) { /* no vertice left */
|
|
|
|
uint32_t *eac;
|
|
|
|
/* connect last edge to first edge */
|
|
|
|
edges[ec].next = &(edges[ec - j]);
|
|
|
|
eac = &(edges[ec].next->vert->eac);
|
|
|
|
|
2014-05-13 18:16:33 +00:00
|
|
|
/* Acceleration struct:
|
|
|
|
* add last edge to edge_array element of first vertice */
|
2014-05-13 18:09:16 +00:00
|
|
|
tmp_ptr = realloc(edges[ec].next->vert->edge_array,
|
|
|
|
sizeof(HE_edge*) * (*eac + 1));
|
|
|
|
CHECK_PTR_VAL(tmp_ptr);
|
|
|
|
edges[ec].next->vert->edge_array = tmp_ptr;
|
|
|
|
edges[ec].next->vert->edge_array[*eac] = &(edges[ec]);
|
|
|
|
(*eac)++;
|
|
|
|
}
|
|
|
|
}
|
2014-05-10 17:33:54 +00:00
|
|
|
|
|
|
|
ec++;
|
|
|
|
j++;
|
|
|
|
}
|
|
|
|
|
|
|
|
faces[i].edge = &(edges[ec - 1]); /* "last" edge */
|
|
|
|
}
|
|
|
|
|
|
|
|
/* find pairs */
|
2014-05-13 18:26:36 +00:00
|
|
|
for (uint32_t i = 0; i < ec; i++) { /* for all edges */
|
2014-05-13 18:09:16 +00:00
|
|
|
uint32_t eac = edges[i].vert->eac;
|
|
|
|
|
2014-05-13 18:26:36 +00:00
|
|
|
for (uint32_t j = 0; j < eac; j++) { /* for all potential pairs */
|
2014-05-13 18:09:16 +00:00
|
|
|
if (edges[i].vert->edge_array[j] &&
|
|
|
|
(edges[i].next->vert ==
|
|
|
|
edges[i].vert->edge_array[j]->vert)) {
|
|
|
|
edges[i].pair = edges[i].vert->edge_array[j];
|
|
|
|
edges[i].vert->edge_array[j] = NULL;
|
|
|
|
break;
|
|
|
|
}
|
|
|
|
}
|
2014-05-10 17:33:54 +00:00
|
|
|
}
|
|
|
|
|
|
|
|
obj = (HE_obj*) malloc(sizeof(HE_obj));
|
|
|
|
CHECK_PTR_VAL(obj);
|
|
|
|
|
|
|
|
obj->edges = edges;
|
2014-05-13 12:05:10 +00:00
|
|
|
obj->vertices = vertices;
|
2014-05-10 17:33:54 +00:00
|
|
|
obj->faces = faces;
|
2014-05-13 12:05:10 +00:00
|
|
|
obj->ec = ec;
|
|
|
|
obj->vc = vc;
|
2014-05-10 17:33:54 +00:00
|
|
|
obj->fc = fc;
|
|
|
|
|
|
|
|
free(string);
|
|
|
|
for (uint32_t i = 0; i < fc; i++)
|
|
|
|
free(face_v[i]);
|
|
|
|
free(face_v);
|
|
|
|
|
|
|
|
return obj;
|
|
|
|
}
|
2014-05-12 21:02:10 +00:00
|
|
|
|
|
|
|
/**
|
|
|
|
* Free the inner structures of an object.
|
|
|
|
*
|
|
|
|
* @param obj the object to free
|
|
|
|
*/
|
|
|
|
void delete_object(HE_obj *obj)
|
|
|
|
{
|
2014-05-13 12:05:02 +00:00
|
|
|
if (!obj)
|
|
|
|
return;
|
|
|
|
|
2014-05-12 21:02:10 +00:00
|
|
|
for (uint32_t i = 0; i < obj->vc; i++)
|
|
|
|
free(obj->vertices[i].vec);
|
|
|
|
free(obj->edges);
|
|
|
|
free(obj->vertices);
|
|
|
|
free(obj->faces);
|
|
|
|
}
|