ootk
Version:
Orbital Object Toolkit including Multiple Propagators, Initial Orbit Determination, and Maneuver Calculations.
84 lines • 3.42 kB
JavaScript
/**
* @author @thkruz Theodore Kruczek
* @license AGPL-3.0-or-later
* @copyright (c) 2025 Kruczek Labs LLC
*
* Orbital Object ToolKit is free software: you can redistribute it and/or modify it under the
* terms of the GNU Affero General Public License as published by the Free Software
* Foundation, either version 3 of the License, or (at your option) any later version.
*
* Orbital Object ToolKit is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY;
* without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.
* See the GNU Affero General Public License for more details.
*
* You should have received a copy of the GNU Affero General Public License along with
* Orbital Object ToolKit. If not, see <http://www.gnu.org/licenses/>.
*/
import { J2000 } from '../main.js';
import { VerletBlendInterpolator } from '../interpolator/VerletBlendInterpolator.js';
import { Propagator } from './Propagator.js';
// / Kepler analytical two-body propagator.
export class KeplerPropagator extends Propagator {
initElements_;
elements_;
cacheState_;
checkpoints_;
constructor(initElements) {
super();
this.initElements_ = initElements;
this.elements_ = initElements;
this.cacheState_ = J2000.fromClassicalElements(initElements);
this.checkpoints_ = [];
}
get state() {
return this.cacheState_;
}
propagate(epoch) {
this.cacheState_ = J2000.fromClassicalElements(this.elements_.propagate(epoch));
return this.cacheState_;
}
reset() {
this.elements_ = this.initElements_;
this.cacheState_ = J2000.fromClassicalElements(this.elements_);
}
// eslint-disable-next-line @typescript-eslint/no-unused-vars
maneuver(maneuver, interval = 60) {
this.cacheState_ = maneuver.apply(this.propagate(maneuver.center));
this.elements_ = this.cacheState_.toClassicalElements();
return [this.cacheState_];
}
ephemerisManeuver(start, finish, maneuvers, interval = 60.0) {
const tMvr = maneuvers.slice(0).filter((mvr) => mvr.center >= start || mvr.center <= finish);
const ephemeris = [];
if (tMvr[0].start > start) {
ephemeris.push(this.propagate(start));
}
for (const mvr of tMvr) {
while (this.cacheState_.epoch < mvr.center) {
const step = Math.min(mvr.center.difference(this.cacheState_.epoch), interval);
this.propagate(this.cacheState_.epoch.roll(step));
if (this.cacheState_.epoch.posix !== mvr.center.posix) {
ephemeris.push(this.cacheState_);
}
}
ephemeris.push(...this.maneuver(mvr, interval));
}
while (this.cacheState_.epoch < finish) {
const step = Math.min(finish.difference(this.cacheState_.epoch), interval);
this.propagate(this.cacheState_.epoch.roll(step));
ephemeris.push(this.cacheState_);
}
return new VerletBlendInterpolator(ephemeris);
}
checkpoint() {
this.checkpoints_.push(this.cacheState_);
return this.checkpoints_.length - 1;
}
clearCheckpoints() {
this.checkpoints_ = [];
}
restore(index) {
this.cacheState_ = this.checkpoints_[index];
}
}
//# sourceMappingURL=KeplerPropagator.js.map