00001 /* 00002 Sheep - A Rigid Body Dynamics Engine 00003 Copyright (C) 2001-2004 Francois Beaune 00004 Contact: http://toxicengine.sourceforge.net/ 00005 00006 This file is part of Sheep. 00007 00008 Sheep is free software; you can redistribute it and/or modify 00009 it under the terms of the GNU General Public License as published by 00010 the Free Software Foundation; either version 2 of the License, or 00011 (at your option) any later version. 00012 00013 Sheep is distributed in the hope that it will be useful, 00014 but WITHOUT ANY WARRANTY; without even the implied warranty of 00015 MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the 00016 GNU General Public License for more details. 00017 00018 You should have received a copy of the GNU General Public License 00019 along with Sheep; if not, write to the Free Software 00020 Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA 00021 */ 00022 00023 #include "arcball.h" // include first 00024 #include "common/math/real.h" 00025 00026 #include <cmath> 00027 00028 using namespace sheep; 00029 00030 Vector3 Arcball::point_on_sphere(const Vector2 &p) const { 00031 Vector3 p1; 00032 00033 p1.m_x = p.m_x; 00034 p1.m_y = p.m_y; 00035 00036 Real m = sq(p.m_x) + sq(p.m_y); 00037 00038 if(m < 1.0) 00039 // Point (x, y) is above the unit sphere. 00040 p1.m_z = sqrt(1.0 - m); 00041 else { 00042 // If the point does not lie on the sphere, push it back 00043 // to the sphere border. 00044 m = sqrt(m); 00045 p1.m_x /= m; 00046 p1.m_y /= m; 00047 p1.m_z = 0.0; 00048 } 00049 00050 assert(feq(p1.Norm(), 1.0)); 00051 00052 return p1; 00053 }
1.3.6