rNd.h

Go to the documentation of this file.

/*
 *
 *      rNd.h
 *      
 *      Routines for fast, geometrically robust clipping operations
 *      and analytic volume/moment computations over 
 *      polytopes of arbitrary dimensionality. 
 *      
 *      Devon Powell
 *      31 August 2015
 *
 *      This program was prepared by Los Alamos National Security, LLC at Los Alamos National
 *      Laboratory (LANL) under contract No. DE-AC52-06NA25396 with the U.S. Department of Energy (DOE). 
 *      All rights in the program are reserved by the DOE and Los Alamos National Security, LLC.  
 *      Permission is granted to the public to copy and use this software without charge, provided that 
 *      this Notice and any statement of authorship are reproduced on all copies.  Neither the U.S. 
 *      Government nor LANS makes any warranty, express or implied, or assumes any liability 
 *      or responsibility for the use of this software.
 *    
 */

#ifndef _RND_H_
#define _RND_H_

#ifdef __cplusplus
extern "C" {
#endif

#include <stdint.h>

/*
 * \brief Dimensionality of the geometric calculations.
 */
#define RND_DIM 4

#ifdef SINGLE_PRECISION
typedef float rNd_real;
#else 
typedef double rNd_real;
#endif

typedef int32_t rNd_int;
typedef int64_t rNd_long;

typedef struct {
    rNd_real xyz[RND_DIM]; 
} rNd_rvec;

typedef struct {
    rNd_int ijk[RND_DIM]; 
} rNd_dvec;

typedef struct {
    rNd_rvec n; 
    rNd_real d; 
} rNd_plane;

typedef struct {
    rNd_int pnbrs[RND_DIM]; 
    rNd_int finds[RND_DIM][RND_DIM]; 
    rNd_rvec pos; 
} rNd_vertex;

typedef struct {
#define RND_MAX_VERTS 1024 
    rNd_vertex verts[RND_MAX_VERTS]; 
    rNd_int nverts, nfaces; 
} rNd_poly;

void rNd_clip(rNd_poly* poly, rNd_plane* planes, rNd_int nplanes);

//#define RND_NUM_MOMENTS(order) ((order+1)*(order+2)*(order+3)/6)
void rNd_reduce(rNd_poly* poly, rNd_real* moments, rNd_int polyorder);

rNd_int rNd_is_good(rNd_poly* poly);

rNd_real rNd_orient(rNd_rvec verts[RND_DIM+1]);

void rNd_print(rNd_poly* poly);

void rNd_rotate(rNd_poly* poly, rNd_real theta, rNd_int ax1, rNd_int ax2);

void rNd_translate(rNd_poly* poly, rNd_rvec shift);

void rNd_scale(rNd_poly* poly, rNd_rvec scale);

void rNd_shear(rNd_poly* poly, rNd_real shear, rNd_int axb, rNd_int axs);

void rNd_affine(rNd_poly* poly, rNd_real mat[RND_DIM+1][RND_DIM+1]);

void rNd_init_simplex(rNd_poly* poly, rNd_rvec verts[RND_DIM+1]);

void rNd_init_box(rNd_poly* poly, rNd_rvec rbounds[2]);

//void rNd_simplex_faces_from_verts(rNd_plane* faces, rNd_rvec3* verts);

//void rNd_box_faces_from_verts(rNd_plane* faces, rNd_rvec3* rbounds);

#ifdef __cplusplus
}
#endif

#endif // _RND_H_