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

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