orbitutils.h

Go to the documentation of this file.

/******************************************************************************
* Copyright (c) ATTX INC 2025. All Rights Reserved.
*
* This software and associated documentation (the "Software") are the 
* proprietary and confidential information of ATTX INC. The Software is 
* furnished under a license agreement between ATTX and the user organization 
* and may be used or copied only in accordance with the terms of the agreement.
* Refer to 'license/attx_license.adoc' for standard license terms.
*
* EXPORT CONTROL NOTICE: THIS SOFTWARE MAY INCLUDE CONTENT CONTROLLED UNDER THE
* INTERNATIONAL TRAFFIC IN ARMS REGULATIONS (ITAR) OR THE EXPORT ADMINISTRATION 
* REGULATIONS (EAR99). No part of the Software may be used, reproduced, or 
* transmitted in any form or by any means, for any purpose, without the express 
* written permission of ATTX INC.
******************************************************************************/
/*
Orbit Utilities header file
---------------------------
This file contains orbit utilities for converting between "classical"
keplerian orbital elements and xyz cartesian position and velocity.
 
It also contains a number of utilities for calculating parameters associated
with keplerian orbits.
 
TODO: Protect edge cases in these files. Right now they only work for 
elliptical orbits
 
The converters in this file are based on those found in Vallado's 
Fundamentals of Astrodynamics and Applications, fourth edition, Chapter
2.5. The test cases are drawn from that text as well.
 
Author: Alex Reynolds
*/
#ifndef UTILS_ORBIT_UTILS_H
#define UTILS_ORBIT_UTILS_H
 
#include "core/Matrix.hpp"
#include "core/CartesianVector.hpp"
#include "frames/Frame.h"
#include "dynamics/DCM.h"
 
namespace warptwin {
    double true2eccentric(double true_anomaly);

    int rv2coe(const CartesianVector3 &pos, 
               const CartesianVector3 &vel,
               double mu,
               CartesianVector6 &elements);

    int coe2rv(CartesianVector6 &elements,
            double mu,
            CartesianVector3 &pos, 
            CartesianVector3 &vel);

    clockwerk::DCM dcmPqr2Xyz(double i, double W, double w);

    int lvlhFrame(const CartesianVector3 &pos, 
                  const CartesianVector3 &vel, 
                  clockwerk::DCM &lvlh_frame, 
                  CartesianVector3 &lvlh_ang_vel__lvlh);
 
}
 
#endif