r3d.h

Go to the documentation of this file.

/*
 *
 *      r3d.h
 *
 *      Routines for fast, geometrically robust clipping operations and analytic
 *      volume/moment computations over polyhedra in 3D.
 *
 *      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 _R3D_H_
#define _R3D_H_

#include <stdint.h>
#include <stdlib.h>

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

typedef int32_t r3d_int;
typedef int64_t r3d_long;

typedef union {
  struct
  {
    r3d_real x, 
        y,      
        z;      
  };
  r3d_real xyz[3]; 
} r3d_rvec3;

typedef union {
  struct
  {
    r3d_int i, 
        j,     
        k;     
  };
  r3d_int ijk[3]; 
} r3d_dvec3;

typedef struct
{
  r3d_rvec3 n; 
  r3d_real d;  
} r3d_plane;

typedef struct
{
  r3d_int pnbrs[3]; 
  r3d_rvec3 pos;    
} r3d_vertex;

typedef struct
{
#define R3D_MAX_VERTS 512
  r3d_vertex verts[R3D_MAX_VERTS]; 
  r3d_int nverts;                  
} r3d_poly;

void r3d_clip(r3d_poly *poly, r3d_plane *planes, r3d_int nplanes);

void r3d_split(r3d_poly *inpolys, r3d_int npolys, r3d_plane plane,
               r3d_poly *out_pos, r3d_poly *out_neg);

#define R3D_NUM_MOMENTS(order) ((order + 1) * (order + 2) * (order + 3) / 6)
void r3d_reduce(r3d_poly *poly, r3d_real *moments, r3d_int polyorder);

r3d_int r3d_is_good(r3d_poly *poly);

r3d_rvec3 r3d_poly_center(r3d_poly* poly);

void r3d_shift_moments(r3d_real* moments, r3d_int polyorder, r3d_rvec3 vc);

r3d_real r3d_orient(r3d_rvec3 *verts);

void r3d_print(r3d_poly *poly);

void r3d_rotate(r3d_poly *poly, r3d_real theta, r3d_int axis);

void r3d_translate(r3d_poly *poly, r3d_rvec3 shift);

void r3d_scale(r3d_poly *poly, r3d_real scale);

void r3d_shear(r3d_poly *poly, r3d_real shear, r3d_int axb, r3d_int axs);

void r3d_affine(r3d_poly *poly, r3d_real mat[4][4]);

void r3d_init_tet(r3d_poly *poly, r3d_rvec3 *verts);

void r3d_init_box(r3d_poly *poly, r3d_rvec3 *rbounds);

void r3d_init_poly(r3d_poly *poly, r3d_rvec3 *vertices, r3d_int numverts,
                   r3d_int **faceinds, r3d_int *numvertsperface,
                   r3d_int numfaces);

void r3d_tet_faces_from_verts(r3d_plane *faces, r3d_rvec3 *verts);

void r3d_box_faces_from_verts(r3d_plane *faces, r3d_rvec3 *rbounds);

void r3d_poly_faces_from_verts(r3d_plane *faces, r3d_rvec3 *vertices,
                               r3d_int numverts, r3d_int **faceinds,
                               r3d_int *numvertsperface, r3d_int numfaces);

typedef struct r3d_brep
{
  r3d_rvec3 *vertices;
  r3d_int numvertices;
  r3d_int **faceinds;
  r3d_int *numvertsperface;
  r3d_int numfaces;
} r3d_brep;

void r3d_init_brep(r3d_poly *poly, r3d_brep **brep, r3d_int *numcomponents);


void r3d_print_brep(r3d_brep **brep, r3d_int numcomponents);

void r3d_free_brep(r3d_brep **brep, r3d_int numcomponents);

#endif // _R3D_H_