geom.h

Go to the documentation of this file.

 
#ifndef __GEOM_H_
#define __GEOM_H_
 
#include "params.h"
 
#ifdef SSE_INTRIN
#define vox_dot_add(a,b,res) _mm_store_ps ((res), _mm_load_ps (a) + _mm_load_ps (b))
#else
#define vox_dot_add(a,b,res) do { \
    res[0] = a[0] + b[0]; \
    res[1] = a[1] + b[1]; \
    res[2] = a[2] + b[2]; \
    }                     \
    while (0);
#endif
 
#ifdef SSE_INTRIN
#define vox_dot_sub(a,b,res) _mm_store_ps ((res), _mm_load_ps (a) - _mm_load_ps (b))
#else
#define vox_dot_sub(a,b,res) do { \
    res[0] = a[0] - b[0]; \
    res[1] = a[1] - b[1]; \
    res[2] = a[2] - b[2]; \
    }                     \
    while (0);
#endif
 
#ifdef SSE_INTRIN
#define vox_dot_scmul(d,sc,res) _mm_store_ps ((res), _mm_load_ps (d) * _mm_set_ps1 (sc))
#else
#define vox_dot_scmul(d,sc,res) do { \
    res[0] = sc * d[0]; \
    res[1] = sc * d[1]; \
    res[2] = sc * d[2]; \
    }                   \
    while (0);
#endif
 
VOX_EXPORT float vox_abs_metric (const vox_dot dot1, const vox_dot dot2);
 
VOX_EXPORT float vox_sqr_metric (const vox_dot dot1, const vox_dot dot2);
 
VOX_EXPORT float vox_sqr_norm (const vox_dot dot);
 
#ifdef VOXTREES_SOURCE
 
VOX_EXPORT int get_subspace_idx (const vox_dot center, const vox_dot dot); // KLUDGE: needed for tests
 
/*
 * Generally this works like get_subspace_idx () but corrects subspace
 * index if center[i] == dot[i]. In this case, direction[i] is used to
 * choose the index at i-th coordinate.
 */
int get_corrected_subspace_idx (const vox_dot center, const vox_dot dot, const vox_dot direction);
 
 
int hit_box (const struct vox_box *box, const vox_dot origin, const vox_dot dir, vox_dot res);
 
int hit_plane_within_box (const vox_dot origin, const vox_dot dir, const vox_dot planedot,
                          int planenum, vox_dot res, const struct vox_box *box);
 
int box_ball_interp (const struct vox_box *box, const vox_dot center, float radius);
int dense_set_p (const struct vox_box *box, size_t n);
int voxel_in_box (const struct vox_box *box, const vox_dot dot);
void closest_vertex (const struct vox_box *box, const vox_dot dot, vox_dot res);
int divide_box (const struct vox_box *box, const vox_dot center, struct vox_box *res, int idx);
void get_dimensions (const struct vox_box *box, size_t dim[]);
int stripep (const struct vox_box *box, int *which);
#endif
 
#endif