Fix merging of audio in various circumstances.
[dcpomatic.git] / src / lib / dcpomatic_time.h
1 /*
2     Copyright (C) 2014-2017 Carl Hetherington <cth@carlh.net>
3
4     This file is part of DCP-o-matic.
5
6     DCP-o-matic is free software; you can redistribute it and/or modify
7     it under the terms of the GNU General Public License as published by
8     the Free Software Foundation; either version 2 of the License, or
9     (at your option) any later version.
10
11     DCP-o-matic is distributed in the hope that it will be useful,
12     but WITHOUT ANY WARRANTY; without even the implied warranty of
13     MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
14     GNU General Public License for more details.
15
16     You should have received a copy of the GNU General Public License
17     along with DCP-o-matic.  If not, see <http://www.gnu.org/licenses/>.
18
19 */
20
21 /** @file  src/lib/dcpomatic_time.h
22  *  @brief Types to describe time.
23  */
24
25 #ifndef DCPOMATIC_TIME_H
26 #define DCPOMATIC_TIME_H
27
28 #include "frame_rate_change.h"
29 #include "dcpomatic_assert.h"
30 #include <boost/optional.hpp>
31 #include <boost/foreach.hpp>
32 #include <stdint.h>
33 #include <cmath>
34 #include <ostream>
35 #include <iomanip>
36 #include <cstdio>
37
38 class dcpomatic_round_up_test;
39
40 /** A time in seconds, expressed as a number scaled up by Time::HZ.  We want two different
41  *  versions of this class, ContentTime and DCPTime, and we want it to be impossible to
42  *  convert implicitly between the two.  Hence there's this template hack.  I'm not
43  *  sure if it's the best way to do it.
44  *
45  *  S is the name of `this' class and O is its opposite (see the typedefs below).
46  */
47 template <class S, class O>
48 class Time
49 {
50 public:
51         Time ()
52                 : _t (0)
53         {}
54
55         typedef int64_t Type;
56
57         explicit Time (Type t)
58                 : _t (t)
59         {}
60
61         explicit Time (Type n, Type d)
62                 : _t (n * HZ / d)
63         {}
64
65         /* Explicit conversion from type O */
66         Time (Time<O, S> d, FrameRateChange f);
67
68         Type get () const {
69                 return _t;
70         }
71
72         bool operator< (Time<S, O> const & o) const {
73                 return _t < o._t;
74         }
75
76         bool operator<= (Time<S, O> const & o) const {
77                 return _t <= o._t;
78         }
79
80         bool operator== (Time<S, O> const & o) const {
81                 return _t == o._t;
82         }
83
84         bool operator!= (Time<S, O> const & o) const {
85                 return _t != o._t;
86         }
87
88         bool operator>= (Time<S, O> const & o) const {
89                 return _t >= o._t;
90         }
91
92         bool operator> (Time<S, O> const & o) const {
93                 return _t > o._t;
94         }
95
96         Time<S, O> operator+ (Time<S, O> const & o) const {
97                 return Time<S, O> (_t + o._t);
98         }
99
100         Time<S, O> & operator+= (Time<S, O> const & o) {
101                 _t += o._t;
102                 return *this;
103         }
104
105         Time<S, O> operator- () const {
106                 return Time<S, O> (-_t);
107         }
108
109         Time<S, O> operator- (Time<S, O> const & o) const {
110                 return Time<S, O> (_t - o._t);
111         }
112
113         Time<S, O> & operator-= (Time<S, O> const & o) {
114                 _t -= o._t;
115                 return *this;
116         }
117
118         /** Round up to the nearest sampling interval
119          *  at some sampling rate.
120          *  @param r Sampling rate.
121          */
122         Time<S, O> ceil (float r) const {
123                 return Time<S, O> (llrint (HZ * frames_ceil(r) / double(r)));
124         }
125
126         Time<S, O> floor (float r) const {
127                 return Time<S, O> (llrint (HZ * frames_floor(r) / double(r)));
128         }
129
130         double seconds () const {
131                 return double (_t) / HZ;
132         }
133
134         Time<S, O> abs () const {
135                 return Time<S, O> (std::abs (_t));
136         }
137
138         template <typename T>
139         int64_t frames_round (T r) const {
140                 /* We must cast to double here otherwise if T is integer
141                    the calculation will round down before we get the chance
142                    to llrint().
143                 */
144                 return llrint (_t * double(r) / HZ);
145         }
146
147         template <typename T>
148         int64_t frames_floor (T r) const {
149                 return ::floor (_t * r / HZ);
150         }
151
152         template <typename T>
153         int64_t frames_ceil (T r) const {
154                 /* We must cast to double here otherwise if T is integer
155                    the calculation will round down before we get the chance
156                    to ceil().
157                 */
158                 return ::ceil (_t * double(r) / HZ);
159         }
160
161         /** @param r Frames per second */
162         template <typename T>
163         void split (T r, int& h, int& m, int& s, int& f) const
164         {
165                 /* Do this calculation with frames so that we can round
166                    to a frame boundary at the start rather than the end.
167                 */
168                 int64_t ff = frames_round (r);
169
170                 h = ff / (3600 * r);
171                 ff -= h * 3600 * r;
172                 m = ff / (60 * r);
173                 ff -= m * 60 * r;
174                 s = ff / r;
175                 ff -= s * r;
176
177                 f = static_cast<int> (ff);
178         }
179
180         template <typename T>
181         std::string timecode (T r) const {
182                 int h;
183                 int m;
184                 int s;
185                 int f;
186                 split (r, h, m, s, f);
187
188                 char buffer[128];
189                 snprintf (buffer, sizeof (buffer), "%02d:%02d:%02d:%02d", h, m, s, f);
190                 return buffer;
191         }
192
193
194         static Time<S, O> from_seconds (double s) {
195                 return Time<S, O> (llrint (s * HZ));
196         }
197
198         template <class T>
199         static Time<S, O> from_frames (int64_t f, T r) {
200                 DCPOMATIC_ASSERT (r > 0);
201                 return Time<S, O> (f * HZ / r);
202         }
203
204         static Time<S, O> delta () {
205                 return Time<S, O> (1);
206         }
207
208         static Time<S, O> min () {
209                 return Time<S, O> (-INT64_MAX);
210         }
211
212         static Time<S, O> max () {
213                 return Time<S, O> (INT64_MAX);
214         }
215
216 private:
217         friend struct dcptime_ceil_test;
218         friend struct dcptime_floor_test;
219
220         Type _t;
221         static const int HZ = 96000;
222 };
223
224 class ContentTimeDifferentiator {};
225 class DCPTimeDifferentiator {};
226
227 /* Specializations for the two allowed explicit conversions */
228
229 template<>
230 Time<ContentTimeDifferentiator, DCPTimeDifferentiator>::Time (Time<DCPTimeDifferentiator, ContentTimeDifferentiator> d, FrameRateChange f);
231
232 template<>
233 Time<DCPTimeDifferentiator, ContentTimeDifferentiator>::Time (Time<ContentTimeDifferentiator, DCPTimeDifferentiator> d, FrameRateChange f);
234
235 /** Time relative to the start or position of a piece of content in its native frame rate */
236 typedef Time<ContentTimeDifferentiator, DCPTimeDifferentiator> ContentTime;
237 /** Time relative to the start of the output DCP in its frame rate */
238 typedef Time<DCPTimeDifferentiator, ContentTimeDifferentiator> DCPTime;
239
240 template <class T>
241 class TimePeriod
242 {
243 public:
244         TimePeriod () {}
245
246         TimePeriod (T f, T t)
247                 : from (f)
248                 , to (t)
249         {}
250
251         /** start time of sampling interval that the period is from */
252         T from;
253         /** start time of next sampling interval after the period */
254         T to;
255
256         T duration () const {
257                 return to - from;
258         }
259
260         TimePeriod<T> operator+ (T const & o) const {
261                 return TimePeriod<T> (from + o, to + o);
262         }
263
264         boost::optional<TimePeriod<T> > overlap (TimePeriod<T> const & other) const {
265                 T const max_from = std::max (from, other.from);
266                 T const min_to = std::min (to, other.to);
267
268                 if (max_from >= min_to) {
269                         return boost::optional<TimePeriod<T> > ();
270                 }
271
272                 return TimePeriod<T> (max_from, min_to);
273         }
274
275         bool contains (T const & other) const {
276                 return (from <= other && other < to);
277         }
278
279         bool operator< (TimePeriod<T> const & o) const {
280                 if (from != o.from) {
281                         return from < o.from;
282                 }
283                 return to < o.to;
284         }
285
286         bool operator== (TimePeriod<T> const & other) const {
287                 return from == other.from && to == other.to;
288         }
289
290         bool operator!= (TimePeriod<T> const & other) const {
291                 return !(*this == other);
292         }
293 };
294
295 /** @param B Periods to subtract from `A', must be in ascending order of start time and must not overlap */
296 template <class T>
297 std::list<TimePeriod<T> > subtract (TimePeriod<T> A, std::list<TimePeriod<T> > const & B)
298 {
299         std::list<TimePeriod<T> > result;
300         result.push_back (A);
301
302         BOOST_FOREACH (TimePeriod<T> i, B) {
303                 std::list<TimePeriod<T> > new_result;
304                 BOOST_FOREACH (TimePeriod<T> j, result) {
305                         boost::optional<TimePeriod<T> > ov = i.overlap (j);
306                         if (ov) {
307                                 if (*ov == i) {
308                                         /* A contains all of B */
309                                         if (i.from != j.from) {
310                                                 new_result.push_back (TimePeriod<T> (j.from, i.from));
311                                         }
312                                         if (i.to != j.to) {
313                                                 new_result.push_back (TimePeriod<T> (i.to, j.to));
314                                         }
315                                 } else if (*ov == j) {
316                                         /* B contains all of A */
317                                 } else if (i.from < j.from) {
318                                         /* B overlaps start of A */
319                                         new_result.push_back (TimePeriod<T> (i.to, j.to));
320                                 } else if (i.to > j.to) {
321                                         /* B overlaps end of A */
322                                         new_result.push_back (TimePeriod<T> (j.from, i.from));
323                                 }
324                         } else {
325                                 new_result.push_back (j);
326                         }
327                 }
328                 result = new_result;
329         }
330
331         return result;
332 }
333
334 typedef TimePeriod<ContentTime> ContentTimePeriod;
335 typedef TimePeriod<DCPTime> DCPTimePeriod;
336
337 DCPTime min (DCPTime a, DCPTime b);
338 DCPTime max (DCPTime a, DCPTime b);
339 ContentTime min (ContentTime a, ContentTime b);
340 ContentTime max (ContentTime a, ContentTime b);
341 std::string to_string (ContentTime t);
342 std::string to_string (DCPTime t);
343 std::string to_string (DCPTimePeriod p);
344
345 #endif