-union elem_val_un { /* holds an element value */
- char c;
- int i;
- long int l;
- float f;
- double d;
-};
-
-typedef union elem_val_un MTX_ELEM_VAL;
-
-
-union elem_ptr_un { /* holds a pointer to a 1D vector of any type */
- char *c;
- int *i;
- long int *l;
- float *f;
- double *d;
-};
-
-typedef union elem_ptr_un MTX_1D;
-typedef MTX_1D *MTX_2D;
-typedef MTX_2D *MTX_3D;
-
-union mtx_val_ptr_un { /* pointer to matrix values */
- MTX_1D m1;
- MTX_2D m2;
- MTX_3D m3;
-};
-
-typedef union mtx_val_ptr_un MTX_PTR;
-
-struct matrix_st {
- unsigned int order; /* order, or dimension, of matrix */
- unsigned int elemtype; /* element type */
- unsigned int elemsize; /* size of element in bytes */
- unsigned int nx, ny, nz; /* size of matrix in each dimension */
- MTX_PTR val; /* pointer to matrix values */
-};
-
-typedef struct matrix_st MTX;
-typedef struct matrix_st *MTXP;
-
-/* DEFINITION IDENTIFICATION
- *
- * Definitions to access a matrix element from an matrix
- *
-
-#define me1(mtx,x)\
- (mtx->elemtype == DT_FLOAT ? mtx->val.m1.f[x] :\
- (mtx->elemtype == DT_DOUBLE ? mtx->val.m1.d[x] :\
- (mtx->elemtype == DT_INT ? mtx->val.m1.i[x] :\
- (mtx->elemtype == DT_LONG ? mtx->val.m1.l[x] :\
- (mtx->elemtype == DT_CHAR ? mtx->val.m1.c[x] :\
- 0\
- )))))
-
-#define me2(mtx,x,y)\
- (mtx->elemtype == DT_FLOAT ? mtx->val.m2[x].f[y] :\
- (mtx->elemtype == DT_DOUBLE ? mtx->val.m2[x].d[y] :\
- (mtx->elemtype == DT_INT ? mtx->val.m2[x].i[y] :\
- (mtx->elemtype == DT_LONG ? mtx->val.m2[x].l[y] :\
- (mtx->elemtype == DT_CHAR ? mtx->val.m2[x].c[y] :\
- 0\
- )))))
-
-
-#define me3(mtx,x,y,z)\
- (mtx->elemtype == DT_FLOAT ? mtx->val.m3[x][y].f[z] :\
- (mtx->elemtype == DT_DOUBLE ? mtx->val.m3[x][y].d[z] :\
- (mtx->elemtype == DT_INT ? mtx->val.m3[x][y].i[z] :\
- (mtx->elemtype == DT_LONG ? mtx->val.m3[x][y].l[z] :\
- (mtx->elemtype == DT_CHAR ? mtx->val.m3[x][y].c[z] :\
- 0\
- )))))
-*/
-
-/* clip.c */
-int clip_rect(double *x1, double *y1, double *x2, double *y2, const double rect[4]);
-int clip_segment(double *x1, double *y1, double *x2, double *y2, const double u, const double v);
-int clip_sector(double *x1, double *y1, double *x2, double *y2, const double u, const double v);
-int clip_circle(double *x1, double *y1, double *x2, double *y2, const double cx, const double cy, const double radius, double t1, double t2);
-int clip_triangle(double *x1, double *y1, double *x2, double *y2, const double u, const double v, const int clip_xaxis);
-
-/* lnearest.c */
-long int lnearest(double x);
-
-/* minmax.c */
-double fmax(const double a, const double b);
-void minmax_dvector(const double array[], const int pts, double *xmin, double *xmax);
-
-/* mtx_disp.c */
-void mtx_show(const MTX *mtx);
-void mtx_prt(const MTX *mtx, FILE *fp);
-int mtx_prt_elem(const MTX *mtx, FILE *fp, unsigned int x, unsigned int y, unsigned int z);
-
-/* mtx_elem.c */
-int mtx_get_elem(const MTX *mtx, MTX_ELEM_VAL *me, const int x, const int y, const int z);
-int mtx_put_elem(MTX *mtx, const MTX_ELEM_VAL *me, unsigned int x, unsigned int y, unsigned int z);
-
-/* mtx_inp.c */
-int mtx_inp_elem(const char *prompt, MTX_ELEM_VAL *mev, const int dtype);
-
-/* mtx_main.c */
-MTX *mtx_init(const unsigned int order, const unsigned int elem_type, const unsigned int nx, const unsigned int ny, const unsigned int nz);
-MTX *mtx_clr(MTX *mtx);
-int mtx_free(MTX *mtx);
-int mtx_elem_size(const int dt);
-int mtx_check(const MTX *mtx, const char *func_name);
-
-/* norm_ang.c */
-double norm_ang(double theta);
-
-/* simpson.c */
-double simpson(const double xmin, const double xmax, const double *y, const int np);
-
-/* xform.c */
+//////////////////////////////////////////////////////////////
+// FUNTION DECLARATIONS
+//////////////////////////////////////////////////////////////
+
+// clip.cpp
+int clip_rect(double& x1, double& y1, double& x2, double& y2, const double rect[4]);
+int clip_segment(double& x1, double& y1, double& x2, double& y2, const double u, const double v);
+int clip_sector(double& x1, double& y1, double& x2, double& y2, const double u, const double v);
+int clip_circle(double& x1, double& y1, double& x2, double& y2, const double cx, const double cy, const double radius, double t1, double t2);
+int clip_triangle(double& x1, double& y1, double& x2, double& y2, const double u, const double v, const int clip_xaxis);
+
+// norm_ang.cpp
+double norm_ang (double theta);
+
+// xform.cpp