+++ /dev/null
-/**
- * \file
- * <!--
- * This file is part of BeRTOS.
- *
- * Bertos is free software; you can redistribute it and/or modify
- * it under the terms of the GNU General Public License as published by
- * the Free Software Foundation; either version 2 of the License, or
- * (at your option) any later version.
- *
- * This program is distributed in the hope that it will be useful,
- * but WITHOUT ANY WARRANTY; without even the implied warranty of
- * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
- * GNU General Public License for more details.
- *
- * You should have received a copy of the GNU General Public License
- * along with this program; if not, write to the Free Software
- * Foundation, Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301 USA
- *
- * As a special exception, you may use this file as part of a free software
- * library without restriction. Specifically, if other files instantiate
- * templates or use macros or inline functions from this file, or you compile
- * this file and link it with other files to produce an executable, this
- * file does not by itself cause the resulting executable to be covered by
- * the GNU General Public License. This exception does not however
- * invalidate any other reasons why the executable file might be covered by
- * the GNU General Public License.
- *
- * Copyright 2010 Develer S.r.l. (http://www.develer.com/)
- *
- * -->
- *
- * \brief DevelGPS: compass routines.
- *
- * \author Andrea Righi <arighi@develer.com>
- */
-
-#include <cfg/compiler.h>
-#include <cfg/debug.h>
-#include <math.h>
-#include "compass.h"
-
-static const char *_compass_heading[] =
-{
- "N", "NNE", "NE", "ENE",
- "E", "ESE", "SE", "SSE",
- "S", "SSW", "SW", "WSW",
- "W", "WNW", "NW", "NNW",
-};
-
-/**
- * Use the Haversine formula to calculate great-circle distances between the
- * two points.
- *
- * The Haversine formula remains particularly well-conditioned for numerical
- * computation even at small distances, unlike calculations based on the
- * spherical law of cosines.
- */
-float distance(float lat1, float lon1, float lat2, float lon2)
-{
- const float PLANET_RADIUS = 6371000;
- float d_lat = deg2rad(lat2 - lat1);
- float d_lon = deg2rad(lon2 - lon1);
-
- float a = sin(d_lat / 2) * sin(d_lat / 2) +
- cos(deg2rad(lat1)) * cos(deg2rad(lat2)) *
- sin(d_lon / 2) * sin(d_lon / 2);
- float c = 2 * atan2(sqrt(a), sqrt(1 - a));
-
- return PLANET_RADIUS * c;
-}
-
-/**
- * Evaluate the bearing (also known as forward azimuth) using spherical law
- * coordinates.
- *
- * The bearing is a straight line along a great-circle arc from the start point
- * to the destination point.
- */
-int bearing(float lat1, float lon1, float lat2, float lon2)
-{
- float res;
-
- res = rad2deg(atan2(sin(deg2rad(lon2 - lon1)) *
- cos(deg2rad(lat2)), cos(deg2rad(lat1)) *
- sin(deg2rad(lat2)) - sin(deg2rad(lat1)) *
- cos(deg2rad(lat2)) * cos(deg2rad(lon2) -
- deg2rad(lon1))));
- return ((int)res + 360) % 360;
-}
-
-const char *compass_heading(int bearing)
-{
- ASSERT(bearing >= 0 && bearing < 360);
- /*
- * bearing / 22.5
- */
- return _compass_heading[(bearing << 4) / 360];
-}