git.sesse.net Git - mlt/blob - src/framework/mlt_frame.c

   1 /*
   2  * mlt_frame.c -- interface for all frame classes
   3  * Copyright (C) 2003-2004 Ushodaya Enterprises Limited
   4  * Author: Charles Yates <charles.yates@pandora.be>
   5  *
   6  * This program 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  * This program 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 this program; if not, write to the Free Software Foundation,
  18  * Inc., 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA.
  19  */
  20
  21 #include "config.h"
  22 #include "mlt_frame.h"
  23 #include "mlt_producer.h"
  24 #include <stdio.h>
  25 #include <stdlib.h>
  26 #include <string.h>
  27
  28 /** Constructor for a frame.
  29 */
  30
  31 mlt_frame mlt_frame_init( )
  32 {
  33         // Allocate a frame
  34         mlt_frame this = calloc( sizeof( struct mlt_frame_s ), 1 );
  35
  36         if ( this != NULL )
  37         {
  38                 // Get the normalisation
  39                 char *normalisation = getenv( "MLT_NORMALISATION" );
  40
  41                 // Initialise the properties
  42                 mlt_properties properties = &this->parent;
  43                 mlt_properties_init( properties, this );
  44
  45                 // Set default properties on the frame
  46                 mlt_properties_set_position( properties, "_position", 0.0 );
  47                 mlt_properties_set_data( properties, "image", NULL, 0, NULL, NULL );
  48
  49                 if ( normalisation == NULL || strcmp( normalisation, "NTSC" ) )
  50                 {
  51                         mlt_properties_set_int( properties, "width", 720 );
  52                         mlt_properties_set_int( properties, "height", 576 );
  53                         mlt_properties_set_int( properties, "normalised_width", 720 );
  54                         mlt_properties_set_int( properties, "normalised_height", 576 );
  55                 }
  56                 else
  57                 {
  58                         mlt_properties_set_int( properties, "width", 720 );
  59                         mlt_properties_set_int( properties, "height", 480 );
  60                         mlt_properties_set_int( properties, "normalised_width", 720 );
  61                         mlt_properties_set_int( properties, "normalised_height", 480 );
  62                 }
  63
  64                 mlt_properties_set_double( properties, "aspect_ratio", 4.0 / 3.0 );
  65                 mlt_properties_set_data( properties, "audio", NULL, 0, NULL, NULL );
  66                 mlt_properties_set_data( properties, "alpha", NULL, 0, NULL, NULL );
  67
  68                 // Construct stacks for frames and methods
  69                 this->stack_get_image = mlt_deque_init( );
  70                 this->stack_frame = mlt_deque_init( );
  71         }
  72
  73         return this;
  74 }
  75
  76 /** Fetch the frames properties.
  77 */
  78
  79 mlt_properties mlt_frame_properties( mlt_frame this )
  80 {
  81         return &this->parent;
  82 }
  83
  84 /** Check if we have a way to derive something other than a test card.
  85 */
  86
  87 int mlt_frame_is_test_card( mlt_frame this )
  88 {
  89         return mlt_properties_get_int( mlt_frame_properties( this ), "test_image" );
  90 }
  91
  92 /** Check if we have a way to derive something than test audio.
  93 */
  94
  95 int mlt_frame_is_test_audio( mlt_frame this )
  96 {
  97         return this->get_audio == NULL || mlt_properties_get_int( mlt_frame_properties( this ), "test_audio" );
  98 }
  99
 100 /** Get the aspect ratio of the frame.
 101 */
 102
 103 double mlt_frame_get_aspect_ratio( mlt_frame this )
 104 {
 105         return mlt_properties_get_double( mlt_frame_properties( this ), "aspect_ratio" );
 106 }
 107
 108 /** Set the aspect ratio of the frame.
 109 */
 110
 111 int mlt_frame_set_aspect_ratio( mlt_frame this, double value )
 112 {
 113         return mlt_properties_set_double( mlt_frame_properties( this ), "aspect_ratio", value );
 114 }
 115
 116 /** Get the position of this frame.
 117 */
 118
 119 mlt_position mlt_frame_get_position( mlt_frame this )
 120 {
 121         return mlt_properties_get_position( mlt_frame_properties( this ), "_position" );
 122 }
 123
 124 /** Set the position of this frame.
 125 */
 126
 127 int mlt_frame_set_position( mlt_frame this, mlt_position value )
 128 {
 129         return mlt_properties_set_position( mlt_frame_properties( this ), "_position", value );
 130 }
 131
 132 /** Stack a get_image callback.
 133 */
 134
 135 int mlt_frame_push_get_image( mlt_frame this, mlt_get_image get_image )
 136 {
 137         return mlt_deque_push_back( this->stack_get_image, get_image );
 138 }
 139
 140 /** Pop a get_image callback.
 141 */
 142
 143 mlt_get_image mlt_frame_pop_get_image( mlt_frame this )
 144 {
 145         return mlt_deque_pop_back( this->stack_get_image );
 146 }
 147
 148 /** Push a frame.
 149 */
 150
 151 int mlt_frame_push_frame( mlt_frame this, mlt_frame that )
 152 {
 153         return mlt_deque_push_back( this->stack_frame, that );
 154 }
 155
 156 /** Pop a frame.
 157 */
 158
 159 mlt_frame mlt_frame_pop_frame( mlt_frame this )
 160 {
 161         return mlt_deque_pop_back( this->stack_frame );
 162 }
 163
 164 int mlt_frame_get_image( mlt_frame this, uint8_t **buffer, mlt_image_format *format, int *width, int *height, int writable )
 165 {
 166         mlt_properties properties = mlt_frame_properties( this );
 167         mlt_get_image get_image = mlt_frame_pop_get_image( this );
 168         mlt_producer producer = mlt_properties_get_data( properties, "test_card_producer", NULL );
 169
 170         if ( get_image != NULL )
 171         {
 172                 return get_image( this, buffer, format, width, height, writable );
 173         }
 174         else if ( mlt_properties_get_data( properties, "image", NULL ) != NULL )
 175         {
 176                 *format = mlt_image_yuv422;
 177                 *buffer = mlt_properties_get_data( properties, "image", NULL );
 178                 *width = mlt_properties_get_int( properties, "width" );
 179                 *height = mlt_properties_get_int( properties, "height" );
 180         }
 181         else if ( producer != NULL )
 182         {
 183                 mlt_frame test_frame = NULL;
 184                 mlt_service_get_frame( mlt_producer_service( producer ), &test_frame, 0 );
 185                 if ( test_frame != NULL )
 186                 {
 187                         mlt_properties test_properties = mlt_frame_properties( test_frame );
 188                         mlt_properties_set_double( test_properties, "consumer_aspect_ratio", mlt_properties_get_double( properties, "consumer_aspect_ratio" ) );
 189                         mlt_properties_set_double( test_properties, "consumer_scale", mlt_properties_get_double( properties, "consumer_scale" ) );
 190                         mlt_properties_set( test_properties, "rescale.interp", "nearest" );
 191                         mlt_frame_get_image( test_frame, buffer, format, width, height, writable );
 192                         mlt_properties_set_data( properties, "test_card_frame", test_frame, 0, ( mlt_destructor )mlt_frame_close, NULL );
 193                         mlt_properties_set_data( properties, "image", *buffer, *width * *height * 2, NULL, NULL );
 194                         mlt_properties_set_int( properties, "width", *width );
 195                         mlt_properties_set_int( properties, "height", *height );
 196                 }
 197                 else
 198                 {
 199                         mlt_properties_set_data( properties, "test_card_producer", NULL, 0, NULL, NULL );
 200                         mlt_frame_get_image( this, buffer, format, width, height, writable );
 201                 }
 202         }
 203         else
 204         {
 205                 uint8_t *p;
 206                 uint8_t *q;
 207                 int size = 0;
 208
 209                 *width = *width == 0 ? 720 : *width;
 210                 *height = *height == 0 ? 576 : *height;
 211                 size = *width * *height;
 212
 213                 mlt_properties_set_int( properties, "width", *width );
 214                 mlt_properties_set_int( properties, "height", *height );
 215
 216                 switch( *format )
 217                 {
 218                         case mlt_image_none:
 219                                 size = 0;
 220                                 *buffer = NULL;
 221                                 break;
 222                         case mlt_image_rgb24:
 223                                 size *= 3;
 224                                 size += *width * 3;
 225                                 *buffer = mlt_pool_alloc( size );
 226                                 if ( *buffer )
 227                                         memset( *buffer, 255, size );
 228                                 break;
 229                         case mlt_image_rgb24a:
 230                                 size *= 4;
 231                                 size += *width * 4;
 232                                 *buffer = mlt_pool_alloc( size );
 233                                 if ( *buffer )
 234                                         memset( *buffer, 255, size );
 235                                 break;
 236                         case mlt_image_yuv422:
 237                                 size *= 2;
 238                                 size += *width * 2;
 239                                 *buffer = mlt_pool_alloc( size );
 240                                 p = *buffer;
 241                                 q = p + size;
 242                                 while ( p != NULL && p != q )
 243                                 {
 244                                         *p ++ = 255;
 245                                         *p ++ = 128;
 246                                 }
 247                                 break;
 248                         case mlt_image_yuv420p:
 249                                 size = size * 3 / 2;
 250                                 *buffer = mlt_pool_alloc( size );
 251                                 if ( *buffer )
 252                                         memset( *buffer, 255, size );
 253                                 break;
 254                 }
 255
 256                 mlt_properties_set_data( properties, "image", *buffer, size, ( mlt_destructor )mlt_pool_release, NULL );
 257                 mlt_properties_set_int( properties, "test_image", 1 );
 258         }
 259
 260         return 0;
 261 }
 262
 263 uint8_t *mlt_frame_get_alpha_mask( mlt_frame this )
 264 {
 265         if ( this->get_alpha_mask != NULL )
 266                 return this->get_alpha_mask( this );
 267         return NULL;
 268 }
 269
 270 int mlt_frame_get_audio( mlt_frame this, int16_t **buffer, mlt_audio_format *format, int *frequency, int *channels, int *samples )
 271 {
 272         mlt_properties properties = mlt_frame_properties( this );
 273
 274         if ( this->get_audio != NULL )
 275         {
 276                 return this->get_audio( this, buffer, format, frequency, channels, samples );
 277         }
 278         else
 279         {
 280                 int size = 0;
 281                 *samples = *samples <= 0 ? 1920 : *samples;
 282                 *channels = *channels <= 0 ? 2 : *channels;
 283                 *frequency = *frequency <= 0 ? 48000 : *frequency;
 284                 size = *samples * *channels * sizeof( int16_t );
 285                 *buffer = mlt_pool_alloc( size );
 286                 if ( *buffer != NULL )
 287                         memset( *buffer, 0, size );
 288                 mlt_properties_set_data( properties, "audio", *buffer, size, ( mlt_destructor )mlt_pool_release, NULL );
 289                 mlt_properties_set_int( properties, "test_audio", 1 );
 290         }
 291         return 0;
 292 }
 293
 294 void mlt_frame_close( mlt_frame this )
 295 {
 296         mlt_deque_close( this->stack_get_image );
 297         mlt_deque_close( this->stack_frame );
 298         mlt_properties_close( &this->parent );
 299         free( this );
 300 }
 301
 302 /***** convenience functions *****/
 303 #define RGB2YUV(r, g, b, y, u, v)\
 304   y = (306*r + 601*g + 117*b)  >> 10;\
 305   u = ((-172*r - 340*g + 512*b) >> 10)  + 128;\
 306   v = ((512*r - 429*g - 83*b) >> 10) + 128;\
 307   y = y < 0 ? 0 : y;\
 308   u = u < 0 ? 0 : u;\
 309   v = v < 0 ? 0 : v;\
 310   y = y > 255 ? 255 : y;\
 311   u = u > 255 ? 255 : u;\
 312   v = v > 255 ? 255 : v
 313
 314 int mlt_convert_rgb24a_to_yuv422( uint8_t *rgba, int width, int height, int stride, uint8_t *yuv, uint8_t *alpha )
 315 {
 316         int ret = 0;
 317         register int y0, y1, u0, u1, v0, v1;
 318         register int r, g, b;
 319         register uint8_t *d = yuv;
 320         register int i, j;
 321
 322         for ( i = 0; i < height; i++ )
 323         {
 324                 register uint8_t *s = rgba + ( stride * i );
 325                 for ( j = 0; j < ( width / 2 ); j++ )
 326                 {
 327                         r = *s++;
 328                         g = *s++;
 329                         b = *s++;
 330                         *alpha++ = *s++;
 331                         RGB2YUV (r, g, b, y0, u0 , v0);
 332                         r = *s++;
 333                         g = *s++;
 334                         b = *s++;
 335                         *alpha++ = *s++;
 336                         RGB2YUV (r, g, b, y1, u1 , v1);
 337                         *d++ = y0;
 338                         *d++ = (u0+u1) >> 1;
 339                         *d++ = y1;
 340                         *d++ = (v0+v1) >> 1;
 341                 }
 342                 if ( width % 2 )
 343                 {
 344                         r = *s++;
 345                         g = *s++;
 346                         b = *s++;
 347                         *alpha++ = *s++;
 348                         RGB2YUV (r, g, b, y0, u0 , v0);
 349                         *d++ = y0;
 350                         *d++ = u0;
 351                 }
 352         }
 353         return ret;
 354 }
 355
 356 int mlt_convert_rgb24_to_yuv422( uint8_t *rgb, int width, int height, int stride, uint8_t *yuv )
 357 {
 358         int ret = 0;
 359         register int y0, y1, u0, u1, v0, v1;
 360         register int r, g, b;
 361         register uint8_t *d = yuv;
 362         register int i, j;
 363
 364         for ( i = 0; i < height; i++ )
 365         {
 366                 register uint8_t *s = rgb + ( stride * i );
 367                 for ( j = 0; j < ( width / 2 ); j++ )
 368                 {
 369                         r = *s++;
 370                         g = *s++;
 371                         b = *s++;
 372                         RGB2YUV (r, g, b, y0, u0 , v0);
 373                         r = *s++;
 374                         g = *s++;
 375                         b = *s++;
 376                         RGB2YUV (r, g, b, y1, u1 , v1);
 377                         *d++ = y0;
 378                         *d++ = (u0+u1) >> 1;
 379                         *d++ = y1;
 380                         *d++ = (v0+v1) >> 1;
 381                 }
 382                 if ( width % 2 )
 383                 {
 384                         r = *s++;
 385                         g = *s++;
 386                         b = *s++;
 387                         RGB2YUV (r, g, b, y0, u0 , v0);
 388                         *d++ = y0;
 389                         *d++ = u0;
 390                 }
 391         }
 392         return ret;
 393 }
 394
 395 int mlt_convert_yuv420p_to_yuv422( uint8_t *yuv420p, int width, int height, int stride, uint8_t *yuv )
 396 {
 397         int ret = 0;
 398         register int i, j;
 399
 400         int half = width >> 1;
 401
 402         uint8_t *Y = yuv420p;
 403         uint8_t *U = Y + width * height;
 404         uint8_t *V = U + width * height / 4;
 405
 406         register uint8_t *d = yuv;
 407
 408         for ( i = 0; i < height; i++ )
 409         {
 410                 register uint8_t *u = U + ( i / 2 ) * ( half );
 411                 register uint8_t *v = V + ( i / 2 ) * ( half );
 412
 413                 for ( j = 0; j < half; j++ )
 414                 {
 415                         *d ++ = *Y ++;
 416                         *d ++ = *u ++;
 417                         *d ++ = *Y ++;
 418                         *d ++ = *v ++;
 419                 }
 420         }
 421         return ret;
 422 }
 423
 424 void mlt_resize_yuv422( uint8_t *output, int owidth, int oheight, uint8_t *input, int iwidth, int iheight )
 425 {
 426         // Calculate strides
 427         int istride = iwidth * 2;
 428         int ostride = owidth * 2;
 429
 430         iwidth = iwidth - ( iwidth % 4 );
 431         owidth = owidth - ( owidth % 4 );
 432         iheight = iheight - ( iheight % 2 );
 433         oheight = oheight - ( oheight % 2 );
 434
 435         // Optimisation point
 436         if ( iwidth == owidth && iheight == oheight )
 437                 memcpy( output, input, iheight * istride );
 438
 439         // Coordinates (0,0 is middle of output)
 440         int y;
 441
 442         // Calculate ranges
 443         int out_x_range = owidth / 2;
 444         int out_y_range = oheight / 2;
 445         int in_x_range = iwidth / 2 < out_x_range ? iwidth / 2 : out_x_range;
 446         int in_y_range = iheight / 2 < out_y_range ? iheight / 2 : out_y_range;
 447
 448         // Output pointers
 449         uint8_t *out_line = output;
 450         uint8_t *out_ptr = out_line;
 451
 452         // Calculate a middle and possibly invalid pointer in the input
 453         uint8_t *in_middle = input + istride * ( iheight / 2 ) + ( iwidth / 2 ) * 2;
 454         int in_line = - in_y_range * istride - in_x_range * 2;
 455
 456         int elements;
 457
 458         // Fill whole section with black
 459         y = out_y_range - ( iheight / 2 );
 460         int blank_elements = ostride * y / 2;
 461         elements = blank_elements;
 462         while ( elements -- )
 463         {
 464                 *out_line ++ = 0;
 465                 *out_line ++ = 128;
 466         }
 467
 468         int active_width = 2 * iwidth;
 469         int inactive_width = out_x_range - in_x_range;
 470
 471         // Loop for the entirety of our output height.
 472         while ( iheight -- )
 473         {
 474         // Start at the beginning of the line
 475         out_ptr = out_line;
 476
 477                 // Fill the outer part with black
 478                 elements = inactive_width;
 479                 while ( elements -- )
 480                 {
 481                         *out_ptr ++ = 0;
 482                         *out_ptr ++ = 128;
 483                 }
 484
 485                 // We're in the input range for this row.
 486                 memcpy( out_ptr, in_middle + in_line, active_width );
 487                 out_ptr += active_width;
 488
 489                 // Fill the outer part with black
 490                 elements = inactive_width;
 491                 while ( elements -- )
 492                 {
 493                         *out_ptr ++ = 0;
 494                         *out_ptr ++ = 128;
 495                 }
 496
 497                 // Move to next input line
 498                 in_line += istride;
 499
 500         // Move to next output line
 501         out_line += ostride;
 502         }
 503
 504         // Fill whole section with black
 505         elements = blank_elements;
 506         while ( elements -- )
 507         {
 508                 *out_line ++ = 0;
 509                 *out_line ++ = 128;
 510         }
 511 }
 512
 513 /** A resizing function for yuv422 frames - this does not rescale, but simply
 514         resizes. It assumes yuv422 images available on the frame so use with care.
 515 */
 516
 517 uint8_t *mlt_frame_resize_yuv422( mlt_frame this, int owidth, int oheight )
 518 {
 519         // Get properties
 520         mlt_properties properties = mlt_frame_properties( this );
 521
 522         // Get the input image, width and height
 523         uint8_t *input = mlt_properties_get_data( properties, "image", NULL );
 524         int iwidth = mlt_properties_get_int( properties, "width" );
 525         int iheight = mlt_properties_get_int( properties, "height" );
 526
 527         // If width and height are correct, don't do anything
 528         if ( iwidth != owidth || iheight != oheight )
 529         {
 530                 // Create the output image
 531                 uint8_t *output = mlt_pool_alloc( owidth * ( oheight + 1 ) * 2 );
 532
 533                 // Call the generic resize
 534                 mlt_resize_yuv422( output, owidth, oheight, input, iwidth, iheight );
 535
 536                 // Now update the frame
 537                 mlt_properties_set_data( properties, "image", output, owidth * ( oheight + 1 ) * 2, ( mlt_destructor )mlt_pool_release, NULL );
 538                 mlt_properties_set_int( properties, "width", owidth );
 539                 mlt_properties_set_int( properties, "height", oheight );
 540
 541                 // Return the output
 542                 return output;
 543         }
 544         // No change, return input
 545         return input;
 546 }
 547
 548 /** A rescaling function for yuv422 frames - low quality, and provided for testing
 549         only. It assumes yuv422 images available on the frame so use with care.
 550 */
 551
 552 uint8_t *mlt_frame_rescale_yuv422( mlt_frame this, int owidth, int oheight )
 553 {
 554         // Get properties
 555         mlt_properties properties = mlt_frame_properties( this );
 556
 557         // Get the input image, width and height
 558         uint8_t *input = mlt_properties_get_data( properties, "image", NULL );
 559         int iwidth = mlt_properties_get_int( properties, "width" );
 560         int iheight = mlt_properties_get_int( properties, "height" );
 561
 562         // If width and height are correct, don't do anything
 563         if ( iwidth != owidth || iheight != oheight )
 564         {
 565                 // Create the output image
 566                 uint8_t *output = mlt_pool_alloc( owidth * ( oheight + 1 ) * 2 );
 567
 568                 // Calculate strides
 569                 int istride = iwidth * 2;
 570                 int ostride = owidth * 2;
 571
 572                 iwidth = iwidth - ( iwidth % 4 );
 573
 574         // Coordinates (0,0 is middle of output)
 575         int y, x;
 576
 577                 // Derived coordinates
 578                 int dy, dx;
 579
 580         // Calculate ranges
 581         int out_x_range = owidth / 2;
 582         int out_y_range = oheight / 2;
 583         int in_x_range = iwidth / 2;
 584         int in_y_range = iheight / 2;
 585
 586         // Output pointers
 587         uint8_t *out_line = output;
 588         uint8_t *out_ptr;
 589
 590         // Calculate a middle pointer
 591         uint8_t *in_middle = input + istride * in_y_range + in_x_range * 2;
 592         uint8_t *in_line;
 593                 uint8_t *in_ptr;
 594
 595                 // Generate the affine transform scaling values
 596                 int scale_width = ( iwidth << 16 ) / owidth;
 597                 int scale_height = ( iheight << 16 ) / oheight;
 598
 599         // Loop for the entirety of our output height.
 600         for ( y = - out_y_range; y < out_y_range ; y ++ )
 601         {
 602                         // Calculate the derived y value
 603                         dy = ( scale_height * y ) >> 16;
 604
 605                 // Start at the beginning of the line
 606                 out_ptr = out_line;
 607
 608                 // Pointer to the middle of the input line
 609                 in_line = in_middle + dy * istride;
 610
 611                 // Loop for the entirety of our output row.
 612                 for ( x = - out_x_range; x < out_x_range; x += 1 )
 613                 {
 614                                 // Calculated the derived x
 615                                 dx = ( scale_width * x ) >> 16;
 616
 617                 // We're in the input range for this row.
 618                                 in_ptr = in_line + ( dx << 1 );
 619                 *out_ptr ++ = *in_ptr ++;
 620                                 in_ptr = in_line + ( ( dx >> 1 ) << 2 ) + ( ( x & 1 ) << 1 ) + 1;
 621                 *out_ptr ++ = *in_ptr;
 622                 }
 623
 624                 // Move to next output line
 625                 out_line += ostride;
 626         }
 627
 628                 // Now update the frame
 629                 mlt_properties_set_data( properties, "image", output, owidth * ( oheight + 1 ) * 2, ( mlt_destructor )mlt_pool_release, NULL );
 630                 mlt_properties_set_int( properties, "width", owidth );
 631                 mlt_properties_set_int( properties, "height", oheight );
 632
 633                 // Return the output
 634                 return output;
 635         }
 636
 637         // No change, return input
 638         return input;
 639 }
 640
 641 int mlt_frame_mix_audio( mlt_frame this, mlt_frame that, float weight_start, float weight_end, int16_t **buffer, mlt_audio_format *format, int *frequency, int *channels, int *samples )
 642 {
 643         int ret = 0;
 644         int16_t *p_src, *p_dest;
 645         int16_t *src, *dest;
 646         //static int16_t *extra_src = NULL, *extra_dest = NULL;
 647         static int extra_src_samples = 0, extra_dest_samples = 0;
 648         int frequency_src = *frequency, frequency_dest = *frequency;
 649         int channels_src = *channels, channels_dest = *channels;
 650         int samples_src = *samples, samples_dest = *samples;
 651         int i, j;
 652         double d = 0, s = 0;
 653
 654         mlt_frame_get_audio( this, &p_dest, format, &frequency_dest, &channels_dest, &samples_dest );
 655         //fprintf( stderr, "mix: frame dest samples %d channels %d position %lld\n", samples_dest, channels_dest, mlt_properties_get_position( mlt_frame_properties( this ), "_position" ) );
 656         mlt_frame_get_audio( that, &p_src, format, &frequency_src, &channels_src, &samples_src );
 657         //fprintf( stderr, "mix: frame src  samples %d channels %d\n", samples_src, channels_src );
 658         src = p_src;
 659         dest = p_dest;
 660         if ( channels_src > 6 )
 661                 channels_src = 0;
 662         if ( channels_dest > 6 )
 663                 channels_dest = 0;
 664         if ( samples_src > 4000 )
 665                 samples_src = 0;
 666         if ( samples_dest > 4000 )
 667                 samples_dest = 0;
 668
 669 #if 0
 670         // Append new samples to leftovers
 671         if ( extra_dest_samples > 0 )
 672         {
 673                 fprintf( stderr, "prepending %d samples to dest\n", extra_dest_samples );
 674                 dest = realloc( extra_dest, ( samples_dest + extra_dest_samples ) * 2 * channels_dest );
 675                 memcpy( &extra_dest[ extra_dest_samples * channels_dest ], p_dest, samples_dest * 2 * channels_dest );
 676         }
 677         else
 678                 dest = p_dest;
 679         if ( extra_src_samples > 0 )
 680         {
 681                 fprintf( stderr, "prepending %d samples to src\n", extra_src_samples );
 682                 src = realloc( extra_src, ( samples_src + extra_src_samples ) * 2 * channels_src );
 683                 memcpy( &extra_src[ extra_src_samples * channels_src ], p_src, samples_src * 2 * channels_src );
 684         }
 685         else
 686                 src = p_src;
 687 #endif
 688
 689         // determine number of samples to process
 690         if ( samples_src + extra_src_samples < samples_dest + extra_dest_samples )
 691                 *samples = samples_src + extra_src_samples;
 692         else if ( samples_dest + extra_dest_samples < samples_src + extra_src_samples )
 693                 *samples = samples_dest + extra_dest_samples;
 694
 695         *channels = channels_src < channels_dest ? channels_src : channels_dest;
 696         *buffer = p_dest;
 697         *frequency = frequency_dest;
 698
 699         // Compute a smooth ramp over start to end
 700         float weight = weight_start;
 701         float weight_step = ( weight_end - weight_start ) / *samples;
 702
 703         // Mixdown
 704         for ( i = 0; i < *samples; i++ )
 705         {
 706                 for ( j = 0; j < *channels; j++ )
 707                 {
 708                         if ( j < channels_dest )
 709                                 d = (double) dest[ i * channels_dest + j ];
 710                         if ( j < channels_src )
 711                                 s = (double) src[ i * channels_src + j ];
 712                         dest[ i * channels_dest + j ] = s * weight + d * ( 1.0 - weight );
 713                 }
 714                 weight += weight_step;
 715         }
 716
 717         // We have to copy --sigh
 718         if ( dest != p_dest )
 719                 memcpy( p_dest, dest, *samples * 2 * *channels );
 720
 721 #if 0
 722         // Store the leftovers
 723         if ( samples_src + extra_src_samples < samples_dest + extra_dest_samples )
 724         {
 725                 extra_dest_samples = ( samples_dest + extra_dest_samples ) - ( samples_src + extra_src_samples );
 726                 size_t size = extra_dest_samples * 2 * channels_dest;
 727                 fprintf( stderr, "storing %d samples from dest\n", extra_dest_samples );
 728                 if ( extra_dest )
 729                         free( extra_dest );
 730                 extra_dest = malloc( size );
 731                 if ( extra_dest )
 732                         memcpy( extra_dest, &p_dest[ ( samples_dest - extra_dest_samples - 1 ) * channels_dest ], size );
 733                 else
 734                         extra_dest_samples = 0;
 735         }
 736         else if ( samples_dest + extra_dest_samples < samples_src + extra_src_samples )
 737         {
 738                 extra_src_samples = ( samples_src + extra_src_samples ) - ( samples_dest + extra_dest_samples );
 739                 size_t size = extra_src_samples * 2 * channels_src;
 740                 fprintf( stderr, "storing %d samples from src\n", extra_dest_samples );
 741                 if ( extra_src )
 742                         free( extra_src );
 743                 extra_src = malloc( size );
 744                 if ( extra_src )
 745                         memcpy( extra_src, &p_src[ ( samples_src - extra_src_samples - 1 ) * channels_src ], size );
 746                 else
 747                         extra_src_samples = 0;
 748         }
 749 #endif
 750
 751         return ret;
 752 }
 753
 754 int mlt_sample_calculator( float fps, int frequency, int64_t position )
 755 {
 756         int samples = 0;
 757
 758         if ( fps > 29 && fps <= 30 )
 759         {
 760                 samples = frequency / 30;
 761
 762                 switch ( frequency )
 763                 {
 764                         case 48000:
 765                                 if ( position % 5 != 0 )
 766                                         samples += 2;
 767                                 break;
 768                         case 44100:
 769                                 if ( position % 300 == 0 )
 770                                         samples = 1471;
 771                                 else if ( position % 30 == 0 )
 772                                         samples = 1470;
 773                                 else if ( position % 2 == 0 )
 774                                         samples = 1472;
 775                                 else
 776                                         samples = 1471;
 777                                 break;
 778                         case 32000:
 779                                 if ( position % 30 == 0 )
 780                                         samples = 1068;
 781                                 else if ( position % 29 == 0 )
 782                                         samples = 1067;
 783                                 else if ( position % 4 == 2 )
 784                                         samples = 1067;
 785                                 else
 786                                         samples = 1068;
 787                                 break;
 788                         default:
 789                                 samples = 0;
 790                 }
 791         }
 792         else if ( fps != 0 )
 793         {
 794                 samples = frequency / fps;
 795         }
 796
 797         return samples;
 798 }
 799