AbstractAccelerationModel
├── AbstractGravityModel{T}
│   ├── Point mass (central body)
│   ├── Third-body gravity
│   ├── Relativistic corrections
│   ├── Spherical harmonics (ICGEM/PDS)
│   ├── Polyhedron
│   ├── Mascons
│   └── PINN (neural network)
└── AbstractSolarPressureModel{T}
    ├── Cannonball SRP
    └── Flat plate SRP

# All models implement this interface
a = compute_acceleration(model, args...)

# Gravity models also implement
U = compute_potential(model, args...)

# Parse coefficient data from ICGEM or PDS format files
data = parse_data(Float64, DataType, filename; kwargs...)
model = parse_model(Float64, ModelType, DataType, args...)