1
0
Fork 0

cleanup of atmosphere patch

Remove some dead code and space changes. Replace boost::tuple with std::pair,
but do use boost::tie.
This commit is contained in:
Tim Moore 2009-09-10 11:42:42 +02:00
parent abc7ae96a5
commit 7b6d15d537
4 changed files with 15 additions and 61 deletions

View file

@ -1,3 +1,5 @@
#include <boost/tuple/tuple.hpp>
#include <simgear/math/SGMath.hxx> #include <simgear/math/SGMath.hxx>
#include <simgear/debug/logstream.hxx> #include <simgear/debug/logstream.hxx>
@ -21,54 +23,6 @@ ISA_layer(7, 80000, 262467, 0.88628, 0.000261718, 196.65, -76.50),
const int ISA_def_size(sizeof(ISA_def) / sizeof(ISA_layer)); const int ISA_def_size(sizeof(ISA_def) / sizeof(ISA_layer));
#if 0
// Pressure within a layer, as a function of height
// and of layer parameters.
// Physics model: standard or nonstandard atmosphere,
// depending on what parameters you pass in.
// Heights in meters, pressures in Pa.
// As always, $lambda is positive in the troposphere,
// and zero in the first part of the stratosphere.
//
double Ph_layer(
const double hh,
const double hb,
const double Pb,
const double Tb,
const double lambda) {
using namespace atmodel;
if (lambda) {
// (1/N) &:=& {g\,m_M \over \lambda\,R}
double N = lambda * Rgas / mm / g;
return Pb * pow((Tb - lambda*(hh - hb)) / Tb, 1/N);
} else {
return Pb * exp(-g * mm / Rgas / Tb * (hh - hb));
}
}
// Pressure as a function of height.
// Valid below 32000 meters,
// i.e. troposphere and first two layers of stratosphere.
// Does not depend on any caching; can be used to
// *construct* caches and interpolation tables.
//
// Height in meters, pressure in pascals.
double old_p_vs_a(const double height) {
using namespace atmodel;
if (height <= 11000.) {
return P_layer(height, 0.0, ISA::P0, ISA::T0, ISA::lam0);
} else if (height <= 20000.) {
return P_layer(height, 11000., 22632.06, 216.65, 0.0);
} else if (height <= 32000.) {
return P_layer(height, 20000., 5474.89, 216.65, -0.001);
}
return 0;
}
#endif
// Pressure within a layer, as a function of height. // Pressure within a layer, as a function of height.
// Physics model: standard or nonstandard atmosphere, // Physics model: standard or nonstandard atmosphere,
// depending on what parameters you pass in. // depending on what parameters you pass in.
@ -112,7 +66,7 @@ double T_layer (
// from the table: the reference height in that layer, // from the table: the reference height in that layer,
// the lapse in that layer, and the cap (if any) for that layer // the lapse in that layer, and the cap (if any) for that layer
// (which we take from the /next/ row of the table, if any). // (which we take from the /next/ row of the table, if any).
tuple<double,double> PT_vs_hpt( pair<double,double> PT_vs_hpt(
const double hh, const double hh,
const double _p0, const double _p0,
const double _t0 const double _t0
@ -145,7 +99,7 @@ tuple<double,double> PT_vs_hpt(
xhgt = (pp+1)->height; xhgt = (pp+1)->height;
} }
if (hh <= xhgt) { if (hh <= xhgt) {
return make_tuple(P_layer(hh, hgt, p0, t0, lapse), return make_pair(P_layer(hh, hgt, p0, t0, lapse),
T_layer(hh, hgt, p0, t0, lapse)); T_layer(hh, hgt, p0, t0, lapse));
} }
p0 = P_layer(xhgt, hgt, p0, t0, lapse); p0 = P_layer(xhgt, hgt, p0, t0, lapse);
@ -155,7 +109,7 @@ tuple<double,double> PT_vs_hpt(
// Should never get here. // Should never get here.
SG_LOG(SG_GENERAL, SG_ALERT, "PT_vs_hpt: ran out of layers"); SG_LOG(SG_GENERAL, SG_ALERT, "PT_vs_hpt: ran out of layers");
return make_tuple(d0,d0); return make_pair(d0, d0);
} }
@ -206,7 +160,7 @@ void FGAtmoCache::tabulate() {
for (double hgt = -1000; hgt <= 32000;) { for (double hgt = -1000; hgt <= 32000;) {
double press,temp; double press,temp;
make_tuple(ref(press), ref(temp)) = PT_vs_hpt(hgt); boost::tie(press, temp) = PT_vs_hpt(hgt);
a_tvs_p->addEntry(press / inHg, hgt / foot); a_tvs_p->addEntry(press / inHg, hgt / foot);
#ifdef DEBUG_EXPORT_P_H #ifdef DEBUG_EXPORT_P_H
@ -258,7 +212,7 @@ void FGAtmoCache::check_model() {
using namespace atmodel; using namespace atmodel;
cache(); cache();
double press,temp; double press,temp;
make_tuple(ref(press), ref(temp)) = PT_vs_hpt(height); boost::tie(press, temp) = PT_vs_hpt(height);
cout << "Height: " << height cout << "Height: " << height
<< " \tpressure: " << press << endl; << " \tpressure: " << press << endl;
cout << "Check: " cout << "Check: "

View file

@ -27,10 +27,9 @@
#include <simgear/compiler.h> #include <simgear/compiler.h>
#include <simgear/math/interpolater.hxx> #include <simgear/math/interpolater.hxx>
#include "boost/tuple/tuple.hpp"
using namespace boost;
#include <cmath> #include <cmath>
#include <utility>
using namespace std; using namespace std;
@ -84,7 +83,7 @@ public:
extern const ISA_layer ISA_def[]; extern const ISA_layer ISA_def[];
tuple<double,double> PT_vs_hpt( std::pair<double,double> PT_vs_hpt(
const double hh, const double hh,
const double _p0 = atmodel::ISA::P0, const double _p0 = atmodel::ISA::P0,
const double _t0 = atmodel::ISA::T0); const double _t0 = atmodel::ISA::T0);

View file

@ -18,7 +18,7 @@
// along with this program; if not, write to the Free Software // along with this program; if not, write to the Free Software
// Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301, USA. // Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301, USA.
// //
// $Id: environment.cxx,v 1.1 2009/01/30 15:07:04 jsd Exp $ // $Id$
#ifdef HAVE_CONFIG_H #ifdef HAVE_CONFIG_H
@ -27,6 +27,8 @@
#include <math.h> #include <math.h>
#include <boost/tuple/tuple.hpp>
#include <plib/sg.h> #include <plib/sg.h>
#include <simgear/constants.h> #include <simgear/constants.h>
@ -634,7 +636,7 @@ FGEnvironment::_recalc_alt_pt ()
} }
#endif #endif
double press, temp; double press, temp;
make_tuple(ref(press), ref(temp)) = PT_vs_hpt(elevation_ft * foot, boost::tie(press, temp) = PT_vs_hpt(elevation_ft * foot,
pressure_sea_level_inhg * inHg, temperature_sea_level_degc + freezing); pressure_sea_level_inhg * inHg, temperature_sea_level_degc + freezing);
temperature_degc = temp - freezing; temperature_degc = temp - freezing;
pressure_inhg = press / inHg; pressure_inhg = press / inHg;

View file

@ -718,7 +718,6 @@ struct GeneralInitOperation : public GraphicsContextOperation
static void fgIdleFunction ( void ) { static void fgIdleFunction ( void ) {
static osg::ref_ptr<GeneralInitOperation> genOp; static osg::ref_ptr<GeneralInitOperation> genOp;
////////cerr << "Idle state: " << idle_state << endl;
if ( idle_state == 0 ) { if ( idle_state == 0 ) {
idle_state++; idle_state++;
// Pick some window on which to do queries. // Pick some window on which to do queries.