Clean-up cmake (linker) dependencies.

[flightgear.git] / src / Time / light.cxx

diff --git a/src/Time/light.cxx b/src/Time/light.cxx
index 2aedce6ede3b4432e7b3ca353e038a57b352f67e..d5e17686aa4782e2f1c23acea2af169d9fe7c5bd 100644 (file)
--- a/src/Time/light.cxx
+++ b/src/Time/light.cxx
@@ -348,11 +348,8 @@ void FGLight::update_adj_fog_color () {
     // first determine the difference between our view angle and local
     // direction to the sun
     //double vert_rotation = pitch + pitch_offset;
-    double hor_rotation = -(_sun_rotation + SGD_PI) - heading + heading_offset;
-    if (hor_rotation < 0 )
-       hor_rotation = fmod(hor_rotation, SGD_2PI) + SGD_2PI;
-    else
-       hor_rotation = fmod(hor_rotation, SGD_2PI);
+
+    double hor_rotation = -_sun_rotation + SGD_PI - heading + heading_offset;
 
     // revert to unmodified values before using them.
     //
@@ -382,12 +379,19 @@ void FGLight::update_adj_fog_color () {
     if (sif < 1e-4)
        sif = 1e-4;
 
-    float rf1 = fabs((hor_rotation - SGD_PI) / SGD_PI);                // 0.0 .. 1.0
+    float rf1 = fabs(fmod(hor_rotation, SGD_2PI) - SGD_PI) / SGD_PI;
     float rf2 = avf * pow(rf1*rf1, 1/sif) * 1.0639 * _saturation * _scattering;
+
+    // HACK: clamp rf2 to 1.0.
+    // Somehow, rf2 is huge at certain sun angles (around midnight), which results in
+    // rf3 being negative and hence negative fog colors and weird display issues...
+    // Something more fundamental may be wrong with the formulas here...
+    if (rf2>1.0) rf2 = 1.0;
+
     float rf3 = 1.0 - rf2;
 
     gamma = system_gamma * (0.9 - sif*avf);
-    _adj_fog_color = rf3 * _fog_color + rf2 * _sun_color;;
+    _adj_fog_color = rf3 * _fog_color + rf2 * _sun_color;
     gamma_correct_rgb( _adj_fog_color.data(), gamma);
 
     // make sure the colors have their original value before they are being
@@ -401,19 +405,18 @@ void FGLight::update_adj_fog_color () {
 void FGLight::updateSunPos()
 {
     SGTime *t = globals->get_time_params();
-    FGViewer *v = globals->get_current_view();
-
     SG_LOG( SG_EVENT, SG_DEBUG, "  Updating Sun position" );
     SG_LOG( SG_EVENT, SG_DEBUG, "  Gst = " << t->getGst() );
 
     fgSunPositionGST(t->getGst(), &_sun_lon, &_sun_lat);
+
     // It might seem that sun_gc_lat needs to be converted to geodetic
     // latitude here, but it doesn't. The sun latitude is the latitude
     // of the point on the earth where the up vector has the same
     // angle from geocentric Z as the sun direction. But geodetic
     // latitude is defined as 90 - angle of up vector from Z!
-    SGVec3d sunpos(SGVec3d::fromGeoc(SGGeoc::fromRadM(_sun_lon, _sun_lat,
-                                                      SGGeodesy::EQURAD)));
+    SGVec3d sunpos = SGVec3d::fromGeoc(SGGeoc::fromRadM(_sun_lon, _sun_lat,
+                                                          SGGeodesy::EQURAD));
 
     SG_LOG( SG_EVENT, SG_DEBUG, "    t->cur_time = " << t->get_cur_time() );
     SG_LOG( SG_EVENT, SG_DEBUG,
@@ -421,28 +424,33 @@ void FGLight::updateSunPos()
             << " Geodcentric lat = " << _sun_lat );
 
     // update the sun light vector
-    sun_vec() = SGVec4f(toVec3f(normalize(sunpos)), 0);
-    sun_vec_inv() = - sun_vec();
+    _sun_vec = SGVec4f(toVec3f(normalize(sunpos)), 0);
+    _sun_vec_inv = - _sun_vec;
 
     // calculate the sun's relative angle to local up
-    SGVec3d viewPos = v->get_view_pos();
-    SGQuatd hlOr = SGQuatd::fromLonLat(SGGeod::fromCart(viewPos));
-    SGVec3d world_up = hlOr.backTransform(-SGVec3d::e3());
-    SGVec3d nsun = normalize(sunpos);
+    FGViewer *v = globals->get_current_view();
+    SGQuatd hlOr =  SGQuatd::fromLonLat( v->getPosition() );
+    SGVec3d world_up = hlOr.backTransform( -SGVec3d::e3() );
     // cout << "nup = " << nup[0] << "," << nup[1] << ","
     //      << nup[2] << endl;
     // cout << "nsun = " << nsun[0] << "," << nsun[1] << ","
     //      << nsun[2] << endl;
 
-    _sun_angle = acos( dot ( world_up, nsun ) );
+    SGVec3d nsun = normalize(sunpos);
+    SGVec3d nup = normalize(world_up);
+    _sun_angle = acos( dot( nup, nsun ) );
+
+    double signedPI = (_sun_angle < 0.0) ? -SGD_PI : SGD_PI;
+    _sun_angle = fmod(_sun_angle+signedPI, SGD_2PI) - signedPI;
+
     SG_LOG( SG_EVENT, SG_DEBUG, "sun angle relative to current location = "
             << get_sun_angle() );
 
     // Get direction to the sun in the local frame.
     SGVec3d local_sun_vec = hlOr.transform(nsun);
 
-    // Angle from south. XXX Is this correct in the southern hemisphere?
-    _sun_rotation = atan2(local_sun_vec.x(), -local_sun_vec.y());
+    // Angle from south. 
+    _sun_rotation = 2.0 * atan2(local_sun_vec.x(), -local_sun_vec.y());
 
     // cout << "  Sky needs to rotate = " << _sun_rotation << " rads = "
     //      << _sun_rotation * SGD_RADIANS_TO_DEGREES << " degrees." << endl;