/Users/jack/Code/basso_dev/inc/geometry.h

#ifndef _GEOMETRY_FUNCTIONS_BASSO_H_
#define _GEOMETRY_FUNCTIONS_BASSO_H_

// std includes

// Gmm++ includes

// Basso includes
#include "basso.h"
#include "Point.h"

namespace basso {

        Numeric area_triangle( const Point &a, const Point &b, const Point &c )
        {
                Numeric x1=a.x(),  y1=a.y(),  z1=a.z();
                Numeric x2=b.x(),  y2=b.y(),  z2=b.z();
                Numeric x3=c.x(),  y3=c.y(),  z3=c.z();
                return 0.5*(-x2*y1 + x3*y1 + x1*y2 - x3*y2 - x1*y3 + x2*y3 + x2*z1 - x3*z1 - 
                                y2*z1 + y3*z1 - x1*z2 + x3*z2 + y1*z2 - y3*z2 + x1*z3 - x2*z3 - 
                                y1*z3 + y2*z3);
        }

        Numeric dist( const Point &p1, const Point &p2 ) 
        { 
                return sqrt( pow(p1.x()-p2.x(),2) + pow(p1.y()-p2.y(),2) + pow(p1.z()-p2.z(),2) );
        }

        class FiniteLine 
        {
        public:
                FiniteLine( const Point &a, const Point &b ) { p0=a; p1=b; }
                
        protected:
                Point p0, p1;
                
        };
        
} // end namespace

#endif

---