2 Copyright (C) 2012 Carl Hetherington <cth@carlh.net>
3 Copyright (C) 2000-2007 Paul Davis
5 This program is free software; you can redistribute it and/or modify
6 it under the terms of the GNU General Public License as published by
7 the Free Software Foundation; either version 2 of the License, or
8 (at your option) any later version.
10 This program is distributed in the hope that it will be useful,
11 but WITHOUT ANY WARRANTY; without even the implied warranty of
12 MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
13 GNU General Public License for more details.
15 You should have received a copy of the GNU General Public License
16 along with this program; if not, write to the Free Software
17 Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA.
21 /** @file src/lib/util.cc
22 * @brief Some utility functions and classes.
29 #ifdef DVDOMATIC_POSIX
33 #include <libssh/libssh.h>
35 #include <boost/algorithm/string.hpp>
36 #include <boost/bind.hpp>
37 #include <boost/lambda/lambda.hpp>
38 #include <boost/lexical_cast.hpp>
39 #include <boost/thread.hpp>
40 #include <boost/filesystem.hpp>
42 #include <openssl/md5.h>
43 #include <magick/MagickCore.h>
44 #include <magick/version.h>
45 #include <libdcp/version.h>
47 #include <libavcodec/avcodec.h>
48 #include <libavformat/avformat.h>
49 #include <libswscale/swscale.h>
50 #include <libavfilter/avfiltergraph.h>
51 #include <libpostproc/postprocess.h>
52 #include <libavutil/pixfmt.h>
55 #include "exceptions.h"
58 #include "dcp_content_type.h"
60 #include "sound_processor.h"
63 using namespace boost;
67 /** Convert some number of seconds to a string representation
68 * in hours, minutes and seconds.
71 * @return String of the form H:M:S (where H is hours, M
72 * is minutes and S is seconds).
75 seconds_to_hms (int s)
85 hms << setfill ('0') << m << ":";
87 hms << setfill ('0') << s;
92 /** @param s Number of seconds.
93 * @return String containing an approximate description of s (e.g. "about 2 hours")
96 seconds_to_approximate_hms (int s)
107 ap << (h + 1) << " hours";
119 ap << m << " minutes";
122 ap << s << " seconds";
128 #ifdef DVDOMATIC_POSIX
129 /** @param l Mangled C++ identifier.
130 * @return Demangled version.
135 string::size_type const b = l.find_first_of ("(");
136 if (b == string::npos) {
140 string::size_type const p = l.find_last_of ("+");
141 if (p == string::npos) {
149 string const fn = l.substr (b + 1, p - b - 1);
154 char* realname = abi::__cxa_demangle (fn.c_str(), 0, 0, &status);
159 } catch (std::exception) {
166 /** Write a stacktrace to an ostream.
167 * @param out Stream to write to.
168 * @param levels Number of levels to go up the call stack.
171 stacktrace (ostream& out, int levels)
178 size = backtrace (array, 200);
179 strings = backtrace_symbols (array, size);
182 for (i = 0; i < size && (levels == 0 || i < size_t(levels)); i++) {
183 out << " " << demangle (strings[i]) << endl;
191 /** @param v Version as used by FFmpeg.
192 * @return A string representation of v.
195 ffmpeg_version_to_string (int v)
198 s << ((v & 0xff0000) >> 16) << "." << ((v & 0xff00) >> 8) << "." << (v & 0xff);
202 /** Return a user-readable string summarising the versions of our dependencies */
204 dependency_version_summary ()
207 s << "libopenjpeg " << opj_version () << ", "
208 << "libavcodec " << ffmpeg_version_to_string (avcodec_version()) << ", "
209 << "libavfilter " << ffmpeg_version_to_string (avfilter_version()) << ", "
210 << "libavformat " << ffmpeg_version_to_string (avformat_version()) << ", "
211 << "libavutil " << ffmpeg_version_to_string (avutil_version()) << ", "
212 << "libpostproc " << ffmpeg_version_to_string (postproc_version()) << ", "
213 << "libswscale " << ffmpeg_version_to_string (swscale_version()) << ", "
214 << MagickVersion << ", "
215 << "libssh " << ssh_version (0) << ", "
216 << "libdcp " << libdcp::version << " git " << libdcp::git_commit;
222 seconds (struct timeval t)
224 return t.tv_sec + (double (t.tv_usec) / 1e6);
227 /** Call the required functions to set up DVD-o-matic's static arrays, etc.
228 * Must be called from the UI thread, if there is one.
233 Format::setup_formats ();
234 DCPContentType::setup_dcp_content_types ();
235 Scaler::setup_scalers ();
236 Filter::setup_filters ();
237 SoundProcessor::setup_sound_processors ();
239 ui_thread = this_thread::get_id ();
242 /** @param start Start position for the crop within the image.
243 * @param size Size of the cropped area.
244 * @return FFmpeg crop filter string.
247 crop_string (Position start, Size size)
250 s << "crop=" << size.width << ":" << size.height << ":" << start.x << ":" << start.y;
254 /** @param s A string.
255 * @return Parts of the string split at spaces, except when a space is within quotation marks.
258 split_at_spaces_considering_quotes (string s)
261 bool in_quotes = false;
263 for (string::size_type i = 0; i < s.length(); ++i) {
264 if (s[i] == ' ' && !in_quotes) {
267 } else if (s[i] == '"') {
268 in_quotes = !in_quotes;
279 md5_digest (void const * data, int size)
282 MD5_Init (&md5_context);
283 MD5_Update (&md5_context, data, size);
284 unsigned char digest[MD5_DIGEST_LENGTH];
285 MD5_Final (digest, &md5_context);
288 for (int i = 0; i < MD5_DIGEST_LENGTH; ++i) {
289 s << hex << setfill('0') << setw(2) << ((int) digest[i]);
295 /** @param file File name.
296 * @return MD5 digest of file's contents.
299 md5_digest (string file)
301 ifstream f (file.c_str(), ios::binary);
303 throw OpenFileError (file);
306 f.seekg (0, ios::end);
307 int bytes = f.tellg ();
308 f.seekg (0, ios::beg);
310 int const buffer_size = 64 * 1024;
311 char buffer[buffer_size];
314 MD5_Init (&md5_context);
316 int const t = min (bytes, buffer_size);
318 MD5_Update (&md5_context, buffer, t);
322 unsigned char digest[MD5_DIGEST_LENGTH];
323 MD5_Final (digest, &md5_context);
326 for (int i = 0; i < MD5_DIGEST_LENGTH; ++i) {
327 s << hex << setfill('0') << setw(2) << ((int) digest[i]);
333 static bool about_equal (float a, float b)
335 /* A film of F seconds at f FPS will be Ff frames;
336 Consider some delta FPS d, so if we run the same
337 film at (f + d) FPS it will last F(f + d) seconds.
339 Hence the difference in length over the length of the film will
340 be F(f + d) - Ff frames
341 = Ff + Fd - Ff frames
345 So if we accept a difference of 1 frame, ie 1/f seconds, we can
352 So for a 3hr film, ie F = 3 * 60 * 60 = 10800, the acceptable
353 FPS error is 1/F ~= 0.0001 ~= 10-e4
356 return (fabs (a - b) < 1e-4);
359 /** @param fps Arbitrary source frames-per-second value */
360 DCPFrameRate::DCPFrameRate (float fps)
361 : frames_per_second (rint (fps))
366 if (about_equal (fps, 50)) {
367 /* XXX: not sure about this; just run at 50?
370 frames_per_second = 25;
372 } else if (fps >= (27.5 / 2) && fps <= (32.5 / 2)) {
373 frames_per_second = 30;
375 } else if (fps >= (24.5 / 2) && fps <= (27.5 / 2)) {
376 frames_per_second = 25;
378 } else if (fps >= (20 / 2) && fps <= (24.5 / 2)) {
379 frames_per_second = 24;
384 /** @param An arbitrary sampling rate.
385 * @return The appropriate DCP-approved sampling rate (48kHz or 96kHz).
388 dcp_audio_sample_rate (int fs)
398 dcp_audio_channels (int f)
401 /* The source is mono, so to put the mono channel into
402 the centre we need to generate a 5.1 soundtrack.
411 bool operator== (Size const & a, Size const & b)
413 return (a.width == b.width && a.height == b.height);
416 bool operator!= (Size const & a, Size const & b)
421 bool operator== (Crop const & a, Crop const & b)
423 return (a.left == b.left && a.right == b.right && a.top == b.top && a.bottom == b.bottom);
426 bool operator!= (Crop const & a, Crop const & b)
431 /** @param index Colour LUT index.
432 * @return Human-readable name.
435 colour_lut_index_to_name (int index)
449 : _deadline (_io_service)
450 , _socket (_io_service)
453 _deadline.expires_at (posix_time::pos_infin);
460 if (_deadline.expires_at() <= asio::deadline_timer::traits_type::now ()) {
462 _deadline.expires_at (posix_time::pos_infin);
465 _deadline.async_wait (boost::bind (&Socket::check, this));
468 /** Blocking connect with timeout.
469 * @param endpoint End-point to connect to.
470 * @param timeout Time-out in seconds.
473 Socket::connect (asio::ip::basic_resolver_entry<asio::ip::tcp> const & endpoint, int timeout)
475 _deadline.expires_from_now (posix_time::seconds (timeout));
476 system::error_code ec = asio::error::would_block;
477 _socket.async_connect (endpoint, lambda::var(ec) = lambda::_1);
479 _io_service.run_one();
480 } while (ec == asio::error::would_block);
482 if (ec || !_socket.is_open ()) {
483 throw NetworkError ("connect timed out");
487 /** Blocking write with timeout.
488 * @param data Buffer to write.
489 * @param size Number of bytes to write.
490 * @param timeout Time-out, in seconds.
493 Socket::write (uint8_t const * data, int size, int timeout)
495 _deadline.expires_from_now (posix_time::seconds (timeout));
496 system::error_code ec = asio::error::would_block;
498 asio::async_write (_socket, asio::buffer (data, size), lambda::var(ec) = lambda::_1);
500 _io_service.run_one ();
501 } while (ec == asio::error::would_block);
504 throw NetworkError ("write timed out");
508 /** Blocking read with timeout.
509 * @param data Buffer to read to.
510 * @param size Number of bytes to read.
511 * @param timeout Time-out, in seconds.
514 Socket::read (uint8_t* data, int size, int timeout)
516 _deadline.expires_from_now (posix_time::seconds (timeout));
517 system::error_code ec = asio::error::would_block;
521 _socket.async_read_some (
522 asio::buffer (data, size),
523 (lambda::var(ec) = lambda::_1, lambda::var(amount_read) = lambda::_2)
527 _io_service.run_one ();
528 } while (ec == asio::error::would_block);
537 /** Mark some data as being `consumed', so that it will not be returned
539 * @param size Amount of data to consume, in bytes.
542 Socket::consume (int size)
544 assert (_buffer_data >= size);
546 _buffer_data -= size;
547 if (_buffer_data > 0) {
548 /* Shift still-valid data to the start of the buffer */
549 memmove (_buffer, _buffer + size, _buffer_data);
553 /** Read a definite amount of data from our socket, and mark
555 * @param data Where to put the data.
556 * @param size Number of bytes to read.
559 Socket::read_definite_and_consume (uint8_t* data, int size, int timeout)
561 int const from_buffer = min (_buffer_data, size);
562 if (from_buffer > 0) {
563 /* Get data from our buffer */
564 memcpy (data, _buffer, from_buffer);
565 consume (from_buffer);
566 /* Update our output state */
571 /* read() the rest */
573 int const n = read (data, size, timeout);
575 throw NetworkError ("could not read");
583 /** Read as much data as is available, up to some limit.
584 * @param data Where to put the data.
585 * @param size Maximum amount of data to read.
588 Socket::read_indefinite (uint8_t* data, int size, int timeout)
590 assert (size < int (sizeof (_buffer)));
592 /* Amount of extra data we need to read () */
593 int to_read = size - _buffer_data;
594 while (to_read > 0) {
595 /* read as much of it as we can (into our buffer) */
596 int const n = read (_buffer + _buffer_data, to_read, timeout);
598 throw NetworkError ("could not read");
605 assert (_buffer_data >= size);
607 /* copy data into the output buffer */
608 assert (size >= _buffer_data);
609 memcpy (data, _buffer, size);
612 /** @param other A Rect.
613 * @return The intersection of this with `other'.
616 Rect::intersection (Rect const & other) const
618 int const tx = max (x, other.x);
619 int const ty = max (y, other.y);
623 min (x + width, other.x + other.width) - tx,
624 min (y + height, other.y + other.height) - ty
628 /** Round a number up to the nearest multiple of another number.
630 * @param s Array of numbers to round, indexed by c.
631 * @param t Multiple to round to.
632 * @return Rounded number.
635 stride_round_up (int c, int const * stride, int t)
637 int const a = stride[c] + (t - 1);
642 stride_lookup (int c, int const * stride)
647 /** Read a sequence of key / value pairs from a text stream;
648 * the keys are the first words on the line, and the values are
649 * the remainder of the line following the key. Lines beginning
650 * with # are ignored.
651 * @param s Stream to read.
652 * @return key/value pairs.
654 multimap<string, string>
655 read_key_value (istream &s)
657 multimap<string, string> kv;
660 while (getline (s, line)) {
665 if (line[0] == '#') {
669 if (line[line.size() - 1] == '\r') {
670 line = line.substr (0, line.size() - 1);
673 size_t const s = line.find (' ');
674 if (s == string::npos) {
678 kv.insert (make_pair (line.substr (0, s), line.substr (s + 1)));
685 get_required_string (multimap<string, string> const & kv, string k)
687 if (kv.count (k) > 1) {
688 throw StringError ("unexpected multiple keys in key-value set");
691 multimap<string, string>::const_iterator i = kv.find (k);
693 if (i == kv.end ()) {
694 throw StringError (String::compose ("missing key %1 in key-value set", k));
701 get_required_int (multimap<string, string> const & kv, string k)
703 string const v = get_required_string (kv, k);
704 return lexical_cast<int> (v);
708 get_required_float (multimap<string, string> const & kv, string k)
710 string const v = get_required_string (kv, k);
711 return lexical_cast<float> (v);
715 get_optional_string (multimap<string, string> const & kv, string k)
717 if (kv.count (k) > 1) {
718 throw StringError ("unexpected multiple keys in key-value set");
721 multimap<string, string>::const_iterator i = kv.find (k);
722 if (i == kv.end ()) {
730 get_optional_int (multimap<string, string> const & kv, string k)
732 if (kv.count (k) > 1) {
733 throw StringError ("unexpected multiple keys in key-value set");
736 multimap<string, string>::const_iterator i = kv.find (k);
737 if (i == kv.end ()) {
741 return lexical_cast<int> (i->second);
744 /** Construct an AudioBuffers. Audio data is undefined after this constructor.
745 * @param channels Number of channels.
746 * @param frames Number of frames to reserve space for.
748 AudioBuffers::AudioBuffers (int channels, int frames)
749 : _channels (channels)
751 , _allocated_frames (frames)
753 _data = new float*[_channels];
754 for (int i = 0; i < _channels; ++i) {
755 _data[i] = new float[frames];
759 /** Copy constructor.
760 * @param other Other AudioBuffers; data is copied.
762 AudioBuffers::AudioBuffers (AudioBuffers const & other)
763 : _channels (other._channels)
764 , _frames (other._frames)
765 , _allocated_frames (other._frames)
767 _data = new float*[_channels];
768 for (int i = 0; i < _channels; ++i) {
769 _data[i] = new float[_frames];
770 memcpy (_data[i], other._data[i], _frames * sizeof (float));
774 /** AudioBuffers destructor */
775 AudioBuffers::~AudioBuffers ()
777 for (int i = 0; i < _channels; ++i) {
784 /** @param c Channel index.
785 * @return Buffer for this channel.
788 AudioBuffers::data (int c) const
790 assert (c >= 0 && c < _channels);
794 /** Set the number of frames that these AudioBuffers will report themselves
796 * @param f Frames; must be less than or equal to the number of allocated frames.
799 AudioBuffers::set_frames (int f)
801 assert (f <= _allocated_frames);
805 /** Make all samples on all channels silent */
807 AudioBuffers::make_silent ()
809 for (int i = 0; i < _channels; ++i) {
814 /** Make all samples on a given channel silent.
818 AudioBuffers::make_silent (int c)
820 assert (c >= 0 && c < _channels);
822 for (int i = 0; i < _frames; ++i) {
827 /** Copy data from another AudioBuffers to this one. All channels are copied.
828 * @param from AudioBuffers to copy from; must have the same number of channels as this.
829 * @param frames_to_copy Number of frames to copy.
830 * @param read_offset Offset to read from in `from'.
831 * @param write_offset Offset to write to in `to'.
834 AudioBuffers::copy_from (AudioBuffers* from, int frames_to_copy, int read_offset, int write_offset)
836 assert (from->channels() == channels());
839 assert (read_offset >= 0 && (read_offset + frames_to_copy) <= from->_allocated_frames);
840 assert (write_offset >= 0 && (write_offset + frames_to_copy) <= _allocated_frames);
842 for (int i = 0; i < _channels; ++i) {
843 memcpy (_data[i] + write_offset, from->_data[i] + read_offset, frames_to_copy * sizeof(float));
847 /** Move audio data around.
848 * @param from Offset to move from.
849 * @param to Offset to move to.
850 * @param frames Number of frames to move.
854 AudioBuffers::move (int from, int to, int frames)
861 assert (from < _frames);
863 assert (to < _frames);
865 assert (frames <= _frames);
866 assert ((from + frames) <= _frames);
867 assert ((to + frames) <= _frames);
869 for (int i = 0; i < _channels; ++i) {
870 memmove (_data[i] + to, _data[i] + from, frames * sizeof(float));
874 /** Trip an assert if the caller is not in the UI thread */
878 assert (this_thread::get_id() == ui_thread);
881 /** @param v Source video frame.
882 * @param audio_sample_rate Source audio sample rate.
883 * @param frames_per_second Number of video frames per second.
884 * @return Equivalent number of audio frames for `v'.
887 video_frames_to_audio_frames (SourceFrame v, float audio_sample_rate, float frames_per_second)
889 return ((int64_t) v * audio_sample_rate / frames_per_second);
892 /** @param f Filename.
893 * @return true if this file is a still image, false if it is something else.
896 still_image_file (string f)
898 #if BOOST_FILESYSTEM_VERSION == 3
899 string ext = boost::filesystem::path(f).extension().string();
901 string ext = boost::filesystem::path(f).extension();
904 transform (ext.begin(), ext.end(), ext.begin(), ::tolower);
906 return (ext == ".tif" || ext == ".tiff" || ext == ".jpg" || ext == ".jpeg" || ext == ".png");
909 /** @return A pair containing CPU model name and the number of processors */
913 pair<string, int> info;
916 #ifdef DVDOMATIC_POSIX
917 ifstream f ("/proc/cpuinfo");
921 if (boost::algorithm::starts_with (l, "model name")) {
922 string::size_type const c = l.find (':');
923 if (c != string::npos) {
924 info.first = l.substr (c + 2);
926 } else if (boost::algorithm::starts_with (l, "processor")) {