You can not select more than 25 topics Topics must start with a letter or number, can include dashes ('-') and can be up to 35 characters long.

1127 lines
40KB

  1. /*
  2. * Copyright (c) 2001-2003 The ffmpeg Project
  3. *
  4. * This file is part of FFmpeg.
  5. *
  6. * FFmpeg is free software; you can redistribute it and/or
  7. * modify it under the terms of the GNU Lesser General Public
  8. * License as published by the Free Software Foundation; either
  9. * version 2.1 of the License, or (at your option) any later version.
  10. *
  11. * FFmpeg 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 GNU
  14. * Lesser General Public License for more details.
  15. *
  16. * You should have received a copy of the GNU Lesser General Public
  17. * License along with FFmpeg; if not, write to the Free Software
  18. * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
  19. */
  20. #include "avcodec.h"
  21. #include "get_bits.h"
  22. #include "put_bits.h"
  23. #include "bytestream.h"
  24. #include "adpcm.h"
  25. #include "adpcm_data.h"
  26. /**
  27. * @file
  28. * ADPCM decoders
  29. * First version by Francois Revol (revol@free.fr)
  30. * Fringe ADPCM codecs (e.g., DK3, DK4, Westwood)
  31. * by Mike Melanson (melanson@pcisys.net)
  32. * CD-ROM XA ADPCM codec by BERO
  33. * EA ADPCM decoder by Robin Kay (komadori@myrealbox.com)
  34. * EA ADPCM R1/R2/R3 decoder by Peter Ross (pross@xvid.org)
  35. * EA IMA EACS decoder by Peter Ross (pross@xvid.org)
  36. * EA IMA SEAD decoder by Peter Ross (pross@xvid.org)
  37. * EA ADPCM XAS decoder by Peter Ross (pross@xvid.org)
  38. * MAXIS EA ADPCM decoder by Robert Marston (rmarston@gmail.com)
  39. * THP ADPCM decoder by Marco Gerards (mgerards@xs4all.nl)
  40. *
  41. * Features and limitations:
  42. *
  43. * Reference documents:
  44. * http://www.pcisys.net/~melanson/codecs/simpleaudio.html
  45. * http://www.geocities.com/SiliconValley/8682/aud3.txt
  46. * http://openquicktime.sourceforge.net/plugins.htm
  47. * XAnim sources (xa_codec.c) http://www.rasnaimaging.com/people/lapus/download.html
  48. * http://www.cs.ucla.edu/~leec/mediabench/applications.html
  49. * SoX source code http://home.sprynet.com/~cbagwell/sox.html
  50. *
  51. * CD-ROM XA:
  52. * http://ku-www.ss.titech.ac.jp/~yatsushi/xaadpcm.html
  53. * vagpack & depack http://homepages.compuserve.de/bITmASTER32/psx-index.html
  54. * readstr http://www.geocities.co.jp/Playtown/2004/
  55. */
  56. /* These are for CD-ROM XA ADPCM */
  57. static const int xa_adpcm_table[5][2] = {
  58. { 0, 0 },
  59. { 60, 0 },
  60. { 115, -52 },
  61. { 98, -55 },
  62. { 122, -60 }
  63. };
  64. static const int ea_adpcm_table[] = {
  65. 0, 240, 460, 392, 0, 0, -208, -220, 0, 1,
  66. 3, 4, 7, 8, 10, 11, 0, -1, -3, -4
  67. };
  68. // padded to zero where table size is less then 16
  69. static const int swf_index_tables[4][16] = {
  70. /*2*/ { -1, 2 },
  71. /*3*/ { -1, -1, 2, 4 },
  72. /*4*/ { -1, -1, -1, -1, 2, 4, 6, 8 },
  73. /*5*/ { -1, -1, -1, -1, -1, -1, -1, -1, 1, 2, 4, 6, 8, 10, 13, 16 }
  74. };
  75. /* end of tables */
  76. typedef struct ADPCMDecodeContext {
  77. ADPCMChannelStatus status[6];
  78. } ADPCMDecodeContext;
  79. static av_cold int adpcm_decode_init(AVCodecContext * avctx)
  80. {
  81. ADPCMDecodeContext *c = avctx->priv_data;
  82. unsigned int max_channels = 2;
  83. switch(avctx->codec->id) {
  84. case CODEC_ID_ADPCM_EA_R1:
  85. case CODEC_ID_ADPCM_EA_R2:
  86. case CODEC_ID_ADPCM_EA_R3:
  87. case CODEC_ID_ADPCM_EA_XAS:
  88. max_channels = 6;
  89. break;
  90. }
  91. if(avctx->channels > max_channels){
  92. return -1;
  93. }
  94. switch(avctx->codec->id) {
  95. case CODEC_ID_ADPCM_CT:
  96. c->status[0].step = c->status[1].step = 511;
  97. break;
  98. case CODEC_ID_ADPCM_IMA_WAV:
  99. if (avctx->bits_per_coded_sample != 4) {
  100. av_log(avctx, AV_LOG_ERROR, "Only 4-bit ADPCM IMA WAV files are supported\n");
  101. return -1;
  102. }
  103. break;
  104. case CODEC_ID_ADPCM_IMA_WS:
  105. if (avctx->extradata && avctx->extradata_size == 2 * 4) {
  106. c->status[0].predictor = AV_RL32(avctx->extradata);
  107. c->status[1].predictor = AV_RL32(avctx->extradata + 4);
  108. }
  109. break;
  110. default:
  111. break;
  112. }
  113. avctx->sample_fmt = AV_SAMPLE_FMT_S16;
  114. return 0;
  115. }
  116. static inline short adpcm_ima_expand_nibble(ADPCMChannelStatus *c, char nibble, int shift)
  117. {
  118. int step_index;
  119. int predictor;
  120. int sign, delta, diff, step;
  121. step = ff_adpcm_step_table[c->step_index];
  122. step_index = c->step_index + ff_adpcm_index_table[(unsigned)nibble];
  123. if (step_index < 0) step_index = 0;
  124. else if (step_index > 88) step_index = 88;
  125. sign = nibble & 8;
  126. delta = nibble & 7;
  127. /* perform direct multiplication instead of series of jumps proposed by
  128. * the reference ADPCM implementation since modern CPUs can do the mults
  129. * quickly enough */
  130. diff = ((2 * delta + 1) * step) >> shift;
  131. predictor = c->predictor;
  132. if (sign) predictor -= diff;
  133. else predictor += diff;
  134. c->predictor = av_clip_int16(predictor);
  135. c->step_index = step_index;
  136. return (short)c->predictor;
  137. }
  138. static inline int adpcm_ima_qt_expand_nibble(ADPCMChannelStatus *c, int nibble, int shift)
  139. {
  140. int step_index;
  141. int predictor;
  142. int diff, step;
  143. step = ff_adpcm_step_table[c->step_index];
  144. step_index = c->step_index + ff_adpcm_index_table[nibble];
  145. step_index = av_clip(step_index, 0, 88);
  146. diff = step >> 3;
  147. if (nibble & 4) diff += step;
  148. if (nibble & 2) diff += step >> 1;
  149. if (nibble & 1) diff += step >> 2;
  150. if (nibble & 8)
  151. predictor = c->predictor - diff;
  152. else
  153. predictor = c->predictor + diff;
  154. c->predictor = av_clip_int16(predictor);
  155. c->step_index = step_index;
  156. return c->predictor;
  157. }
  158. static inline short adpcm_ms_expand_nibble(ADPCMChannelStatus *c, char nibble)
  159. {
  160. int predictor;
  161. predictor = (((c->sample1) * (c->coeff1)) + ((c->sample2) * (c->coeff2))) / 64;
  162. predictor += (signed)((nibble & 0x08)?(nibble - 0x10):(nibble)) * c->idelta;
  163. c->sample2 = c->sample1;
  164. c->sample1 = av_clip_int16(predictor);
  165. c->idelta = (ff_adpcm_AdaptationTable[(int)nibble] * c->idelta) >> 8;
  166. if (c->idelta < 16) c->idelta = 16;
  167. return c->sample1;
  168. }
  169. static inline short adpcm_ct_expand_nibble(ADPCMChannelStatus *c, char nibble)
  170. {
  171. int sign, delta, diff;
  172. int new_step;
  173. sign = nibble & 8;
  174. delta = nibble & 7;
  175. /* perform direct multiplication instead of series of jumps proposed by
  176. * the reference ADPCM implementation since modern CPUs can do the mults
  177. * quickly enough */
  178. diff = ((2 * delta + 1) * c->step) >> 3;
  179. /* predictor update is not so trivial: predictor is multiplied on 254/256 before updating */
  180. c->predictor = ((c->predictor * 254) >> 8) + (sign ? -diff : diff);
  181. c->predictor = av_clip_int16(c->predictor);
  182. /* calculate new step and clamp it to range 511..32767 */
  183. new_step = (ff_adpcm_AdaptationTable[nibble & 7] * c->step) >> 8;
  184. c->step = av_clip(new_step, 511, 32767);
  185. return (short)c->predictor;
  186. }
  187. static inline short adpcm_sbpro_expand_nibble(ADPCMChannelStatus *c, char nibble, int size, int shift)
  188. {
  189. int sign, delta, diff;
  190. sign = nibble & (1<<(size-1));
  191. delta = nibble & ((1<<(size-1))-1);
  192. diff = delta << (7 + c->step + shift);
  193. /* clamp result */
  194. c->predictor = av_clip(c->predictor + (sign ? -diff : diff), -16384,16256);
  195. /* calculate new step */
  196. if (delta >= (2*size - 3) && c->step < 3)
  197. c->step++;
  198. else if (delta == 0 && c->step > 0)
  199. c->step--;
  200. return (short) c->predictor;
  201. }
  202. static inline short adpcm_yamaha_expand_nibble(ADPCMChannelStatus *c, unsigned char nibble)
  203. {
  204. if(!c->step) {
  205. c->predictor = 0;
  206. c->step = 127;
  207. }
  208. c->predictor += (c->step * ff_adpcm_yamaha_difflookup[nibble]) / 8;
  209. c->predictor = av_clip_int16(c->predictor);
  210. c->step = (c->step * ff_adpcm_yamaha_indexscale[nibble]) >> 8;
  211. c->step = av_clip(c->step, 127, 24567);
  212. return c->predictor;
  213. }
  214. static void xa_decode(short *out, const unsigned char *in,
  215. ADPCMChannelStatus *left, ADPCMChannelStatus *right, int inc)
  216. {
  217. int i, j;
  218. int shift,filter,f0,f1;
  219. int s_1,s_2;
  220. int d,s,t;
  221. for(i=0;i<4;i++) {
  222. shift = 12 - (in[4+i*2] & 15);
  223. filter = in[4+i*2] >> 4;
  224. f0 = xa_adpcm_table[filter][0];
  225. f1 = xa_adpcm_table[filter][1];
  226. s_1 = left->sample1;
  227. s_2 = left->sample2;
  228. for(j=0;j<28;j++) {
  229. d = in[16+i+j*4];
  230. t = (signed char)(d<<4)>>4;
  231. s = ( t<<shift ) + ((s_1*f0 + s_2*f1+32)>>6);
  232. s_2 = s_1;
  233. s_1 = av_clip_int16(s);
  234. *out = s_1;
  235. out += inc;
  236. }
  237. if (inc==2) { /* stereo */
  238. left->sample1 = s_1;
  239. left->sample2 = s_2;
  240. s_1 = right->sample1;
  241. s_2 = right->sample2;
  242. out = out + 1 - 28*2;
  243. }
  244. shift = 12 - (in[5+i*2] & 15);
  245. filter = in[5+i*2] >> 4;
  246. f0 = xa_adpcm_table[filter][0];
  247. f1 = xa_adpcm_table[filter][1];
  248. for(j=0;j<28;j++) {
  249. d = in[16+i+j*4];
  250. t = (signed char)d >> 4;
  251. s = ( t<<shift ) + ((s_1*f0 + s_2*f1+32)>>6);
  252. s_2 = s_1;
  253. s_1 = av_clip_int16(s);
  254. *out = s_1;
  255. out += inc;
  256. }
  257. if (inc==2) { /* stereo */
  258. right->sample1 = s_1;
  259. right->sample2 = s_2;
  260. out -= 1;
  261. } else {
  262. left->sample1 = s_1;
  263. left->sample2 = s_2;
  264. }
  265. }
  266. }
  267. /* DK3 ADPCM support macro */
  268. #define DK3_GET_NEXT_NIBBLE() \
  269. if (decode_top_nibble_next) \
  270. { \
  271. nibble = last_byte >> 4; \
  272. decode_top_nibble_next = 0; \
  273. } \
  274. else \
  275. { \
  276. last_byte = *src++; \
  277. if (src >= buf + buf_size) break; \
  278. nibble = last_byte & 0x0F; \
  279. decode_top_nibble_next = 1; \
  280. }
  281. static int adpcm_decode_frame(AVCodecContext *avctx,
  282. void *data, int *data_size,
  283. AVPacket *avpkt)
  284. {
  285. const uint8_t *buf = avpkt->data;
  286. int buf_size = avpkt->size;
  287. ADPCMDecodeContext *c = avctx->priv_data;
  288. ADPCMChannelStatus *cs;
  289. int n, m, channel, i;
  290. int block_predictor[2];
  291. short *samples;
  292. short *samples_end;
  293. const uint8_t *src;
  294. int st; /* stereo */
  295. /* DK3 ADPCM accounting variables */
  296. unsigned char last_byte = 0;
  297. unsigned char nibble;
  298. int decode_top_nibble_next = 0;
  299. int diff_channel;
  300. /* EA ADPCM state variables */
  301. uint32_t samples_in_chunk;
  302. int32_t previous_left_sample, previous_right_sample;
  303. int32_t current_left_sample, current_right_sample;
  304. int32_t next_left_sample, next_right_sample;
  305. int32_t coeff1l, coeff2l, coeff1r, coeff2r;
  306. uint8_t shift_left, shift_right;
  307. int count1, count2;
  308. int coeff[2][2], shift[2];//used in EA MAXIS ADPCM
  309. if (!buf_size)
  310. return 0;
  311. //should protect all 4bit ADPCM variants
  312. //8 is needed for CODEC_ID_ADPCM_IMA_WAV with 2 channels
  313. //
  314. if(*data_size/4 < buf_size + 8)
  315. return -1;
  316. samples = data;
  317. samples_end= samples + *data_size/2;
  318. *data_size= 0;
  319. src = buf;
  320. st = avctx->channels == 2 ? 1 : 0;
  321. switch(avctx->codec->id) {
  322. case CODEC_ID_ADPCM_IMA_QT:
  323. n = buf_size - 2*avctx->channels;
  324. for (channel = 0; channel < avctx->channels; channel++) {
  325. int16_t predictor;
  326. int step_index;
  327. cs = &(c->status[channel]);
  328. /* (pppppp) (piiiiiii) */
  329. /* Bits 15-7 are the _top_ 9 bits of the 16-bit initial predictor value */
  330. predictor = AV_RB16(src);
  331. step_index = predictor & 0x7F;
  332. predictor &= 0xFF80;
  333. src += 2;
  334. if (cs->step_index == step_index) {
  335. int diff = (int)predictor - cs->predictor;
  336. if (diff < 0)
  337. diff = - diff;
  338. if (diff > 0x7f)
  339. goto update;
  340. } else {
  341. update:
  342. cs->step_index = step_index;
  343. cs->predictor = predictor;
  344. }
  345. if (cs->step_index > 88){
  346. av_log(avctx, AV_LOG_ERROR, "ERROR: step_index = %i\n", cs->step_index);
  347. cs->step_index = 88;
  348. }
  349. samples = (short*)data + channel;
  350. for(m=32; n>0 && m>0; n--, m--) { /* in QuickTime, IMA is encoded by chuncks of 34 bytes (=64 samples) */
  351. *samples = adpcm_ima_qt_expand_nibble(cs, src[0] & 0x0F, 3);
  352. samples += avctx->channels;
  353. *samples = adpcm_ima_qt_expand_nibble(cs, src[0] >> 4 , 3);
  354. samples += avctx->channels;
  355. src ++;
  356. }
  357. }
  358. if (st)
  359. samples--;
  360. break;
  361. case CODEC_ID_ADPCM_IMA_WAV:
  362. if (avctx->block_align != 0 && buf_size > avctx->block_align)
  363. buf_size = avctx->block_align;
  364. // samples_per_block= (block_align-4*chanels)*8 / (bits_per_sample * chanels) + 1;
  365. for(i=0; i<avctx->channels; i++){
  366. cs = &(c->status[i]);
  367. cs->predictor = *samples++ = (int16_t)bytestream_get_le16(&src);
  368. cs->step_index = *src++;
  369. if (cs->step_index > 88){
  370. av_log(avctx, AV_LOG_ERROR, "ERROR: step_index = %i\n", cs->step_index);
  371. cs->step_index = 88;
  372. }
  373. if (*src++) av_log(avctx, AV_LOG_ERROR, "unused byte should be null but is %d!!\n", src[-1]); /* unused */
  374. }
  375. while(src < buf + buf_size){
  376. for(m=0; m<4; m++){
  377. for(i=0; i<=st; i++)
  378. *samples++ = adpcm_ima_expand_nibble(&c->status[i], src[4*i] & 0x0F, 3);
  379. for(i=0; i<=st; i++)
  380. *samples++ = adpcm_ima_expand_nibble(&c->status[i], src[4*i] >> 4 , 3);
  381. src++;
  382. }
  383. src += 4*st;
  384. }
  385. break;
  386. case CODEC_ID_ADPCM_4XM:
  387. cs = &(c->status[0]);
  388. c->status[0].predictor= (int16_t)bytestream_get_le16(&src);
  389. if(st){
  390. c->status[1].predictor= (int16_t)bytestream_get_le16(&src);
  391. }
  392. c->status[0].step_index= (int16_t)bytestream_get_le16(&src);
  393. if(st){
  394. c->status[1].step_index= (int16_t)bytestream_get_le16(&src);
  395. }
  396. if (cs->step_index < 0) cs->step_index = 0;
  397. if (cs->step_index > 88) cs->step_index = 88;
  398. m= (buf_size - (src - buf))>>st;
  399. for(i=0; i<m; i++) {
  400. *samples++ = adpcm_ima_expand_nibble(&c->status[0], src[i] & 0x0F, 4);
  401. if (st)
  402. *samples++ = adpcm_ima_expand_nibble(&c->status[1], src[i+m] & 0x0F, 4);
  403. *samples++ = adpcm_ima_expand_nibble(&c->status[0], src[i] >> 4, 4);
  404. if (st)
  405. *samples++ = adpcm_ima_expand_nibble(&c->status[1], src[i+m] >> 4, 4);
  406. }
  407. src += m<<st;
  408. break;
  409. case CODEC_ID_ADPCM_MS:
  410. if (avctx->block_align != 0 && buf_size > avctx->block_align)
  411. buf_size = avctx->block_align;
  412. n = buf_size - 7 * avctx->channels;
  413. if (n < 0)
  414. return -1;
  415. block_predictor[0] = av_clip(*src++, 0, 6);
  416. block_predictor[1] = 0;
  417. if (st)
  418. block_predictor[1] = av_clip(*src++, 0, 6);
  419. c->status[0].idelta = (int16_t)bytestream_get_le16(&src);
  420. if (st){
  421. c->status[1].idelta = (int16_t)bytestream_get_le16(&src);
  422. }
  423. c->status[0].coeff1 = ff_adpcm_AdaptCoeff1[block_predictor[0]];
  424. c->status[0].coeff2 = ff_adpcm_AdaptCoeff2[block_predictor[0]];
  425. c->status[1].coeff1 = ff_adpcm_AdaptCoeff1[block_predictor[1]];
  426. c->status[1].coeff2 = ff_adpcm_AdaptCoeff2[block_predictor[1]];
  427. c->status[0].sample1 = bytestream_get_le16(&src);
  428. if (st) c->status[1].sample1 = bytestream_get_le16(&src);
  429. c->status[0].sample2 = bytestream_get_le16(&src);
  430. if (st) c->status[1].sample2 = bytestream_get_le16(&src);
  431. *samples++ = c->status[0].sample2;
  432. if (st) *samples++ = c->status[1].sample2;
  433. *samples++ = c->status[0].sample1;
  434. if (st) *samples++ = c->status[1].sample1;
  435. for(;n>0;n--) {
  436. *samples++ = adpcm_ms_expand_nibble(&c->status[0 ], src[0] >> 4 );
  437. *samples++ = adpcm_ms_expand_nibble(&c->status[st], src[0] & 0x0F);
  438. src ++;
  439. }
  440. break;
  441. case CODEC_ID_ADPCM_IMA_DK4:
  442. if (avctx->block_align != 0 && buf_size > avctx->block_align)
  443. buf_size = avctx->block_align;
  444. c->status[0].predictor = (int16_t)bytestream_get_le16(&src);
  445. c->status[0].step_index = *src++;
  446. src++;
  447. *samples++ = c->status[0].predictor;
  448. if (st) {
  449. c->status[1].predictor = (int16_t)bytestream_get_le16(&src);
  450. c->status[1].step_index = *src++;
  451. src++;
  452. *samples++ = c->status[1].predictor;
  453. }
  454. while (src < buf + buf_size) {
  455. /* take care of the top nibble (always left or mono channel) */
  456. *samples++ = adpcm_ima_expand_nibble(&c->status[0],
  457. src[0] >> 4, 3);
  458. /* take care of the bottom nibble, which is right sample for
  459. * stereo, or another mono sample */
  460. if (st)
  461. *samples++ = adpcm_ima_expand_nibble(&c->status[1],
  462. src[0] & 0x0F, 3);
  463. else
  464. *samples++ = adpcm_ima_expand_nibble(&c->status[0],
  465. src[0] & 0x0F, 3);
  466. src++;
  467. }
  468. break;
  469. case CODEC_ID_ADPCM_IMA_DK3:
  470. if (avctx->block_align != 0 && buf_size > avctx->block_align)
  471. buf_size = avctx->block_align;
  472. if(buf_size + 16 > (samples_end - samples)*3/8)
  473. return -1;
  474. c->status[0].predictor = (int16_t)AV_RL16(src + 10);
  475. c->status[1].predictor = (int16_t)AV_RL16(src + 12);
  476. c->status[0].step_index = src[14];
  477. c->status[1].step_index = src[15];
  478. /* sign extend the predictors */
  479. src += 16;
  480. diff_channel = c->status[1].predictor;
  481. /* the DK3_GET_NEXT_NIBBLE macro issues the break statement when
  482. * the buffer is consumed */
  483. while (1) {
  484. /* for this algorithm, c->status[0] is the sum channel and
  485. * c->status[1] is the diff channel */
  486. /* process the first predictor of the sum channel */
  487. DK3_GET_NEXT_NIBBLE();
  488. adpcm_ima_expand_nibble(&c->status[0], nibble, 3);
  489. /* process the diff channel predictor */
  490. DK3_GET_NEXT_NIBBLE();
  491. adpcm_ima_expand_nibble(&c->status[1], nibble, 3);
  492. /* process the first pair of stereo PCM samples */
  493. diff_channel = (diff_channel + c->status[1].predictor) / 2;
  494. *samples++ = c->status[0].predictor + c->status[1].predictor;
  495. *samples++ = c->status[0].predictor - c->status[1].predictor;
  496. /* process the second predictor of the sum channel */
  497. DK3_GET_NEXT_NIBBLE();
  498. adpcm_ima_expand_nibble(&c->status[0], nibble, 3);
  499. /* process the second pair of stereo PCM samples */
  500. diff_channel = (diff_channel + c->status[1].predictor) / 2;
  501. *samples++ = c->status[0].predictor + c->status[1].predictor;
  502. *samples++ = c->status[0].predictor - c->status[1].predictor;
  503. }
  504. break;
  505. case CODEC_ID_ADPCM_IMA_ISS:
  506. c->status[0].predictor = (int16_t)AV_RL16(src + 0);
  507. c->status[0].step_index = src[2];
  508. src += 4;
  509. if(st) {
  510. c->status[1].predictor = (int16_t)AV_RL16(src + 0);
  511. c->status[1].step_index = src[2];
  512. src += 4;
  513. }
  514. while (src < buf + buf_size) {
  515. if (st) {
  516. *samples++ = adpcm_ima_expand_nibble(&c->status[0],
  517. src[0] >> 4 , 3);
  518. *samples++ = adpcm_ima_expand_nibble(&c->status[1],
  519. src[0] & 0x0F, 3);
  520. } else {
  521. *samples++ = adpcm_ima_expand_nibble(&c->status[0],
  522. src[0] & 0x0F, 3);
  523. *samples++ = adpcm_ima_expand_nibble(&c->status[0],
  524. src[0] >> 4 , 3);
  525. }
  526. src++;
  527. }
  528. break;
  529. case CODEC_ID_ADPCM_IMA_WS:
  530. /* no per-block initialization; just start decoding the data */
  531. while (src < buf + buf_size) {
  532. if (st) {
  533. *samples++ = adpcm_ima_expand_nibble(&c->status[0],
  534. src[0] >> 4 , 3);
  535. *samples++ = adpcm_ima_expand_nibble(&c->status[1],
  536. src[0] & 0x0F, 3);
  537. } else {
  538. *samples++ = adpcm_ima_expand_nibble(&c->status[0],
  539. src[0] >> 4 , 3);
  540. *samples++ = adpcm_ima_expand_nibble(&c->status[0],
  541. src[0] & 0x0F, 3);
  542. }
  543. src++;
  544. }
  545. break;
  546. case CODEC_ID_ADPCM_XA:
  547. while (buf_size >= 128) {
  548. xa_decode(samples, src, &c->status[0], &c->status[1],
  549. avctx->channels);
  550. src += 128;
  551. samples += 28 * 8;
  552. buf_size -= 128;
  553. }
  554. break;
  555. case CODEC_ID_ADPCM_IMA_EA_EACS:
  556. samples_in_chunk = bytestream_get_le32(&src) >> (1-st);
  557. if (samples_in_chunk > buf_size-4-(8<<st)) {
  558. src += buf_size - 4;
  559. break;
  560. }
  561. for (i=0; i<=st; i++)
  562. c->status[i].step_index = bytestream_get_le32(&src);
  563. for (i=0; i<=st; i++)
  564. c->status[i].predictor = bytestream_get_le32(&src);
  565. for (; samples_in_chunk; samples_in_chunk--, src++) {
  566. *samples++ = adpcm_ima_expand_nibble(&c->status[0], *src>>4, 3);
  567. *samples++ = adpcm_ima_expand_nibble(&c->status[st], *src&0x0F, 3);
  568. }
  569. break;
  570. case CODEC_ID_ADPCM_IMA_EA_SEAD:
  571. for (; src < buf+buf_size; src++) {
  572. *samples++ = adpcm_ima_expand_nibble(&c->status[0], src[0] >> 4, 6);
  573. *samples++ = adpcm_ima_expand_nibble(&c->status[st],src[0]&0x0F, 6);
  574. }
  575. break;
  576. case CODEC_ID_ADPCM_EA:
  577. /* Each EA ADPCM frame has a 12-byte header followed by 30-byte pieces,
  578. each coding 28 stereo samples. */
  579. if (buf_size < 12) {
  580. av_log(avctx, AV_LOG_ERROR, "frame too small\n");
  581. return AVERROR(EINVAL);
  582. }
  583. samples_in_chunk = AV_RL32(src);
  584. if (samples_in_chunk / 28 > (buf_size - 12) / 30) {
  585. av_log(avctx, AV_LOG_ERROR, "invalid frame\n");
  586. return AVERROR(EINVAL);
  587. }
  588. src += 4;
  589. current_left_sample = (int16_t)bytestream_get_le16(&src);
  590. previous_left_sample = (int16_t)bytestream_get_le16(&src);
  591. current_right_sample = (int16_t)bytestream_get_le16(&src);
  592. previous_right_sample = (int16_t)bytestream_get_le16(&src);
  593. for (count1 = 0; count1 < samples_in_chunk/28;count1++) {
  594. coeff1l = ea_adpcm_table[ *src >> 4 ];
  595. coeff2l = ea_adpcm_table[(*src >> 4 ) + 4];
  596. coeff1r = ea_adpcm_table[*src & 0x0F];
  597. coeff2r = ea_adpcm_table[(*src & 0x0F) + 4];
  598. src++;
  599. shift_left = (*src >> 4 ) + 8;
  600. shift_right = (*src & 0x0F) + 8;
  601. src++;
  602. for (count2 = 0; count2 < 28; count2++) {
  603. next_left_sample = (int32_t)((*src & 0xF0) << 24) >> shift_left;
  604. next_right_sample = (int32_t)((*src & 0x0F) << 28) >> shift_right;
  605. src++;
  606. next_left_sample = (next_left_sample +
  607. (current_left_sample * coeff1l) +
  608. (previous_left_sample * coeff2l) + 0x80) >> 8;
  609. next_right_sample = (next_right_sample +
  610. (current_right_sample * coeff1r) +
  611. (previous_right_sample * coeff2r) + 0x80) >> 8;
  612. previous_left_sample = current_left_sample;
  613. current_left_sample = av_clip_int16(next_left_sample);
  614. previous_right_sample = current_right_sample;
  615. current_right_sample = av_clip_int16(next_right_sample);
  616. *samples++ = (unsigned short)current_left_sample;
  617. *samples++ = (unsigned short)current_right_sample;
  618. }
  619. }
  620. if (src - buf == buf_size - 2)
  621. src += 2; // Skip terminating 0x0000
  622. break;
  623. case CODEC_ID_ADPCM_EA_MAXIS_XA:
  624. for(channel = 0; channel < avctx->channels; channel++) {
  625. for (i=0; i<2; i++)
  626. coeff[channel][i] = ea_adpcm_table[(*src >> 4) + 4*i];
  627. shift[channel] = (*src & 0x0F) + 8;
  628. src++;
  629. }
  630. for (count1 = 0; count1 < (buf_size - avctx->channels) / avctx->channels; count1++) {
  631. for(i = 4; i >= 0; i-=4) { /* Pairwise samples LL RR (st) or LL LL (mono) */
  632. for(channel = 0; channel < avctx->channels; channel++) {
  633. int32_t sample = (int32_t)(((*(src+channel) >> i) & 0x0F) << 0x1C) >> shift[channel];
  634. sample = (sample +
  635. c->status[channel].sample1 * coeff[channel][0] +
  636. c->status[channel].sample2 * coeff[channel][1] + 0x80) >> 8;
  637. c->status[channel].sample2 = c->status[channel].sample1;
  638. c->status[channel].sample1 = av_clip_int16(sample);
  639. *samples++ = c->status[channel].sample1;
  640. }
  641. }
  642. src+=avctx->channels;
  643. }
  644. break;
  645. case CODEC_ID_ADPCM_EA_R1:
  646. case CODEC_ID_ADPCM_EA_R2:
  647. case CODEC_ID_ADPCM_EA_R3: {
  648. /* channel numbering
  649. 2chan: 0=fl, 1=fr
  650. 4chan: 0=fl, 1=rl, 2=fr, 3=rr
  651. 6chan: 0=fl, 1=c, 2=fr, 3=rl, 4=rr, 5=sub */
  652. const int big_endian = avctx->codec->id == CODEC_ID_ADPCM_EA_R3;
  653. int32_t previous_sample, current_sample, next_sample;
  654. int32_t coeff1, coeff2;
  655. uint8_t shift;
  656. unsigned int channel;
  657. uint16_t *samplesC;
  658. const uint8_t *srcC;
  659. const uint8_t *src_end = buf + buf_size;
  660. samples_in_chunk = (big_endian ? bytestream_get_be32(&src)
  661. : bytestream_get_le32(&src)) / 28;
  662. if (samples_in_chunk > UINT32_MAX/(28*avctx->channels) ||
  663. 28*samples_in_chunk*avctx->channels > samples_end-samples) {
  664. src += buf_size - 4;
  665. break;
  666. }
  667. for (channel=0; channel<avctx->channels; channel++) {
  668. int32_t offset = (big_endian ? bytestream_get_be32(&src)
  669. : bytestream_get_le32(&src))
  670. + (avctx->channels-channel-1) * 4;
  671. if ((offset < 0) || (offset >= src_end - src - 4)) break;
  672. srcC = src + offset;
  673. samplesC = samples + channel;
  674. if (avctx->codec->id == CODEC_ID_ADPCM_EA_R1) {
  675. current_sample = (int16_t)bytestream_get_le16(&srcC);
  676. previous_sample = (int16_t)bytestream_get_le16(&srcC);
  677. } else {
  678. current_sample = c->status[channel].predictor;
  679. previous_sample = c->status[channel].prev_sample;
  680. }
  681. for (count1=0; count1<samples_in_chunk; count1++) {
  682. if (*srcC == 0xEE) { /* only seen in R2 and R3 */
  683. srcC++;
  684. if (srcC > src_end - 30*2) break;
  685. current_sample = (int16_t)bytestream_get_be16(&srcC);
  686. previous_sample = (int16_t)bytestream_get_be16(&srcC);
  687. for (count2=0; count2<28; count2++) {
  688. *samplesC = (int16_t)bytestream_get_be16(&srcC);
  689. samplesC += avctx->channels;
  690. }
  691. } else {
  692. coeff1 = ea_adpcm_table[ *srcC>>4 ];
  693. coeff2 = ea_adpcm_table[(*srcC>>4) + 4];
  694. shift = (*srcC++ & 0x0F) + 8;
  695. if (srcC > src_end - 14) break;
  696. for (count2=0; count2<28; count2++) {
  697. if (count2 & 1)
  698. next_sample = (int32_t)((*srcC++ & 0x0F) << 28) >> shift;
  699. else
  700. next_sample = (int32_t)((*srcC & 0xF0) << 24) >> shift;
  701. next_sample += (current_sample * coeff1) +
  702. (previous_sample * coeff2);
  703. next_sample = av_clip_int16(next_sample >> 8);
  704. previous_sample = current_sample;
  705. current_sample = next_sample;
  706. *samplesC = current_sample;
  707. samplesC += avctx->channels;
  708. }
  709. }
  710. }
  711. if (avctx->codec->id != CODEC_ID_ADPCM_EA_R1) {
  712. c->status[channel].predictor = current_sample;
  713. c->status[channel].prev_sample = previous_sample;
  714. }
  715. }
  716. src = src + buf_size - (4 + 4*avctx->channels);
  717. samples += 28 * samples_in_chunk * avctx->channels;
  718. break;
  719. }
  720. case CODEC_ID_ADPCM_EA_XAS:
  721. if (samples_end-samples < 32*4*avctx->channels
  722. || buf_size < (4+15)*4*avctx->channels) {
  723. src += buf_size;
  724. break;
  725. }
  726. for (channel=0; channel<avctx->channels; channel++) {
  727. int coeff[2][4], shift[4];
  728. short *s2, *s = &samples[channel];
  729. for (n=0; n<4; n++, s+=32*avctx->channels) {
  730. for (i=0; i<2; i++)
  731. coeff[i][n] = ea_adpcm_table[(src[0]&0x0F)+4*i];
  732. shift[n] = (src[2]&0x0F) + 8;
  733. for (s2=s, i=0; i<2; i++, src+=2, s2+=avctx->channels)
  734. s2[0] = (src[0]&0xF0) + (src[1]<<8);
  735. }
  736. for (m=2; m<32; m+=2) {
  737. s = &samples[m*avctx->channels + channel];
  738. for (n=0; n<4; n++, src++, s+=32*avctx->channels) {
  739. for (s2=s, i=0; i<8; i+=4, s2+=avctx->channels) {
  740. int level = (int32_t)((*src & (0xF0>>i)) << (24+i)) >> shift[n];
  741. int pred = s2[-1*avctx->channels] * coeff[0][n]
  742. + s2[-2*avctx->channels] * coeff[1][n];
  743. s2[0] = av_clip_int16((level + pred + 0x80) >> 8);
  744. }
  745. }
  746. }
  747. }
  748. samples += 32*4*avctx->channels;
  749. break;
  750. case CODEC_ID_ADPCM_IMA_AMV:
  751. case CODEC_ID_ADPCM_IMA_SMJPEG:
  752. c->status[0].predictor = (int16_t)bytestream_get_le16(&src);
  753. c->status[0].step_index = bytestream_get_le16(&src);
  754. if (avctx->codec->id == CODEC_ID_ADPCM_IMA_AMV)
  755. src+=4;
  756. while (src < buf + buf_size) {
  757. char hi, lo;
  758. lo = *src & 0x0F;
  759. hi = *src >> 4;
  760. if (avctx->codec->id == CODEC_ID_ADPCM_IMA_AMV)
  761. FFSWAP(char, hi, lo);
  762. *samples++ = adpcm_ima_expand_nibble(&c->status[0],
  763. lo, 3);
  764. *samples++ = adpcm_ima_expand_nibble(&c->status[0],
  765. hi, 3);
  766. src++;
  767. }
  768. break;
  769. case CODEC_ID_ADPCM_CT:
  770. while (src < buf + buf_size) {
  771. if (st) {
  772. *samples++ = adpcm_ct_expand_nibble(&c->status[0],
  773. src[0] >> 4);
  774. *samples++ = adpcm_ct_expand_nibble(&c->status[1],
  775. src[0] & 0x0F);
  776. } else {
  777. *samples++ = adpcm_ct_expand_nibble(&c->status[0],
  778. src[0] >> 4);
  779. *samples++ = adpcm_ct_expand_nibble(&c->status[0],
  780. src[0] & 0x0F);
  781. }
  782. src++;
  783. }
  784. break;
  785. case CODEC_ID_ADPCM_SBPRO_4:
  786. case CODEC_ID_ADPCM_SBPRO_3:
  787. case CODEC_ID_ADPCM_SBPRO_2:
  788. if (!c->status[0].step_index) {
  789. /* the first byte is a raw sample */
  790. *samples++ = 128 * (*src++ - 0x80);
  791. if (st)
  792. *samples++ = 128 * (*src++ - 0x80);
  793. c->status[0].step_index = 1;
  794. }
  795. if (avctx->codec->id == CODEC_ID_ADPCM_SBPRO_4) {
  796. while (src < buf + buf_size) {
  797. *samples++ = adpcm_sbpro_expand_nibble(&c->status[0],
  798. src[0] >> 4, 4, 0);
  799. *samples++ = adpcm_sbpro_expand_nibble(&c->status[st],
  800. src[0] & 0x0F, 4, 0);
  801. src++;
  802. }
  803. } else if (avctx->codec->id == CODEC_ID_ADPCM_SBPRO_3) {
  804. while (src < buf + buf_size && samples + 2 < samples_end) {
  805. *samples++ = adpcm_sbpro_expand_nibble(&c->status[0],
  806. src[0] >> 5 , 3, 0);
  807. *samples++ = adpcm_sbpro_expand_nibble(&c->status[0],
  808. (src[0] >> 2) & 0x07, 3, 0);
  809. *samples++ = adpcm_sbpro_expand_nibble(&c->status[0],
  810. src[0] & 0x03, 2, 0);
  811. src++;
  812. }
  813. } else {
  814. while (src < buf + buf_size && samples + 3 < samples_end) {
  815. *samples++ = adpcm_sbpro_expand_nibble(&c->status[0],
  816. src[0] >> 6 , 2, 2);
  817. *samples++ = adpcm_sbpro_expand_nibble(&c->status[st],
  818. (src[0] >> 4) & 0x03, 2, 2);
  819. *samples++ = adpcm_sbpro_expand_nibble(&c->status[0],
  820. (src[0] >> 2) & 0x03, 2, 2);
  821. *samples++ = adpcm_sbpro_expand_nibble(&c->status[st],
  822. src[0] & 0x03, 2, 2);
  823. src++;
  824. }
  825. }
  826. break;
  827. case CODEC_ID_ADPCM_SWF:
  828. {
  829. GetBitContext gb;
  830. const int *table;
  831. int k0, signmask, nb_bits, count;
  832. int size = buf_size*8;
  833. init_get_bits(&gb, buf, size);
  834. //read bits & initial values
  835. nb_bits = get_bits(&gb, 2)+2;
  836. //av_log(NULL,AV_LOG_INFO,"nb_bits: %d\n", nb_bits);
  837. table = swf_index_tables[nb_bits-2];
  838. k0 = 1 << (nb_bits-2);
  839. signmask = 1 << (nb_bits-1);
  840. while (get_bits_count(&gb) <= size - 22*avctx->channels) {
  841. for (i = 0; i < avctx->channels; i++) {
  842. *samples++ = c->status[i].predictor = get_sbits(&gb, 16);
  843. c->status[i].step_index = get_bits(&gb, 6);
  844. }
  845. for (count = 0; get_bits_count(&gb) <= size - nb_bits*avctx->channels && count < 4095; count++) {
  846. int i;
  847. for (i = 0; i < avctx->channels; i++) {
  848. // similar to IMA adpcm
  849. int delta = get_bits(&gb, nb_bits);
  850. int step = ff_adpcm_step_table[c->status[i].step_index];
  851. long vpdiff = 0; // vpdiff = (delta+0.5)*step/4
  852. int k = k0;
  853. do {
  854. if (delta & k)
  855. vpdiff += step;
  856. step >>= 1;
  857. k >>= 1;
  858. } while(k);
  859. vpdiff += step;
  860. if (delta & signmask)
  861. c->status[i].predictor -= vpdiff;
  862. else
  863. c->status[i].predictor += vpdiff;
  864. c->status[i].step_index += table[delta & (~signmask)];
  865. c->status[i].step_index = av_clip(c->status[i].step_index, 0, 88);
  866. c->status[i].predictor = av_clip_int16(c->status[i].predictor);
  867. *samples++ = c->status[i].predictor;
  868. if (samples >= samples_end) {
  869. av_log(avctx, AV_LOG_ERROR, "allocated output buffer is too small\n");
  870. return -1;
  871. }
  872. }
  873. }
  874. }
  875. src += buf_size;
  876. break;
  877. }
  878. case CODEC_ID_ADPCM_YAMAHA:
  879. while (src < buf + buf_size) {
  880. if (st) {
  881. *samples++ = adpcm_yamaha_expand_nibble(&c->status[0],
  882. src[0] & 0x0F);
  883. *samples++ = adpcm_yamaha_expand_nibble(&c->status[1],
  884. src[0] >> 4 );
  885. } else {
  886. *samples++ = adpcm_yamaha_expand_nibble(&c->status[0],
  887. src[0] & 0x0F);
  888. *samples++ = adpcm_yamaha_expand_nibble(&c->status[0],
  889. src[0] >> 4 );
  890. }
  891. src++;
  892. }
  893. break;
  894. case CODEC_ID_ADPCM_THP:
  895. {
  896. int table[2][16];
  897. unsigned int samplecnt;
  898. int prev[2][2];
  899. int ch;
  900. if (buf_size < 80) {
  901. av_log(avctx, AV_LOG_ERROR, "frame too small\n");
  902. return -1;
  903. }
  904. src+=4;
  905. samplecnt = bytestream_get_be32(&src);
  906. for (i = 0; i < 32; i++)
  907. table[0][i] = (int16_t)bytestream_get_be16(&src);
  908. /* Initialize the previous sample. */
  909. for (i = 0; i < 4; i++)
  910. prev[0][i] = (int16_t)bytestream_get_be16(&src);
  911. if (samplecnt >= (samples_end - samples) / (st + 1)) {
  912. av_log(avctx, AV_LOG_ERROR, "allocated output buffer is too small\n");
  913. return -1;
  914. }
  915. for (ch = 0; ch <= st; ch++) {
  916. samples = (unsigned short *) data + ch;
  917. /* Read in every sample for this channel. */
  918. for (i = 0; i < samplecnt / 14; i++) {
  919. int index = (*src >> 4) & 7;
  920. unsigned int exp = 28 - (*src++ & 15);
  921. int factor1 = table[ch][index * 2];
  922. int factor2 = table[ch][index * 2 + 1];
  923. /* Decode 14 samples. */
  924. for (n = 0; n < 14; n++) {
  925. int32_t sampledat;
  926. if(n&1) sampledat= *src++ <<28;
  927. else sampledat= (*src&0xF0)<<24;
  928. sampledat = ((prev[ch][0]*factor1
  929. + prev[ch][1]*factor2) >> 11) + (sampledat>>exp);
  930. *samples = av_clip_int16(sampledat);
  931. prev[ch][1] = prev[ch][0];
  932. prev[ch][0] = *samples++;
  933. /* In case of stereo, skip one sample, this sample
  934. is for the other channel. */
  935. samples += st;
  936. }
  937. }
  938. }
  939. /* In the previous loop, in case stereo is used, samples is
  940. increased exactly one time too often. */
  941. samples -= st;
  942. break;
  943. }
  944. default:
  945. return -1;
  946. }
  947. *data_size = (uint8_t *)samples - (uint8_t *)data;
  948. return src - buf;
  949. }
  950. #define ADPCM_DECODER(id_, name_, long_name_) \
  951. AVCodec ff_ ## name_ ## _decoder = { \
  952. .name = #name_, \
  953. .type = AVMEDIA_TYPE_AUDIO, \
  954. .id = id_, \
  955. .priv_data_size = sizeof(ADPCMDecodeContext), \
  956. .init = adpcm_decode_init, \
  957. .decode = adpcm_decode_frame, \
  958. .long_name = NULL_IF_CONFIG_SMALL(long_name_), \
  959. }
  960. /* Note: Do not forget to add new entries to the Makefile as well. */
  961. ADPCM_DECODER(CODEC_ID_ADPCM_4XM, adpcm_4xm, "ADPCM 4X Movie");
  962. ADPCM_DECODER(CODEC_ID_ADPCM_CT, adpcm_ct, "ADPCM Creative Technology");
  963. ADPCM_DECODER(CODEC_ID_ADPCM_EA, adpcm_ea, "ADPCM Electronic Arts");
  964. ADPCM_DECODER(CODEC_ID_ADPCM_EA_MAXIS_XA, adpcm_ea_maxis_xa, "ADPCM Electronic Arts Maxis CDROM XA");
  965. ADPCM_DECODER(CODEC_ID_ADPCM_EA_R1, adpcm_ea_r1, "ADPCM Electronic Arts R1");
  966. ADPCM_DECODER(CODEC_ID_ADPCM_EA_R2, adpcm_ea_r2, "ADPCM Electronic Arts R2");
  967. ADPCM_DECODER(CODEC_ID_ADPCM_EA_R3, adpcm_ea_r3, "ADPCM Electronic Arts R3");
  968. ADPCM_DECODER(CODEC_ID_ADPCM_EA_XAS, adpcm_ea_xas, "ADPCM Electronic Arts XAS");
  969. ADPCM_DECODER(CODEC_ID_ADPCM_IMA_AMV, adpcm_ima_amv, "ADPCM IMA AMV");
  970. ADPCM_DECODER(CODEC_ID_ADPCM_IMA_DK3, adpcm_ima_dk3, "ADPCM IMA Duck DK3");
  971. ADPCM_DECODER(CODEC_ID_ADPCM_IMA_DK4, adpcm_ima_dk4, "ADPCM IMA Duck DK4");
  972. ADPCM_DECODER(CODEC_ID_ADPCM_IMA_EA_EACS, adpcm_ima_ea_eacs, "ADPCM IMA Electronic Arts EACS");
  973. ADPCM_DECODER(CODEC_ID_ADPCM_IMA_EA_SEAD, adpcm_ima_ea_sead, "ADPCM IMA Electronic Arts SEAD");
  974. ADPCM_DECODER(CODEC_ID_ADPCM_IMA_ISS, adpcm_ima_iss, "ADPCM IMA Funcom ISS");
  975. ADPCM_DECODER(CODEC_ID_ADPCM_IMA_QT, adpcm_ima_qt, "ADPCM IMA QuickTime");
  976. ADPCM_DECODER(CODEC_ID_ADPCM_IMA_SMJPEG, adpcm_ima_smjpeg, "ADPCM IMA Loki SDL MJPEG");
  977. ADPCM_DECODER(CODEC_ID_ADPCM_IMA_WAV, adpcm_ima_wav, "ADPCM IMA WAV");
  978. ADPCM_DECODER(CODEC_ID_ADPCM_IMA_WS, adpcm_ima_ws, "ADPCM IMA Westwood");
  979. ADPCM_DECODER(CODEC_ID_ADPCM_MS, adpcm_ms, "ADPCM Microsoft");
  980. ADPCM_DECODER(CODEC_ID_ADPCM_SBPRO_2, adpcm_sbpro_2, "ADPCM Sound Blaster Pro 2-bit");
  981. ADPCM_DECODER(CODEC_ID_ADPCM_SBPRO_3, adpcm_sbpro_3, "ADPCM Sound Blaster Pro 2.6-bit");
  982. ADPCM_DECODER(CODEC_ID_ADPCM_SBPRO_4, adpcm_sbpro_4, "ADPCM Sound Blaster Pro 4-bit");
  983. ADPCM_DECODER(CODEC_ID_ADPCM_SWF, adpcm_swf, "ADPCM Shockwave Flash");
  984. ADPCM_DECODER(CODEC_ID_ADPCM_THP, adpcm_thp, "ADPCM Nintendo Gamecube THP");
  985. ADPCM_DECODER(CODEC_ID_ADPCM_XA, adpcm_xa, "ADPCM CDROM XA");
  986. ADPCM_DECODER(CODEC_ID_ADPCM_YAMAHA, adpcm_yamaha, "ADPCM Yamaha");