geo.h

Go to the documentation of this file.

/*
**      PFSGeo.h - header file for the computational geometry package
**
**      Version:        1.0
**
**      History:        Created 28-Mar-88       Luiz F. Martha
**                      Modified 7-jun-88       PWW     version 1.0
**                      Modified 24-Nov-94      LFM
**                         Changed float to real.
**      Modified 11-May-95      Kubrusly
**         Changed real to double
**      Modified 08-Feb-96      LFM
**         Modified GeoPoly structure, replancing 
**         fields a, u, v, n by a GeoNormPlane 
**         structure field.
**      Modified 17-Apr-06      Alexandre A. Del Savio
**         Modified the name of this package from "Geo"
**         to "PFSGeo".
*/

#ifndef _PFSGeo_h
#define _PFSGeo_h

/*      --------------------------------------------------
**      computational geometry data types :     
*/


typedef struct _PFSGeoPoint { double x, y, z ; }            PFSGeoPoint ;

typedef struct _PFSGeoLine { PFSGeoPoint d, e ; }         PFSGeoLine ;

typedef struct _PFSGeoPlane { PFSGeoPoint a, b, c ; }     PFSGeoPlane ;

typedef struct _PFSGeoNormPlane { PFSGeoPoint a, u, 
                                        v, n ; }    PFSGeoNormPlane ;

typedef struct _PFSGeoPoly { PFSGeoNormPlane nplane ;
                  int n_verts; PFSGeoPoint *verts ; }  PFSGeoPoly ;

typedef struct _PFSGeoPlanEqn { double a, b, c, d ; }   PFSGeoPlanEqn ;

typedef struct _PFSGeoSurfPoint { double x, y, z ;
                               double xn, yn, zn ; } PFSGeoSurfPoint ;

typedef double                                      PFSGeoPointPar ;

typedef struct _PFSGeoSurfPar { double u, v ; }         PFSGeoSurfPar ;

typedef enum _PFSGeoStatus { PFSGeoInvalid, PFSGeoValid }    PFSGeoStatus ;

typedef enum _PFSGeoPredicate { PFSGeoFalse, PFSGeoTrue }    PFSGeoPredicate ;

typedef enum _PFSGeoPtTest { PFSGeoAbove,
                          PFSGeoOn,
                          PFSGeoBelow         }     PFSGeoPtTest ;

/*      --------------------------------------------------
**      function prototypes :   
*/

