geometry.h

Go to the documentation of this file.

/* *******************************************************************
 * Rocstar Simulation Suite                                          *
 * Copyright@2015, Illinois Rocstar LLC. All rights reserved.        *
 *                                                                   *
 * Illinois Rocstar LLC                                              *
 * Champaign, IL                                                     *
 * www.illinoisrocstar.com                                           *
 * sales@illinoisrocstar.com                                         *
 *                                                                   *
 * License: See LICENSE file in top level of distribution package or *
 * http://opensource.org/licenses/NCSA                               *
 *********************************************************************/
/* *******************************************************************
 * THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND,   *
 * EXPRESS OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES   *
 * OF MERCHANTABILITY, FITNESS FOR A PARTICULAR PURPOSE AND          *
 * NONINFRINGEMENT.  IN NO EVENT SHALL THE CONTRIBUTORS OR           *
 * COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER       *
 * LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE,   *
 * ARISING FROM, OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE    *
 * USE OR OTHER DEALINGS WITH THE SOFTWARE.                          *
 *********************************************************************/
// $Id: geometry.h,v 1.3 2008/12/06 08:45:24 mtcampbe Exp $

#ifndef __GEOMETRY_H__
#define __GEOMETRY_H__

#include "mopbasic.h"
#include <cmath>
#include <vector>

MOP_BEGIN_NAMESPACE

void min_max3(double a1, double a2, double a3, double &min, 
              double &max);

void min_max4(double a1, double a2, double a3, double a4, 
              double &min, double &max);

void min_max6(double a1, double a2, double a3, double a4, 
              double a5, double a6, double &min, double &max);

void min_max12(double a1, double a2, double a3, double a4, 
               double a5, double a6, double a7, double a8,
               double a9, double a10, double a11, double a12, 
               double &min, double &max);

double angle(Vector_3<double> v1, 
             Vector_3<double> v2);

double f_angle(Vector_3<double> v1, Vector_3<double>v2);

inline double dotP(const Vector_3<double> & v1,
                   const Vector_3<double> & v2){
  return (v1[0]*v2[0] + v1[1]*v2[1] + v1[2]*v2[2]);
}

inline void crossP(const Vector_3<double> & v1, 
            const Vector_3<double> & v2, 
            Vector_3<double> & v3){
  v3[0]= (v1[1]*v2[2]-v2[1]*v1[2]); v3[1]= (v2[0]*v1[2]-v1[0]*v2[2]);
  v3[2]= (v1[0]*v2[1]-v2[0]*v1[1]);  
}

void unit_normal(const Vector_3<double> & v1, 
                 const Vector_3<double> & v2,
                 Vector_3<double> & v3);

double edge_length(const Vector_3<double> & v);

double tri_area(const Vector_3<double> & v1,
                const Vector_3<double> & v2,
                const Vector_3<double> & v3);


double tet_vol(const Vector_3<double> & a, 
               const Vector_3<double> & b, 
               const Vector_3<double> & c);

void printv(const Vector_3<double> & v);

MOP_END_NAMESPACE

#endif