]> git.mxchange.org Git - flightgear.git/blobdiff - src/AIModel/AIBase.cxx
Make the scenerio's work again (now for real) and a small number of updates.
[flightgear.git] / src / AIModel / AIBase.cxx
index 182c3e11459b5918ea14889c5864ad1f3d2e103a..935d5c660ac9ee3c0d310f70b69e435da6d53b9e 100644 (file)
 #include "AIBase.hxx"
 #include "AIManager.hxx"
 
-FGAIBase::FGAIBase() {
+
+const double FGAIBase::e = 2.71828183;
+const double FGAIBase::lbs_to_slugs = 0.031080950172;   //conversion factor
+
+
+FGAIBase::FGAIBase()
+ :  fp( NULL ),
+    model( NULL ),
+    props( NULL ),
+    manager( NULL )
+{
     _type_str = "model";
     tgt_roll = roll = tgt_pitch = tgt_yaw = tgt_vs = vs = pitch = 0.0;
     bearing = elevation = range = rdot = 0.0;
     x_shift = y_shift = rotation = 0.0;
+    in_range = false;
     invisible = true;
     no_roll = true;
+    life = 900;
     model_path = "";
-    model = 0;
     _otype = otNull;
     index = 0;
+    delete_me = false;
 }
 
 FGAIBase::~FGAIBase() {
@@ -61,9 +73,31 @@ FGAIBase::~FGAIBase() {
     // unbind();
     SGPropertyNode *root = globals->get_props()->getNode("ai/models", true);
     root->removeChild(_type_str.c_str(), index);
+    if (fp) delete fp;
+    fp = NULL;
 }
 
 void FGAIBase::update(double dt) {
+    ft_per_deg_lat = 366468.96 - 3717.12 * cos(pos.lat()/SG_RADIANS_TO_DEGREES);
+    ft_per_deg_lon = 365228.16 * cos(pos.lat() / SG_RADIANS_TO_DEGREES);
+
+    // Calculate rho at altitude, using standard atmosphere
+    // For the temperature T and the pressure p,
+
+    if (altitude < 36152) {            // curve fits for the troposphere
+      T = 59 - 0.00356 * altitude;
+      p = 2116 * pow( ((T + 459.7) / 518.6) , 5.256);
+
+    } else if ( 36152 < altitude && altitude < 82345 ) {    // lower stratosphere
+      T = -70;
+      p = 473.1 * pow( e , 1.73 - (0.000048 * altitude) );
+
+    } else {                                    //  upper stratosphere
+      T = -205.05 + (0.00164 * altitude);
+      p = 51.97 * pow( ((T + 459.7) / 389.98) , -11.388);
+    }
+
+    rho = p / (1718 * (T + 459.7));
 }
 
 
@@ -83,8 +117,10 @@ void FGAIBase::Transform() {
 bool FGAIBase::init() {
 
    SGPropertyNode *root = globals->get_props()->getNode("ai/models", true);
+
    index = manager->getNum(_otype) - 1;
    props = root->getNode(_type_str.c_str(), index, true);
+
    if (model_path != "") {
       model = sgLoad3DModel( globals->get_fg_root(),
                             model_path.c_str(),
@@ -114,30 +150,32 @@ bool FGAIBase::isa( object_type otype ) {
 
 
 void FGAIBase::bind() {
-   props->tie("id", SGRawValuePointer<int>(&id));
+   props->tie("id", SGRawValueMethods<FGAIBase,int>(*this,
+                                         &FGAIBase::_getID));
    props->tie("velocities/true-airspeed-kt",  SGRawValuePointer<double>(&speed));
    props->tie("velocities/vertical-speed-fps",
                SGRawValueMethods<FGAIBase,double>(*this,
-                                         FGAIBase::_getVS_fps,
-                                         FGAIBase::_setVS_fps));
+                                         &FGAIBase::_getVS_fps,
+                                         &FGAIBase::_setVS_fps));
 
    props->tie("position/altitude-ft",
                SGRawValueMethods<FGAIBase,double>(*this,
-                                         FGAIBase::_getAltitude,
-                                         FGAIBase::_setAltitude));
+                                         &FGAIBase::_getAltitude,
+                                         &FGAIBase::_setAltitude));
    props->tie("position/latitude-deg",
                SGRawValueMethods<FGAIBase,double>(*this,
-                                         FGAIBase::_getLatitude,
-                                         FGAIBase::_setLatitude));
+                                         &FGAIBase::_getLatitude,
+                                         &FGAIBase::_setLatitude));
    props->tie("position/longitude-deg",
                SGRawValueMethods<FGAIBase,double>(*this,
-                                         FGAIBase::_getLongitude,
-                                         FGAIBase::_setLongitude));
+                                         &FGAIBase::_getLongitude,
+                                         &FGAIBase::_setLongitude));
 
    props->tie("orientation/pitch-deg",   SGRawValuePointer<double>(&pitch));
    props->tie("orientation/roll-deg",    SGRawValuePointer<double>(&roll));
    props->tie("orientation/true-heading-deg", SGRawValuePointer<double>(&hdg));
 
+   props->tie("radar/in-range", SGRawValuePointer<bool>(&in_range));
    props->tie("radar/bearing-deg",   SGRawValuePointer<double>(&bearing));
    props->tie("radar/elevation-deg", SGRawValuePointer<double>(&elevation));
    props->tie("radar/range-nm", SGRawValuePointer<double>(&range));
