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.
30 #ifdef DVDOMATIC_POSIX
34 #include <libssh/libssh.h>
36 #include <boost/algorithm/string.hpp>
37 #include <boost/bind.hpp>
38 #include <boost/lambda/lambda.hpp>
39 #include <boost/lexical_cast.hpp>
40 #include <boost/thread.hpp>
41 #include <boost/filesystem.hpp>
43 #include <openssl/md5.h>
44 #include <magick/MagickCore.h>
45 #include <magick/version.h>
46 #include <libdcp/version.h>
48 #include <libavcodec/avcodec.h>
49 #include <libavformat/avformat.h>
50 #include <libswscale/swscale.h>
51 #include <libavfilter/avfiltergraph.h>
52 #include <libpostproc/postprocess.h>
53 #include <libavutil/pixfmt.h>
56 #include "exceptions.h"
59 #include "dcp_content_type.h"
61 #include "sound_processor.h"
67 using namespace boost;
72 /** Convert some number of seconds to a string representation
73 * in hours, minutes and seconds.
76 * @return String of the form H:M:S (where H is hours, M
77 * is minutes and S is seconds).
80 seconds_to_hms (int s)
90 hms << setfill ('0') << m << N_(":");
92 hms << setfill ('0') << s;
97 /** @param s Number of seconds.
98 * @return String containing an approximate description of s (e.g. "about 2 hours")
101 seconds_to_approximate_hms (int s)
112 ap << (h + 1) << N_(" ") << _("hours");
115 ap << N_("1 ") << _("hour");
117 ap << h << N_(" ") << _("hours");
122 ap << N_("1 ") << _("minute");
124 ap << m << N_(" ") << _("minutes");
127 ap << s << N_(" ") << _("seconds");
133 #ifdef DVDOMATIC_POSIX
134 /** @param l Mangled C++ identifier.
135 * @return Demangled version.
140 string::size_type const b = l.find_first_of (N_("("));
141 if (b == string::npos) {
145 string::size_type const p = l.find_last_of (N_("+"));
146 if (p == string::npos) {
154 string const fn = l.substr (b + 1, p - b - 1);
159 char* realname = abi::__cxa_demangle (fn.c_str(), 0, 0, &status);
164 } catch (std::exception) {
171 /** Write a stacktrace to an ostream.
172 * @param out Stream to write to.
173 * @param levels Number of levels to go up the call stack.
176 stacktrace (ostream& out, int levels)
183 size = backtrace (array, 200);
184 strings = backtrace_symbols (array, size);
187 for (i = 0; i < size && (levels == 0 || i < size_t(levels)); i++) {
188 out << N_(" ") << demangle (strings[i]) << endl;
196 /** @param v Version as used by FFmpeg.
197 * @return A string representation of v.
200 ffmpeg_version_to_string (int v)
203 s << ((v & 0xff0000) >> 16) << N_(".") << ((v & 0xff00) >> 8) << N_(".") << (v & 0xff);
207 /** Return a user-readable string summarising the versions of our dependencies */
209 dependency_version_summary ()
212 s << N_("libopenjpeg ") << opj_version () << N_(", ")
213 << N_("libavcodec ") << ffmpeg_version_to_string (avcodec_version()) << N_(", ")
214 << N_("libavfilter ") << ffmpeg_version_to_string (avfilter_version()) << N_(", ")
215 << N_("libavformat ") << ffmpeg_version_to_string (avformat_version()) << N_(", ")
216 << N_("libavutil ") << ffmpeg_version_to_string (avutil_version()) << N_(", ")
217 << N_("libpostproc ") << ffmpeg_version_to_string (postproc_version()) << N_(", ")
218 << N_("libswscale ") << ffmpeg_version_to_string (swscale_version()) << N_(", ")
219 << MagickVersion << N_(", ")
220 << N_("libssh ") << ssh_version (0) << N_(", ")
221 << N_("libdcp ") << libdcp::version << N_(" git ") << libdcp::git_commit;
227 seconds (struct timeval t)
229 return t.tv_sec + (double (t.tv_usec) / 1e6);
232 /** Call the required functions to set up DVD-o-matic's static arrays, etc.
233 * Must be called from the UI thread, if there is one.
238 avfilter_register_all ();
240 Format::setup_formats ();
241 DCPContentType::setup_dcp_content_types ();
242 Scaler::setup_scalers ();
243 Filter::setup_filters ();
244 SoundProcessor::setup_sound_processors ();
246 ui_thread = this_thread::get_id ();
249 #ifdef DVDOMATIC_WINDOWS
250 boost::filesystem::path
254 GetModuleFileName (0, buffer, 512 * sizeof(wchar_t));
255 boost::filesystem::path p (buffer);
256 p = p.parent_path ();
257 p = p.parent_path ();
264 dvdomatic_setup_i18n (string lang)
266 #ifdef DVDOMATIC_POSIX
270 if (!lang.empty ()) {
271 /* Override our environment language; this is essential on
275 snprintf (cmd, sizeof(cmd), "LANGUAGE=%s", lang.c_str ());
277 snprintf (cmd, sizeof(cmd), "LANG=%s", lang.c_str ());
281 setlocale (LC_ALL, "");
282 textdomain ("libdvdomatic");
284 #ifdef DVDOMATIC_WINDOWS
285 bindtextdomain ("libdvdomatic", mo_path().string().c_str());
288 #ifdef DVDOMATIC_POSIX
289 bindtextdomain ("libdvdomatic", POSIX_LOCALE_PREFIX);
293 /** @param start Start position for the crop within the image.
294 * @param size Size of the cropped area.
295 * @return FFmpeg crop filter string.
298 crop_string (Position start, libdcp::Size size)
301 s << N_("crop=") << size.width << N_(":") << size.height << N_(":") << start.x << N_(":") << start.y;
305 /** @param s A string.
306 * @return Parts of the string split at spaces, except when a space is within quotation marks.
309 split_at_spaces_considering_quotes (string s)
312 bool in_quotes = false;
314 for (string::size_type i = 0; i < s.length(); ++i) {
315 if (s[i] == ' ' && !in_quotes) {
318 } else if (s[i] == '"') {
319 in_quotes = !in_quotes;
330 md5_digest (void const * data, int size)
333 MD5_Init (&md5_context);
334 MD5_Update (&md5_context, data, size);
335 unsigned char digest[MD5_DIGEST_LENGTH];
336 MD5_Final (digest, &md5_context);
339 for (int i = 0; i < MD5_DIGEST_LENGTH; ++i) {
340 s << hex << setfill('0') << setw(2) << ((int) digest[i]);
346 /** @param file File name.
347 * @return MD5 digest of file's contents.
350 md5_digest (boost::filesystem::path file)
352 ifstream f (file.string().c_str(), ios::binary);
354 throw OpenFileError (file.string());
357 f.seekg (0, ios::end);
358 int bytes = f.tellg ();
359 f.seekg (0, ios::beg);
361 int const buffer_size = 64 * 1024;
362 char buffer[buffer_size];
365 MD5_Init (&md5_context);
367 int const t = min (bytes, buffer_size);
369 MD5_Update (&md5_context, buffer, t);
373 unsigned char digest[MD5_DIGEST_LENGTH];
374 MD5_Final (digest, &md5_context);
377 for (int i = 0; i < MD5_DIGEST_LENGTH; ++i) {
378 s << hex << setfill('0') << setw(2) << ((int) digest[i]);
385 about_equal (float a, float b)
387 /* A film of F seconds at f FPS will be Ff frames;
388 Consider some delta FPS d, so if we run the same
389 film at (f + d) FPS it will last F(f + d) seconds.
391 Hence the difference in length over the length of the film will
392 be F(f + d) - Ff frames
393 = Ff + Fd - Ff frames
397 So if we accept a difference of 1 frame, ie 1/f seconds, we can
404 So for a 3hr film, ie F = 3 * 60 * 60 = 10800, the acceptable
405 FPS error is 1/F ~= 0.0001 ~= 10-e4
408 return (fabs (a - b) < 1e-4);
411 class FrameRateCandidate
414 FrameRateCandidate (float source_, int dcp_)
424 best_dcp_frame_rate (float source_fps)
426 list<int> const allowed_dcp_frame_rates = Config::instance()->allowed_dcp_frame_rates ();
428 /* Work out what rates we could manage, including those achieved by using skip / repeat. */
429 list<FrameRateCandidate> candidates;
431 /* Start with the ones without skip / repeat so they will get matched in preference to skipped/repeated ones */
432 for (list<int>::const_iterator i = allowed_dcp_frame_rates.begin(); i != allowed_dcp_frame_rates.end(); ++i) {
433 candidates.push_back (FrameRateCandidate (*i, *i));
436 /* Then the skip/repeat ones */
437 for (list<int>::const_iterator i = allowed_dcp_frame_rates.begin(); i != allowed_dcp_frame_rates.end(); ++i) {
438 candidates.push_back (FrameRateCandidate (float (*i) / 2, *i));
439 candidates.push_back (FrameRateCandidate (float (*i) * 2, *i));
442 /* Pick the best one, bailing early if we hit an exact match */
443 float error = numeric_limits<float>::max ();
444 boost::optional<FrameRateCandidate> best;
445 list<FrameRateCandidate>::iterator i = candidates.begin();
446 while (i != candidates.end()) {
448 if (about_equal (i->source, source_fps)) {
453 float const e = fabs (i->source - source_fps);
466 /** @param An arbitrary sampling rate.
467 * @return The appropriate DCP-approved sampling rate (48kHz or 96kHz).
470 dcp_audio_sample_rate (int fs)
479 /** @param index Colour LUT index.
480 * @return Human-readable name.
483 colour_lut_index_to_name (int index)
496 Socket::Socket (int timeout)
497 : _deadline (_io_service)
498 , _socket (_io_service)
501 _deadline.expires_at (posix_time::pos_infin);
508 if (_deadline.expires_at() <= asio::deadline_timer::traits_type::now ()) {
510 _deadline.expires_at (posix_time::pos_infin);
513 _deadline.async_wait (boost::bind (&Socket::check, this));
516 /** Blocking connect.
517 * @param endpoint End-point to connect to.
520 Socket::connect (asio::ip::basic_resolver_entry<asio::ip::tcp> const & endpoint)
522 _deadline.expires_from_now (posix_time::seconds (_timeout));
523 system::error_code ec = asio::error::would_block;
524 _socket.async_connect (endpoint, lambda::var(ec) = lambda::_1);
526 _io_service.run_one();
527 } while (ec == asio::error::would_block);
529 if (ec || !_socket.is_open ()) {
530 throw NetworkError (_("connect timed out"));
535 * @param data Buffer to write.
536 * @param size Number of bytes to write.
539 Socket::write (uint8_t const * data, int size)
541 _deadline.expires_from_now (posix_time::seconds (_timeout));
542 system::error_code ec = asio::error::would_block;
544 asio::async_write (_socket, asio::buffer (data, size), lambda::var(ec) = lambda::_1);
547 _io_service.run_one ();
548 } while (ec == asio::error::would_block);
551 throw NetworkError (ec.message ());
556 Socket::write (uint32_t v)
559 write (reinterpret_cast<uint8_t*> (&v), 4);
563 * @param data Buffer to read to.
564 * @param size Number of bytes to read.
567 Socket::read (uint8_t* data, int size)
569 _deadline.expires_from_now (posix_time::seconds (_timeout));
570 system::error_code ec = asio::error::would_block;
572 asio::async_read (_socket, asio::buffer (data, size), lambda::var(ec) = lambda::_1);
575 _io_service.run_one ();
576 } while (ec == asio::error::would_block);
579 throw NetworkError (ec.message ());
584 Socket::read_uint32 ()
587 read (reinterpret_cast<uint8_t *> (&v), 4);
591 /** @param other A Rect.
592 * @return The intersection of this with `other'.
595 Rect::intersection (Rect const & other) const
597 int const tx = max (x, other.x);
598 int const ty = max (y, other.y);
602 min (x + width, other.x + other.width) - tx,
603 min (y + height, other.y + other.height) - ty
607 /** Round a number up to the nearest multiple of another number.
609 * @param s Array of numbers to round, indexed by c.
610 * @param t Multiple to round to.
611 * @return Rounded number.
614 stride_round_up (int c, int const * stride, int t)
616 int const a = stride[c] + (t - 1);
621 stride_lookup (int c, int const * stride)
626 /** Read a sequence of key / value pairs from a text stream;
627 * the keys are the first words on the line, and the values are
628 * the remainder of the line following the key. Lines beginning
629 * with # are ignored.
630 * @param s Stream to read.
631 * @return key/value pairs.
633 multimap<string, string>
634 read_key_value (istream &s)
636 multimap<string, string> kv;
639 while (getline (s, line)) {
644 if (line[0] == '#') {
648 if (line[line.size() - 1] == '\r') {
649 line = line.substr (0, line.size() - 1);
652 size_t const s = line.find (' ');
653 if (s == string::npos) {
657 kv.insert (make_pair (line.substr (0, s), line.substr (s + 1)));
664 get_required_string (multimap<string, string> const & kv, string k)
666 if (kv.count (k) > 1) {
667 throw StringError (N_("unexpected multiple keys in key-value set"));
670 multimap<string, string>::const_iterator i = kv.find (k);
672 if (i == kv.end ()) {
673 throw StringError (String::compose (_("missing key %1 in key-value set"), k));
680 get_required_int (multimap<string, string> const & kv, string k)
682 string const v = get_required_string (kv, k);
683 return lexical_cast<int> (v);
687 get_required_float (multimap<string, string> const & kv, string k)
689 string const v = get_required_string (kv, k);
690 return lexical_cast<float> (v);
694 get_optional_string (multimap<string, string> const & kv, string k)
696 if (kv.count (k) > 1) {
697 throw StringError (N_("unexpected multiple keys in key-value set"));
700 multimap<string, string>::const_iterator i = kv.find (k);
701 if (i == kv.end ()) {
709 get_optional_int (multimap<string, string> const & kv, string k)
711 if (kv.count (k) > 1) {
712 throw StringError (N_("unexpected multiple keys in key-value set"));
715 multimap<string, string>::const_iterator i = kv.find (k);
716 if (i == kv.end ()) {
720 return lexical_cast<int> (i->second);
723 /** Construct an AudioBuffers. Audio data is undefined after this constructor.
724 * @param channels Number of channels.
725 * @param frames Number of frames to reserve space for.
727 AudioBuffers::AudioBuffers (int channels, int frames)
728 : _channels (channels)
730 , _allocated_frames (frames)
732 _data = new float*[_channels];
733 for (int i = 0; i < _channels; ++i) {
734 _data[i] = new float[frames];
738 /** Copy constructor.
739 * @param other Other AudioBuffers; data is copied.
741 AudioBuffers::AudioBuffers (AudioBuffers const & other)
742 : _channels (other._channels)
743 , _frames (other._frames)
744 , _allocated_frames (other._frames)
746 _data = new float*[_channels];
747 for (int i = 0; i < _channels; ++i) {
748 _data[i] = new float[_frames];
749 memcpy (_data[i], other._data[i], _frames * sizeof (float));
753 /** AudioBuffers destructor */
754 AudioBuffers::~AudioBuffers ()
756 for (int i = 0; i < _channels; ++i) {
763 /** @param c Channel index.
764 * @return Buffer for this channel.
767 AudioBuffers::data (int c) const
769 assert (c >= 0 && c < _channels);
773 /** Set the number of frames that these AudioBuffers will report themselves
775 * @param f Frames; must be less than or equal to the number of allocated frames.
778 AudioBuffers::set_frames (int f)
780 assert (f <= _allocated_frames);
784 /** Make all samples on all channels silent */
786 AudioBuffers::make_silent ()
788 for (int i = 0; i < _channels; ++i) {
793 /** Make all samples on a given channel silent.
797 AudioBuffers::make_silent (int c)
799 assert (c >= 0 && c < _channels);
801 for (int i = 0; i < _frames; ++i) {
806 /** Copy data from another AudioBuffers to this one. All channels are copied.
807 * @param from AudioBuffers to copy from; must have the same number of channels as this.
808 * @param frames_to_copy Number of frames to copy.
809 * @param read_offset Offset to read from in `from'.
810 * @param write_offset Offset to write to in `to'.
813 AudioBuffers::copy_from (AudioBuffers* from, int frames_to_copy, int read_offset, int write_offset)
815 assert (from->channels() == channels());
818 assert (read_offset >= 0 && (read_offset + frames_to_copy) <= from->_allocated_frames);
819 assert (write_offset >= 0 && (write_offset + frames_to_copy) <= _allocated_frames);
821 for (int i = 0; i < _channels; ++i) {
822 memcpy (_data[i] + write_offset, from->_data[i] + read_offset, frames_to_copy * sizeof(float));
826 /** Move audio data around.
827 * @param from Offset to move from.
828 * @param to Offset to move to.
829 * @param frames Number of frames to move.
833 AudioBuffers::move (int from, int to, int frames)
840 assert (from < _frames);
842 assert (to < _frames);
844 assert (frames <= _frames);
845 assert ((from + frames) <= _frames);
846 assert ((to + frames) <= _frames);
848 for (int i = 0; i < _channels; ++i) {
849 memmove (_data[i] + to, _data[i] + from, frames * sizeof(float));
853 /** Trip an assert if the caller is not in the UI thread */
857 assert (this_thread::get_id() == ui_thread);
860 /** @param v Content video frame.
861 * @param audio_sample_rate Source audio sample rate.
862 * @param frames_per_second Number of video frames per second.
863 * @return Equivalent number of audio frames for `v'.
866 video_frames_to_audio_frames (ContentVideoFrame v, float audio_sample_rate, float frames_per_second)
868 return ((int64_t) v * audio_sample_rate / frames_per_second);
871 /** @return A pair containing CPU model name and the number of processors */
875 pair<string, int> info;
878 #ifdef DVDOMATIC_POSIX
879 ifstream f (N_("/proc/cpuinfo"));
883 if (boost::algorithm::starts_with (l, N_("model name"))) {
884 string::size_type const c = l.find (':');
885 if (c != string::npos) {
886 info.first = l.substr (c + 2);
888 } else if (boost::algorithm::starts_with (l, N_("processor"))) {
898 audio_channel_name (int c)
900 assert (MAX_AUDIO_CHANNELS == 6);
902 /* TRANSLATORS: these are the names of audio channels; Lfe (sub) is the low-frequency
903 enhancement channel (sub-woofer)./
905 string const channels[] = {
917 AudioMapping::AudioMapping (int c)
918 : _source_channels (c)
923 optional<libdcp::Channel>
924 AudioMapping::source_to_dcp (int c) const
926 if (c >= _source_channels) {
927 return optional<libdcp::Channel> ();
930 if (_source_channels == 1) {
931 /* mono sources to centre */
932 return libdcp::CENTRE;
935 return static_cast<libdcp::Channel> (c);
939 AudioMapping::dcp_to_source (libdcp::Channel c) const
941 if (_source_channels == 1) {
942 if (c == libdcp::CENTRE) {
945 return optional<int> ();
949 if (static_cast<int> (c) >= _source_channels) {
950 return optional<int> ();
953 return static_cast<int> (c);
957 AudioMapping::dcp_channels () const
959 if (_source_channels == 1) {
960 /* The source is mono, so to put the mono channel into
961 the centre we need to generate a 5.1 soundtrack.
966 return _source_channels;
969 FrameRateConversion::FrameRateConversion (float source, int dcp)
972 , change_speed (false)
974 if (fabs (source / 2.0 - dcp) < (fabs (source - dcp))) {
976 } else if (fabs (source * 2 - dcp) < fabs (source - dcp)) {
980 change_speed = !about_equal (source * factor(), dcp);
982 if (!skip && !repeat && !change_speed) {
983 description = _("DCP and source have the same rate.\n");
986 description = _("DCP will use every other frame of the source.\n");
988 description = _("Each source frame will be doubled in the DCP.\n");
992 float const pc = dcp * 100 / (source * factor());
993 description += String::compose (_("DCP will run at %1%% of the source speed."), pc);