#ifdef __cplusplus
extern "C" {
#endif

void PFSGeoSetTolerance( double tol );
double PFSGeoGetTolerance( void );

PFSGeoPredicate PFSGeoIntscPolyLine ( PFSGeoSurfPar *line, PFSGeoSurfPar *spar ) ;
PFSGeoPredicate PFSGeoIsInsidePoly  ( PFSGeoPoly *poly, PFSGeoSurfPar *spar ) ;
PFSGeoPredicate    PFSGeoIsOnPolyBdry  ( PFSGeoPoly *poly, PFSGeoPoint *pt );

void      PFSGeoGenLine  ( PFSGeoPoint *pt0, PFSGeoPoint *pt1, PFSGeoLine *line ) ;
void      PFSGeoGenPlane ( PFSGeoPoint *pt0, PFSGeoPoint *pt1, PFSGeoPoint *pt2,
                        PFSGeoPlane *plane ) ;
void      PFSGeoGenPlanEqn( PFSGeoPoint *pt0, PFSGeoPoint *pt1, PFSGeoPoint *pt2,
                         PFSGeoPlanEqn *plan_eqn ) ;

PFSGeoStatus PFSGeoGenPoly  ( int *n_verts, PFSGeoPoint *verts, PFSGeoNormPlane *plane ) ;

void      PFSGeoEvalLine    ( PFSGeoPointPar *par, PFSGeoLine *line, PFSGeoPoint *pt ) ;
void      PFSGeoEvalPlane   ( PFSGeoSurfPar  *par, PFSGeoPlane *pln, PFSGeoPoint *pt ) ;

double        PFSGeoDotProd    ( PFSGeoPoint *in0, PFSGeoPoint *in1 ) ;
void         PFSGeoCrossProd  ( PFSGeoPoint *in0, PFSGeoPoint *in1, PFSGeoPoint *out ) ;
double        PFSGeoTripleProd ( PFSGeoPoint *in0, PFSGeoPoint *in1, PFSGeoPoint *in2 ) ;
double       PFSGeoVecLen     ( PFSGeoPoint *vec ) ;
void         PFSGeoVecNormalize ( PFSGeoPoint *vec, PFSGeoPoint *n ) ;
double       PFSGeoDist    ( PFSGeoPoint *in0, PFSGeoPoint *in1 ) ;
PFSGeoPredicate PFSGeoNullVec    ( PFSGeoPoint *vec, double *tolerance ) ;

void    PFSGeoDiffVec( PFSGeoPoint *in0, PFSGeoPoint *in1, PFSGeoPoint *out );
void    PFSGeoSumVec( PFSGeoPoint *in0, PFSGeoPoint *in1, PFSGeoPoint *out );
void    PFSGeoProdScalVec( PFSGeoPoint *in, double *csi, PFSGeoPoint *out );
void    PFSGeoLinCombVec( PFSGeoPoint *in0, double *csi0, PFSGeoPoint *in1, double *csi1,
                       PFSGeoPoint *out );

void      PFSGeoRotateX ( double *angle, PFSGeoPoint *in, PFSGeoPoint *out ) ;
void      PFSGeoRotateY ( double *angle, PFSGeoPoint *in, PFSGeoPoint *out ) ;
void      PFSGeoRotateZ ( double *angle, PFSGeoPoint *in, PFSGeoPoint *out ) ;
void      PFSGeoRotateAxis ( PFSGeoPoint *axis, double *angle, 
                          PFSGeoPoint *in, PFSGeoPoint *out ) ;

PFSGeoStatus PFSGeoIntscLinePlane ( PFSGeoLine *line, PFSGeoPlane *plane, 
                              PFSGeoPoint *pt, PFSGeoPointPar *ptpar,
                              PFSGeoSurfPar *spar ) ;
PFSGeoStatus PFSGeoIntscLinePoly ( PFSGeoLine *line, PFSGeoPoly *poly, 
                              PFSGeoPoint *pt, PFSGeoPointPar *ptpar,
                              PFSGeoSurfPar *spar ) ;

PFSGeoStatus PFSGeoPolyNormal  ( int *n_verts, PFSGeoPoint *verts,  PFSGeoPoint *norm ) ;

PFSGeoStatus  PFSGeoLstSqrPlanEqn( int *num_pts, PFSGeoPoint *pts, 
                             PFSGeoPlanEqn *plan_eqn ) ;

void      PFSGeoClstPtLine  ( PFSGeoLine *line, PFSGeoPoint *pt, PFSGeoPoint *clst,
                           PFSGeoPointPar *par ) ;
void      PFSGeoClstPtNormPlane  ( PFSGeoNormPlane *plane, PFSGeoPoint *pt, 
                                PFSGeoPoint *clst, PFSGeoSurfPar *spar ) ;
PFSGeoStatus PFSGeoClstPtPoly  ( PFSGeoPoly *poly, PFSGeoPoint *pt, 
                           PFSGeoPoint *clst, PFSGeoSurfPar *spar ) ;

PFSGeoPredicate PFSGeoIsOnPlane( PFSGeoPoint  *pt, PFSGeoNormPlane  *pla ) ;

PFSGeoStatus PFSGeoInvert2x2 ( double a[2][2], double b[2][2] ) ;
PFSGeoStatus PFSGeoInvert3x3 ( double a[3][3], double b[3][3] ) ;
PFSGeoStatus PFSGeodInvert3x3 ( double a[3][3], double b[3][3] ) ;

void    PFSGeoCentroid( int *num_pts,  PFSGeoPoint *pts,  PFSGeoPoint *cntrd ) ;

PFSGeoStatus PFSGeoPerpLine ( PFSGeoLine *in0, PFSGeoLine *in1, PFSGeoLine *out,
                        PFSGeoPointPar *par0, PFSGeoPointPar *par1,
                        PFSGeoPointPar *parout ) ;

void PFSGeoPtNormToParPlan( PFSGeoPoint *pt, PFSGeoPoint *normal, PFSGeoNormPlane *plane ) ;

void PFSGeoPlanToNormPlan( PFSGeoPlane *plane, PFSGeoNormPlane *normplane );
PFSGeoStatus  PFSGeo3PointsToNormPlane( PFSGeoPoint *in0, PFSGeoPoint *in1, 
                                  PFSGeoPoint *in2, PFSGeoNormPlane *plane );

void PFSGeoPlanEqnToParPlan( PFSGeoPlanEqn *plan_eqn, PFSGeoNormPlane *plane ) ;

PFSGeoPtTest PFSGeoSideofPlan( PFSGeoPlanEqn *plan_eqn, PFSGeoPoint *pt ) ;

void PFSGeoSurfInterp( int *num_pts, double *map_func,PFSGeoSurfPar *pts_val, 
                        PFSGeoSurfPar *interp_val ) ;

void PFSGeoSurfJacob( int *num_pts, double *map_deriv_s, double *map_deriv_t, 
                        PFSGeoSurfPar *pts_val, double jac[2][2] ) ;

double PFSGeoAngVecAboutLine( PFSGeoLine *l, PFSGeoPoint *v0, PFSGeoPoint *v1 ) ;
double PFSGeoAngPtAboutLine ( PFSGeoLine *l, PFSGeoPoint *p0, PFSGeoPoint *p1 ) ;
double PFSGeoAngPtAboutPt( PFSGeoPoint *pivot, PFSGeoPoint *p0, PFSGeoPoint *p1 ) ;

PFSGeoStatus PFSGeoProjVecPlane( PFSGeoPoint *vec, PFSGeoPlane *plane, PFSGeoPoint *proj ) ;
double PFSGeoDistPtNormPlane( PFSGeoNormPlane *plane, PFSGeoPoint *pt );

#ifdef __cplusplus
}
#endif

#endif

---