geometrytools.h

   1 //
   2 // ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
   3 // +                                                                      +
   4 // + This file is part of enGrid.                                         +
   5 // +                                                                      +
   6 // + Copyright 2008,2009 Oliver Gloth                                     +
   7 // +                                                                      +
   8 // + enGrid is free software: you can redistribute it and/or modify       +
   9 // + it under the terms of the GNU General Public License as published by +
  10 // + the Free Software Foundation, either version 3 of the License, or    +
  11 // + (at your option) any later version.                                  +
  12 // +                                                                      +
  13 // + enGrid is distributed in the hope that it will be useful,            +
  14 // + but WITHOUT ANY WARRANTY; without even the implied warranty of       +
  15 // + MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the        +
  16 // + GNU General Public License for more details.                         +
  17 // +                                                                      +
  18 // + You should have received a copy of the GNU General Public License    +
  19 // + along with enGrid. If not, see <http://www.gnu.org/licenses/>.       +
  20 // +                                                                      +
  21 // ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
  22 //
  23 #ifndef geometrytools_H
  24 #define geometrytools_H
  25
  26 #include "math/mathvector.h"
  27 #include "math/smallsquarematrix.h"
  28
  29 #include <vtkUnstructuredGrid.h>
  30
  31 namespace GeometryTools
  32 {
  33
  34 /** Converts radians to degrees */
  35 double rad2deg( double rad );
  36
  37 /** Converts degrees to radians */
  38 double deg2rad( double deg );
  39
  40 void rotate(vec3_t g1, vec3_t g2, vec3_t g3, vec3_t &b, double theta);
  41
  42 /** Rotates vector v around vector axis by an angle theta */
  43 vec3_t rotate(vec3_t v, vec3_t axis, double theta);
  44
  45 vec3_t orthogonalVector(vec3_t v);
  46
  47 double intersection(vec3_t x_straight, vec3_t v_straight,
  48                     vec3_t x_plane, vec3_t n_plane);
  49
  50 double intersection(vec3_t x_straight, vec3_t v_straight,
  51                     vec3_t x_plane, vec3_t u_plane, vec3_t v_plane);
  52
  53 bool intersection (double &k1, double &k2, vec2_t r1, vec2_t u1, vec2_t r2, vec2_t u2);
  54
  55 void sliceTriangle(const vector<vec3_t> &Tin, vec3_t x, vec3_t n, vector<vector<vec3_t> > &Tout);
  56
  57 double tetraVol(vec3_t x1, vec3_t x2, vec3_t x3, vec3_t x4, bool neg = false);
  58
  59 double pyraVol(vec3_t x1, vec3_t x2, vec3_t x3, vec3_t x4, vec3_t x5, bool neg = false);
  60
  61 double prismVol(vec3_t x1, vec3_t x2, vec3_t x3, vec3_t x4, vec3_t x5, vec3_t x6, bool neg = false);
  62
  63 double hexaVol(vec3_t x1, vec3_t x2, vec3_t x3, vec3_t x4, vec3_t x5, vec3_t x6, vec3_t x7, vec3_t x8,  bool neg = false);
  64
  65 double triArea(vec3_t x1, vec3_t x2, vec3_t x3);
  66
  67 double quadArea(vec3_t x1, vec3_t x2, vec3_t x3, vec3_t x4);
  68
  69 vec3_t triNormal(vec3_t x1, vec3_t x2, vec3_t x3);
  70
  71 vec3_t quadNormal(vec3_t x1, vec3_t x2, vec3_t x3, vec3_t x4);
  72
  73 vec3_t triNormal(vtkUnstructuredGrid *grid, vtkIdType p1, vtkIdType p2, vtkIdType p3);
  74
  75 vec3_t quadNormal(vtkUnstructuredGrid *grid, vtkIdType p1, vtkIdType p2, vtkIdType p3, vtkIdType p4);
  76
  77 vec3_t cellNormal(vtkUnstructuredGrid *grid, vtkIdType i);
  78
  79 double cellVA(vtkUnstructuredGrid *grid, vtkIdType cellId, bool neg = false);
  80
  81 inline vec2_t turnRight(const vec2_t &v)
  82 {
  83   vec2_t u;
  84   u[0] =  v[1];
  85   u[1] = -v[0];
  86   return u;
  87 };
  88
  89 inline vec2_t turnLeft(const vec2_t &v)
  90 {
  91   vec2_t u;
  92   u[0] = -v[1];
  93   u[1] =  v[0];
  94   return u;
  95 };
  96
  97 /** return the angle w.r.t. another 3-vector */
  98 double angle(const vec3_t & u, const vec3_t & v);
  99
 100 /** return the deviation p1->p2->p3 (angle(p2-p1,p3-p2)) */
 101 double deviation(vtkUnstructuredGrid *grid, vtkIdType p1, vtkIdType p2, vtkIdType p3);
 102
 103 /** return the angle p1,p2,p3 (angle(p1-p2,p3-p2)) */
 104 double angle(vtkUnstructuredGrid *grid, vtkIdType p1, vtkIdType p2, vtkIdType p3);
 105
 106 };
 107
 108 #endif