Alc/mixer_inc.c

   1 #include "config.h"
   2
   3 #include "alMain.h"
   4 #include "alSource.h"
   5 #include "mixer_defs.h"
   6
   7 #ifdef __GNUC__
   8 #define LIKELY(x) __builtin_expect(!!(x), 1)
   9 #define UNLIKELY(x) __builtin_expect(!!(x), 0)
  10 #else
  11 #define LIKELY(x) (x)
  12 #define UNLIKELY(x) (x)
  13 #endif
  14
  15 #define REAL_MERGE2(a,b) a##b
  16 #define MERGE2(a,b) REAL_MERGE2(a,b)
  17
  18 #define MixDirect_Hrtf MERGE2(MixDirect_Hrtf_,SUFFIX)
  19
  20
  21 static inline void ApplyCoeffsStep(const ALuint irSize,
  22                                    ALfloat (*restrict Coeffs)[2],
  23                                    const ALfloat (*restrict CoeffStep)[2]);
  24 static inline void ApplyCoeffs(ALuint Offset, ALfloat (*restrict Values)[2],
  25                                const ALuint irSize,
  26                                ALfloat (*restrict Coeffs)[2],
  27                                ALfloat left, ALfloat right);
  28
  29
  30 void MixDirect_Hrtf(const DirectParams *params, const ALfloat *restrict data, ALuint srcchan,
  31   ALuint OutPos, ALuint SamplesToDo, ALuint BufferSize)
  32 {
  33     ALfloat (*restrict DryBuffer)[BUFFERSIZE] = params->OutBuffer;
  34     ALfloat *restrict ClickRemoval = params->ClickRemoval;
  35     ALfloat *restrict PendingClicks = params->PendingClicks;
  36     const ALuint IrSize = params->Hrtf.Params.IrSize;
  37     const ALint *restrict DelayStep = params->Hrtf.Params.DelayStep;
  38     const ALfloat (*restrict CoeffStep)[2] = params->Hrtf.Params.CoeffStep;
  39     const ALfloat (*restrict TargetCoeffs)[2] = params->Hrtf.Params.Coeffs[srcchan];
  40     const ALuint *restrict TargetDelay = params->Hrtf.Params.Delay[srcchan];
  41     ALfloat *restrict History = params->Hrtf.State->History[srcchan];
  42     ALfloat (*restrict Values)[2] = params->Hrtf.State->Values[srcchan];
  43     ALint Counter = maxu(params->Hrtf.State->Counter, OutPos) - OutPos;
  44     ALuint Offset = params->Hrtf.State->Offset + OutPos;
  45     ALIGN(16) ALfloat Coeffs[HRIR_LENGTH][2];
  46     ALuint Delay[2];
  47     ALfloat left, right;
  48     ALuint pos;
  49     ALuint c;
  50
  51     pos = 0;
  52     for(c = 0;c < IrSize;c++)
  53     {
  54         Coeffs[c][0] = TargetCoeffs[c][0] - (CoeffStep[c][0]*Counter);
  55         Coeffs[c][1] = TargetCoeffs[c][1] - (CoeffStep[c][1]*Counter);
  56     }
  57
  58     Delay[0] = TargetDelay[0] - (DelayStep[0]*Counter);
  59     Delay[1] = TargetDelay[1] - (DelayStep[1]*Counter);
  60
  61     if(LIKELY(OutPos == 0))
  62     {
  63         History[Offset&SRC_HISTORY_MASK] = data[pos];
  64         left  = lerp(History[(Offset-(Delay[0]>>HRTFDELAY_BITS))&SRC_HISTORY_MASK],
  65                      History[(Offset-(Delay[0]>>HRTFDELAY_BITS)-1)&SRC_HISTORY_MASK],
  66                      (Delay[0]&HRTFDELAY_MASK)*(1.0f/HRTFDELAY_FRACONE));
  67         right = lerp(History[(Offset-(Delay[1]>>HRTFDELAY_BITS))&SRC_HISTORY_MASK],
  68                      History[(Offset-(Delay[1]>>HRTFDELAY_BITS)-1)&SRC_HISTORY_MASK],
  69                      (Delay[1]&HRTFDELAY_MASK)*(1.0f/HRTFDELAY_FRACONE));
  70
  71         ClickRemoval[FrontLeft]  -= Values[(Offset+1)&HRIR_MASK][0] +
  72                                     Coeffs[0][0] * left;
  73         ClickRemoval[FrontRight] -= Values[(Offset+1)&HRIR_MASK][1] +
  74                                     Coeffs[0][1] * right;
  75     }
  76     for(pos = 0;pos < BufferSize && Counter > 0;pos++)
  77     {
  78         History[Offset&SRC_HISTORY_MASK] = data[pos];
  79         left  = lerp(History[(Offset-(Delay[0]>>HRTFDELAY_BITS))&SRC_HISTORY_MASK],
  80                      History[(Offset-(Delay[0]>>HRTFDELAY_BITS)-1)&SRC_HISTORY_MASK],
  81                      (Delay[0]&HRTFDELAY_MASK)*(1.0f/HRTFDELAY_FRACONE));
  82         right = lerp(History[(Offset-(Delay[1]>>HRTFDELAY_BITS))&SRC_HISTORY_MASK],
  83                      History[(Offset-(Delay[1]>>HRTFDELAY_BITS)-1)&SRC_HISTORY_MASK],
  84                      (Delay[1]&HRTFDELAY_MASK)*(1.0f/HRTFDELAY_FRACONE));
  85
  86         Delay[0] += DelayStep[0];
  87         Delay[1] += DelayStep[1];
  88
  89         Values[(Offset+IrSize)&HRIR_MASK][0] = 0.0f;
  90         Values[(Offset+IrSize)&HRIR_MASK][1] = 0.0f;
  91         Offset++;
  92
  93         ApplyCoeffs(Offset, Values, IrSize, Coeffs, left, right);
  94         DryBuffer[FrontLeft][OutPos]  += Values[Offset&HRIR_MASK][0];
  95         DryBuffer[FrontRight][OutPos] += Values[Offset&HRIR_MASK][1];
  96         ApplyCoeffsStep(IrSize, Coeffs, CoeffStep);
  97
  98         OutPos++;
  99         Counter--;
 100     }
 101
 102     Delay[0] >>= HRTFDELAY_BITS;
 103     Delay[1] >>= HRTFDELAY_BITS;
 104     for(;pos < BufferSize;pos++)
 105     {
 106         History[Offset&SRC_HISTORY_MASK] = data[pos];
 107         left = History[(Offset-Delay[0])&SRC_HISTORY_MASK];
 108         right = History[(Offset-Delay[1])&SRC_HISTORY_MASK];
 109
 110         Values[(Offset+IrSize)&HRIR_MASK][0] = 0.0f;
 111         Values[(Offset+IrSize)&HRIR_MASK][1] = 0.0f;
 112         Offset++;
 113
 114         ApplyCoeffs(Offset, Values, IrSize, Coeffs, left, right);
 115         DryBuffer[FrontLeft][OutPos]  += Values[Offset&HRIR_MASK][0];
 116         DryBuffer[FrontRight][OutPos] += Values[Offset&HRIR_MASK][1];
 117
 118         OutPos++;
 119     }
 120     if(LIKELY(OutPos == SamplesToDo))
 121     {
 122         History[Offset&SRC_HISTORY_MASK] = data[pos];
 123         left = History[(Offset-Delay[0])&SRC_HISTORY_MASK];
 124         right = History[(Offset-Delay[1])&SRC_HISTORY_MASK];
 125
 126         PendingClicks[FrontLeft]  += Values[(Offset+1)&HRIR_MASK][0] +
 127                                      Coeffs[0][0] * left;
 128         PendingClicks[FrontRight] += Values[(Offset+1)&HRIR_MASK][1] +
 129                                      Coeffs[0][1] * right;
 130     }
 131 }
 132
 133
 134 #undef MixDirect_Hrtf
 135
 136 #undef MERGE2
 137 #undef REAL_MERGE2
 138
 139 #undef UNLIKELY
 140 #undef LIKELY