Math/vector.hxx

   1 // vector.hxx -- additional vector routines
   2 //
   3 // Written by Curtis Olson, started December 1997.
   4 //
   5 // Copyright (C) 1997  Curtis L. Olson  - curt@infoplane.com
   6 //
   7 // This program is free software; you can redistribute it and/or
   8 // modify it under the terms of the GNU General Public License as
   9 // published by the Free Software Foundation; either version 2 of the
  10 // License, or (at your option) any later version.
  11 //
  12 // This program is distributed in the hope that it will be useful, but
  13 // WITHOUT ANY WARRANTY; without even the implied warranty of
  14 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  15 // General Public License for more details.
  16 //
  17 // You should have received a copy of the GNU General Public License
  18 // along with this program; if not, write to the Free Software
  19 // Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA.
  20 //
  21 // $Id$
  22 // (Log is kept at end of this file)
  23
  24
  25 #ifndef _VECTOR_HXX
  26 #define _VECTOR_HXX
  27
  28
  29 #ifndef __cplusplus
  30 # error This library requires C++
  31 #endif
  32
  33
  34 #include "mat3.h"
  35
  36
  37 // Map a vector onto the plane specified by normal
  38 #if defined( USE_XTRA_MAT3_INLINES )
  39 #  define map_vec_onto_cur_surface_plane(normal, v0, vec, result) { \
  40         double scale = ((normal[0]*vec[0]+normal[1]*vec[1]+normal[2]*vec[2]) / \
  41                (normal[0]*normal[0]+normal[1]*normal[1]+normal[2]*normal[2])); \
  42         result[0] = vec[0]-normal[0]*scale; \
  43         result[1] = vec[1]-normal[1]*scale; \
  44         result[2] = vec[2]-normal[2]*scale; \
  45   }
  46 #else
  47   void map_vec_onto_cur_surface_plane(MAT3vec normal, MAT3vec v0, MAT3vec vec,
  48                                     MAT3vec result);
  49 #endif //defined( USE_XTRA_MAT3_INLINES )
  50
  51
  52 // Given a point p, and a line through p0 with direction vector d,
  53 // find the shortest distance from the point to the line
  54 double fgPointLine(MAT3vec p, MAT3vec p0, MAT3vec d);
  55
  56
  57 // Given a point p, and a line through p0 with direction vector d,
  58 // find the shortest distance (squared) from the point to the line
  59 double fgPointLineSquared(MAT3vec p, MAT3vec p0, MAT3vec d);
  60
  61
  62 #endif // _VECTOR_HXX
  63
  64
  65 // $Log$
  66 // Revision 1.4  1998/10/16 23:36:39  curt
  67 // c++-ifying.
  68 //
  69 // Revision 1.3  1998/08/24 20:04:13  curt
  70 // Various "inline" code optimizations contributed by Norman Vine.
  71 //
  72 // Revision 1.2  1998/07/24 21:34:38  curt
  73 // fgPointLine() rewritten into fgPointLineSquared() ... this ultimately saves
  74 // us from doing a sqrt().
  75 //
  76 // Revision 1.1  1998/07/08 14:40:10  curt
  77 // polar3d.[ch] renamed to polar3d.[ch]xx, vector.[ch] renamed to vector.[ch]xx
  78 // Updated fg_geodesy comments to reflect that routines expect and produce
  79 //   meters.
  80 //
  81 // Revision 1.4  1998/04/21 17:03:51  curt
  82 // Prepairing for C++ integration.
  83 //
  84 // Revision 1.3  1998/01/22 02:59:39  curt
  85 // Changed #ifdef FILE_H to #ifdef _FILE_H
  86 //
  87 // Revision 1.2  1998/01/19 19:27:14  curt
  88 // Merged in make system changes from Bob Kuehne <rpk@sgi.com>
  89 // This should simplify things tremendously.
  90 //
  91 // Revision 1.1  1997/12/22 04:13:18  curt
  92 // Initial revision.
  93 //
  94