simgear/math/project.cxx

   1 // Copyright (C) 2010  Tim Moore moore@bricoworks.com
   2 //
   3 // This library is free software; you can redistribute it and/or
   4 // modify it under the terms of the GNU Library General Public
   5 // License as published by the Free Software Foundation; either
   6 // version 2 of the License, or (at your option) any later version.
   7 //
   8 // This library is distributed in the hope that it will be useful,
   9 // but WITHOUT ANY WARRANTY; without even the implied warranty of
  10 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  11 // Library General Public License for more details.
  12 //
  13 // You should have received a copy of the GNU General Public License
  14 // along with this program; if not, write to the Free Software
  15 // Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA  02110-1301, USA.
  16 //
  17
  18 #include "project.hxx"
  19
  20 #include <osg/Math>
  21 #include <osg/Matrixd>
  22
  23 namespace simgear
  24 {
  25 GLint project(GLdouble objX, GLdouble objY, GLdouble objZ,
  26               const GLdouble *model, const GLdouble *proj, const GLint *view,
  27               GLdouble* winX, GLdouble* winY, GLdouble* winZ)
  28 {
  29     using namespace osg;
  30     Vec4d obj(objX, objY, objZ, 1.0);
  31     Matrixd Mmodel(model), Mproj(proj);
  32     Matrixd Mwin = (Matrixd::translate(1.0, 1.0, 1.0)
  33                     * Matrixd::scale(0.5 * view[2], 0.5 * view[3], 0.5)
  34                     * Matrixd::translate(view[0], view[1], 0.0));
  35     Vec4d result = obj * Mmodel * Mproj * Mwin;
  36     if (equivalent(result.w(), 0.0))
  37         return GL_FALSE;
  38     result = result / result.w();
  39     *winX = result.x();  *winY = result.y();  *winZ = result.z();
  40     return GL_TRUE;
  41 }
  42
  43 }