git.sesse.net Git - casparcg/blob - accelerator/cpu/image/image_mixer.cpp

   1 /*
   2 * Copyright (c) 2011 Sveriges Television AB <info@casparcg.com>
   3 *
   4 * This file is part of CasparCG (www.casparcg.com).
   5 *
   6 * CasparCG 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 3 of the License, or
   9 * (at your option) any later version.
  10 *
  11 * CasparCG 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 CasparCG. If not, see <http://www.gnu.org/licenses/>.
  18 *
  19 * Author: Robert Nagy, ronag89@gmail.com
  20 */
  21
  22 #include "../../StdAfx.h"
  23
  24 #include "image_mixer.h"
  25
  26 #include "../util/xmm.h"
  27
  28 #include <common/assert.h>
  29 #include <common/gl/gl_check.h>
  30 #include <common/future.h>
  31 #include <common/array.h>
  32
  33 #include <core/frame/frame.h>
  34 #include <core/frame/frame_transform.h>
  35 #include <core/frame/pixel_format.h>
  36 #include <core/video_format.h>
  37
  38 #include <modules/ffmpeg/producer/util/util.h>
  39
  40 #include <GL/glew.h>
  41
  42 #include <tbb/parallel_for.h>
  43 #include <tbb/parallel_for_each.h>
  44 #include <tbb/concurrent_queue.h>
  45
  46 #include <boost/range/algorithm_ext/erase.hpp>
  47 #include <boost/thread/future.hpp>
  48
  49 #include <algorithm>
  50 #include <cstdint>
  51 #include <vector>
  52 #include <set>
  53 #include <array>
  54
  55 #if defined(_MSC_VER)
  56 #pragma warning (push)
  57 #pragma warning (disable : 4244)
  58 #endif
  59 extern "C"
  60 {
  61         #include <libswscale/swscale.h>
  62         #include <libavcodec/avcodec.h>
  63         #include <libavformat/avformat.h>
  64 }
  65 #if defined(_MSC_VER)
  66 #pragma warning (pop)
  67 #endif
  68
  69 namespace caspar { namespace accelerator { namespace cpu {
  70
  71 struct item
  72 {
  73         core::pixel_format_desc                 pix_desc        = core::pixel_format::invalid;
  74         std::array<const uint8_t*, 4>   data;
  75         core::image_transform                   transform;
  76
  77         item()
  78         {
  79                 data.fill(0);
  80         }
  81 };
  82
  83 bool operator==(const item& lhs, const item& rhs)
  84 {
  85         return lhs.data == rhs.data && lhs.transform == rhs.transform;
  86 }
  87
  88 bool operator!=(const item& lhs, const item& rhs)
  89 {
  90         return !(lhs == rhs);
  91 }
  92
  93 // 100% accurate blending with correct rounding.
  94 inline xmm::s8_x blend(xmm::s8_x d, xmm::s8_x s)
  95 {
  96         using namespace xmm;
  97
  98         // C(S, D) = S + D - (((T >> 8) + T) >> 8);
  99         // T(S, D) = S * D[A] + 0x80
 100
 101         auto aaaa   = s8_x::shuffle(d, s8_x(15, 15, 15, 15, 11, 11, 11, 11, 7, 7, 7, 7, 3, 3, 3, 3));
 102         d                       = s8_x(u8_x::min(u8_x(d), u8_x(aaaa))); // Overflow guard. Some source files have color values which incorrectly exceed pre-multiplied alpha values, e.g. red(255) > alpha(254).
 103
 104         auto xaxa       = s16_x(aaaa) >> 8;
 105
 106         auto t1         = s16_x::multiply_low(s16_x(s) & 0x00FF, xaxa) + 0x80;
 107         auto t2         = s16_x::multiply_low(s16_x(s) >> 8    , xaxa) + 0x80;
 108
 109         auto xyxy       = s8_x(((t1 >> 8) + t1) >> 8);
 110         auto yxyx       = s8_x((t2 >> 8) + t2);
 111         auto argb   = s8_x::blend(xyxy, yxyx, s8_x(-1, 0, -1, 0));
 112
 113         return s8_x(s) + (d - argb);
 114 }
 115
 116 template<typename temporal, typename alignment>
 117 static void kernel(uint8_t* dest, const uint8_t* source, size_t count)
 118 {
 119         using namespace xmm;
 120
 121         for(auto n = 0; n < count; n += 32)
 122         {
 123                 auto s0 = s8_x::load<temporal_tag, alignment>(dest+n+0);
 124                 auto s1 = s8_x::load<temporal_tag, alignment>(dest+n+16);
 125
 126                 auto d0 = s8_x::load<temporal_tag, alignment>(source+n+0);
 127                 auto d1 = s8_x::load<temporal_tag, alignment>(source+n+16);
 128
 129                 auto argb0 = blend(d0, s0);
 130                 auto argb1 = blend(d1, s1);
 131
 132                 s8_x::store<temporal, alignment>(argb0, dest+n+0 );
 133                 s8_x::store<temporal, alignment>(argb1, dest+n+16);
 134         }
 135 }
 136
 137 template<typename temporal>
 138 static void kernel(uint8_t* dest, const uint8_t* source, size_t count)
 139 {
 140         using namespace xmm;
 141
 142         if(reinterpret_cast<std::uint64_t>(dest) % 16 != 0 || reinterpret_cast<std::uint64_t>(source) % 16 != 0)
 143                 kernel<temporal_tag, unaligned_tag>(dest, source, count);
 144         else
 145                 kernel<temporal_tag, aligned_tag>(dest, source, count);
 146 }
 147
 148 class image_renderer
 149 {
 150         tbb::concurrent_unordered_map<int64_t, tbb::concurrent_bounded_queue<std::shared_ptr<SwsContext>>>      sws_devices_;
 151         tbb::concurrent_bounded_queue<spl::shared_ptr<buffer>>                                                                                          temp_buffers_;
 152         core::video_format_desc                                                                                                                                                         format_desc_;
 153 public:
 154         std::future<array<const std::uint8_t>> operator()(std::vector<item> items, const core::video_format_desc& format_desc)
 155         {
 156                 if (format_desc != format_desc_)
 157                 {
 158                         format_desc_ = format_desc;
 159                         sws_devices_.clear();
 160                 }
 161
 162                 convert(items, format_desc.width, format_desc.height);
 163
 164                 // Remove first field stills.
 165                 boost::range::remove_erase_if(items, [&](const item& item)
 166                 {
 167                         return item.transform.is_still && item.transform.field_mode == format_desc.field_mode; // only us last field for stills.
 168                 });
 169
 170                 // Stills are progressive
 171                 for (auto& item : items)
 172                 {
 173                         if(item.transform.is_still)
 174                                 item.transform.field_mode = core::field_mode::progressive;
 175                 }
 176
 177                 auto result = spl::make_shared<buffer>(format_desc.size, 0);
 178                 if(format_desc.field_mode != core::field_mode::progressive)
 179                 {
 180                         draw(items, result->data(), format_desc.width, format_desc.height, core::field_mode::upper);
 181                         draw(items, result->data(), format_desc.width, format_desc.height, core::field_mode::lower);
 182                 }
 183                 else
 184                 {
 185                         draw(items, result->data(), format_desc.width, format_desc.height,  core::field_mode::progressive);
 186                 }
 187
 188                 temp_buffers_.clear();
 189
 190                 return make_ready_future(array<const std::uint8_t>(result->data(), format_desc.size, true, result));
 191         }
 192
 193 private:
 194
 195         void draw(std::vector<item> items, uint8_t* dest, std::size_t width, std::size_t height, core::field_mode field_mode)
 196         {
 197                 for (auto& item : items)
 198                         item.transform.field_mode &= field_mode;
 199
 200                 // Remove empty items.
 201                 boost::range::remove_erase_if(items, [&](const item& item)
 202                 {
 203                         return item.transform.field_mode == core::field_mode::empty;
 204                 });
 205
 206                 if(items.empty())
 207                         return;
 208
 209                 auto start = field_mode == core::field_mode::lower ? 1 : 0;
 210                 auto step  = field_mode == core::field_mode::progressive ? 1 : 2;
 211
 212                 // TODO: Add support for fill translations.
 213                 // TODO: Add support for mask rect.
 214                 // TODO: Add support for opacity.
 215                 // TODO: Add support for mix transition.
 216                 // TODO: Add support for push transition.
 217                 // TODO: Add support for wipe transition.
 218                 // TODO: Add support for slide transition.
 219                 tbb::parallel_for(tbb::blocked_range<std::size_t>(0, height/step), [&](const tbb::blocked_range<std::size_t>& r)
 220                 {
 221                         for(auto i = r.begin(); i != r.end(); ++i)
 222                         {
 223                                 auto y = i*step+start;
 224
 225                                 for(std::size_t n = 0; n < items.size()-1; ++n)
 226                                         kernel<xmm::temporal_tag>(dest + y*width*4, items[n].data.at(0) + y*width*4, width*4);
 227
 228                                 std::size_t n = items.size()-1;
 229                                 kernel<xmm::nontemporal_tag>(dest + y*width*4, items[n].data.at(0) + y*width*4, width*4);
 230                         }
 231
 232                         _mm_mfence();
 233                 });
 234         }
 235
 236         void convert(std::vector<item>& source_items, int width, int height)
 237         {
 238                 std::set<std::array<const uint8_t*, 4>> buffers;
 239
 240                 for (auto& item : source_items)
 241                         buffers.insert(item.data);
 242
 243                 auto dest_items = source_items;
 244
 245                 tbb::parallel_for_each(buffers.begin(), buffers.end(), [&](const std::array<const uint8_t*, 4>& data)
 246                 {
 247                         auto pix_desc = std::find_if(source_items.begin(), source_items.end(), [&](const item& item){return item.data == data;})->pix_desc;
 248
 249                         if(pix_desc.format == core::pixel_format::bgra &&
 250                                 pix_desc.planes.at(0).width == width &&
 251                                 pix_desc.planes.at(0).height == height)
 252                                 return;
 253
 254                         std::array<uint8_t*, 4> data2 = {};
 255                         for(std::size_t n = 0; n < data.size(); ++n)
 256                                 data2.at(n) = const_cast<uint8_t*>(data[n]);
 257
 258                         auto input_av_frame = ffmpeg::make_av_frame(data2, pix_desc);
 259
 260
 261                         int64_t key = ((static_cast<int64_t>(input_av_frame->width)      << 32) & 0xFFFF00000000) |
 262                                                   ((static_cast<int64_t>(input_av_frame->height) << 16) & 0xFFFF0000) |
 263                                                   ((static_cast<int64_t>(input_av_frame->format) <<  8) & 0xFF00);
 264
 265                         auto& pool = sws_devices_[key];
 266
 267                         std::shared_ptr<SwsContext> sws_device;
 268                         if(!pool.try_pop(sws_device))
 269                         {
 270                                 double param;
 271                                 sws_device.reset(sws_getContext(input_av_frame->width, input_av_frame->height, static_cast<AVPixelFormat>(input_av_frame->format), width, height, AVPixelFormat::AV_PIX_FMT_BGRA, SWS_BILINEAR, nullptr, nullptr, &param), sws_freeContext);
 272                         }
 273
 274                         if(!sws_device)
 275                                 CASPAR_THROW_EXCEPTION(operation_failed() << msg_info("Could not create software scaling device.") << boost::errinfo_api_function("sws_getContext"));
 276
 277                         auto dest_frame = spl::make_shared<buffer>(width*height*4);
 278                         temp_buffers_.push(dest_frame);
 279
 280                         {
 281                                 auto dest_av_frame = ffmpeg::create_frame();
 282                                 avpicture_fill(reinterpret_cast<AVPicture*>(dest_av_frame.get()), dest_frame->data(), AVPixelFormat::AV_PIX_FMT_BGRA, width, height);
 283
 284                                 sws_scale(sws_device.get(), input_av_frame->data, input_av_frame->linesize, 0, input_av_frame->height, dest_av_frame->data, dest_av_frame->linesize);
 285                                 pool.push(sws_device);
 286                         }
 287
 288                         for(std::size_t n = 0; n < source_items.size(); ++n)
 289                         {
 290                                 if(source_items[n].data == data)
 291                                 {
 292                                         dest_items[n].data.fill(0);
 293                                         dest_items[n].data[0]                   = dest_frame->data();
 294                                         dest_items[n].pix_desc                  = core::pixel_format_desc(core::pixel_format::bgra);
 295                                         dest_items[n].pix_desc.planes   = { core::pixel_format_desc::plane(width, height, 4) };
 296                                         dest_items[n].transform                 = source_items[n].transform;
 297                                 }
 298                         }
 299                 });
 300
 301                 source_items = std::move(dest_items);
 302         }
 303 };
 304
 305 struct image_mixer::impl : boost::noncopyable
 306 {
 307         image_renderer                                          renderer_;
 308         std::vector<core::image_transform>      transform_stack_;
 309         std::vector<item>                                       items_; // layer/stream/items
 310 public:
 311         impl(int channel_id)
 312                 : transform_stack_(1)
 313         {
 314                 CASPAR_LOG(info) << L"Initialized Streaming SIMD Extensions Accelerated CPU Image Mixer for channel " << channel_id;
 315         }
 316
 317         void push(const core::frame_transform& transform)
 318         {
 319                 transform_stack_.push_back(transform_stack_.back()*transform.image_transform);
 320         }
 321
 322         void visit(const core::const_frame& frame)
 323         {
 324                 if(frame.pixel_format_desc().format == core::pixel_format::invalid)
 325                         return;
 326
 327                 if(frame.pixel_format_desc().planes.empty())
 328                         return;
 329
 330                 if(frame.pixel_format_desc().planes.at(0).size < 16)
 331                         return;
 332
 333                 if(transform_stack_.back().field_mode == core::field_mode::empty)
 334                         return;
 335
 336                 item item;
 337                 item.pix_desc   = frame.pixel_format_desc();
 338                 item.transform  = transform_stack_.back();
 339                 for(int n = 0; n < item.pix_desc.planes.size(); ++n)
 340                         item.data.at(n) = frame.image_data(n).begin();
 341
 342                 items_.push_back(item);
 343         }
 344
 345         void pop()
 346         {
 347                 transform_stack_.pop_back();
 348         }
 349
 350         std::future<array<const std::uint8_t>> render(const core::video_format_desc& format_desc)
 351         {
 352                 return renderer_(std::move(items_), format_desc);
 353         }
 354
 355         core::mutable_frame create_frame(const void* tag, const core::pixel_format_desc& desc, const core::audio_channel_layout& channel_layout)
 356         {
 357                 std::vector<array<std::uint8_t>> buffers;
 358                 for (auto& plane : desc.planes)
 359                 {
 360                         auto buf = spl::make_shared<buffer>(plane.size);
 361                         buffers.push_back(array<std::uint8_t>(buf->data(), plane.size, true, buf));
 362                 }
 363                 return core::mutable_frame(std::move(buffers), core::mutable_audio_buffer(), tag, desc, channel_layout);
 364         }
 365 };
 366
 367 image_mixer::image_mixer(int channel_id) : impl_(new impl(channel_id)){}
 368 image_mixer::~image_mixer(){}
 369 void image_mixer::push(const core::frame_transform& transform){impl_->push(transform);}
 370 void image_mixer::visit(const core::const_frame& frame){impl_->visit(frame);}
 371 void image_mixer::pop(){impl_->pop();}
 372 int image_mixer::get_max_frame_size() { return std::numeric_limits<int>::max(); }
 373 std::future<array<const std::uint8_t>> image_mixer::operator()(const core::video_format_desc& format_desc, bool /* straighten_alpha */){return impl_->render(format_desc);}
 374 core::mutable_frame image_mixer::create_frame(const void* tag, const core::pixel_format_desc& desc, const core::audio_channel_layout& channel_layout) {return impl_->create_frame(tag, desc, channel_layout);}
 375
 376 }}}