Fastest body flattening retrieval method

[luigigisolfi]

What is the fastest way to retrieve information on the flattening of a given body?
If there is no current option to get the flattening from BodyShapeModel, I'd propose to have something like:

body.get("BODY_NAME").shape_model.flattening,

which makes sense to me, since we have we have

body.get("BODY_NAME").shape_model.average_radius (coming from exposing the getEquatorialRadius function)

Given that we have already implemented a getPolarRadius function in C++ (not exposed), we could in principle use it to get the flatting via something like

double getFlattening( const std::string &body )
{
    setSpiceErrorHandling( );

    double radii[ 3 ];
    SpiceInt n;

    // 1. Retrieve the 3 radii (a, b, c)
    bodvrd_c( body.c_str( ), "RADII", 3, &n, radii );

    if( failed_c( ) )
    {
        handleSpiceException( );
        return 0.0; // Return safe default or handle error
    }

    // 2. Compute f = (Re - Rp) / Re
    // radii[0] = Equatorial Radius (Re)
    // radii[2] = Polar Radius (Rp)
    
    if (radii[0] == 0.0) return 0.0; // Prevent division by zero

    return ( radii[ 0 ] - radii[ 2 ] ) / radii[ 0 ];
}

and then expose it.
Opinions/faster solutions?

[DominicDirkx]

The average_radius is an attribute of the shape model base class, the flattening only of the derived oblate shape model (since it is not a defined quantity for all shape models).

I'm not sure if the flattening of the oblate shape model is exposed as an attribute in Python, though. The link to spice should not be made inside the shape model, but only during its creation (where it already is/can be).

[luigigisolfi]

Thanks, I solved it via

f = body_settings.get('Earth').shape_settings.flattening