2 * SBG (SBaGen) file format decoder
3 * Copyright (c) 2011 Nicolas George
5 * This file is part of FFmpeg.
7 * FFmpeg is free software; you can redistribute it and/or
8 * modify it under the terms of the GNU Lesser General Public
9 * License as published by the Free Software Foundation; either
10 * version 2.1 of the License, or (at your option) any later version.
12 * FFmpeg is distributed in the hope that it will be useful,
13 * but WITHOUT ANY WARRANTY; without even the implied warranty of
14 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
15 * Lesser General Public License for more details.
17 * You should have received a copy of the GNU Lesser General Public
18 * License along with FFmpeg; if not, write to the Free Software
19 * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
25 #include "libavutil/bprint.h"
26 #include "libavutil/channel_layout.h"
27 #include "libavutil/intreadwrite.h"
28 #include "libavutil/log.h"
29 #include "libavutil/mem.h"
30 #include "libavutil/opt.h"
31 #include "libavutil/time_internal.h"
36 #define SBG_SCALE (1 << 16)
37 #define DAY (24 * 60 * 60)
38 #define DAY_TS ((int64_t)DAY * AV_TIME_BASE)
59 int8_t in
, out
, slide
;
71 /* bell: freq constant, ampl decreases exponentially, can be approx lin */
73 struct sbg_timestamp
{
75 char type
; /* 0 for relative, 'N' for now, 'T' for absolute */
78 struct sbg_script_definition
{
81 int elements
, nb_elements
;
82 char type
; /* 'S' or 'B' */
85 struct sbg_script_synth
{
89 enum sbg_synth_type type
;
95 struct sbg_script_tseq
{
96 struct sbg_timestamp ts
;
100 struct sbg_fade fade
;
103 struct sbg_script_event
{
105 int64_t ts_int
, ts_trans
, ts_next
;
106 int elements
, nb_elements
;
107 struct sbg_fade fade
;
111 struct sbg_script_definition
*def
;
112 struct sbg_script_synth
*synth
;
113 struct sbg_script_tseq
*tseq
;
114 struct sbg_script_tseq
*block_tseq
;
115 struct sbg_script_event
*events
;
122 int64_t opt_fade_time
;
123 int64_t opt_duration
;
126 uint8_t opt_start_at_first
;
127 uint8_t opt_end_at_last
;
134 struct sbg_script scs
;
135 struct sbg_timestamp current_time
;
137 int nb_def_max
, nb_synth_max
, nb_tseq_max
, nb_block_tseq_max
;
142 enum ws_interval_type
{
143 WS_SINE
= MKTAG('S','I','N','E'),
144 WS_NOISE
= MKTAG('N','O','I','S'),
149 enum ws_interval_type type
;
156 struct ws_intervals
{
157 struct ws_interval
*inter
;
162 static void *alloc_array_elem(void **array
, size_t elsize
,
163 int *size
, int *max_size
)
167 if (*size
== *max_size
) {
168 int m
= FFMAX(32, FFMIN(*max_size
, INT_MAX
/ 2) * 2);
171 *array
= av_realloc_f(*array
, m
, elsize
);
176 ret
= (char *)*array
+ elsize
* *size
;
177 memset(ret
, 0, elsize
);
182 static int str_to_time(const char *str
, int64_t *rtime
)
184 const char *cur
= str
;
190 if (*cur
< '0' || *cur
> '9')
192 hours
= strtol(cur
, &end
, 10);
193 if (end
== cur
|| *end
!= ':' || end
[1] < '0' || end
[1] > '9')
196 minutes
= strtol(cur
, &end
, 10);
201 seconds
= strtod(cur
+ 1, &end
);
204 ts
= av_clipd(seconds
* AV_TIME_BASE
, INT64_MIN
/2, INT64_MAX
/2);
206 *rtime
= av_sat_add64((hours
* 3600LL + minutes
* 60LL) * AV_TIME_BASE
, ts
);
210 static inline int is_space(char c
)
212 return c
== ' ' || c
== '\t' || c
== '\r';
215 static inline int scale_double(void *log
, double d
, double m
, int *r
)
218 if (m
< INT_MIN
|| m
>= INT_MAX
) {
220 av_log(log
, AV_LOG_ERROR
, "%g is too large\n", d
);
221 return AVERROR(EDOM
);
227 static int lex_space(struct sbg_parser
*p
)
231 while (p
->cursor
< p
->end
&& is_space(*p
->cursor
))
233 return p
->cursor
> c
;
236 static int lex_char(struct sbg_parser
*p
, char c
)
238 int r
= p
->cursor
< p
->end
&& *p
->cursor
== c
;
244 static int lex_double(struct sbg_parser
*p
, double *r
)
249 if (p
->cursor
== p
->end
|| is_space(*p
->cursor
) || *p
->cursor
== '\n')
251 d
= strtod(p
->cursor
, &end
);
252 if (end
> p
->cursor
) {
260 static int lex_fixed(struct sbg_parser
*p
, const char *t
, int l
)
262 if (p
->end
- p
->cursor
< l
|| memcmp(p
->cursor
, t
, l
))
268 static int lex_line_end(struct sbg_parser
*p
)
270 if (p
->cursor
< p
->end
&& *p
->cursor
== '#') {
272 while (p
->cursor
< p
->end
&& *p
->cursor
!= '\n')
275 if (p
->cursor
== p
->end
)
276 /* simulate final LF for files lacking it */
278 if (*p
->cursor
!= '\n')
286 static int lex_wsword(struct sbg_parser
*p
, struct sbg_string
*rs
)
288 char *s
= p
->cursor
, *c
= s
;
290 if (s
== p
->end
|| *s
== '\n')
292 while (c
< p
->end
&& *c
!= '\n' && !is_space(*c
))
295 rs
->e
= p
->cursor
= c
;
300 static int lex_name(struct sbg_parser
*p
, struct sbg_string
*rs
)
302 char *s
= p
->cursor
, *c
= s
;
304 while (c
< p
->end
&& ((*c
>= 'a' && *c
<= 'z') || (*c
>= 'A' && *c
<= 'Z')
305 || (*c
>= '0' && *c
<= '9') || *c
== '_' || *c
== '-'))
310 rs
->e
= p
->cursor
= c
;
314 static int lex_time(struct sbg_parser
*p
, int64_t *rt
)
316 int r
= str_to_time(p
->cursor
, rt
);
321 #define FORWARD_ERROR(c) \
325 return errcode ? errcode : AVERROR_INVALIDDATA; \
328 static int parse_immediate(struct sbg_parser
*p
)
330 snprintf(p
->err_msg
, sizeof(p
->err_msg
),
331 "immediate sequences not yet implemented");
332 return AVERROR_PATCHWELCOME
;
335 static int parse_preprogrammed(struct sbg_parser
*p
)
337 snprintf(p
->err_msg
, sizeof(p
->err_msg
),
338 "preprogrammed sequences not yet implemented");
339 return AVERROR_PATCHWELCOME
;
342 static int parse_optarg(struct sbg_parser
*p
, char o
, struct sbg_string
*r
)
344 if (!lex_wsword(p
, r
)) {
345 snprintf(p
->err_msg
, sizeof(p
->err_msg
),
346 "option '%c' requires an argument", o
);
347 return AVERROR_INVALIDDATA
;
352 static int parse_options(struct sbg_parser
*p
)
354 struct sbg_string ostr
, oarg
;
360 if (p
->cursor
== p
->end
|| *p
->cursor
!= '-')
362 while (lex_char(p
, '-') && lex_wsword(p
, &ostr
)) {
363 for (; ostr
.s
< ostr
.e
; ostr
.s
++) {
367 p
->scs
.opt_start_at_first
= 1;
370 p
->scs
.opt_end_at_last
= 1;
379 FORWARD_ERROR(parse_optarg(p
, opt
, &oarg
));
380 v
= strtod(oarg
.s
, &tptr
);
381 if (oarg
.e
!= tptr
) {
382 snprintf(p
->err_msg
, sizeof(p
->err_msg
),
383 "syntax error for option -F");
384 return AVERROR_INVALIDDATA
;
386 p
->scs
.opt_fade_time
= v
* AV_TIME_BASE
/ 1000;
389 FORWARD_ERROR(parse_optarg(p
, opt
, &oarg
));
390 r
= str_to_time(oarg
.s
, &p
->scs
.opt_duration
);
391 if (oarg
.e
!= oarg
.s
+ r
|| p
->scs
.opt_duration
< 0) {
392 snprintf(p
->err_msg
, sizeof(p
->err_msg
),
393 "syntax error for option -L");
394 return AVERROR_INVALIDDATA
;
398 FORWARD_ERROR(parse_optarg(p
, opt
, &oarg
));
399 r
= str_to_time(oarg
.s
, &p
->scs
.start_ts
);
400 if (oarg
.e
!= oarg
.s
+ r
) {
401 snprintf(p
->err_msg
, sizeof(p
->err_msg
),
402 "syntax error for option -T");
403 return AVERROR_INVALIDDATA
;
407 FORWARD_ERROR(parse_optarg(p
, opt
, &oarg
));
408 tptr
= av_malloc(oarg
.e
- oarg
.s
+ 1);
410 return AVERROR(ENOMEM
);
411 memcpy(tptr
, oarg
.s
, oarg
.e
- oarg
.s
);
412 tptr
[oarg
.e
- oarg
.s
] = 0;
413 av_free(p
->scs
.opt_mix
);
414 p
->scs
.opt_mix
= tptr
;
417 FORWARD_ERROR(parse_optarg(p
, opt
, &oarg
));
418 v
= strtod(oarg
.s
, &tptr
);
419 if (oarg
.e
!= tptr
) {
420 snprintf(p
->err_msg
, sizeof(p
->err_msg
),
421 "syntax error for option -q");
422 return AVERROR_INVALIDDATA
;
425 snprintf(p
->err_msg
, sizeof(p
->err_msg
),
426 "speed factor other than 1 not supported");
427 return AVERROR_PATCHWELCOME
;
431 FORWARD_ERROR(parse_optarg(p
, opt
, &oarg
));
432 r
= strtol(oarg
.s
, &tptr
, 10);
433 if (oarg
.e
!= tptr
) {
434 snprintf(p
->err_msg
, sizeof(p
->err_msg
),
435 "syntax error for option -r");
436 return AVERROR_INVALIDDATA
;
439 snprintf(p
->err_msg
, sizeof(p
->err_msg
),
440 "invalid sample rate");
441 return AVERROR_PATCHWELCOME
;
443 p
->scs
.sample_rate
= r
;
446 snprintf(p
->err_msg
, sizeof(p
->err_msg
),
447 "unknown option: '%c'", *ostr
.s
);
448 return AVERROR_INVALIDDATA
;
454 return parse_immediate(p
);
456 return parse_preprogrammed(p
);
458 if (!lex_line_end(p
))
459 return AVERROR_INVALIDDATA
;
465 static int parse_timestamp(struct sbg_parser
*p
,
466 struct sbg_timestamp
*rts
, int64_t *rrel
)
468 int64_t abs
= 0, rel
= 0, dt
;
472 if (lex_fixed(p
, "NOW", 3)) {
476 r
= lex_time(p
, &abs
);
480 while (lex_char(p
, '+')) {
481 if (!lex_time(p
, &dt
))
482 return AVERROR_INVALIDDATA
;
483 if (av_sat_add64(rel
, dt
) - dt
!= rel
)
484 return AVERROR_INVALIDDATA
;
490 return AVERROR_INVALIDDATA
;
498 static int parse_fade(struct sbg_parser
*p
, struct sbg_fade
*fr
)
500 struct sbg_fade f
= {0};
502 if (lex_char(p
, '<'))
503 f
.in
= SBG_FADE_SILENCE
;
504 else if (lex_char(p
, '-'))
505 f
.in
= SBG_FADE_SAME
;
506 else if (lex_char(p
, '='))
507 f
.in
= SBG_FADE_ADAPT
;
510 if (lex_char(p
, '>'))
511 f
.out
= SBG_FADE_SILENCE
;
512 else if (lex_char(p
, '-'))
513 f
.out
= SBG_FADE_SAME
;
514 else if (lex_char(p
, '='))
515 f
.out
= SBG_FADE_ADAPT
;
517 return AVERROR_INVALIDDATA
;
522 static int parse_time_sequence(struct sbg_parser
*p
, int inblock
)
524 struct sbg_timestamp ts
;
527 struct sbg_fade fade
= { SBG_FADE_SAME
, SBG_FADE_SAME
, 0 };
528 struct sbg_string name
;
529 struct sbg_script_tseq
*tseq
;
531 r
= parse_timestamp(p
, &ts
, &rel_ts
);
538 return AVERROR_INVALIDDATA
;
539 p
->current_time
.type
= ts
.type
;
540 p
->current_time
.t
= ts
.t
;
541 } else if(!inblock
&& !p
->current_time
.type
) {
542 snprintf(p
->err_msg
, sizeof(p
->err_msg
),
543 "relative time without previous absolute time");
544 return AVERROR_INVALIDDATA
;
546 ts
.type
= p
->current_time
.type
;
548 if (av_sat_add64(p
->current_time
.t
, rel_ts
) != p
->current_time
.t
+ (uint64_t)rel_ts
)
549 return AVERROR_INVALIDDATA
;
550 ts
.t
= p
->current_time
.t
+ rel_ts
;
551 r
= parse_fade(p
, &fade
);
555 if (!lex_name(p
, &name
))
556 return AVERROR_INVALIDDATA
;
558 if (lex_fixed(p
, "->", 2)) {
559 fade
.slide
= SBG_FADE_ADAPT
;
562 if (!lex_line_end(p
))
563 return AVERROR_INVALIDDATA
;
565 alloc_array_elem((void **)&p
->scs
.block_tseq
, sizeof(*tseq
),
566 &p
->nb_block_tseq
, &p
->nb_block_tseq_max
) :
567 alloc_array_elem((void **)&p
->scs
.tseq
, sizeof(*tseq
),
568 &p
->scs
.nb_tseq
, &p
->nb_tseq_max
);
570 return AVERROR(ENOMEM
);
573 tseq
->name_len
= name
.e
- name
.s
;
578 static int parse_wave_def(struct sbg_parser
*p
, int wavenum
)
580 snprintf(p
->err_msg
, sizeof(p
->err_msg
),
581 "waveform definitions not yet implemented");
582 return AVERROR_PATCHWELCOME
;
585 static int parse_block_def(struct sbg_parser
*p
,
586 struct sbg_script_definition
*def
)
591 if (!lex_line_end(p
))
592 return AVERROR_INVALIDDATA
;
593 tseq
= p
->nb_block_tseq
;
595 r
= parse_time_sequence(p
, 1);
601 if (!lex_char(p
, '}'))
602 return AVERROR_INVALIDDATA
;
604 if (!lex_line_end(p
))
605 return AVERROR_INVALIDDATA
;
607 def
->elements
= tseq
;
608 def
->nb_elements
= p
->nb_block_tseq
- tseq
;
609 if (!def
->nb_elements
)
610 return AVERROR_INVALIDDATA
;
614 static int parse_volume(struct sbg_parser
*p
, int *vol
)
618 if (!lex_char(p
, '/'))
620 if (!lex_double(p
, &v
))
621 return AVERROR_INVALIDDATA
;
622 if (scale_double(p
->log
, v
, 0.01, vol
))
623 return AVERROR(ERANGE
);
627 static int parse_synth_channel_sine(struct sbg_parser
*p
,
628 struct sbg_script_synth
*synth
)
630 double carrierf
, beatf
;
631 int carrier
, beat
, vol
;
633 if (!lex_double(p
, &carrierf
))
635 if (!lex_double(p
, &beatf
))
637 FORWARD_ERROR(parse_volume(p
, &vol
));
638 if (scale_double(p
->log
, carrierf
, 1, &carrier
) < 0 ||
639 scale_double(p
->log
, beatf
, 1, &beat
) < 0)
640 return AVERROR(EDOM
);
641 synth
->type
= SBG_TYPE_SINE
;
642 synth
->carrier
= carrier
;
648 static int parse_synth_channel_pink(struct sbg_parser
*p
,
649 struct sbg_script_synth
*synth
)
653 if (!lex_fixed(p
, "pink", 4))
655 FORWARD_ERROR(parse_volume(p
, &vol
));
656 synth
->type
= SBG_TYPE_NOISE
;
661 static int parse_synth_channel_bell(struct sbg_parser
*p
,
662 struct sbg_script_synth
*synth
)
667 if (!lex_fixed(p
, "bell", 4))
669 if (!lex_double(p
, &carrierf
))
670 return AVERROR_INVALIDDATA
;
671 FORWARD_ERROR(parse_volume(p
, &vol
));
672 if (scale_double(p
->log
, carrierf
, 1, &carrier
) < 0)
673 return AVERROR(EDOM
);
674 synth
->type
= SBG_TYPE_BELL
;
675 synth
->carrier
= carrier
;
680 static int parse_synth_channel_mix(struct sbg_parser
*p
,
681 struct sbg_script_synth
*synth
)
685 if (!lex_fixed(p
, "mix", 3))
687 FORWARD_ERROR(parse_volume(p
, &vol
));
688 synth
->type
= SBG_TYPE_MIX
;
693 static int parse_synth_channel_spin(struct sbg_parser
*p
,
694 struct sbg_script_synth
*synth
)
696 double carrierf
, beatf
;
697 int carrier
, beat
, vol
;
699 if (!lex_fixed(p
, "spin:", 5))
701 if (!lex_double(p
, &carrierf
))
702 return AVERROR_INVALIDDATA
;
703 if (!lex_double(p
, &beatf
))
704 return AVERROR_INVALIDDATA
;
705 FORWARD_ERROR(parse_volume(p
, &vol
));
706 if (scale_double(p
->log
, carrierf
, 1, &carrier
) < 0 ||
707 scale_double(p
->log
, beatf
, 1, &beat
) < 0)
708 return AVERROR(EDOM
);
709 synth
->type
= SBG_TYPE_SPIN
;
710 synth
->carrier
= carrier
;
716 static int parse_synth_channel(struct sbg_parser
*p
)
719 struct sbg_script_synth
*synth
;
721 synth
= alloc_array_elem((void **)&p
->scs
.synth
, sizeof(*synth
),
722 &p
->scs
.nb_synth
, &p
->nb_synth_max
);
724 return AVERROR(ENOMEM
);
725 r
= lex_char(p
, '-');
727 r
= parse_synth_channel_pink(p
, synth
);
729 r
= parse_synth_channel_bell(p
, synth
);
731 r
= parse_synth_channel_mix(p
, synth
);
733 r
= parse_synth_channel_spin(p
, synth
);
734 /* Unimplemented: wave%d:%f%f/vol (carrier, beat) */
736 r
= parse_synth_channel_sine(p
, synth
);
742 static int parse_synth_def(struct sbg_parser
*p
,
743 struct sbg_script_definition
*def
)
747 synth
= p
->scs
.nb_synth
;
749 r
= parse_synth_channel(p
);
752 if (!r
|| !lex_space(p
))
756 if (synth
== p
->scs
.nb_synth
)
757 return AVERROR_INVALIDDATA
;
758 if (!lex_line_end(p
))
759 return AVERROR_INVALIDDATA
;
761 def
->elements
= synth
;
762 def
->nb_elements
= p
->scs
.nb_synth
- synth
;
766 static int parse_named_def(struct sbg_parser
*p
)
768 char *cursor_save
= p
->cursor
;
769 struct sbg_string name
;
770 struct sbg_script_definition
*def
;
772 if (!lex_name(p
, &name
) || !lex_char(p
, ':') || !lex_space(p
)) {
773 p
->cursor
= cursor_save
;
776 if (name
.e
- name
.s
== 6 && !memcmp(name
.s
, "wave", 4) &&
777 name
.s
[4] >= '0' && name
.s
[4] <= '9' &&
778 name
.s
[5] >= '0' && name
.s
[5] <= '9') {
779 int wavenum
= (name
.s
[4] - '0') * 10 + (name
.s
[5] - '0');
780 return parse_wave_def(p
, wavenum
);
782 def
= alloc_array_elem((void **)&p
->scs
.def
, sizeof(*def
),
783 &p
->scs
.nb_def
, &p
->nb_def_max
);
785 return AVERROR(ENOMEM
);
787 def
->name_len
= name
.e
- name
.s
;
788 if (lex_char(p
, '{'))
789 return parse_block_def(p
, def
);
790 return parse_synth_def(p
, def
);
793 static void free_script(struct sbg_script
*s
)
798 av_freep(&s
->block_tseq
);
799 av_freep(&s
->events
);
800 av_freep(&s
->opt_mix
);
803 static int parse_script(void *log
, char *script
, int script_len
,
804 struct sbg_script
*rscript
)
806 struct sbg_parser sp
= {
809 .end
= script
+ script_len
,
815 .start_ts
= AV_NOPTS_VALUE
,
816 .sample_rate
= 44100,
817 .opt_fade_time
= 60 * AV_TIME_BASE
,
823 while (sp
.cursor
< sp
.end
) {
824 r
= parse_options(&sp
);
827 if (!r
&& !lex_line_end(&sp
))
830 while (sp
.cursor
< sp
.end
) {
831 r
= parse_named_def(&sp
);
833 r
= parse_time_sequence(&sp
, 0);
835 r
= lex_line_end(&sp
) ? 1 : AVERROR_INVALIDDATA
;
842 free_script(&sp
.scs
);
844 if (r
== AVERROR_INVALIDDATA
)
845 snprintf(sp
.err_msg
, sizeof(sp
.err_msg
), "syntax error");
846 if (log
&& *sp
.err_msg
) {
847 const char *ctx
= sp
.cursor
;
848 const char *ectx
= av_x_if_null(memchr(ctx
, '\n', sp
.end
- sp
.cursor
),
850 int lctx
= ectx
- ctx
;
851 const char *quote
= "\"";
852 if (lctx
> 0 && ctx
[lctx
- 1] == '\r')
855 ctx
= "the end of line";
859 av_log(log
, AV_LOG_ERROR
, "Error line %d: %s near %s%.*s%s.\n",
860 sp
.line_no
, sp
.err_msg
, quote
, lctx
, ctx
, quote
);
865 static int read_whole_file(AVIOContext
*io
, int max_size
, AVBPrint
*rbuf
)
867 int ret
= avio_read_to_bprint(io
, rbuf
, max_size
);
870 if (!av_bprint_is_complete(rbuf
))
871 return AVERROR(ENOMEM
);
872 /* Check if we have read the whole file. AVIOContext.eof_reached is only
873 * set after a read failed due to EOF, so this check is incorrect in case
874 * max_size equals the actual file size, but checking for that would
875 * require attempting to read beyond max_size. */
876 if (!io
->eof_reached
)
877 return AVERROR(EFBIG
);
881 static int expand_timestamps(void *log
, struct sbg_script
*s
)
884 int64_t now
, cur_ts
, delta
= 0;
886 for (i
= 0; i
< s
->nb_tseq
; i
++)
887 nb_rel
+= s
->tseq
[i
].ts
.type
== 'N';
888 if (nb_rel
== s
->nb_tseq
) {
889 /* All ts are relative to NOW: consider NOW = 0 */
891 if (s
->start_ts
!= AV_NOPTS_VALUE
)
892 av_log(log
, AV_LOG_WARNING
,
893 "Start time ignored in a purely relative script.\n");
894 } else if (nb_rel
== 0 && s
->start_ts
!= AV_NOPTS_VALUE
||
895 s
->opt_start_at_first
) {
896 /* All ts are absolute and start time is specified */
897 if (s
->start_ts
== AV_NOPTS_VALUE
)
898 s
->start_ts
= s
->tseq
[0].ts
.t
;
901 /* Mixed relative/absolute ts: expand */
903 struct tm
*tm
, tmpbuf
;
905 av_log(log
, AV_LOG_WARNING
,
906 "Scripts with mixed absolute and relative timestamps can give "
907 "unexpected results (pause, seeking, time zone change).\n");
910 tm
= localtime_r(&now0
, &tmpbuf
);
911 now
= tm
? tm
->tm_hour
* 3600 + tm
->tm_min
* 60 + tm
->tm_sec
:
913 av_log(log
, AV_LOG_INFO
, "Using %02d:%02d:%02d as NOW.\n",
914 (int)(now
/ 3600), (int)(now
/ 60) % 60, (int)now
% 60);
916 for (i
= 0; i
< s
->nb_tseq
; i
++) {
917 if (s
->tseq
[i
].ts
.type
== 'N') {
918 s
->tseq
[i
].ts
.t
+= now
;
919 s
->tseq
[i
].ts
.type
= 'T'; /* not necessary */
923 if (s
->start_ts
== AV_NOPTS_VALUE
)
924 s
->start_ts
= (s
->opt_start_at_first
&& s
->tseq
) ? s
->tseq
[0].ts
.t
: now
;
925 if (s
->start_ts
> INT64_MAX
- s
->opt_duration
)
926 return AVERROR_INVALIDDATA
;
928 s
->end_ts
= s
->opt_duration
? s
->start_ts
+ s
->opt_duration
:
929 AV_NOPTS_VALUE
; /* may be overridden later by -E option */
931 for (i
= 0; i
< s
->nb_tseq
; i
++) {
932 if (av_sat_add64(s
->tseq
[i
].ts
.t
, delta
) != s
->tseq
[i
].ts
.t
+ (uint64_t)delta
)
933 return AVERROR_INVALIDDATA
;
934 if (s
->tseq
[i
].ts
.t
+ delta
< cur_ts
)
936 cur_ts
= s
->tseq
[i
].ts
.t
+= delta
;
941 static int expand_tseq(void *log
, struct sbg_script
*s
, int *nb_ev_max
,
942 int64_t t0
, struct sbg_script_tseq
*tseq
)
945 struct sbg_script_definition
*def
;
946 struct sbg_script_tseq
*be
;
947 struct sbg_script_event
*ev
;
950 av_log(log
, AV_LOG_ERROR
, "Recursion loop on \"%.*s\"\n",
951 tseq
->name_len
, tseq
->name
);
952 return AVERROR(EINVAL
);
954 if (t0
+ (uint64_t)tseq
->ts
.t
!= av_sat_add64(t0
, tseq
->ts
.t
))
955 return AVERROR(EINVAL
);
958 for (i
= 0; i
< s
->nb_def
; i
++) {
959 if (s
->def
[i
].name_len
== tseq
->name_len
&&
960 !memcmp(s
->def
[i
].name
, tseq
->name
, tseq
->name_len
))
963 if (i
>= s
->nb_def
) {
964 av_log(log
, AV_LOG_ERROR
, "Tone-set \"%.*s\" not defined\n",
965 tseq
->name_len
, tseq
->name
);
966 return AVERROR(EINVAL
);
969 if (def
->type
== 'B') {
970 be
= s
->block_tseq
+ def
->elements
;
971 for (i
= 0; i
< def
->nb_elements
; i
++) {
972 r
= expand_tseq(log
, s
, nb_ev_max
, t0
, &be
[i
]);
977 ev
= alloc_array_elem((void **)&s
->events
, sizeof(*ev
),
978 &s
->nb_events
, nb_ev_max
);
980 return AVERROR(ENOMEM
);
982 ev
->elements
= def
->elements
;
983 ev
->nb_elements
= def
->nb_elements
;
984 ev
->fade
= tseq
->fade
;
990 static int expand_script(void *log
, struct sbg_script
*s
)
992 int i
, r
, nb_events_max
= 0;
994 r
= expand_timestamps(log
, s
);
997 for (i
= 0; i
< s
->nb_tseq
; i
++) {
998 r
= expand_tseq(log
, s
, &nb_events_max
, 0, &s
->tseq
[i
]);
1002 if (!s
->nb_events
) {
1003 av_log(log
, AV_LOG_ERROR
, "No events in script\n");
1004 return AVERROR_INVALIDDATA
;
1006 if (s
->opt_end_at_last
)
1007 s
->end_ts
= s
->events
[s
->nb_events
- 1].ts
;
1011 static int add_interval(struct ws_intervals
*inter
,
1012 enum ws_interval_type type
, uint32_t channels
, int ref
,
1013 int64_t ts1
, int32_t f1
, int32_t a1
,
1014 int64_t ts2
, int32_t f2
, int32_t a2
)
1016 struct ws_interval
*i
, *ri
;
1019 ri
= &inter
->inter
[ref
];
1020 /* ref and new intervals are constant, identical and adjacent */
1021 if (ri
->type
== type
&& ri
->channels
== channels
&&
1022 ri
->f1
== ri
->f2
&& ri
->f2
== f1
&& f1
== f2
&&
1023 ri
->a1
== ri
->a2
&& ri
->a2
== a1
&& a1
== a2
&&
1029 i
= alloc_array_elem((void **)&inter
->inter
, sizeof(*i
),
1030 &inter
->nb_inter
, &inter
->max_inter
);
1032 return AVERROR(ENOMEM
);
1036 i
->channels
= channels
;
1041 i
->phi
= ref
>= 0 ? ref
| 0x80000000 : 0;
1042 return i
- inter
->inter
;
1045 static int add_bell(struct ws_intervals
*inter
, struct sbg_script
*s
,
1046 int64_t ts1
, int64_t ts2
, int32_t f
, int32_t a
)
1048 /* SBaGen uses an exponential decrease every 50ms.
1049 We approximate it with piecewise affine segments. */
1050 int32_t cpoints
[][2] = {
1060 int64_t dt
= s
->sample_rate
/ 20, ts3
= ts1
, ts4
;
1061 for (i
= 0; i
< FF_ARRAY_ELEMS(cpoints
); i
++) {
1062 ts4
= FFMIN(ts2
, ts1
+ cpoints
[i
][0] * dt
);
1063 r
= add_interval(inter
, WS_SINE
, 3, -1,
1064 ts3
, f
, a
, ts4
, f
, cpoints
[i
][1]);
1073 static int generate_interval(void *log
, struct sbg_script
*s
,
1074 struct ws_intervals
*inter
,
1075 int64_t ts1
, int64_t ts2
,
1076 struct sbg_script_synth
*s1
,
1077 struct sbg_script_synth
*s2
,
1082 if (ts2
<= ts1
|| (s1
->vol
== 0 && s2
->vol
== 0))
1088 if (s1
->beat
== 0 && s2
->beat
== 0) {
1089 r
= add_interval(inter
, WS_SINE
, 3, s1
->ref
.l
,
1090 ts1
, s1
->carrier
, s1
->vol
,
1091 ts2
, s2
->carrier
, s2
->vol
);
1094 s2
->ref
.l
= s2
->ref
.r
= r
;
1096 r
= add_interval(inter
, WS_SINE
, 1, s1
->ref
.l
,
1097 ts1
, s1
->carrier
+ s1
->beat
/ 2, s1
->vol
,
1098 ts2
, s2
->carrier
+ s2
->beat
/ 2, s2
->vol
);
1102 r
= add_interval(inter
, WS_SINE
, 2, s1
->ref
.r
,
1103 ts1
, s1
->carrier
- s1
->beat
/ 2, s1
->vol
,
1104 ts2
, s2
->carrier
- s2
->beat
/ 2, s2
->vol
);
1112 if (transition
== 2) {
1113 r
= add_bell(inter
, s
, ts1
, ts2
, s1
->carrier
, s2
->vol
);
1120 av_log(log
, AV_LOG_WARNING
, "Spinning noise not implemented, "
1121 "using pink noise instead.\n");
1123 case SBG_TYPE_NOISE
:
1124 /* SBaGen's pink noise generator uses:
1125 - 1 band of white noise, mean square: 1/3;
1126 - 9 bands of subsampled white noise with linear
1127 interpolation, mean square: 2/3 each;
1128 with 1/10 weight each: the total mean square is 7/300.
1129 Our pink noise generator uses 8 bands of white noise with
1130 rectangular subsampling: the total mean square is 1/24.
1131 Therefore, to match SBaGen's volume, we must multiply vol by
1132 sqrt((7/300) / (1/24)) = sqrt(14/25) =~ 0.748
1134 r
= add_interval(inter
, WS_NOISE
, 3, s1
->ref
.l
,
1135 ts1
, 0, s1
->vol
- s1
->vol
/ 4,
1136 ts2
, 0, s2
->vol
- s2
->vol
/ 4);
1139 s2
->ref
.l
= s2
->ref
.r
= r
;
1143 /* Unimplemented: silence; warning present elsewhere */
1145 av_log(log
, AV_LOG_ERROR
,
1146 "Type %d is not implemented\n", s1
->type
);
1147 return AVERROR_PATCHWELCOME
;
1152 static int generate_plateau(void *log
, struct sbg_script
*s
,
1153 struct ws_intervals
*inter
,
1154 struct sbg_script_event
*ev1
)
1156 int64_t ts1
= ev1
->ts_int
, ts2
= ev1
->ts_trans
;
1158 struct sbg_script_synth
*s1
;
1160 for (i
= 0; i
< ev1
->nb_elements
; i
++) {
1161 s1
= &s
->synth
[ev1
->elements
+ i
];
1162 r
= generate_interval(log
, s
, inter
, ts1
, ts2
, s1
, s1
, 0);
1171 ts1 ts2 ts1 tsmid ts2
1176 ''''....____ ''....''
1178 compatible transition incompatible transition
1181 static int generate_transition(void *log
, struct sbg_script
*s
,
1182 struct ws_intervals
*inter
,
1183 struct sbg_script_event
*ev1
,
1184 struct sbg_script_event
*ev2
)
1186 int64_t ts1
= ev1
->ts_trans
, ts2
= ev1
->ts_next
;
1187 /* (ts1 + ts2) / 2 without overflow */
1188 int64_t tsmid
= (ts1
>> 1) + (ts2
>> 1) + (ts1
& ts2
& 1);
1189 enum sbg_fade_type type
= ev1
->fade
.slide
| (ev1
->fade
.out
& ev2
->fade
.in
);
1190 int nb_elements
= FFMAX(ev1
->nb_elements
, ev2
->nb_elements
);
1191 struct sbg_script_synth
*s1
, *s2
, s1mod
, s2mod
, smid
;
1194 for (pass
= 0; pass
< 2; pass
++) {
1195 /* pass = 0 -> compatible and first half of incompatible
1196 pass = 1 -> second half of incompatible
1197 Using two passes like that ensures that the intervals are generated
1198 in increasing order according to their start timestamp.
1199 Otherwise it would be necessary to sort them
1200 while keeping the mutual references.
1202 for (i
= 0; i
< nb_elements
; i
++) {
1203 s1
= i
< ev1
->nb_elements
? &s
->synth
[ev1
->elements
+ i
] : &s1mod
;
1204 s2
= i
< ev2
->nb_elements
? &s
->synth
[ev2
->elements
+ i
] : &s2mod
;
1205 s1mod
= s1
!= &s1mod
? *s1
: (struct sbg_script_synth
){ 0 };
1206 s2mod
= s2
!= &s2mod
? *s2
: (struct sbg_script_synth
){ 0 };
1207 if (ev1
->fade
.slide
) {
1208 /* for slides, and only for slides, silence ("-") is equivalent
1209 to anything with volume 0 */
1210 if (s1mod
.type
== SBG_TYPE_NONE
) {
1213 } else if (s2mod
.type
== SBG_TYPE_NONE
) {
1218 if (s1mod
.type
== s2mod
.type
&&
1219 s1mod
.type
!= SBG_TYPE_BELL
&&
1220 (type
== SBG_FADE_ADAPT
||
1221 (s1mod
.carrier
== s2mod
.carrier
&&
1222 s1mod
.beat
== s2mod
.beat
))) {
1223 /* compatible: single transition */
1225 r
= generate_interval(log
, s
, inter
,
1226 ts1
, ts2
, &s1mod
, &s2mod
, 3);
1229 s2
->ref
= s2mod
.ref
;
1232 /* incompatible: silence at midpoint */
1236 r
= generate_interval(log
, s
, inter
,
1237 ts1
, tsmid
, &s1mod
, &smid
, 1);
1243 r
= generate_interval(log
, s
, inter
,
1244 tsmid
, ts2
, &smid
, &s2mod
, 2);
1247 s2
->ref
= s2mod
.ref
;
1256 ev1 trats ev2 intts endts ev3
1261 '''....________________....''' '''...._______________
1263 \_________/\______________/\_________/\______________/\_________/\_____________/
1264 tr x->1 int1 tr 1->2 int2 tr 2->3 int3
1267 static int generate_intervals(void *log
, struct sbg_script
*s
, int sample_rate
,
1268 struct ws_intervals
*inter
)
1270 int64_t trans_time
= s
->opt_fade_time
/ 2;
1271 struct sbg_script_event ev0
, *ev1
, *ev2
;
1275 /* SBaGen handles the time before and after the extremal events,
1276 and the corresponding transitions, as if the sequence were cyclic
1277 with a 24-hours period. */
1278 period
= s
->events
[s
->nb_events
- 1].ts
- (uint64_t)s
->events
[0].ts
;
1280 return AVERROR_INVALIDDATA
;
1282 period
= (period
+ (DAY_TS
- 1)) / DAY_TS
* DAY_TS
;
1283 period
= FFMAX(period
, DAY_TS
);
1285 /* Prepare timestamps for transitions */
1286 for (i
= 0; i
< s
->nb_events
; i
++) {
1287 ev1
= &s
->events
[i
];
1288 ev2
= &s
->events
[(i
+ 1) % s
->nb_events
];
1289 ev1
->ts_int
= ev1
->ts
;
1291 if (!ev1
->fade
.slide
&& ev1
>= ev2
&& ev2
->ts
> INT64_MAX
- period
)
1292 return AVERROR_INVALIDDATA
;
1294 ev1
->ts_trans
= ev1
->fade
.slide
? ev1
->ts
1295 : ev2
->ts
+ (ev1
< ev2
? 0 : period
);
1297 for (i
= 0; i
< s
->nb_events
; i
++) {
1298 ev1
= &s
->events
[i
];
1299 ev2
= &s
->events
[(i
+ 1) % s
->nb_events
];
1300 if (!ev1
->fade
.slide
) {
1301 ev1
->ts_trans
= FFMAX(ev1
->ts_int
, ev1
->ts_trans
- trans_time
);
1302 ev2
->ts_int
= FFMIN(ev2
->ts_trans
, ev2
->ts_int
+ trans_time
);
1304 ev1
->ts_next
= ev2
->ts_int
+ (ev1
< ev2
? 0 : period
);
1307 /* Pseudo event before the first one */
1308 ev0
= s
->events
[s
->nb_events
- 1];
1309 if (av_sat_sub64(ev0
.ts_int
, period
) != (uint64_t)ev0
.ts_int
- period
)
1310 return AVERROR_INVALIDDATA
;
1311 ev0
.ts_int
-= period
;
1312 ev0
.ts_trans
-= period
;
1313 ev0
.ts_next
-= period
;
1315 /* Convert timestamps */
1316 for (i
= -1; i
< s
->nb_events
; i
++) {
1317 ev1
= i
< 0 ? &ev0
: &s
->events
[i
];
1318 ev1
->ts_int
= av_rescale(ev1
->ts_int
, sample_rate
, AV_TIME_BASE
);
1319 ev1
->ts_trans
= av_rescale(ev1
->ts_trans
, sample_rate
, AV_TIME_BASE
);
1320 ev1
->ts_next
= av_rescale(ev1
->ts_next
, sample_rate
, AV_TIME_BASE
);
1323 /* Generate intervals */
1324 for (i
= 0; i
< s
->nb_synth
; i
++)
1325 s
->synth
[i
].ref
.l
= s
->synth
[i
].ref
.r
= -1;
1326 for (i
= -1; i
< s
->nb_events
; i
++) {
1327 ev1
= i
< 0 ? &ev0
: &s
->events
[i
];
1328 ev2
= &s
->events
[(i
+ 1) % s
->nb_events
];
1329 r
= generate_plateau(log
, s
, inter
, ev1
);
1332 r
= generate_transition(log
, s
, inter
, ev1
, ev2
);
1336 if (!inter
->nb_inter
)
1337 av_log(log
, AV_LOG_WARNING
, "Completely silent script.\n");
1341 static int encode_intervals(struct sbg_script
*s
, AVCodecParameters
*par
,
1342 struct ws_intervals
*inter
)
1344 int i
, edata_size
= 4, ret
;
1347 for (i
= 0; i
< inter
->nb_inter
; i
++) {
1348 edata_size
+= inter
->inter
[i
].type
== WS_SINE
? 44 :
1349 inter
->inter
[i
].type
== WS_NOISE
? 32 : 0;
1351 return AVERROR(ENOMEM
);
1353 if ((ret
= ff_alloc_extradata(par
, edata_size
)) < 0)
1355 edata
= par
->extradata
;
1357 #define ADD_EDATA32(v) do { AV_WL32(edata, (v)); edata += 4; } while(0)
1358 #define ADD_EDATA64(v) do { AV_WL64(edata, (v)); edata += 8; } while(0)
1359 ADD_EDATA32(inter
->nb_inter
);
1360 for (i
= 0; i
< inter
->nb_inter
; i
++) {
1361 ADD_EDATA64(inter
->inter
[i
].ts1
);
1362 ADD_EDATA64(inter
->inter
[i
].ts2
);
1363 ADD_EDATA32(inter
->inter
[i
].type
);
1364 ADD_EDATA32(inter
->inter
[i
].channels
);
1365 switch (inter
->inter
[i
].type
) {
1367 ADD_EDATA32(inter
->inter
[i
].f1
);
1368 ADD_EDATA32(inter
->inter
[i
].f2
);
1369 ADD_EDATA32(inter
->inter
[i
].a1
);
1370 ADD_EDATA32(inter
->inter
[i
].a2
);
1371 ADD_EDATA32(inter
->inter
[i
].phi
);
1374 ADD_EDATA32(inter
->inter
[i
].a1
);
1375 ADD_EDATA32(inter
->inter
[i
].a2
);
1379 if (edata
!= par
->extradata
+ edata_size
)
1384 static av_cold
int sbg_read_probe(const AVProbeData
*p
)
1387 struct sbg_script script
= { 0 };
1389 r
= parse_script(NULL
, p
->buf
, p
->buf_size
, &script
);
1390 score
= r
< 0 || !script
.nb_def
|| !script
.nb_tseq
? 0 :
1391 AVPROBE_SCORE_MAX
/ 3;
1392 free_script(&script
);
1396 static av_cold
int sbg_read_header(AVFormatContext
*avf
)
1398 struct sbg_demuxer
*sbg
= avf
->priv_data
;
1401 struct sbg_script script
= { 0 };
1404 struct ws_intervals inter
= { 0 };
1406 av_bprint_init(&bprint
, 0, sbg
->max_file_size
+ 1U);
1407 r
= read_whole_file(avf
->pb
, sbg
->max_file_size
, &bprint
);
1411 r
= parse_script(avf
, bprint
.str
, bprint
.len
, &script
);
1414 if (!sbg
->sample_rate
)
1415 sbg
->sample_rate
= script
.sample_rate
;
1417 script
.sample_rate
= sbg
->sample_rate
;
1418 if (!sbg
->frame_size
)
1419 sbg
->frame_size
= FFMAX(1, sbg
->sample_rate
/ 10);
1421 av_log(avf
, AV_LOG_WARNING
, "Mix feature not implemented: "
1422 "-m is ignored and mix channels will be silent.\n");
1423 r
= expand_script(avf
, &script
);
1426 av_bprint_finalize(&bprint
, NULL
);
1427 r
= generate_intervals(avf
, &script
, sbg
->sample_rate
, &inter
);
1431 if (script
.end_ts
!= AV_NOPTS_VALUE
&& script
.end_ts
< script
.start_ts
) {
1432 r
= AVERROR_INVALIDDATA
;
1436 st
= avformat_new_stream(avf
, NULL
);
1438 return AVERROR(ENOMEM
);
1440 st
->codecpar
->codec_type
= AVMEDIA_TYPE_AUDIO
;
1441 st
->codecpar
->codec_id
= AV_CODEC_ID_FFWAVESYNTH
;
1442 st
->codecpar
->ch_layout
= (AVChannelLayout
)AV_CHANNEL_LAYOUT_STEREO
;
1443 st
->codecpar
->sample_rate
= sbg
->sample_rate
;
1444 st
->codecpar
->frame_size
= sbg
->frame_size
;
1445 avpriv_set_pts_info(st
, 64, 1, st
->codecpar
->sample_rate
);
1446 sti
->probe_packets
= 0;
1447 st
->start_time
= av_rescale(script
.start_ts
,
1448 sbg
->sample_rate
, AV_TIME_BASE
);
1449 st
->duration
= script
.end_ts
== AV_NOPTS_VALUE
? AV_NOPTS_VALUE
:
1450 av_rescale(script
.end_ts
- script
.start_ts
,
1451 sbg
->sample_rate
, AV_TIME_BASE
);
1453 if (st
->duration
!= AV_NOPTS_VALUE
&& (
1454 st
->duration
< 0 || st
->start_time
> INT64_MAX
- st
->duration
)) {
1455 r
= AVERROR_INVALIDDATA
;
1459 sti
->cur_dts
= st
->start_time
;
1460 r
= encode_intervals(&script
, st
->codecpar
, &inter
);
1464 av_free(inter
.inter
);
1465 free_script(&script
);
1469 av_bprint_finalize(&bprint
, NULL
);
1471 av_free(inter
.inter
);
1472 free_script(&script
);
1476 static int sbg_read_packet(AVFormatContext
*avf
, AVPacket
*packet
)
1481 ts
= ffstream(avf
->streams
[0])->cur_dts
;
1482 end_ts
= av_sat_add64(ts
, avf
->streams
[0]->codecpar
->frame_size
);
1483 if (avf
->streams
[0]->duration
!= AV_NOPTS_VALUE
)
1484 end_ts
= FFMIN(avf
->streams
[0]->start_time
+ avf
->streams
[0]->duration
,
1488 if ((ret
= av_new_packet(packet
, 12)) < 0)
1490 packet
->dts
= packet
->pts
= ts
;
1491 packet
->duration
= end_ts
- ts
;
1492 AV_WL64(packet
->data
+ 0, ts
);
1493 AV_WL32(packet
->data
+ 8, packet
->duration
);
1494 return packet
->size
;
1497 static int sbg_read_seek2(AVFormatContext
*avf
, int stream_index
,
1498 int64_t min_ts
, int64_t ts
, int64_t max_ts
, int flags
)
1500 if (flags
|| stream_index
> 0)
1501 return AVERROR(EINVAL
);
1502 if (stream_index
< 0)
1503 ts
= av_rescale_q(ts
, AV_TIME_BASE_Q
, avf
->streams
[0]->time_base
);
1504 ffstream(avf
->streams
[0])->cur_dts
= ts
;
1508 static int sbg_read_seek(AVFormatContext
*avf
, int stream_index
,
1509 int64_t ts
, int flags
)
1511 return sbg_read_seek2(avf
, stream_index
, ts
, ts
, ts
, 0);
1514 static const AVOption sbg_options
[] = {
1515 { "sample_rate", "", offsetof(struct sbg_demuxer
, sample_rate
),
1516 AV_OPT_TYPE_INT
, { .i64
= 0 }, 0, INT_MAX
,
1517 AV_OPT_FLAG_DECODING_PARAM
},
1518 { "frame_size", "", offsetof(struct sbg_demuxer
, frame_size
),
1519 AV_OPT_TYPE_INT
, { .i64
= 0 }, 0, INT_MAX
,
1520 AV_OPT_FLAG_DECODING_PARAM
},
1521 { "max_file_size", "", offsetof(struct sbg_demuxer
, max_file_size
),
1522 AV_OPT_TYPE_INT
, { .i64
= 5000000 }, 0, INT_MAX
,
1523 AV_OPT_FLAG_DECODING_PARAM
},
1527 static const AVClass sbg_demuxer_class
= {
1528 .class_name
= "sbg_demuxer",
1529 .item_name
= av_default_item_name
,
1530 .option
= sbg_options
,
1531 .version
= LIBAVUTIL_VERSION_INT
,
1534 const FFInputFormat ff_sbg_demuxer
= {
1536 .p
.long_name
= NULL_IF_CONFIG_SMALL("SBaGen binaural beats script"),
1537 .p
.extensions
= "sbg",
1538 .p
.priv_class
= &sbg_demuxer_class
,
1539 .priv_data_size
= sizeof(struct sbg_demuxer
),
1540 .read_probe
= sbg_read_probe
,
1541 .read_header
= sbg_read_header
,
1542 .read_packet
= sbg_read_packet
,
1543 .read_seek
= sbg_read_seek
,
1544 .read_seek2
= sbg_read_seek2
,