123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210211212213214215216217218219220221222223224225226227228229230231232233234235236237238239240241242243244245246247248249250251252253254255256257258259260261262263264265266267268269270271272273274275276277278279280281282283284285286287288289290291292293294295296297298299300301302303304305306307308309310311312313314315316317318319320321322323324325326327328329330331332333334335336337338339340341342343344345346347348349350351352353354355356357358359360361362363364365366367368369370371372373374375376377378379380381382383384385386387388389390391392393394395396397398399400401402403404405406407408409410411412413414415416417418419420421422423424425426427428429430431432433434435436437438439440441442443444445446447448449450451452453454455456457458459460461462463464465466467468469470471472473474475476477478479480481482483484485486487488489490491492493494495496497498499500501502503504505506507508509510511512513514515516517518519520521522523524525526527528529530531532533534535536537538539540541542543544545546547548549550551552553554555556557558559560561562563564565566567568569570571572573574575576577578579580581582583584585586587588589 |
-
- #include "coder.h"
- #include "assembly.h"
- #define RND_VAL (1 << (FBITS_OUT_IMDCT-1))
- #ifndef AAC_ENABLE_SBR
- static void DecWindowOverlap(int *buf0, int *over0, short *pcm0, int nChans, int winTypeCurr, int winTypePrev)
- {
- int in, w0, w1, f0, f1;
- int *buf1, *over1;
- short *pcm1;
- const int *wndPrev, *wndCurr;
- buf0 += (1024 >> 1);
- buf1 = buf0 - 1;
- pcm1 = pcm0 + (1024 - 1) * nChans;
- over1 = over0 + 1024 - 1;
- wndPrev = (winTypePrev == 1 ? kbdWindow + kbdWindowOffset[1] : sinWindow + sinWindowOffset[1]);
- if (winTypeCurr == winTypePrev) {
-
- do {
- w0 = *wndPrev++;
- w1 = *wndPrev++;
- in = *buf0++;
- f0 = MULSHIFT32(w0, in);
- f1 = MULSHIFT32(w1, in);
- in = *over0;
- *pcm0 = CLIPTOSHORT( (in - f0 + RND_VAL) >> FBITS_OUT_IMDCT );
- pcm0 += nChans;
- in = *over1;
- *pcm1 = CLIPTOSHORT( (in + f1 + RND_VAL) >> FBITS_OUT_IMDCT );
- pcm1 -= nChans;
- in = *buf1--;
- *over1-- = MULSHIFT32(w0, in);
- *over0++ = MULSHIFT32(w1, in);
- } while (over0 < over1);
- } else {
-
- wndCurr = (winTypeCurr == 1 ? kbdWindow + kbdWindowOffset[1] : sinWindow + sinWindowOffset[1]);
- do {
- w0 = *wndPrev++;
- w1 = *wndPrev++;
- in = *buf0++;
- f0 = MULSHIFT32(w0, in);
- f1 = MULSHIFT32(w1, in);
- in = *over0;
- *pcm0 = CLIPTOSHORT( (in - f0 + RND_VAL) >> FBITS_OUT_IMDCT );
- pcm0 += nChans;
- in = *over1;
- *pcm1 = CLIPTOSHORT( (in + f1 + RND_VAL) >> FBITS_OUT_IMDCT );
- pcm1 -= nChans;
- w0 = *wndCurr++;
- w1 = *wndCurr++;
- in = *buf1--;
- *over1-- = MULSHIFT32(w0, in);
- *over0++ = MULSHIFT32(w1, in);
- } while (over0 < over1);
- }
- }
- static void DecWindowOverlapLongStart(int *buf0, int *over0, short *pcm0, int nChans, int winTypeCurr, int winTypePrev)
- {
- int i, in, w0, w1, f0, f1;
- int *buf1, *over1;
- short *pcm1;
- const int *wndPrev, *wndCurr;
- buf0 += (1024 >> 1);
- buf1 = buf0 - 1;
- pcm1 = pcm0 + (1024 - 1) * nChans;
- over1 = over0 + 1024 - 1;
- wndPrev = (winTypePrev == 1 ? kbdWindow + kbdWindowOffset[1] : sinWindow + sinWindowOffset[1]);
- i = 448;
- do {
- w0 = *wndPrev++;
- w1 = *wndPrev++;
- in = *buf0++;
- f0 = MULSHIFT32(w0, in);
- f1 = MULSHIFT32(w1, in);
- in = *over0;
- *pcm0 = CLIPTOSHORT( (in - f0 + RND_VAL) >> FBITS_OUT_IMDCT );
- pcm0 += nChans;
- in = *over1;
- *pcm1 = CLIPTOSHORT( (in + f1 + RND_VAL) >> FBITS_OUT_IMDCT );
- pcm1 -= nChans;
- in = *buf1--;
- *over1-- = 0;
- *over0++ = in >> 1;
- } while (--i);
- wndCurr = (winTypeCurr == 1 ? kbdWindow + kbdWindowOffset[0] : sinWindow + sinWindowOffset[0]);
-
- do {
- w0 = *wndPrev++;
- w1 = *wndPrev++;
- in = *buf0++;
- f0 = MULSHIFT32(w0, in);
- f1 = MULSHIFT32(w1, in);
- in = *over0;
- *pcm0 = CLIPTOSHORT( (in - f0 + RND_VAL) >> FBITS_OUT_IMDCT );
- pcm0 += nChans;
- in = *over1;
- *pcm1 = CLIPTOSHORT( (in + f1 + RND_VAL) >> FBITS_OUT_IMDCT );
- pcm1 -= nChans;
- w0 = *wndCurr++;
- w1 = *wndCurr++;
- in = *buf1--;
- *over1-- = MULSHIFT32(w0, in);
- *over0++ = MULSHIFT32(w1, in);
- } while (over0 < over1);
- }
- static void DecWindowOverlapLongStop(int *buf0, int *over0, short *pcm0, int nChans, int winTypeCurr, int winTypePrev)
- {
- int i, in, w0, w1, f0, f1;
- int *buf1, *over1;
- short *pcm1;
- const int *wndPrev, *wndCurr;
- buf0 += (1024 >> 1);
- buf1 = buf0 - 1;
- pcm1 = pcm0 + (1024 - 1) * nChans;
- over1 = over0 + 1024 - 1;
- wndPrev = (winTypePrev == 1 ? kbdWindow + kbdWindowOffset[0] : sinWindow + sinWindowOffset[0]);
- wndCurr = (winTypeCurr == 1 ? kbdWindow + kbdWindowOffset[1] : sinWindow + sinWindowOffset[1]);
- i = 448;
- do {
-
-
- in = *buf0++;
- f1 = in >> 1;
- in = *over0;
- *pcm0 = CLIPTOSHORT( (in + RND_VAL) >> FBITS_OUT_IMDCT );
- pcm0 += nChans;
- in = *over1;
- *pcm1 = CLIPTOSHORT( (in + f1 + RND_VAL) >> FBITS_OUT_IMDCT );
- pcm1 -= nChans;
- w0 = *wndCurr++;
- w1 = *wndCurr++;
- in = *buf1--;
- *over1-- = MULSHIFT32(w0, in);
- *over0++ = MULSHIFT32(w1, in);
- } while (--i);
-
- do {
- w0 = *wndPrev++;
- w1 = *wndPrev++;
- in = *buf0++;
- f0 = MULSHIFT32(w0, in);
- f1 = MULSHIFT32(w1, in);
- in = *over0;
- *pcm0 = CLIPTOSHORT( (in - f0 + RND_VAL) >> FBITS_OUT_IMDCT );
- pcm0 += nChans;
- in = *over1;
- *pcm1 = CLIPTOSHORT( (in + f1 + RND_VAL) >> FBITS_OUT_IMDCT );
- pcm1 -= nChans;
- w0 = *wndCurr++;
- w1 = *wndCurr++;
- in = *buf1--;
- *over1-- = MULSHIFT32(w0, in);
- *over0++ = MULSHIFT32(w1, in);
- } while (over0 < over1);
- }
- static void DecWindowOverlapShort(int *buf0, int *over0, short *pcm0, int nChans, int winTypeCurr, int winTypePrev)
- {
- int i, in, w0, w1, f0, f1;
- int *buf1, *over1;
- short *pcm1;
- const int *wndPrev, *wndCurr;
- wndPrev = (winTypePrev == 1 ? kbdWindow + kbdWindowOffset[0] : sinWindow + sinWindowOffset[0]);
- wndCurr = (winTypeCurr == 1 ? kbdWindow + kbdWindowOffset[0] : sinWindow + sinWindowOffset[0]);
-
- i = 448;
- do {
- f0 = *over0++;
- f1 = *over0++;
- *pcm0 = CLIPTOSHORT( (f0 + RND_VAL) >> FBITS_OUT_IMDCT ); pcm0 += nChans;
- *pcm0 = CLIPTOSHORT( (f1 + RND_VAL) >> FBITS_OUT_IMDCT ); pcm0 += nChans;
- i -= 2;
- } while (i);
-
- pcm1 = pcm0 + (128 - 1) * nChans;
- over1 = over0 + 128 - 1;
- buf0 += 64;
- buf1 = buf0 - 1;
- do {
- w0 = *wndPrev++;
- w1 = *wndPrev++;
- in = *buf0++;
- f0 = MULSHIFT32(w0, in);
- f1 = MULSHIFT32(w1, in);
- in = *over0;
- *pcm0 = CLIPTOSHORT( (in - f0 + RND_VAL) >> FBITS_OUT_IMDCT );
- pcm0 += nChans;
- in = *over1;
- *pcm1 = CLIPTOSHORT( (in + f1 + RND_VAL) >> FBITS_OUT_IMDCT );
- pcm1 -= nChans;
- w0 = *wndCurr++;
- w1 = *wndCurr++;
- in = *buf1--;
-
- *over1-- = MULSHIFT32(w0, in);
- *over0++ = MULSHIFT32(w1, in);
- } while (over0 < over1);
-
- for (i = 0; i < 3; i++) {
- pcm0 += 64 * nChans;
- pcm1 = pcm0 + (128 - 1) * nChans;
- over0 += 64;
- over1 = over0 + 128 - 1;
- buf0 += 64;
- buf1 = buf0 - 1;
- wndCurr -= 128;
- do {
- w0 = *wndCurr++;
- w1 = *wndCurr++;
- in = *buf0++;
- f0 = MULSHIFT32(w0, in);
- f1 = MULSHIFT32(w1, in);
- in = *(over0 - 128);
- in += *(over0 + 0);
- *pcm0 = CLIPTOSHORT( (in - f0 + RND_VAL) >> FBITS_OUT_IMDCT );
- pcm0 += nChans;
- in = *(over1 - 128);
- in += *(over1 + 0);
- *pcm1 = CLIPTOSHORT( (in + f1 + RND_VAL) >> FBITS_OUT_IMDCT );
- pcm1 -= nChans;
-
- in = *buf1--;
- *over1-- = MULSHIFT32(w0, in);
- *over0++ = MULSHIFT32(w1, in);
- } while (over0 < over1);
- }
-
- pcm0 += 64 * nChans;
- over0 -= 832;
- over1 = over0 + 128 - 1;
- buf0 += 64;
- buf1 = buf0 - 1;
- wndCurr -= 128;
- do {
- w0 = *wndCurr++;
- w1 = *wndCurr++;
- in = *buf0++;
- f0 = MULSHIFT32(w0, in);
- f1 = MULSHIFT32(w1, in);
- in = *(over0 + 768);
- in += *(over0 + 896);
- *pcm0 = CLIPTOSHORT( (in - f0 + RND_VAL) >> FBITS_OUT_IMDCT );
- pcm0 += nChans;
- in = *(over1 + 768);
- *(over1 - 128) = in + f1;
- in = *buf1--;
- *over1-- = MULSHIFT32(w0, in);
- *over0++ = MULSHIFT32(w1, in);
- } while (over0 < over1);
-
-
-
-
- for (i = 0; i < 3; i++) {
- over0 += 64;
- over1 = over0 + 128 - 1;
- buf0 += 64;
- buf1 = buf0 - 1;
- wndCurr -= 128;
- do {
- w0 = *wndCurr++;
- w1 = *wndCurr++;
- in = *buf0++;
- f0 = MULSHIFT32(w0, in);
- f1 = MULSHIFT32(w1, in);
-
- *(over0 - 128) -= f0;
- *(over1 - 128)+= f1;
- in = *buf1--;
- *over1-- = MULSHIFT32(w0, in);
- *over0++ = MULSHIFT32(w1, in);
- } while (over0 < over1);
- }
-
- i = 448;
- over0 += 64;
- do {
- *over0++ = 0;
- *over0++ = 0;
- *over0++ = 0;
- *over0++ = 0;
- i -= 4;
- } while (i);
- }
- #endif
- int IMDCT(AACDecInfo *aacDecInfo, int ch, int chOut, short *outbuf)
- {
- int i;
- PSInfoBase *psi;
- ICSInfo *icsInfo;
-
- if (!aacDecInfo || !aacDecInfo->psInfoBase)
- return -1;
- psi = (PSInfoBase *)(aacDecInfo->psInfoBase);
- icsInfo = (ch == 1 && psi->commonWin == 1) ? &(psi->icsInfo[0]) : &(psi->icsInfo[ch]);
- outbuf += chOut;
-
- if (icsInfo->winSequence == 2) {
-
- for (i = 0; i < 8; i++)
- DCT4(0, psi->coef[ch] + i*128, psi->gbCurrent[ch]);
- } else {
-
- DCT4(1, psi->coef[ch], psi->gbCurrent[ch]);
- }
- #ifdef AAC_ENABLE_SBR
-
- if (icsInfo->winSequence == 0)
- DecWindowOverlapNoClip(psi->coef[ch], psi->overlap[chOut], psi->sbrWorkBuf[ch], icsInfo->winShape, psi->prevWinShape[chOut]);
- else if (icsInfo->winSequence == 1)
- DecWindowOverlapLongStartNoClip(psi->coef[ch], psi->overlap[chOut], psi->sbrWorkBuf[ch], icsInfo->winShape, psi->prevWinShape[chOut]);
- else if (icsInfo->winSequence == 2)
- DecWindowOverlapShortNoClip(psi->coef[ch], psi->overlap[chOut], psi->sbrWorkBuf[ch], icsInfo->winShape, psi->prevWinShape[chOut]);
- else if (icsInfo->winSequence == 3)
- DecWindowOverlapLongStopNoClip(psi->coef[ch], psi->overlap[chOut], psi->sbrWorkBuf[ch], icsInfo->winShape, psi->prevWinShape[chOut]);
- if (!aacDecInfo->sbrEnabled) {
- for (i = 0; i < AAC_MAX_NSAMPS; i++) {
- *outbuf = CLIPTOSHORT((psi->sbrWorkBuf[ch][i] + RND_VAL) >> FBITS_OUT_IMDCT);
- outbuf += aacDecInfo->nChans;
- }
- }
- aacDecInfo->rawSampleBuf[ch] = psi->sbrWorkBuf[ch];
- aacDecInfo->rawSampleBytes = sizeof(int);
- aacDecInfo->rawSampleFBits = FBITS_OUT_IMDCT;
- #else
-
- if (icsInfo->winSequence == 0)
- DecWindowOverlap(psi->coef[ch], psi->overlap[chOut], outbuf, aacDecInfo->nChans, icsInfo->winShape, psi->prevWinShape[chOut]);
- else if (icsInfo->winSequence == 1)
- DecWindowOverlapLongStart(psi->coef[ch], psi->overlap[chOut], outbuf, aacDecInfo->nChans, icsInfo->winShape, psi->prevWinShape[chOut]);
- else if (icsInfo->winSequence == 2)
- DecWindowOverlapShort(psi->coef[ch], psi->overlap[chOut], outbuf, aacDecInfo->nChans, icsInfo->winShape, psi->prevWinShape[chOut]);
- else if (icsInfo->winSequence == 3)
- DecWindowOverlapLongStop(psi->coef[ch], psi->overlap[chOut], outbuf, aacDecInfo->nChans, icsInfo->winShape, psi->prevWinShape[chOut]);
- aacDecInfo->rawSampleBuf[ch] = 0;
- aacDecInfo->rawSampleBytes = 0;
- aacDecInfo->rawSampleFBits = 0;
- #endif
- psi->prevWinShape[chOut] = icsInfo->winShape;
- return 0;
- }
|