int synth_1to1(struct mpstr *mp, real *bandPtr,int channel,unsigned char *out,int *pnt) { static const int step = 2; int bo; short *samples = (short *) (out + *pnt); real *b0,(*buf)[0x110]; int clip = 0; int bo1; bo = mp->synth_bo; if(!channel) { bo--; bo &= 0xf; buf = mp->synth_buffs[0]; } else { samples++; buf = mp->synth_buffs[1]; } if(bo & 0x1) { b0 = buf[0]; bo1 = bo; dct64(buf[1]+((bo+1)&0xf),buf[0]+bo,bandPtr); } else { b0 = buf[1]; bo1 = bo+1; dct64(buf[0]+bo,buf[1]+bo+1,bandPtr); } mp->synth_bo = bo; { register int j; real *window = decwin + 16 - bo1; for (j=16;j;j--,b0+=0x10,window+=0x20,samples+=step) { real sum; sum = window[0x0] * b0[0x0]; sum -= window[0x1] * b0[0x1]; sum += window[0x2] * b0[0x2]; sum -= window[0x3] * b0[0x3]; sum += window[0x4] * b0[0x4]; sum -= window[0x5] * b0[0x5]; sum += window[0x6] * b0[0x6]; sum -= window[0x7] * b0[0x7]; sum += window[0x8] * b0[0x8]; sum -= window[0x9] * b0[0x9]; sum += window[0xA] * b0[0xA]; sum -= window[0xB] * b0[0xB]; sum += window[0xC] * b0[0xC]; sum -= window[0xD] * b0[0xD]; sum += window[0xE] * b0[0xE]; sum -= window[0xF] * b0[0xF]; WRITE_SAMPLE(samples,sum,clip); } { real sum; sum = window[0x0] * b0[0x0]; sum += window[0x2] * b0[0x2]; sum += window[0x4] * b0[0x4]; sum += window[0x6] * b0[0x6]; sum += window[0x8] * b0[0x8]; sum += window[0xA] * b0[0xA]; sum += window[0xC] * b0[0xC]; sum += window[0xE] * b0[0xE]; WRITE_SAMPLE(samples,sum,clip); b0-=0x10,window-=0x20,samples+=step; } window += bo1<<1; for (j=15;j;j--,b0-=0x10,window-=0x20,samples+=step) { real sum; sum = -window[-0x1] * b0[0x0]; sum -= window[-0x2] * b0[0x1]; sum -= window[-0x3] * b0[0x2]; sum -= window[-0x4] * b0[0x3]; sum -= window[-0x5] * b0[0x4]; sum -= window[-0x6] * b0[0x5]; sum -= window[-0x7] * b0[0x6]; sum -= window[-0x8] * b0[0x7]; sum -= window[-0x9] * b0[0x8]; sum -= window[-0xA] * b0[0x9]; sum -= window[-0xB] * b0[0xA]; sum -= window[-0xC] * b0[0xB]; sum -= window[-0xD] * b0[0xC]; sum -= window[-0xE] * b0[0xD]; sum -= window[-0xF] * b0[0xE]; sum -= window[-0x0] * b0[0xF]; WRITE_SAMPLE(samples,sum,clip); } } *pnt += 128; return clip; }
static int synth_1to1(real *bandPtr,int channel,unsigned char *out,int *pnt) { static real buffs[2][2][0x110]; static const int step = 2; static int bo = 1; short *samples = (short *) (out + *pnt); real *b0,(*buf)[0x110]; int clip = 0; int bo1; *pnt += 128; /* optimized for x86 */ #if defined(CAN_COMPILE_X86_ASM) if ( synth_func ) { // printf("Calling %p, bandPtr=%p channel=%d samples=%p\n",synth_func,bandPtr,channel,samples); // FIXME: synth_func() may destroy EBP, don't rely on stack contents!!! return (*synth_func)( bandPtr,channel,samples); } #endif if(!channel) { /* channel=0 */ bo--; bo &= 0xf; buf = buffs[0]; } else { samples++; buf = buffs[1]; } if(bo & 0x1) { b0 = buf[0]; bo1 = bo; dct64_base(buf[1]+((bo+1)&0xf),buf[0]+bo,bandPtr); } else { b0 = buf[1]; bo1 = bo+1; dct64_base(buf[0]+bo,buf[1]+bo+1,bandPtr); } { register int j; real *window = mp3lib_decwin + 16 - bo1; for (j=16; j; j--,b0+=0x10,window+=0x20,samples+=step) { real sum; sum = window[0x0] * b0[0x0]; sum -= window[0x1] * b0[0x1]; sum += window[0x2] * b0[0x2]; sum -= window[0x3] * b0[0x3]; sum += window[0x4] * b0[0x4]; sum -= window[0x5] * b0[0x5]; sum += window[0x6] * b0[0x6]; sum -= window[0x7] * b0[0x7]; sum += window[0x8] * b0[0x8]; sum -= window[0x9] * b0[0x9]; sum += window[0xA] * b0[0xA]; sum -= window[0xB] * b0[0xB]; sum += window[0xC] * b0[0xC]; sum -= window[0xD] * b0[0xD]; sum += window[0xE] * b0[0xE]; sum -= window[0xF] * b0[0xF]; WRITE_SAMPLE(samples,sum,clip); } { real sum; sum = window[0x0] * b0[0x0]; sum += window[0x2] * b0[0x2]; sum += window[0x4] * b0[0x4]; sum += window[0x6] * b0[0x6]; sum += window[0x8] * b0[0x8]; sum += window[0xA] * b0[0xA]; sum += window[0xC] * b0[0xC]; sum += window[0xE] * b0[0xE]; WRITE_SAMPLE(samples,sum,clip); b0-=0x10,window-=0x20,samples+=step; } window += bo1<<1; for (j=15; j; j--,b0-=0x10,window-=0x20,samples+=step) { real sum; sum = -window[-0x1] * b0[0x0]; sum -= window[-0x2] * b0[0x1]; sum -= window[-0x3] * b0[0x2]; sum -= window[-0x4] * b0[0x3]; sum -= window[-0x5] * b0[0x4]; sum -= window[-0x6] * b0[0x5]; sum -= window[-0x7] * b0[0x6]; sum -= window[-0x8] * b0[0x7]; sum -= window[-0x9] * b0[0x8]; sum -= window[-0xA] * b0[0x9]; sum -= window[-0xB] * b0[0xA]; sum -= window[-0xC] * b0[0xB]; sum -= window[-0xD] * b0[0xC]; sum -= window[-0xE] * b0[0xD]; sum -= window[-0xF] * b0[0xE]; sum -= window[-0x0] * b0[0xF]; WRITE_SAMPLE(samples,sum,clip); } } return clip; }
int synth_ntom(real *bandPtr,int channel,unsigned char *out,int *pnt) { static real buffs[2][2][0x110]; static const int step = 2; static int bo = 1; short *samples = (short *) (out + *pnt); real *b0,(*buf)[0x110]; int clip = 0; int bo1; int ntom; #ifndef NO_EQUALIZER if(param.enable_equalizer) do_equalizer(bandPtr,channel); #endif if(!channel) { bo--; bo &= 0xf; buf = buffs[0]; ntom = ntom_val[1] = ntom_val[0]; } else { samples++; out += 2; /* to compute the right *pnt value */ buf = buffs[1]; ntom = ntom_val[1]; } if(bo & 0x1) { b0 = buf[0]; bo1 = bo; dct64(buf[1]+((bo+1)&0xf),buf[0]+bo,bandPtr); } else { b0 = buf[1]; bo1 = bo+1; dct64(buf[0]+bo,buf[1]+bo+1,bandPtr); } { register int j; real *window = decwin + 16 - bo1; for (j=16;j;j--,window+=0x10) { real sum; ntom += ntom_step; if(ntom < NTOM_MUL) { window += 16; b0 += 16; continue; } sum = *window++ * *b0++; sum -= *window++ * *b0++; sum += *window++ * *b0++; sum -= *window++ * *b0++; sum += *window++ * *b0++; sum -= *window++ * *b0++; sum += *window++ * *b0++; sum -= *window++ * *b0++; sum += *window++ * *b0++; sum -= *window++ * *b0++; sum += *window++ * *b0++; sum -= *window++ * *b0++; sum += *window++ * *b0++; sum -= *window++ * *b0++; sum += *window++ * *b0++; sum -= *window++ * *b0++; while(ntom >= NTOM_MUL) { WRITE_SAMPLE(samples,sum,clip); samples += step; ntom -= NTOM_MUL; } } ntom += ntom_step; if(ntom >= NTOM_MUL) { real sum; sum = window[0x0] * b0[0x0]; sum += window[0x2] * b0[0x2]; sum += window[0x4] * b0[0x4]; sum += window[0x6] * b0[0x6]; sum += window[0x8] * b0[0x8]; sum += window[0xA] * b0[0xA]; sum += window[0xC] * b0[0xC]; sum += window[0xE] * b0[0xE]; while(ntom >= NTOM_MUL) { WRITE_SAMPLE(samples,sum,clip); samples += step; ntom -= NTOM_MUL; } } b0-=0x10,window-=0x20; window += bo1<<1; for (j=15;j;j--,b0-=0x20,window-=0x10) { real sum; ntom += ntom_step; if(ntom < NTOM_MUL) { window -= 16; b0 += 16; continue; } sum = -*(--window) * *b0++; sum -= *(--window) * *b0++; sum -= *(--window) * *b0++; sum -= *(--window) * *b0++; sum -= *(--window) * *b0++; sum -= *(--window) * *b0++; sum -= *(--window) * *b0++; sum -= *(--window) * *b0++; sum -= *(--window) * *b0++; sum -= *(--window) * *b0++; sum -= *(--window) * *b0++; sum -= *(--window) * *b0++; sum -= *(--window) * *b0++; sum -= *(--window) * *b0++; sum -= *(--window) * *b0++; sum -= *(--window) * *b0++; while(ntom >= NTOM_MUL) { WRITE_SAMPLE(samples,sum,clip); samples += step; ntom -= NTOM_MUL; } } } ntom_val[channel] = ntom; *pnt = ((unsigned char *) samples - out); return clip; }
int mpeg3audio_synth_stereo(mpeg3audio_t *audio, float *bandPtr, int channel, float *out, int *pnt) { const int step = 2; float *samples = out + *pnt; register float sum; float *b0, (*buf)[0x110]; int bo1; if(!channel) { audio->bo--; audio->bo &= 0xf; buf = audio->synth_stereo_buffs[0]; } else { samples++; buf = audio->synth_stereo_buffs[1]; } if(audio->bo & 0x1) { b0 = buf[0]; bo1 = audio->bo; mpeg3audio_dct64(buf[1] + ((audio->bo + 1) & 0xf), buf[0] + audio->bo, bandPtr); } else { b0 = buf[1]; bo1 = audio->bo + 1; mpeg3audio_dct64(buf[0] + audio->bo, buf[1] + audio->bo + 1, bandPtr); } /*printf("%f %f %f\n", buf[0][0], buf[1][0], bandPtr[0]); */ { register int j; float *window = mpeg3_decwin + 16 - bo1; for(j = 16; j; j--, b0 += 0x10, window += 0x20, samples += step) { sum = window[0x0] * b0[0x0]; sum -= window[0x1] * b0[0x1]; sum += window[0x2] * b0[0x2]; sum -= window[0x3] * b0[0x3]; sum += window[0x4] * b0[0x4]; sum -= window[0x5] * b0[0x5]; sum += window[0x6] * b0[0x6]; sum -= window[0x7] * b0[0x7]; sum += window[0x8] * b0[0x8]; sum -= window[0x9] * b0[0x9]; sum += window[0xA] * b0[0xA]; sum -= window[0xB] * b0[0xB]; sum += window[0xC] * b0[0xC]; sum -= window[0xD] * b0[0xD]; sum += window[0xE] * b0[0xE]; sum -= window[0xF] * b0[0xF]; WRITE_SAMPLE(samples, sum); } sum = window[0x0] * b0[0x0]; sum += window[0x2] * b0[0x2]; sum += window[0x4] * b0[0x4]; sum += window[0x6] * b0[0x6]; sum += window[0x8] * b0[0x8]; sum += window[0xA] * b0[0xA]; sum += window[0xC] * b0[0xC]; sum += window[0xE] * b0[0xE]; WRITE_SAMPLE(samples, sum); b0 -= 0x10; window -= 0x20; samples += step; window += bo1 << 1; for(j = 15; j; j--, b0 -= 0x10, window -= 0x20, samples += step) { sum = -window[-0x1] * b0[0x0]; sum -= window[-0x2] * b0[0x1]; sum -= window[-0x3] * b0[0x2]; sum -= window[-0x4] * b0[0x3]; sum -= window[-0x5] * b0[0x4]; sum -= window[-0x6] * b0[0x5]; sum -= window[-0x7] * b0[0x6]; sum -= window[-0x8] * b0[0x7]; sum -= window[-0x9] * b0[0x8]; sum -= window[-0xA] * b0[0x9]; sum -= window[-0xB] * b0[0xA]; sum -= window[-0xC] * b0[0xB]; sum -= window[-0xD] * b0[0xC]; sum -= window[-0xE] * b0[0xD]; sum -= window[-0xF] * b0[0xE]; sum -= window[-0x0] * b0[0xF]; WRITE_SAMPLE(samples, sum); } } *pnt += 64; return 0; }
int synth_1to1(real *bandPtr,int channel,unsigned char *out,int *pnt) { #ifndef PENTIUM_OPT static real buffs[2][2][0x110]; static const int step = 2; static int bo = 1; short *samples = (short *) (out + *pnt); real *b0,(*buf)[0x110]; int clip = 0; int bo1; #endif #ifndef NO_EQUALIZER if(param.enable_equalizer) do_equalizer(bandPtr,channel); #endif #ifndef PENTIUM_OPT if(!channel) { bo--; bo &= 0xf; buf = buffs[0]; } else { samples++; buf = buffs[1]; } if(bo & 0x1) { b0 = buf[0]; bo1 = bo; dct64(buf[1]+((bo+1)&0xf),buf[0]+bo,bandPtr); } else { b0 = buf[1]; bo1 = bo+1; dct64(buf[0]+bo,buf[1]+bo+1,bandPtr); } { register int j; real *window = decwin + 16 - bo1; for (j=16;j;j--,b0+=0x10,window+=0x20,samples+=step) { real sum; sum = window[0x0] * b0[0x0]; sum -= window[0x1] * b0[0x1]; sum += window[0x2] * b0[0x2]; sum -= window[0x3] * b0[0x3]; sum += window[0x4] * b0[0x4]; sum -= window[0x5] * b0[0x5]; sum += window[0x6] * b0[0x6]; sum -= window[0x7] * b0[0x7]; sum += window[0x8] * b0[0x8]; sum -= window[0x9] * b0[0x9]; sum += window[0xA] * b0[0xA]; sum -= window[0xB] * b0[0xB]; sum += window[0xC] * b0[0xC]; sum -= window[0xD] * b0[0xD]; sum += window[0xE] * b0[0xE]; sum -= window[0xF] * b0[0xF]; WRITE_SAMPLE(samples,sum,clip); } { real sum; sum = window[0x0] * b0[0x0]; sum += window[0x2] * b0[0x2]; sum += window[0x4] * b0[0x4]; sum += window[0x6] * b0[0x6]; sum += window[0x8] * b0[0x8]; sum += window[0xA] * b0[0xA]; sum += window[0xC] * b0[0xC]; sum += window[0xE] * b0[0xE]; WRITE_SAMPLE(samples,sum,clip); b0-=0x10,window-=0x20,samples+=step; } window += bo1<<1; for (j=15;j;j--,b0-=0x10,window-=0x20,samples+=step) { real sum; sum = -window[-0x1] * b0[0x0]; sum -= window[-0x2] * b0[0x1]; sum -= window[-0x3] * b0[0x2]; sum -= window[-0x4] * b0[0x3]; sum -= window[-0x5] * b0[0x4]; sum -= window[-0x6] * b0[0x5]; sum -= window[-0x7] * b0[0x6]; sum -= window[-0x8] * b0[0x7]; sum -= window[-0x9] * b0[0x8]; sum -= window[-0xA] * b0[0x9]; sum -= window[-0xB] * b0[0xA]; sum -= window[-0xC] * b0[0xB]; sum -= window[-0xD] * b0[0xC]; sum -= window[-0xE] * b0[0xD]; sum -= window[-0xF] * b0[0xE]; sum -= window[-0x0] * b0[0xF]; WRITE_SAMPLE(samples,sum,clip); } } *pnt += 128; return clip; #elif defined(USE_MMX) { static short buffs[2][2][0x110]; static int bo = 1; short *samples = (short *) (out + *pnt); synth_1to1_MMX(bandPtr, channel, samples, (short *) buffs, &bo); *pnt += 128; return 0; } #else { int ret; ret = synth_1to1_pent(bandPtr,channel,out+*pnt); *pnt += 128; return ret; } #endif }
int synth_2to1(real *bandPtr,int channel,unsigned char *out) { static real buffs[2][2][0x110]; static const int step = 2; static int bo = 1; short *samples = (short *) out; real *b0,(*buf)[0x110]; int clip = 0; int bo1; if(flags.equalizer) do_equalizer(bandPtr,channel); if(!channel) { bo--; bo &= 0xf; buf = buffs[0]; } else { samples++; buf = buffs[1]; } if(bo & 0x1) { b0 = buf[0]; bo1 = bo; dct64(buf[1]+((bo+1)&0xf),buf[0]+bo,bandPtr); } else { b0 = buf[1]; bo1 = bo+1; dct64(buf[0]+bo,buf[1]+bo+1,bandPtr); } { register int j; real *window = decwin + 16 - bo1; for (j=8;j;j--,b0+=0x10,window+=0x30) { real sum; sum = *window++ * *b0++; sum -= *window++ * *b0++; sum += *window++ * *b0++; sum -= *window++ * *b0++; sum += *window++ * *b0++; sum -= *window++ * *b0++; sum += *window++ * *b0++; sum -= *window++ * *b0++; sum += *window++ * *b0++; sum -= *window++ * *b0++; sum += *window++ * *b0++; sum -= *window++ * *b0++; sum += *window++ * *b0++; sum -= *window++ * *b0++; sum += *window++ * *b0++; sum -= *window++ * *b0++; WRITE_SAMPLE(samples,sum,clip); samples += step; #if 0 WRITE_SAMPLE(samples,sum,clip); samples += step; #endif } { real sum; sum = window[0x0] * b0[0x0]; sum += window[0x2] * b0[0x2]; sum += window[0x4] * b0[0x4]; sum += window[0x6] * b0[0x6]; sum += window[0x8] * b0[0x8]; sum += window[0xA] * b0[0xA]; sum += window[0xC] * b0[0xC]; sum += window[0xE] * b0[0xE]; WRITE_SAMPLE(samples,sum,clip); samples += step; #if 0 WRITE_SAMPLE(samples,sum,clip); samples += step; #endif b0-=0x20,window-=0x40; } window += bo1<<1; for (j=7;j;j--,b0-=0x30,window-=0x30) { real sum; sum = -*(--window) * *b0++; sum -= *(--window) * *b0++; sum -= *(--window) * *b0++; sum -= *(--window) * *b0++; sum -= *(--window) * *b0++; sum -= *(--window) * *b0++; sum -= *(--window) * *b0++; sum -= *(--window) * *b0++; sum -= *(--window) * *b0++; sum -= *(--window) * *b0++; sum -= *(--window) * *b0++; sum -= *(--window) * *b0++; sum -= *(--window) * *b0++; sum -= *(--window) * *b0++; sum -= *(--window) * *b0++; sum -= *(--window) * *b0++; WRITE_SAMPLE(samples,sum,clip); samples += step; #if 0 WRITE_SAMPLE(samples,sum,clip); samples += step; #endif } } return clip; }