1 // replay.cxx - a system to record and replay FlightGear flights
3 // Written by Curtis Olson, started Juley 2003.
5 // Copyright (C) 2003 Curtis L. Olson - curt@flightgear.org
7 // This program is free software; you can redistribute it and/or
8 // modify it under the terms of the GNU General Public License as
9 // published by the Free Software Foundation; either version 2 of the
10 // License, or (at your option) any later version.
12 // This program is distributed in the hope that it will be useful, but
13 // WITHOUT ANY WARRANTY; without even the implied warranty of
14 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
15 // General Public License for more details.
17 // You should have received a copy of the GNU General Public License
18 // along with this program; if not, write to the Free Software
19 // Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA.
24 #include <simgear/constants.h>
26 #include <FDM/flight.hxx>
27 #include <Main/fg_props.hxx>
28 #include <Network/native_ctrls.hxx>
29 #include <Network/native_fdm.hxx>
30 #include <Network/net_ctrls.hxx>
31 #include <Network/net_fdm.hxx>
35 const double FGReplay::st_list_time = 60.0; // 60 secs of high res data
36 const double FGReplay::mt_list_time = 600.0; // 10 mins of 1 fps data
37 const double FGReplay::lt_list_time = 3600.0; // 1 hr of 10 spf data
39 // short term sample rate is as every frame
40 const double FGReplay::mt_dt = 0.5; // medium term sample rate (sec)
41 const double FGReplay::lt_dt = 5.0; // long term sample rate (sec)
47 FGReplay::FGReplay() {
55 FGReplay::~FGReplay() {
56 // no dynamically allocated memory to free
61 * Initialize the data structures
64 void FGReplay::init() {
69 // Make sure all queues are flushed
70 while ( !short_term.empty() ) {
71 short_term.pop_front();
73 while ( !medium_term.empty() ) {
74 medium_term.pop_front();
76 while ( !medium_term.empty() ) {
77 medium_term.pop_front();
83 * Bind to the property tree
86 void FGReplay::bind() {
92 * Unbind from the property tree
95 void FGReplay::unbind() {
101 * Update the saved data
104 void FGReplay::update( double dt ) {
105 static SGPropertyNode *replay_master
106 = fgGetNode( "/sim/freeze/replay", true );
108 if ( replay_master->getBoolValue() ) {
109 // don't record the replay session
115 // build the replay record
117 FGProps2NetFDM( &f, false );
119 // sanity check, don't collect data if FDM data isn't good
120 if ( !cur_fdm_state->get_inited() ) {
125 FGProps2NetCtrls( &c, false, false );
128 r.sim_time = sim_time;
132 // update the short term list
133 short_term.push_back( r );
135 FGReplayData st_front = short_term.front();
136 if ( sim_time - st_front.sim_time > st_list_time ) {
137 while ( sim_time - st_front.sim_time > st_list_time ) {
138 st_front = short_term.front();
139 short_term.pop_front();
142 // update the medium term list
143 if ( sim_time - last_mt_time > mt_dt ) {
144 last_mt_time = sim_time;
145 medium_term.push_back( st_front );
147 FGReplayData mt_front = medium_term.front();
148 if ( sim_time - mt_front.sim_time > mt_list_time ) {
149 while ( sim_time - mt_front.sim_time > mt_list_time ) {
150 mt_front = medium_term.front();
151 medium_term.pop_front();
154 // update the long term list
155 if ( sim_time - last_lt_time > lt_dt ) {
156 last_lt_time = sim_time;
157 long_term.push_back( mt_front );
159 FGReplayData lt_front = long_term.front();
160 if ( sim_time - lt_front.sim_time > lt_list_time ) {
161 while ( sim_time - lt_front.sim_time > lt_list_time ) {
162 lt_front = long_term.front();
163 long_term.pop_front();
172 cout << "short term size = " << short_term.size()
173 << " time = " << sim_time - short_term.front().sim_time
175 cout << "medium term size = " << medium_term.size()
176 << " time = " << sim_time - medium_term.front().sim_time
178 cout << "long term size = " << long_term.size()
179 << " time = " << sim_time - long_term.front().sim_time
185 static double weight( double data1, double data2, double ratio,
186 bool rotational = false ) {
188 // special handling of rotational data
189 double tmp = data2 - data1;
190 if ( tmp > SGD_PI ) {
192 } else if ( tmp < -SGD_PI ) {
195 return data1 + tmp * ratio;
197 // normal "linear" data
198 return data1 + ( data2 - data1 ) * ratio;
203 * given two FGReplayData elements and a time, interpolate between them
205 static void update_fdm( FGReplayData frame ) {
206 FGNetFDM2Props( &frame.fdm, false );
207 FGNetCtrls2Props( &frame.ctrls, false, false );
211 * given two FGReplayData elements and a time, interpolate between them
213 static FGReplayData interpolate( double time, FGReplayData f1, FGReplayData f2 )
215 FGReplayData result = f1;
217 FGNetFDM fdm1 = f1.fdm;
218 FGNetFDM fdm2 = f2.fdm;
220 FGNetCtrls ctrls1 = f1.ctrls;
221 FGNetCtrls ctrls2 = f2.ctrls;
223 double ratio = (time - f1.sim_time) / (f2.sim_time - f1.sim_time);
225 // Interpolate FDM data
228 result.fdm.longitude = weight( fdm1.longitude, fdm2.longitude, ratio );
229 result.fdm.latitude = weight( fdm1.latitude, fdm2.latitude, ratio );
230 result.fdm.altitude = weight( fdm1.altitude, fdm2.altitude, ratio );
231 result.fdm.agl = weight( fdm1.agl, fdm2.agl, ratio );
232 result.fdm.phi = weight( fdm1.phi, fdm2.phi, ratio, true );
233 result.fdm.theta = weight( fdm1.theta, fdm2.theta, ratio, true );
234 result.fdm.psi = weight( fdm1.psi, fdm2.psi, ratio, true );
237 result.fdm.phidot = weight( fdm1.phidot, fdm2.phidot, ratio, true );
238 result.fdm.thetadot = weight( fdm1.thetadot, fdm2.thetadot, ratio, true );
239 result.fdm.psidot = weight( fdm1.psidot, fdm2.psidot, ratio, true );
240 result.fdm.vcas = weight( fdm1.vcas, fdm2.vcas, ratio );
241 result.fdm.climb_rate = weight( fdm1.climb_rate, fdm2.climb_rate, ratio );
242 result.fdm.v_north = weight( fdm1.v_north, fdm2.v_north, ratio );
243 result.fdm.v_east = weight( fdm1.v_east, fdm2.v_east, ratio );
244 result.fdm.v_down = weight( fdm1.v_down, fdm2.v_down, ratio );
246 result.fdm.v_wind_body_north
247 = weight( fdm1.v_wind_body_north, fdm2.v_wind_body_north, ratio );
248 result.fdm.v_wind_body_east
249 = weight( fdm1.v_wind_body_east, fdm2.v_wind_body_east, ratio );
250 result.fdm.v_wind_body_down
251 = weight( fdm1.v_wind_body_down, fdm2.v_wind_body_down, ratio );
254 result.fdm.stall_warning
255 = weight( fdm1.stall_warning, fdm2.stall_warning, ratio );
258 result.fdm.A_X_pilot = weight( fdm1.A_X_pilot, fdm2.A_X_pilot, ratio );
259 result.fdm.A_Y_pilot = weight( fdm1.A_Y_pilot, fdm2.A_Y_pilot, ratio );
260 result.fdm.A_Z_pilot = weight( fdm1.A_Z_pilot, fdm2.A_Z_pilot, ratio );
265 for ( i = 0; i < fdm1.num_engines; ++i ) {
266 result.fdm.eng_state[i] = fdm1.eng_state[i];
267 result.fdm.rpm[i] = weight( fdm1.rpm[i], fdm2.rpm[i], ratio );
268 result.fdm.fuel_flow[i]
269 = weight( fdm1.fuel_flow[i], fdm2.fuel_flow[i], ratio );
270 result.fdm.egt[i] = weight( fdm1.egt[i], fdm2.egt[i], ratio );
271 result.fdm.mp_osi[i] = weight( fdm1.mp_osi[i], fdm2.mp_osi[i], ratio );
272 result.fdm.oil_temp[i]
273 = weight( fdm1.oil_temp[i], fdm2.oil_temp[i], ratio );
274 result.fdm.oil_px[i] = weight( fdm1.oil_px[i], fdm2.oil_px[i], ratio );
278 for ( i = 0; i < fdm1.num_tanks; ++i ) {
279 result.fdm.fuel_quantity[i]
280 = weight( fdm1.fuel_quantity[i], fdm2.fuel_quantity[i], ratio );
284 for ( i = 0; i < fdm1.num_wheels; ++i ) {
285 result.fdm.wow[i] = weight( fdm1.wow[i], fdm2.wow[i], ratio );
286 result.fdm.gear_pos[i]
287 = weight( fdm1.gear_pos[i], fdm2.gear_pos[i], ratio );
288 result.fdm.gear_steer[i]
289 = weight( fdm1.gear_steer[i], fdm2.gear_steer[i], ratio );
290 result.fdm.gear_compression[i]
291 = weight( fdm1.gear_compression[i], fdm2.gear_compression[i],
296 result.fdm.cur_time = fdm1.cur_time;
297 result.fdm.warp = fdm1.warp;
298 result.fdm.visibility = weight( fdm1.visibility, fdm2.visibility, ratio );
300 // Control surface positions (normalized values)
301 result.fdm.elevator = weight( fdm1.elevator, fdm2.elevator, ratio );
302 result.fdm.flaps = weight( fdm1.flaps, fdm2.flaps, ratio );
303 result.fdm.left_aileron
304 = weight( fdm1.left_aileron, fdm2.left_aileron, ratio );
305 result.fdm.right_aileron
306 = weight( fdm1.right_aileron, fdm2.right_aileron, ratio );
307 result.fdm.rudder = weight( fdm1.rudder, fdm2.rudder, ratio );
308 result.fdm.speedbrake = weight( fdm1.speedbrake, fdm2.speedbrake, ratio );
309 result.fdm.spoilers = weight( fdm1.spoilers, fdm2.spoilers, ratio );
311 // Interpolate Control input data
314 result.ctrls.aileron = weight( ctrls1.aileron, ctrls2.aileron, ratio );
315 result.ctrls.elevator = weight( ctrls1.elevator, ctrls2.elevator, ratio );
316 result.ctrls.elevator_trim
317 = weight( ctrls1.elevator_trim, ctrls2.elevator_trim, ratio );
318 result.ctrls.rudder = weight( ctrls1.rudder, ctrls2.rudder, ratio );
319 result.ctrls.flaps = weight( ctrls1.flaps, ctrls2.flaps, ratio );
320 result.ctrls.flaps_power = ctrls1.flaps_power;
321 result.ctrls.flap_motor_ok = ctrls1.flap_motor_ok;
324 for ( i = 0; i < ctrls1.num_engines; ++i ) {
325 result.ctrls.master_bat[i] = ctrls1.master_bat[i];
326 result.ctrls.master_alt[i] = ctrls1.master_alt[i];
327 result.ctrls.magnetos[i] = ctrls1.magnetos[i];
328 result.ctrls.starter_power[i] = ctrls1.starter_power[i];
329 result.ctrls.throttle[i]
330 = weight( ctrls1.throttle[i], ctrls2.throttle[i], ratio );
331 result.ctrls.mixture[i]
332 = weight( ctrls1.mixture[i], ctrls2.mixture[i], ratio );
333 result.ctrls.fuel_pump_power[i] = ctrls1.fuel_pump_power[i];
334 result.ctrls.prop_advance[i]
335 = weight( ctrls1.prop_advance[i], ctrls2.prop_advance[i], ratio );
336 result.ctrls.engine_ok[i] = ctrls1.engine_ok[i];
337 result.ctrls.mag_left_ok[i] = ctrls1.mag_left_ok[i];
338 result.ctrls.mag_right_ok[i] = ctrls1.mag_right_ok[i];
339 result.ctrls.spark_plugs_ok[i] = ctrls1.spark_plugs_ok[i];
340 result.ctrls.oil_press_status[i] = ctrls1.oil_press_status[i];
341 result.ctrls.fuel_pump_ok[i] = ctrls1.fuel_pump_ok[i];
345 for ( i = 0; i < ctrls1.num_tanks; ++i ) {
346 result.ctrls.fuel_selector[i] = ctrls1.fuel_selector[i];
350 result.ctrls.brake_left
351 = weight( ctrls1.brake_left, ctrls2.brake_right, ratio );
352 result.ctrls.brake_right
353 = weight( ctrls1.brake_right, ctrls2.brake_right, ratio );
354 result.ctrls.brake_parking
355 = weight( ctrls1.brake_parking, ctrls2.brake_parking, ratio );
358 result.ctrls.gear_handle = ctrls1.gear_handle;
361 result.ctrls.turbulence_norm = ctrls1.turbulence_norm;
363 // wind and turbulance
364 result.ctrls.wind_speed_kt
365 = weight( ctrls1.wind_speed_kt, ctrls2.wind_speed_kt, ratio );
366 result.ctrls.wind_dir_deg
367 = weight( ctrls1.wind_dir_deg, ctrls2.wind_dir_deg, ratio );
368 result.ctrls.turbulence_norm
369 = weight( ctrls1.turbulence_norm, ctrls2.turbulence_norm, ratio );
371 // other information about environment
372 result.ctrls.hground = weight( ctrls1.hground, ctrls2.hground, ratio );
373 result.ctrls.magvar = weight( ctrls1.magvar, ctrls2.magvar, ratio );
375 // simulation control
376 result.ctrls.speedup = ctrls1.speedup;
377 result.ctrls.freeze = ctrls1.freeze;
383 * interpolate a specific time from a specific list
385 static void interpolate( double time, const replay_list_type &list ) {
387 if ( list.size() == 0 ) {
390 } else if ( list.size() == 1 ) {
391 // handle list size == 1
392 update_fdm( list[0] );
396 unsigned int last = list.size() - 1;
397 unsigned int first = 0;
398 unsigned int mid = ( last + first ) / 2;
403 // cout << " " << first << " <=> " << last << endl;
404 if ( last == first ) {
406 } else if ( list[mid].sim_time < time && list[mid+1].sim_time < time ) {
409 mid = ( last + first ) / 2;
410 } else if ( list[mid].sim_time > time && list[mid+1].sim_time > time ) {
413 mid = ( last + first ) / 2;
419 FGReplayData result = interpolate( time, list[mid], list[mid+1] );
421 update_fdm( result );
426 * Replay a saved frame based on time, interpolate from the two
427 * nearest saved frames.
430 void FGReplay::replay( double time ) {
431 // cout << "replay: " << time << " ";
432 // find the two frames to interpolate between
435 if ( short_term.size() > 0 ) {
436 t1 = short_term.back().sim_time;
437 t2 = short_term.front().sim_time;
439 // replay the most recent frame
440 update_fdm( short_term.back() );
441 // cout << "first frame" << endl;
442 } else if ( time <= t1 && time >= t2 ) {
443 interpolate( time, short_term );
444 // cout << "from short term" << endl;
445 } else if ( medium_term.size() > 0 ) {
446 t1 = short_term.front().sim_time;
447 t2 = medium_term.back().sim_time;
448 if ( time <= t1 && time >= t2 ) {
449 FGReplayData result = interpolate( time,
451 short_term.front() );
452 update_fdm( result );
453 // cout << "from short/medium term" << endl;
455 t1 = medium_term.back().sim_time;
456 t2 = medium_term.front().sim_time;
457 if ( time <= t1 && time >= t2 ) {
458 interpolate( time, medium_term );
459 // cout << "from medium term" << endl;
460 } else if ( long_term.size() > 0 ) {
461 t1 = medium_term.front().sim_time;
462 t2 = long_term.back().sim_time;
463 if ( time <= t1 && time >= t2 ) {
464 FGReplayData result = interpolate( time,
466 medium_term.front());
467 update_fdm( result );
468 // cout << "from medium/long term" << endl;
470 t1 = long_term.back().sim_time;
471 t2 = long_term.front().sim_time;
472 if ( time <= t1 && time >= t2 ) {
473 interpolate( time, long_term );
474 // cout << "from long term" << endl;
476 // replay the oldest long term frame
477 update_fdm( long_term.front() );
478 // cout << "oldest long term frame" << endl;
482 // replay the oldest medium term frame
483 update_fdm( medium_term.front() );
484 // cout << "oldest medium term frame" << endl;
488 // replay the oldest short term frame
489 update_fdm( short_term.front() );
490 // cout << "oldest short term frame" << endl;
498 double FGReplay::get_start_time() {
499 if ( long_term.size() > 0 ) {
500 return long_term.front().sim_time;
501 } else if ( medium_term.size() > 0 ) {
502 return medium_term.front().sim_time;
503 } else if ( short_term.size() ) {
504 return short_term.front().sim_time;
510 double FGReplay::get_end_time() {
511 if ( short_term.size() ) {
512 return short_term.back().sim_time;