r2d.h

Go to the documentation of this file.

/*
 *
 *      r2d.h
 *      
 *      Routines for fast, geometrically robust clipping operations
 *      and analytic area/moment computations over polygons in 2D. 
 *      
 *      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 _R2D_H_
#define _R2D_H_

#ifdef __cplusplus
extern "C" {
#endif

#include <stdint.h>

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

typedef int32_t r2d_int;
typedef int64_t r2d_long;

typedef union {
    struct {
        r2d_real x, 
                 y; 
    };
    r2d_real xy[2]; 
} r2d_rvec2;

typedef union {
    struct {
        r2d_int i, 
                j; 
    };
    r2d_int ij[2]; 
} r2d_dvec2;

typedef struct {
    r2d_rvec2 n; 
    r2d_real d; 
} r2d_plane;

typedef struct {
    r2d_int pnbrs[2]; 
    r2d_rvec2 pos; 
} r2d_vertex;

typedef struct {
#define R2D_MAX_VERTS 256 
    r2d_vertex verts[R2D_MAX_VERTS]; 
    r2d_int nverts; 
} r2d_poly;

void r2d_clip(r2d_poly* poly, r2d_plane* planes, r2d_int nplanes);

void r2d_split(r2d_poly* inpolys, r2d_int npolys, r2d_plane plane, r2d_poly* out_pos, r2d_poly* out_neg);

#define R2D_NUM_MOMENTS(order) ((order+1)*(order+2)/2)
void r2d_reduce(r2d_poly* poly, r2d_real* moments, r2d_int polyorder);

r2d_int r2d_is_good(r2d_poly* poly);

r2d_rvec2 r2d_poly_center(r2d_poly* poly);

void r2d_shift_moments(r2d_real* moments, r2d_int polyorder, r2d_rvec2 vc);

r2d_real r2d_orient(r2d_rvec2 pa, r2d_rvec2 pb, r2d_rvec2 pc);

void r2d_print(r2d_poly* poly);

void r2d_rotate(r2d_poly* poly, r2d_real theta);

void r2d_translate(r2d_poly* poly, r2d_rvec2 shift);

void r2d_scale(r2d_poly* poly, r2d_real scale);

void r2d_shear(r2d_poly* poly, r2d_real shear, r2d_int axb, r2d_int axs);

void r2d_affine(r2d_poly* poly, r2d_real mat[3][3]);

void r2d_init_box(r2d_poly* poly, r2d_rvec2* rbounds);

void r2d_init_poly(r2d_poly* poly, r2d_rvec2* vertices, r2d_int numverts);

void r2d_box_faces_from_verts(r2d_plane* faces, r2d_rvec2* rbounds);

void r2d_poly_faces_from_verts(r2d_plane* faces, r2d_rvec2* vertices, r2d_int numverts);

#ifdef __cplusplus
}
#endif

#endif // _R2D_H_