src/FDM/YASim/yasim-test.cpp

   1 #include <stdio.h>
   2
   3 #include <simgear/props/props.hxx>
   4 #include <simgear/xml/easyxml.hxx>
   5
   6 #include "FGFDM.hpp"
   7 #include "Airplane.hpp"
   8
   9 using namespace yasim;
  10
  11 // Stubs.  Not needed by a batch program, but required to link.
  12 bool fgSetFloat (const char * name, float val) { return false; }
  13 bool fgSetBool(char const * name, bool val) { return false; }
  14 bool fgGetBool(char const * name, bool def) { return false; }
  15 SGPropertyNode* fgGetNode (const char * path, bool create) { return 0; }
  16 float fgGetFloat (const char * name, float defaultValue) { return 0; }
  17 float fgGetDouble (const char * name, double defaultValue) { return 0; }
  18 float fgSetDouble (const char * name, double defaultValue) { return 0; }
  19
  20 static const float RAD2DEG = 57.2957795131;
  21
  22 int main(int argc, char** argv)
  23 {
  24     FGFDM fdm;
  25     Airplane* a = fdm.getAirplane();
  26
  27     // Read
  28     try {
  29         readXML(argv[1], fdm);
  30     } catch (const sg_exception &e) {
  31         printf("XML parse error: %s (%s)\n",
  32                e.getFormattedMessage().c_str(), e.getOrigin().c_str());
  33     }
  34
  35     // ... and run
  36     a->compile();
  37
  38     float aoa = a->getCruiseAoA() * RAD2DEG;
  39     float tail = -1 * a->getTailIncidence() * RAD2DEG;
  40     float drag = 1000 * a->getDragCoefficient();
  41     float cg[3];
  42     a->getModel()->getBody()->getCG(cg);
  43
  44     printf("Solution results:");
  45     printf("       Iterations: %d\n", a->getSolutionIterations());
  46     printf(" Drag Coefficient: %f\n", drag);
  47     printf("       Lift Ratio: %f\n", a->getLiftRatio());
  48     printf("       Cruise AoA: %f\n", aoa);
  49     printf("   Tail Incidence: %f\n", tail);
  50     printf("Approach Elevator: %f\n", a->getApproachElevator());
  51     printf("               CG: %.3f, %.3f, %.3f\n", cg[0], cg[1], cg[2]);
  52
  53     if(a->getFailureMsg())
  54         printf("SOLUTION FAILURE: %s\n", a->getFailureMsg());
  55 }