universe.cpp File Reference

#include <celmath/mathlib.h>
#include <celmath/vecmath.h>
#include <celmath/intersect.h>
#include <celutil/utf8.h>
#include <cassert>
#include "astro.h"
#include "3dsmesh.h"
#include "meshmanager.h"
#include "universe.h"

Include dependency graph for universe.cpp:

Go to the source code of this file.

Defines
#define	ANGULAR_RES   3.5e-6
Functions
static bool	ApproxPlanetPickTraversal (Body *body, void *info)
static bool	ExactPlanetPickTraversal (Body *body, void *info)

---

Define Documentation

#define ANGULAR_RES   3.5e-6

Definition at line 22 of file universe.cpp. Referenced by ApproxPlanetPickTraversal(), CloseStarPicker::CloseStarPicker(), DSOPicker::DSOPicker(), Universe::pickDeepSkyObject(), Universe::pickPlanet(), Universe::pickStar(), DSOPicker::process(), CloseStarPicker::process(), StarPicker::process(), and StarPicker::StarPicker().

---

Function Documentation

static bool ApproxPlanetPickTraversal (  Body *  body, void *  info )  [static]

Definition at line 303 of file universe.cpp. References ANGULAR_RES, PlanetPickInfo::atanTolerance, PlanetPickInfo::closestApproxDistance, PlanetPickInfo::closestBody, Ray3< T >::direction, distance(), PlanetPickInfo::jd, Vector3< T >::length(), Vector3< T >::normalize(), Ray3< T >::origin, PlanetPickInfo::pickRay, PlanetPickInfo::sinAngle2Closest, and sqrt(). Referenced by Universe::pickPlanet(). { PlanetPickInfo* pickInfo = (PlanetPickInfo*) info; // Reject invisible bodies and bodies that don't exist at the current time if (body->getClassification() == Body::Invisible || !body->extant(pickInfo->jd)) return true; Point3d bpos = body->getHeliocentricPosition(pickInfo->jd); Vec3d bodyDir = bpos - pickInfo->pickRay.origin; double distance = bodyDir.length(); // Check the apparent radius of the orbit against our tolerance factor. // This check exists to make sure than when picking a distant, we select // the planet rather than one of its satellites. float appOrbitRadius = (float) (body->getOrbit()->getBoundingRadius() / distance); if ((pickInfo->atanTolerance > ANGULAR_RES ? pickInfo->atanTolerance: ANGULAR_RES) > appOrbitRadius) { return true; } bodyDir.normalize(); Vec3d bodyMiss = bodyDir - pickInfo->pickRay.direction; double sinAngle2 = sqrt(bodyMiss * bodyMiss)/2.0; if (sinAngle2 <= pickInfo->sinAngle2Closest) { pickInfo->sinAngle2Closest = sinAngle2 > ANGULAR_RES ? sinAngle2 : ANGULAR_RES; pickInfo->closestBody = body; pickInfo->closestApproxDistance = distance; } return true; }

static bool ExactPlanetPickTraversal (  Body *  body, void *  info )  [static]

Definition at line 344 of file universe.cpp. References PlanetPickInfo::closestBody, PlanetPickInfo::closestDistance, conjugate(), Ray3< T >::direction, distance(), ResourceManager< T >::find(), GetModelManager(), InvalidResource, PlanetPickInfo::jd, Vector3< T >::normalize(), Ray3< T >::origin, Model::pick(), PlanetPickInfo::pickRay, sqrt(), testIntersection(), Quaternion< T >::w, Quaternion< T >::x, Quaternion< T >::y, and Quaternion< T >::z. Referenced by Universe::pickPlanet(). { PlanetPickInfo* pickInfo = reinterpret_cast<PlanetPickInfo*>(info); Point3d bpos = body->getHeliocentricPosition(pickInfo->jd); float radius = body->getRadius(); double distance = -1.0; // Test for intersection with the bounding sphere if (body->getClassification() != Body::Invisible && body->extant(pickInfo->jd) && testIntersection(pickInfo->pickRay, Sphered(bpos, radius), distance)) { if (body->getModel() == InvalidResource) { // There's no mesh, so the object is an ellipsoid. If it's // oblate, do a ray intersection test to see if the object was // picked. Otherwise, the object is spherical and we've already // done all the work we need to. if (body->getOblateness() != 0.0f) { // Oblate sphere; use ray ellipsoid intersection calculation Vec3d ellipsoidAxes(radius, radius * (1.0 - body->getOblateness()), radius); // Transform rotate the pick ray into object coordinates Mat3d m = conjugate(body->getEclipticalToEquatorial(pickInfo->jd)).toMatrix3(); Ray3d r(Point3d(0, 0, 0) + (pickInfo->pickRay.origin - bpos), pickInfo->pickRay.direction); r = r * m; if (!testIntersection(r, Ellipsoidd(ellipsoidAxes), distance)) distance = -1.0; } } else { // Transform rotate the pick ray into object coordinates Quatf qf = body->getOrientation(); Quatd qd(qf.w, qf.x, qf.y, qf.z); Mat3d m = conjugate(qd * body->getEclipticalToGeographic(pickInfo->jd)).toMatrix3(); Ray3d r(Point3d(0, 0, 0) + (pickInfo->pickRay.origin - bpos), pickInfo->pickRay.direction); r = r * m; // The mesh vertices are normalized, then multiplied by a scale // factor. Thus, the ray needs to be multiplied by the inverse of // the mesh scale factor. double s = 1.0 / radius; r.origin.x *= s; r.origin.y *= s; r.origin.z *= s; r.direction *= s; Model* model = GetModelManager()->find(body->getModel()); if (model != NULL) { if (!model->pick(r, distance)) distance = -1.0; } } // Make also sure that the pickRay does not intersect the body in the // opposite hemisphere! Hence, need again the "bodyMiss" angle Vec3d bodyDir = bpos - pickInfo->pickRay.origin; bodyDir.normalize(); Vec3d bodyMiss = bodyDir - pickInfo->pickRay.direction; double sinAngle2 = sqrt(bodyMiss * bodyMiss)/2.0; if (sinAngle2 < sin(PI/4.0) && distance > 0.0 && distance <= pickInfo->closestDistance) { pickInfo->closestDistance = distance; pickInfo->closestBody = body; } } return true; }

---