@@ -148,7 +186,7 @@ void FGAIBase::bind() {
    props->tie("radar/rotation", SGRawValuePointer<double>(&rotation));
 
    props->tie("controls/lighting/nav-lights",
-               SGRawValueMethods<FGAIBase,bool>(*this, _isNight));
+               SGRawValueFunctions<bool>(_isNight));
    props->setBoolValue("controls/lighting/beacon", true);
    props->setBoolValue("controls/lighting/strobe", true);
 }
@@ -166,6 +204,7 @@ void FGAIBase::unbind() {
     props->untie("orientation/roll-deg");
     props->untie("orientation/true-heading-deg");
 
+    props->untie("radar/in-range");
     props->untie("radar/bearing-deg");
     props->untie("radar/elevation-deg");
     props->untie("radar/range-nm");
@@ -178,6 +217,91 @@ void FGAIBase::unbind() {
     props->untie("controls/lighting/nav-lights");
 }
 
+double FGAIBase::UpdateRadar(FGAIManager* manager)
+{
+   double radar_range_ft2 = fgGetDouble("/instrumentation/radar/range");
+   radar_range_ft2 *= SG_NM_TO_METER * SG_METER_TO_FEET * 1.1; // + 10%
+   radar_range_ft2 *= radar_range_ft2;
+
+   double user_latitude  = manager->get_user_latitude();
+   double user_longitude = manager->get_user_longitude();
+   double lat_range = fabs(pos.lat() - user_latitude) * ft_per_deg_lat;
+   double lon_range = fabs(pos.lon() - user_longitude) * ft_per_deg_lon;
+   double range_ft2 = lat_range*lat_range + lon_range*lon_range;
+
+   //
+   // Test whether the target is within radar range.
+   //
+   in_range = (range_ft2 && (range_ft2 <= radar_range_ft2));
+   if ( in_range )
+   {
+     props->setBoolValue("radar/in-range", true);
+
+     // copy values from the AIManager
+     double user_altitude  = manager->get_user_altitude();
+     double user_heading   = manager->get_user_heading();
+     double user_pitch     = manager->get_user_pitch();
+     double user_yaw       = manager->get_user_yaw();
+     double user_speed     = manager->get_user_speed();
+
+     // calculate range to target in feet and nautical miles
+     double range_ft = sqrt( range_ft2 );
+     range = range_ft / 6076.11549;
+
+     // calculate bearing to target
+     if (pos.lat() >= user_latitude) {
+        bearing = atan2(lat_range, lon_range) * SG_RADIANS_TO_DEGREES;
+        if (pos.lon() >= user_longitude) {
+           bearing = 90.0 - bearing;
+        } else {
+           bearing = 270.0 + bearing;
+        }
+     } else {
+        bearing = atan2(lon_range, lat_range) * SG_RADIANS_TO_DEGREES;
+        if (pos.lon() >= user_longitude) {
+           bearing = 180.0 - bearing;
+        } else {
+           bearing = 180.0 + bearing;
+        }
+     }
+
+     // calculate look left/right to target, without yaw correction
+     horiz_offset = bearing - user_heading;
+     if (horiz_offset > 180.0) horiz_offset -= 360.0;
+     if (horiz_offset < -180.0) horiz_offset += 360.0;
+
+     // calculate elevation to target
+     elevation = atan2( altitude * SG_METER_TO_FEET - user_altitude, range_ft )
+                        * SG_RADIANS_TO_DEGREES;
+
+     // calculate look up/down to target
+     vert_offset = elevation + user_pitch;
+
+     /* this calculation needs to be fixed, but it isn't important anyway
+     // calculate range rate
+     double recip_bearing = bearing + 180.0;
+     if (recip_bearing > 360.0) recip_bearing -= 360.0;
+     double my_horiz_offset = recip_bearing - hdg;
+     if (my_horiz_offset > 180.0) my_horiz_offset -= 360.0;
+     if (my_horiz_offset < -180.0) my_horiz_offset += 360.0;
+     rdot = (-user_speed * cos( horiz_offset * SG_DEGREES_TO_RADIANS ))
+             +(-speed * 1.686 * cos( my_horiz_offset * SG_DEGREES_TO_RADIANS ));
+*/
+
+     // now correct look left/right for yaw
+     horiz_offset += user_yaw;
+
+     // calculate values for radar display
+     y_shift = range * cos( horiz_offset * SG_DEGREES_TO_RADIANS);
+     x_shift = range * sin( horiz_offset * SG_DEGREES_TO_RADIANS);
+     rotation = hdg - user_heading;
+     if (rotation < 0.0) rotation += 360.0;
+
+   }
+
+   return range_ft2;
+}
+
 
 /*
  * getters and Setters
@@ -212,7 +336,10 @@ void FGAIBase::_setAltitude( double _alt ) {
     setAltitude( _alt );
 }
 
-bool FGAIBase::_isNight() const {
+bool FGAIBase::_isNight() {
     return (fgGetFloat("/sim/time/sun-angle-rad") > 1.57);
 }
 
+int FGAIBase::_getID() const {
+    return (int)(this);
+}