test/src/RunCirclingWind.cpp

   1 /*
   2 Copyright_License {
   3
   4   XCSoar Glide Computer - http://www.xcsoar.org/
   5   Copyright (C) 2000-2013 The XCSoar Project
   6   A detailed list of copyright holders can be found in the file "AUTHORS".
   7
   8   This program is free software; you can redistribute it and/or
   9   modify it under the terms of the GNU General Public License
  10   as published by the Free Software Foundation; either version 2
  11   of the License, or (at your option) any later version.
  12
  13   This program is distributed in the hope that it will be useful,
  14   but WITHOUT ANY WARRANTY; without even the implied warranty of
  15   MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
  16   GNU General Public License for more details.
  17
  18   You should have received a copy of the GNU General Public License
  19   along with this program; if not, write to the Free Software
  20   Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA  02111-1307, USA.
  21 }
  22 */
  23
  24 #include "Computer/CirclingComputer.hpp"
  25 #include "Computer/Wind/CirclingWind.hpp"
  26 #include "OS/Args.hpp"
  27 #include "DebugReplay.hpp"
  28 #include "Formatter/TimeFormatter.hpp"
  29 #include "Computer/Settings.hpp"
  30
  31 #include <stdio.h>
  32 #include <memory>
  33
  34 int main(int argc, char **argv)
  35 {
  36   Args args(argc, argv, "DRIVER FILE");
  37   std::unique_ptr<DebugReplay> replay(CreateDebugReplay(args));
  38   if (!replay)
  39     return EXIT_FAILURE;
  40
  41   args.ExpectEnd();
  42
  43   printf("# time quality wind_bearing (deg) wind_speed (m/s)\n");
  44
  45   CirclingSettings circling_settings;
  46   circling_settings.SetDefaults();
  47
  48   CirclingComputer circling_computer;
  49   circling_computer.Reset();
  50
  51   CirclingWind circling_wind;
  52   circling_wind.Reset();
  53
  54   while (replay->Next()) {
  55     circling_computer.TurnRate(replay->SetCalculated(),
  56                                replay->Basic(),
  57                                replay->Calculated().flight);
  58     circling_computer.Turning(replay->SetCalculated(),
  59                               replay->Basic(),
  60                               replay->Calculated().flight,
  61                               circling_settings);
  62
  63     CirclingWind::Result result = circling_wind.NewSample(replay->Basic(),
  64                                                           replay->Calculated());
  65     if (result.quality > 0) {
  66       TCHAR time_buffer[32];
  67       FormatTime(time_buffer, replay->Basic().time);
  68
  69       _tprintf(_T("%s %d %d %g\n"),
  70                time_buffer, result.quality,
  71                (int)result.wind.bearing.Degrees(),
  72                (double)result.wind.norm);
  73     }
  74   }
  75 }
  76