spiceutils.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.
******************************************************************************/
#ifndef UTILS_SPICE_H
#define UTILS_SPICE_H
 
#include <vector>
#include <string>
 
#include "core/CartesianVector.hpp"
#include "dynamics/DCM.h"
#include "frames/Frame.h"
#include "core/mathmacros.h"
 
namespace warptwin {
    void loadKernel(std::string kernel);

    class SpiceManager {
    public:
        SpiceManager() {}
        ~SpiceManager() {}

        void loadKernels(std::vector<std::string> kernel_list);

        void registerSpiceBody(const std::string &body_name);

        std::string getBodyFrameName(const std::string &body_name);

        void setCentralBody(const std::string &central_body);

        std::string centralBody() {return _central_body;}

        double getGravitationalParameter(const std::string &body_name);

        CartesianVector3 getSpicePosition(const std::string &body_name, double tdb_time);

        CartesianVector3 getSpiceVelocity(const std::string &body_name, double tdb_time);

        clockwerk::DCM getSpiceAttitude(const std::string &body_name, double tdb_time);

        CartesianVector3 getSpiceAngularVelocity(const std::string &body_name, double tdb_time);
    protected:
        // Record of all SPICE bodies tracked by the SPICE manager
        std::vector<std::string> _spice_bodies;
 
        // String indicating the central body for the SPICE manager
        std::string _central_body; 
 
        // Boolean tracking whether central body was already set manually
        bool _manual_set = false;
 
        // Boolean to protect from duplicate spice kernel loading
        bool _kernels_loaded = false;
 
        // Temporary variables for calculation
        double _tmp_state[6];
        double _lt = 0.0;
        double _tmp_dbl[3][3];
        clockwerk::DCM _tmp_dcm_p_pci;
        double _tmp_vec[3];
        CartesianVector3 _tmp_omega_p_pci__pci;
        double _ftmtrx[6][6];
    };
}
 
#endif