Albedo/IR numerical partials also working

tudat-team · Oct 23, 2023 · 5bb8b05 · 5bb8b05
1 parent f781174
commit 5bb8b05
Show file tree

Hide file tree

Showing 4 changed files with 1,540 additions and 1,532 deletions.
diff --git a/include/tudat/simulation/estimation_setup/createAccelerationPartials.h b/include/tudat/simulation/estimation_setup/createAccelerationPartials.h
@@ -179,6 +179,10 @@ std::vector< std::shared_ptr< orbit_determination::TidalLoveNumberPartialInterfa
         const std::string& acceleratingBodyName );
 
 
+std::shared_ptr< estimatable_parameters::NumericalAccelerationPartialSettings > getDefaultPanelledSurfaceRadiationPressurePartialSettings(
+    const std::string bodyUndergoingAcceleration,
+    const std::string bodyExertingAcceleration );
+
 template< typename InitialStateParameterType = double >
 std::shared_ptr< acceleration_partials::AccelerationPartial > createRadiationPressureAccelerationPartial(
     std::shared_ptr< electromagnetism::RadiationPressureAcceleration > radiationPressureAccelerationModel,
@@ -216,16 +220,16 @@ std::shared_ptr< acceleration_partials::AccelerationPartial > createRadiationPre
             throw std::runtime_error( "Error when creating custom partials for non-isotropic source radiation pressure; no estimated parameters found" );
         }
 
-        customPartialCalculator = createCustomAccelerationPartial< InitialStateParameterType >(
-            parametersToEstimate, radiationPressureAccelerationModel, acceleratedBody, acceleratingBody, bodies );
         if( customPartialCalculator == nullptr )
         {
             for( unsigned  int i = 0; i < parametersToEstimate->getEstimatedInitialStateParameters( ).size( ); i++ )
             {
-                parametersToEstimate->getEstimatedInitialStateParameters( ).at( i )->setCustomPartialSettings(
-                    std::make_shared< estimatable_parameters::NumericalAccelerationPartialSettings >(
-                        ( Eigen::VectorXd( 6 ) << 10.0, 10.0, 10.0, 1.0E-3, 1.0E-3, 1.0E-3 ).finished( ),
-                        acceleratedBody.first, acceleratingBody.first, basic_astrodynamics::radiation_pressure ) );
+                if( parametersToEstimate->getEstimatedInitialStateParameters( ).at( i )->getCustomPartialSettings( ) == nullptr )
+                {
+                    parametersToEstimate->getEstimatedInitialStateParameters( ).at( i )->setCustomPartialSettings(
+                        getDefaultPanelledSurfaceRadiationPressurePartialSettings( acceleratedBody.first, acceleratingBody.first ) );
+                }
+
             }
             customPartialCalculator = createCustomAccelerationPartial< InitialStateParameterType >(
                 parametersToEstimate, radiationPressureAccelerationModel, acceleratedBody, acceleratingBody, bodies );

diff --git a/...ro/orbit_determination/acceleration_partials/fullRadiationPressureAccelerationPartial.cpp b/...ro/orbit_determination/acceleration_partials/fullRadiationPressureAccelerationPartial.cpp
@@ -28,11 +28,6 @@ RadiationPressureAccelerationPartial::RadiationPressureAccelerationPartial(
     currentPartialWrtUndergoingState_.setZero( );
     currentPartialWrtExertingState_.setZero( );
 
-    std::cout<<"Initial state partials "<<std::endl;
-    for( auto it : customAccelerationPartialSet->customInitialStatePartials_ )
-    {
-        std::cout<<"Entries: "<<it.first.first<<" "<<it.first.second.first<<" "<<it.first.second.second<<std::endl;
-    }
     estimatable_parameters::EstimatebleParameterIdentifier undergoingBodyIdentifier =
         std::make_pair( estimatable_parameters::initial_body_state, std::make_pair( acceleratedBody, "" ) );
     if( customAccelerationPartialSet->customInitialStatePartials_.count( undergoingBodyIdentifier ) > 0 )
@@ -52,10 +47,8 @@ void RadiationPressureAccelerationPartial::update( const double currentTime )
 {
     radiationPressureAcceleration_->updateMembers( currentTime );
 
-    std::cout<<currentTime<<" "<<currentTime_<<std::endl;
     if( !( currentTime_ == currentTime ) )
     {
-        std::cout<<"Partials "<<bodyUndergoingPositionPartial_<<" "<<bodyExertingPositionPartial_<<std::endl;
         if( bodyUndergoingPositionPartial_ != nullptr )
         {
             currentPartialWrtUndergoingState_ = bodyUndergoingPositionPartial_->computePartial(

diff --git a/src/simulation/estimation_setup/createAccelerationPartials.cpp b/src/simulation/estimation_setup/createAccelerationPartials.cpp
@@ -75,6 +75,15 @@ std::vector< std::shared_ptr< orbit_determination::TidalLoveNumberPartialInterfa
 }
 
 
+std::shared_ptr< estimatable_parameters::NumericalAccelerationPartialSettings > getDefaultPanelledSurfaceRadiationPressurePartialSettings(
+    const std::string bodyUndergoingAcceleration,
+    const std::string bodyExertingAcceleration )
+{
+    return std::make_shared< estimatable_parameters::NumericalAccelerationPartialSettings >(
+        ( Eigen::VectorXd( 6 ) << 100.0, 100.0, 100.0, 1.0E-3, 1.0E-3, 1.0E-3 ).finished( ),
+        bodyUndergoingAcceleration, bodyExertingAcceleration, basic_astrodynamics::radiation_pressure );
+}
+
 }
 
 }