?  search:
summary | log | graphiclog1 | graphiclog2 | commit | commitdiff | tree | refs | edit | fork
blame | history | raw | HEAD
Increase the band-split IR for decoding ambisonics to HRTF
[openal-soft.git] / Alc / hrtf.c
blob7e5f6e9673203454c5fe71743a711ad7423d6e85
1 /**
2 * OpenAL cross platform audio library
3 * Copyright (C) 2011 by Chris Robinson
4 * This library is free software; you can redistribute it and/or
5 * modify it under the terms of the GNU Library General Public
6 * License as published by the Free Software Foundation; either
7 * version 2 of the License, or (at your option) any later version.
8 *
9 * This library is distributed in the hope that it will be useful,
10 * but WITHOUT ANY WARRANTY; without even the implied warranty of
11 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
12 * Library General Public License for more details.
13 *
14 * You should have received a copy of the GNU Library General Public
15 * License along with this library; if not, write to the
16 * Free Software Foundation, Inc.,
17 * 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA.
18 * Or go to http://www.gnu.org/copyleft/lgpl.html
19 */
20
21 #include "config.h"
22
23 #include <stdlib.h>
24 #include <ctype.h>
25
26 #include "AL/al.h"
27 #include "AL/alc.h"
28 #include "alMain.h"
29 #include "alSource.h"
30 #include "alu.h"
31 #include "hrtf.h"
32 #include "alconfig.h"
33 #include "filters/splitter.h"
34
35 #include "compat.h"
36 #include "almalloc.h"
37
38
39 /* Current data set limits defined by the makehrtf utility. */
40 #define MIN_IR_SIZE (8)
41 #define MAX_IR_SIZE (512)
42 #define MOD_IR_SIZE (8)
43
44 #define MIN_FD_COUNT (1)
45 #define MAX_FD_COUNT (16)
46
47 #define MIN_FD_DISTANCE (50)
48 #define MAX_FD_DISTANCE (2500)
49
50 #define MIN_EV_COUNT (5)
51 #define MAX_EV_COUNT (128)
52
53 #define MIN_AZ_COUNT (1)
54 #define MAX_AZ_COUNT (128)
55
56 #define MAX_HRIR_DELAY (HRTF_HISTORY_LENGTH-1)
57
58 struct HrtfEntry {
59 struct HrtfEntry *next;
60 struct Hrtf *handle;
61 char filename[];
62 };
63
64 static const ALchar magicMarker00[8] = "MinPHR00";
65 static const ALchar magicMarker01[8] = "MinPHR01";
66 static const ALchar magicMarker02[8] = "MinPHR02";
67
68 /* First value for pass-through coefficients (remaining are 0), used for omni-
69 * directional sounds. */
70 static const ALfloat PassthruCoeff = 0.707106781187f/*sqrt(0.5)*/;
71
72 static ATOMIC_FLAG LoadedHrtfLock = ATOMIC_FLAG_INIT;
73 static struct HrtfEntry *LoadedHrtfs = NULL;
74
75
76 /* Calculate the elevation index given the polar elevation in radians. This
77 * will return an index between 0 and (evcount - 1). Assumes the FPU is in
78 * round-to-zero mode.
79 */
80 static ALsizei CalcEvIndex(ALsizei evcount, ALfloat ev, ALfloat *mu)
81 {
82 ALsizei idx;
83 ev = (F_PI_2+ev) * (evcount-1) / F_PI;
84 idx = mini(fastf2i(ev), evcount-1);
85
86 *mu = ev - idx;
87 return idx;
88 }
89
90 /* Calculate the azimuth index given the polar azimuth in radians. This will
91 * return an index between 0 and (azcount - 1). Assumes the FPU is in round-to-
92 * zero mode.
93 */
94 static ALsizei CalcAzIndex(ALsizei azcount, ALfloat az, ALfloat *mu)
95 {
96 ALsizei idx;
97 az = (F_TAU+az) * azcount / F_TAU;
98
99 idx = fastf2i(az) % azcount;
100 *mu = az - floorf(az);
101 return idx;
102 }
103
104 /* Calculates static HRIR coefficients and delays for the given polar elevation
105 * and azimuth in radians. The coefficients are normalized.
106 */
107 void GetHrtfCoeffs(const struct Hrtf *Hrtf, ALfloat elevation, ALfloat azimuth, ALfloat spread,
108 ALfloat (*restrict coeffs)[2], ALsizei *delays)
109 {
110 ALsizei evidx, azidx, idx[4];
111 ALsizei evoffset;
112 ALfloat emu, amu[2];
113 ALfloat blend[4];
114 ALfloat dirfact;
115 ALsizei i, c;
116
117 dirfact = 1.0f - (spread / F_TAU);
118
119 /* Claculate the lower elevation index. */
120 evidx = CalcEvIndex(Hrtf->evCount, elevation, &emu);
121 evoffset = Hrtf->evOffset[evidx];
122
123 /* Calculate lower azimuth index. */
124 azidx= CalcAzIndex(Hrtf->azCount[evidx], azimuth, &amu[0]);
125
126 /* Calculate the lower HRIR indices. */
127 idx[0] = evoffset + azidx;
128 idx[1] = evoffset + ((azidx+1) % Hrtf->azCount[evidx]);
129 if(evidx < Hrtf->evCount-1)
130 {
131 /* Increment elevation to the next (upper) index. */
132 evidx++;
133 evoffset = Hrtf->evOffset[evidx];
134
135 /* Calculate upper azimuth index. */
136 azidx = CalcAzIndex(Hrtf->azCount[evidx], azimuth, &amu[1]);
137
138 /* Calculate the upper HRIR indices. */
139 idx[2] = evoffset + azidx;
140 idx[3] = evoffset + ((azidx+1) % Hrtf->azCount[evidx]);
141 }
142 else
143 {
144 /* If the lower elevation is the top index, the upper elevation is the
145 * same as the lower.
146 */
147 amu[1] = amu[0];
148 idx[2] = idx[0];
149 idx[3] = idx[1];
150 }
151
152 /* Calculate bilinear blending weights, attenuated according to the
153 * directional panning factor.
154 */
155 blend[0] = (1.0f-emu) * (1.0f-amu[0]) * dirfact;
156 blend[1] = (1.0f-emu) * ( amu[0]) * dirfact;
157 blend[2] = ( emu) * (1.0f-amu[1]) * dirfact;
158 blend[3] = ( emu) * ( amu[1]) * dirfact;
159
160 /* Calculate the blended HRIR delays. */
161 delays[0] = fastf2i(
162 Hrtf->delays[idx[0]][0]*blend[0] + Hrtf->delays[idx[1]][0]*blend[1] +
163 Hrtf->delays[idx[2]][0]*blend[2] + Hrtf->delays[idx[3]][0]*blend[3] + 0.5f
164 );
165 delays[1] = fastf2i(
166 Hrtf->delays[idx[0]][1]*blend[0] + Hrtf->delays[idx[1]][1]*blend[1] +
167 Hrtf->delays[idx[2]][1]*blend[2] + Hrtf->delays[idx[3]][1]*blend[3] + 0.5f
168 );
169
170 /* Calculate the sample offsets for the HRIR indices. */
171 idx[0] *= Hrtf->irSize;
172 idx[1] *= Hrtf->irSize;
173 idx[2] *= Hrtf->irSize;
174 idx[3] *= Hrtf->irSize;
175
176 ASSUME(Hrtf->irSize >= MIN_IR_SIZE && (Hrtf->irSize%MOD_IR_SIZE) == 0);
177 coeffs = ASSUME_ALIGNED(coeffs, 16);
178 /* Calculate the blended HRIR coefficients. */
179 coeffs[0][0] = PassthruCoeff * (1.0f-dirfact);
180 coeffs[0][1] = PassthruCoeff * (1.0f-dirfact);
181 for(i = 1;i < Hrtf->irSize;i++)
182 {
183 coeffs[i][0] = 0.0f;
184 coeffs[i][1] = 0.0f;
185 }
186 for(c = 0;c < 4;c++)
187 {
188 const ALfloat (*restrict srccoeffs)[2] = ASSUME_ALIGNED(Hrtf->coeffs+idx[c], 16);
189 for(i = 0;i < Hrtf->irSize;i++)
190 {
191 coeffs[i][0] += srccoeffs[i][0] * blend[c];
192 coeffs[i][1] += srccoeffs[i][1] * blend[c];
193 }
194 }
195 }
196
197
198 void BuildBFormatHrtf(const struct Hrtf *Hrtf, DirectHrtfState *state, ALsizei NumChannels, const struct AngularPoint *AmbiPoints, const ALfloat (*restrict AmbiMatrix)[MAX_AMBI_COEFFS], ALsizei AmbiCount, const ALfloat *restrict AmbiOrderHFGain)
199 {
200 /* Set this to 2 for dual-band HRTF processing. May require a higher quality
201 * band-splitter, or better calculation of the new IR length to deal with the
202 * tail generated by the filter.
203 */
204 #define NUM_BANDS 2
205 BandSplitter splitter;
206 ALdouble (*tmpres)[HRIR_LENGTH][2];
207 ALsizei idx[HRTF_AMBI_MAX_CHANNELS];
208 ALsizei min_delay = HRTF_HISTORY_LENGTH;
209 ALfloat temps[3][HRIR_LENGTH];
210 ALsizei max_length = 0;
211 ALsizei i, c, b;
212
213 for(c = 0;c < AmbiCount;c++)
214 {
215 ALuint evidx, azidx;
216 ALuint evoffset;
217 ALuint azcount;
218
219 /* Calculate elevation index. */
220 evidx = (ALsizei)floorf((F_PI_2 + AmbiPoints[c].Elev) *
221 (Hrtf->evCount-1)/F_PI + 0.5f);
222 evidx = clampi(evidx, 0, Hrtf->evCount-1);
223
224 azcount = Hrtf->azCount[evidx];
225 evoffset = Hrtf->evOffset[evidx];
226
227 /* Calculate azimuth index for this elevation. */
228 azidx = (ALsizei)floorf((F_TAU+AmbiPoints[c].Azim) *
229 azcount/F_TAU + 0.5f) % azcount;
230
231 /* Calculate indices for left and right channels. */
232 idx[c] = evoffset + azidx;
233
234 min_delay = mini(min_delay, mini(Hrtf->delays[idx[c]][0], Hrtf->delays[idx[c]][1]));
235 }
236
237 tmpres = al_calloc(16, NumChannels * sizeof(*tmpres));
238
239 memset(temps, 0, sizeof(temps));
240 bandsplit_init(&splitter, 400.0f / (ALfloat)Hrtf->sampleRate);
241 for(c = 0;c < AmbiCount;c++)
242 {
243 const ALfloat (*fir)[2] = &Hrtf->coeffs[idx[c] * Hrtf->irSize];
244 ALsizei ldelay = Hrtf->delays[idx[c]][0] - min_delay;
245 ALsizei rdelay = Hrtf->delays[idx[c]][1] - min_delay;
246
247 if(NUM_BANDS == 1)
248 {
249 max_length = maxi(max_length,
250 mini(maxi(ldelay, rdelay) + Hrtf->irSize, HRIR_LENGTH)
251 );
252
253 for(i = 0;i < NumChannels;++i)
254 {
255 ALdouble mult = (ALdouble)AmbiOrderHFGain[(ALsizei)sqrt(i)] * AmbiMatrix[c][i];
256 ALsizei lidx = ldelay, ridx = rdelay;
257 ALsizei j = 0;
258 while(lidx < HRIR_LENGTH && ridx < HRIR_LENGTH && j < Hrtf->irSize)
259 {
260 tmpres[i][lidx++][0] += fir[j][0] * mult;
261 tmpres[i][ridx++][1] += fir[j][1] * mult;
262 j++;
263 }
264 }
265 }
266 else
267 {
268 /* Increase the IR size by 2/3rds to account for the tail generated
269 * by the band-split filter.
270 */
271 const ALsizei irsize = mini(Hrtf->irSize*5/3, HRIR_LENGTH);
272
273 max_length = maxi(max_length,
274 mini(maxi(ldelay, rdelay) + irsize, HRIR_LENGTH)
275 );
276
277 /* Band-split left HRIR into low and high frequency responses. */
278 bandsplit_clear(&splitter);
279 for(i = 0;i < Hrtf->irSize;i++)
280 temps[2][i] = fir[i][0];
281 bandsplit_process(&splitter, temps[0], temps[1], temps[2], HRIR_LENGTH);
282
283 /* Apply left ear response with delay. */
284 for(i = 0;i < NumChannels;++i)
285 {
286 ALfloat hfgain = AmbiOrderHFGain[(ALsizei)sqrt(i)];
287 for(b = 0;b < NUM_BANDS;b++)
288 {
289 ALdouble mult = AmbiMatrix[c][i] * (ALdouble)((b==0) ? hfgain : 1.0);
290 ALsizei lidx = ldelay;
291 ALsizei j = 0;
292 while(lidx < HRIR_LENGTH)
293 tmpres[i][lidx++][0] += temps[b][j++] * mult;
294 }
295 }
296
297 /* Band-split right HRIR into low and high frequency responses. */
298 bandsplit_clear(&splitter);
299 for(i = 0;i < Hrtf->irSize;i++)
300 temps[2][i] = fir[i][1];
301 bandsplit_process(&splitter, temps[0], temps[1], temps[2], HRIR_LENGTH);
302
303 /* Apply right ear response with delay. */
304 for(i = 0;i < NumChannels;++i)
305 {
306 ALfloat hfgain = AmbiOrderHFGain[(ALsizei)sqrt(i)];
307 for(b = 0;b < NUM_BANDS;b++)
308 {
309 ALdouble mult = AmbiMatrix[c][i] * (ALdouble)((b==0) ? hfgain : 1.0);
310 ALsizei ridx = rdelay;
311 ALsizei j = 0;
312 while(ridx < HRIR_LENGTH)
313 tmpres[i][ridx++][1] += temps[b][j++] * mult;
314 }
315 }
316 }
317 }
318 /* Round up to the next IR size multiple. */
319 max_length += MOD_IR_SIZE-1;
320 max_length -= max_length%MOD_IR_SIZE;
321
322 for(i = 0;i < NumChannels;++i)
323 {
324 int idx;
325 for(idx = 0;idx < HRIR_LENGTH;idx++)
326 {
327 state->Chan[i].Coeffs[idx][0] = (ALfloat)tmpres[i][idx][0];
328 state->Chan[i].Coeffs[idx][1] = (ALfloat)tmpres[i][idx][1];
329 }
330 }
331
332 al_free(tmpres);
333 tmpres = NULL;
334
335 TRACE("Skipped delay: %d, new FIR length: %d\n", min_delay, max_length);
336 state->IrSize = max_length;
337 #undef NUM_BANDS
338 }
339
340
341 static struct Hrtf *CreateHrtfStore(ALuint rate, ALsizei irSize,
342 ALfloat distance, ALsizei evCount, ALsizei irCount, const ALubyte *azCount,
343 const ALushort *evOffset, const ALfloat (*coeffs)[2], const ALubyte (*delays)[2],
344 const char *filename)
345 {
346 struct Hrtf *Hrtf;
347 size_t total;
348
349 total = sizeof(struct Hrtf);
350 total += sizeof(Hrtf->azCount[0])*evCount;
351 total = RoundUp(total, sizeof(ALushort)); /* Align for ushort fields */
352 total += sizeof(Hrtf->evOffset[0])*evCount;
353 total = RoundUp(total, 16); /* Align for coefficients using SIMD */
354 total += sizeof(Hrtf->coeffs[0])*irSize*irCount;
355 total += sizeof(Hrtf->delays[0])*irCount;
356
357 Hrtf = al_calloc(16, total);
358 if(Hrtf == NULL)
359 ERR("Out of memory allocating storage for %s.\n", filename);
360 else
361 {
362 uintptr_t offset = sizeof(struct Hrtf);
363 char *base = (char*)Hrtf;
364 ALushort *_evOffset;
365 ALubyte *_azCount;
366 ALubyte (*_delays)[2];
367 ALfloat (*_coeffs)[2];
368 ALsizei i;
369
370 InitRef(&Hrtf->ref, 0);
371 Hrtf->sampleRate = rate;
372 Hrtf->irSize = irSize;
373 Hrtf->distance = distance;
374 Hrtf->evCount = evCount;
375
376 /* Set up pointers to storage following the main HRTF struct. */
377 _azCount = (ALubyte*)(base + offset);
378 offset += sizeof(_azCount[0])*evCount;
379
380 offset = RoundUp(offset, sizeof(ALushort)); /* Align for ushort fields */
381 _evOffset = (ALushort*)(base + offset);
382 offset += sizeof(_evOffset[0])*evCount;
383
384 offset = RoundUp(offset, 16); /* Align for coefficients using SIMD */
385 _coeffs = (ALfloat(*)[2])(base + offset);
386 offset += sizeof(_coeffs[0])*irSize*irCount;
387
388 _delays = (ALubyte(*)[2])(base + offset);
389 offset += sizeof(_delays[0])*irCount;
390
391 assert(offset == total);
392
393 /* Copy input data to storage. */
394 for(i = 0;i < evCount;i++) _azCount[i] = azCount[i];
395 for(i = 0;i < evCount;i++) _evOffset[i] = evOffset[i];
396 for(i = 0;i < irSize*irCount;i++)
397 {
398 _coeffs[i][0] = coeffs[i][0];
399 _coeffs[i][1] = coeffs[i][1];
400 }
401 for(i = 0;i < irCount;i++)
402 {
403 _delays[i][0] = delays[i][0];
404 _delays[i][1] = delays[i][1];
405 }
406
407 /* Finally, assign the storage pointers. */
408 Hrtf->azCount = _azCount;
409 Hrtf->evOffset = _evOffset;
410 Hrtf->coeffs = _coeffs;
411 Hrtf->delays = _delays;
412 }
413
414 return Hrtf;
415 }
416
417 static ALubyte GetLE_ALubyte(const ALubyte **data, size_t *len)
418 {
419 ALubyte ret = (*data)[0];
420 *data += 1; *len -= 1;
421 return ret;
422 }
423
424 static ALshort GetLE_ALshort(const ALubyte **data, size_t *len)
425 {
426 ALshort ret = (*data)[0] | ((*data)[1]<<8);
427 *data += 2; *len -= 2;
428 return ret;
429 }
430
431 static ALushort GetLE_ALushort(const ALubyte **data, size_t *len)
432 {
433 ALushort ret = (*data)[0] | ((*data)[1]<<8);
434 *data += 2; *len -= 2;
435 return ret;
436 }
437
438 static ALint GetLE_ALint24(const ALubyte **data, size_t *len)
439 {
440 ALint ret = (*data)[0] | ((*data)[1]<<8) | ((*data)[2]<<16);
441 *data += 3; *len -= 3;
442 return (ret^0x800000) - 0x800000;
443 }
444
445 static ALuint GetLE_ALuint(const ALubyte **data, size_t *len)
446 {
447 ALuint ret = (*data)[0] | ((*data)[1]<<8) | ((*data)[2]<<16) | ((*data)[3]<<24);
448 *data += 4; *len -= 4;
449 return ret;
450 }
451
452 static const ALubyte *Get_ALubytePtr(const ALubyte **data, size_t *len, size_t size)
453 {
454 const ALubyte *ret = *data;
455 *data += size; *len -= size;
456 return ret;
457 }
458
459 static struct Hrtf *LoadHrtf00(const ALubyte *data, size_t datalen, const char *filename)
460 {
461 struct Hrtf *Hrtf = NULL;
462 ALboolean failed = AL_FALSE;
463 ALuint rate = 0;
464 ALushort irCount = 0;
465 ALushort irSize = 0;
466 ALubyte evCount = 0;
467 ALubyte *azCount = NULL;
468 ALushort *evOffset = NULL;
469 ALfloat (*coeffs)[2] = NULL;
470 ALubyte (*delays)[2] = NULL;
471 ALsizei i, j;
472
473 if(datalen < 9)
474 {
475 ERR("Unexpected end of %s data (req %d, rem "SZFMT")\n", filename, 9, datalen);
476 return NULL;
477 }
478
479 rate = GetLE_ALuint(&data, &datalen);
480
481 irCount = GetLE_ALushort(&data, &datalen);
482
483 irSize = GetLE_ALushort(&data, &datalen);
484
485 evCount = GetLE_ALubyte(&data, &datalen);
486
487 if(irSize < MIN_IR_SIZE || irSize > MAX_IR_SIZE || (irSize%MOD_IR_SIZE))
488 {
489 ERR("Unsupported HRIR size: irSize=%d (%d to %d by %d)\n",
490 irSize, MIN_IR_SIZE, MAX_IR_SIZE, MOD_IR_SIZE);
491 failed = AL_TRUE;
492 }
493 if(evCount < MIN_EV_COUNT || evCount > MAX_EV_COUNT)
494 {
495 ERR("Unsupported elevation count: evCount=%d (%d to %d)\n",
496 evCount, MIN_EV_COUNT, MAX_EV_COUNT);
497 failed = AL_TRUE;
498 }
499 if(failed)
500 return NULL;
501
502 if(datalen < evCount*2u)
503 {
504 ERR("Unexpected end of %s data (req %d, rem "SZFMT")\n", filename, evCount*2, datalen);
505 return NULL;
506 }
507
508 azCount = malloc(sizeof(azCount[0])*evCount);
509 evOffset = malloc(sizeof(evOffset[0])*evCount);
510 if(azCount == NULL || evOffset == NULL)
511 {
512 ERR("Out of memory.\n");
513 failed = AL_TRUE;
514 }
515
516 if(!failed)
517 {
518 evOffset[0] = GetLE_ALushort(&data, &datalen);
519 for(i = 1;i < evCount;i++)
520 {
521 evOffset[i] = GetLE_ALushort(&data, &datalen);
522 if(evOffset[i] <= evOffset[i-1])
523 {
524 ERR("Invalid evOffset: evOffset[%d]=%d (last=%d)\n",
525 i, evOffset[i], evOffset[i-1]);
526 failed = AL_TRUE;
527 }
528
529 azCount[i-1] = evOffset[i] - evOffset[i-1];
530 if(azCount[i-1] < MIN_AZ_COUNT || azCount[i-1] > MAX_AZ_COUNT)
531 {
532 ERR("Unsupported azimuth count: azCount[%d]=%d (%d to %d)\n",
533 i-1, azCount[i-1], MIN_AZ_COUNT, MAX_AZ_COUNT);
534 failed = AL_TRUE;
535 }
536 }
537 if(irCount <= evOffset[i-1])
538 {
539 ERR("Invalid evOffset: evOffset[%d]=%d (irCount=%d)\n",
540 i-1, evOffset[i-1], irCount);
541 failed = AL_TRUE;
542 }
543
544 azCount[i-1] = irCount - evOffset[i-1];
545 if(azCount[i-1] < MIN_AZ_COUNT || azCount[i-1] > MAX_AZ_COUNT)
546 {
547 ERR("Unsupported azimuth count: azCount[%d]=%d (%d to %d)\n",
548 i-1, azCount[i-1], MIN_AZ_COUNT, MAX_AZ_COUNT);
549 failed = AL_TRUE;
550 }
551 }
552
553 if(!failed)
554 {
555 coeffs = malloc(sizeof(coeffs[0])*irSize*irCount);
556 delays = malloc(sizeof(delays[0])*irCount);
557 if(coeffs == NULL || delays == NULL)
558 {
559 ERR("Out of memory.\n");
560 failed = AL_TRUE;
561 }
562 }
563
564 if(!failed)
565 {
566 size_t reqsize = 2*irSize*irCount + irCount;
567 if(datalen < reqsize)
568 {
569 ERR("Unexpected end of %s data (req "SZFMT", rem "SZFMT")\n",
570 filename, reqsize, datalen);
571 failed = AL_TRUE;
572 }
573 }
574
575 if(!failed)
576 {
577 for(i = 0;i < irCount;i++)
578 {
579 for(j = 0;j < irSize;j++)
580 coeffs[i*irSize + j][0] = GetLE_ALshort(&data, &datalen) / 32768.0f;
581 }
582
583 for(i = 0;i < irCount;i++)
584 {
585 delays[i][0] = GetLE_ALubyte(&data, &datalen);
586 if(delays[i][0] > MAX_HRIR_DELAY)
587 {
588 ERR("Invalid delays[%d]: %d (%d)\n", i, delays[i][0], MAX_HRIR_DELAY);
589 failed = AL_TRUE;
590 }
591 }
592 }
593
594 if(!failed)
595 {
596 /* Mirror the left ear responses to the right ear. */
597 for(i = 0;i < evCount;i++)
598 {
599 ALushort evoffset = evOffset[i];
600 ALubyte azcount = azCount[i];
601 for(j = 0;j < azcount;j++)
602 {
603 ALsizei lidx = evoffset + j;
604 ALsizei ridx = evoffset + ((azcount-j) % azcount);
605 ALsizei k;
606
607 for(k = 0;k < irSize;k++)
608 coeffs[ridx*irSize + k][1] = coeffs[lidx*irSize + k][0];
609 delays[ridx][1] = delays[lidx][0];
610 }
611 }
612
613 Hrtf = CreateHrtfStore(rate, irSize, 0.0f, evCount, irCount, azCount,
614 evOffset, coeffs, delays, filename);
615 }
616
617 free(azCount);
618 free(evOffset);
619 free(coeffs);
620 free(delays);
621 return Hrtf;
622 }
623
624 static struct Hrtf *LoadHrtf01(const ALubyte *data, size_t datalen, const char *filename)
625 {
626 struct Hrtf *Hrtf = NULL;
627 ALboolean failed = AL_FALSE;
628 ALuint rate = 0;
629 ALushort irCount = 0;
630 ALushort irSize = 0;
631 ALubyte evCount = 0;
632 const ALubyte *azCount = NULL;
633 ALushort *evOffset = NULL;
634 ALfloat (*coeffs)[2] = NULL;
635 ALubyte (*delays)[2] = NULL;
636 ALsizei i, j;
637
638 if(datalen < 6)
639 {
640 ERR("Unexpected end of %s data (req %d, rem "SZFMT"\n", filename, 6, datalen);
641 return NULL;
642 }
643
644 rate = GetLE_ALuint(&data, &datalen);
645
646 irSize = GetLE_ALubyte(&data, &datalen);
647
648 evCount = GetLE_ALubyte(&data, &datalen);
649
650 if(irSize < MIN_IR_SIZE || irSize > MAX_IR_SIZE || (irSize%MOD_IR_SIZE))
651 {
652 ERR("Unsupported HRIR size: irSize=%d (%d to %d by %d)\n",
653 irSize, MIN_IR_SIZE, MAX_IR_SIZE, MOD_IR_SIZE);
654 failed = AL_TRUE;
655 }
656 if(evCount < MIN_EV_COUNT || evCount > MAX_EV_COUNT)
657 {
658 ERR("Unsupported elevation count: evCount=%d (%d to %d)\n",
659 evCount, MIN_EV_COUNT, MAX_EV_COUNT);
660 failed = AL_TRUE;
661 }
662 if(failed)
663 return NULL;
664
665 if(datalen < evCount)
666 {
667 ERR("Unexpected end of %s data (req %d, rem "SZFMT"\n", filename, evCount, datalen);
668 return NULL;
669 }
670
671 azCount = Get_ALubytePtr(&data, &datalen, evCount);
672
673 evOffset = malloc(sizeof(evOffset[0])*evCount);
674 if(azCount == NULL || evOffset == NULL)
675 {
676 ERR("Out of memory.\n");
677 failed = AL_TRUE;
678 }
679
680 if(!failed)
681 {
682 for(i = 0;i < evCount;i++)
683 {
684 if(azCount[i] < MIN_AZ_COUNT || azCount[i] > MAX_AZ_COUNT)
685 {
686 ERR("Unsupported azimuth count: azCount[%d]=%d (%d to %d)\n",
687 i, azCount[i], MIN_AZ_COUNT, MAX_AZ_COUNT);
688 failed = AL_TRUE;
689 }
690 }
691 }
692
693 if(!failed)
694 {
695 evOffset[0] = 0;
696 irCount = azCount[0];
697 for(i = 1;i < evCount;i++)
698 {
699 evOffset[i] = evOffset[i-1] + azCount[i-1];
700 irCount += azCount[i];
701 }
702
703 coeffs = malloc(sizeof(coeffs[0])*irSize*irCount);
704 delays = malloc(sizeof(delays[0])*irCount);
705 if(coeffs == NULL || delays == NULL)
706 {
707 ERR("Out of memory.\n");
708 failed = AL_TRUE;
709 }
710 }
711
712 if(!failed)
713 {
714 size_t reqsize = 2*irSize*irCount + irCount;
715 if(datalen < reqsize)
716 {
717 ERR("Unexpected end of %s data (req "SZFMT", rem "SZFMT"\n",
718 filename, reqsize, datalen);
719 failed = AL_TRUE;
720 }
721 }
722
723 if(!failed)
724 {
725 for(i = 0;i < irCount;i++)
726 {
727 for(j = 0;j < irSize;j++)
728 coeffs[i*irSize + j][0] = GetLE_ALshort(&data, &datalen) / 32768.0f;
729 }
730
731 for(i = 0;i < irCount;i++)
732 {
733 delays[i][0] = GetLE_ALubyte(&data, &datalen);
734 if(delays[i][0] > MAX_HRIR_DELAY)
735 {
736 ERR("Invalid delays[%d]: %d (%d)\n", i, delays[i][0], MAX_HRIR_DELAY);
737 failed = AL_TRUE;
738 }
739 }
740 }
741
742 if(!failed)
743 {
744 /* Mirror the left ear responses to the right ear. */
745 for(i = 0;i < evCount;i++)
746 {
747 ALushort evoffset = evOffset[i];
748 ALubyte azcount = azCount[i];
749 for(j = 0;j < azcount;j++)
750 {
751 ALsizei lidx = evoffset + j;
752 ALsizei ridx = evoffset + ((azcount-j) % azcount);
753 ALsizei k;
754
755 for(k = 0;k < irSize;k++)
756 coeffs[ridx*irSize + k][1] = coeffs[lidx*irSize + k][0];
757 delays[ridx][1] = delays[lidx][0];
758 }
759 }
760
761 Hrtf = CreateHrtfStore(rate, irSize, 0.0f, evCount, irCount, azCount,
762 evOffset, coeffs, delays, filename);
763 }
764
765 free(evOffset);
766 free(coeffs);
767 free(delays);
768 return Hrtf;
769 }
770
771 #define SAMPLETYPE_S16 0
772 #define SAMPLETYPE_S24 1
773
774 #define CHANTYPE_LEFTONLY 0
775 #define CHANTYPE_LEFTRIGHT 1
776
777 static struct Hrtf *LoadHrtf02(const ALubyte *data, size_t datalen, const char *filename)
778 {
779 struct Hrtf *Hrtf = NULL;
780 ALboolean failed = AL_FALSE;
781 ALuint rate = 0;
782 ALubyte sampleType;
783 ALubyte channelType;
784 ALushort irCount = 0;
785 ALushort irSize = 0;
786 ALubyte fdCount = 0;
787 ALushort distance = 0;
788 ALubyte evCount = 0;
789 const ALubyte *azCount = NULL;
790 ALushort *evOffset = NULL;
791 ALfloat (*coeffs)[2] = NULL;
792 ALubyte (*delays)[2] = NULL;
793 ALsizei i, j;
794
795 if(datalen < 8)
796 {
797 ERR("Unexpected end of %s data (req %d, rem "SZFMT"\n", filename, 8, datalen);
798 return NULL;
799 }
800
801 rate = GetLE_ALuint(&data, &datalen);
802 sampleType = GetLE_ALubyte(&data, &datalen);
803 channelType = GetLE_ALubyte(&data, &datalen);
804
805 irSize = GetLE_ALubyte(&data, &datalen);
806
807 fdCount = GetLE_ALubyte(&data, &datalen);
808
809 if(sampleType > SAMPLETYPE_S24)
810 {
811 ERR("Unsupported sample type: %d\n", sampleType);
812 failed = AL_TRUE;
813 }
814 if(channelType > CHANTYPE_LEFTRIGHT)
815 {
816 ERR("Unsupported channel type: %d\n", channelType);
817 failed = AL_TRUE;
818 }
819
820 if(irSize < MIN_IR_SIZE || irSize > MAX_IR_SIZE || (irSize%MOD_IR_SIZE))
821 {
822 ERR("Unsupported HRIR size: irSize=%d (%d to %d by %d)\n",
823 irSize, MIN_IR_SIZE, MAX_IR_SIZE, MOD_IR_SIZE);
824 failed = AL_TRUE;
825 }
826 if(fdCount != 1)
827 {
828 ERR("Multiple field-depths not supported: fdCount=%d (%d to %d)\n",
829 evCount, MIN_FD_COUNT, MAX_FD_COUNT);
830 failed = AL_TRUE;
831 }
832 if(failed)
833 return NULL;
834
835 for(i = 0;i < fdCount;i++)
836 {
837 if(datalen < 3)
838 {
839 ERR("Unexpected end of %s data (req %d, rem "SZFMT"\n", filename, 3, datalen);
840 return NULL;
841 }
842
843 distance = GetLE_ALushort(&data, &datalen);
844 if(distance < MIN_FD_DISTANCE || distance > MAX_FD_DISTANCE)
845 {
846 ERR("Unsupported field distance: distance=%d (%dmm to %dmm)\n",
847 distance, MIN_FD_DISTANCE, MAX_FD_DISTANCE);
848 failed = AL_TRUE;
849 }
850
851 evCount = GetLE_ALubyte(&data, &datalen);
852 if(evCount < MIN_EV_COUNT || evCount > MAX_EV_COUNT)
853 {
854 ERR("Unsupported elevation count: evCount=%d (%d to %d)\n",
855 evCount, MIN_EV_COUNT, MAX_EV_COUNT);
856 failed = AL_TRUE;
857 }
858 if(failed)
859 return NULL;
860
861 if(datalen < evCount)
862 {
863 ERR("Unexpected end of %s data (req %d, rem "SZFMT"\n", filename, evCount, datalen);
864 return NULL;
865 }
866
867 azCount = Get_ALubytePtr(&data, &datalen, evCount);
868 for(j = 0;j < evCount;j++)
869 {
870 if(azCount[j] < MIN_AZ_COUNT || azCount[j] > MAX_AZ_COUNT)
871 {
872 ERR("Unsupported azimuth count: azCount[%d]=%d (%d to %d)\n",
873 j, azCount[j], MIN_AZ_COUNT, MAX_AZ_COUNT);
874 failed = AL_TRUE;
875 }
876 }
877 }
878 if(failed)
879 return NULL;
880
881 evOffset = malloc(sizeof(evOffset[0])*evCount);
882 if(azCount == NULL || evOffset == NULL)
883 {
884 ERR("Out of memory.\n");
885 failed = AL_TRUE;
886 }
887
888 if(!failed)
889 {
890 evOffset[0] = 0;
891 irCount = azCount[0];
892 for(i = 1;i < evCount;i++)
893 {
894 evOffset[i] = evOffset[i-1] + azCount[i-1];
895 irCount += azCount[i];
896 }
897
898 coeffs = malloc(sizeof(coeffs[0])*irSize*irCount);
899 delays = malloc(sizeof(delays[0])*irCount);
900 if(coeffs == NULL || delays == NULL)
901 {
902 ERR("Out of memory.\n");
903 failed = AL_TRUE;
904 }
905 }
906
907 if(!failed)
908 {
909 size_t reqsize = 2*irSize*irCount + irCount;
910 if(datalen < reqsize)
911 {
912 ERR("Unexpected end of %s data (req "SZFMT", rem "SZFMT"\n",
913 filename, reqsize, datalen);
914 failed = AL_TRUE;
915 }
916 }
917
918 if(!failed)
919 {
920 if(channelType == CHANTYPE_LEFTONLY)
921 {
922 if(sampleType == SAMPLETYPE_S16)
923 for(i = 0;i < irCount;i++)
924 {
925 for(j = 0;j < irSize;j++)
926 coeffs[i*irSize + j][0] = GetLE_ALshort(&data, &datalen) / 32768.0f;
927 }
928 else if(sampleType == SAMPLETYPE_S24)
929 for(i = 0;i < irCount;i++)
930 {
931 for(j = 0;j < irSize;j++)
932 coeffs[i*irSize + j][0] = GetLE_ALint24(&data, &datalen) / 8388608.0f;
933 }
934
935 for(i = 0;i < irCount;i++)
936 {
937 delays[i][0] = GetLE_ALubyte(&data, &datalen);
938 if(delays[i][0] > MAX_HRIR_DELAY)
939 {
940 ERR("Invalid delays[%d][0]: %d (%d)\n", i, delays[i][0], MAX_HRIR_DELAY);
941 failed = AL_TRUE;
942 }
943 }
944 }
945 else if(channelType == CHANTYPE_LEFTRIGHT)
946 {
947 if(sampleType == SAMPLETYPE_S16)
948 for(i = 0;i < irCount;i++)
949 {
950 for(j = 0;j < irSize;j++)
951 {
952 coeffs[i*irSize + j][0] = GetLE_ALshort(&data, &datalen) / 32768.0f;
953 coeffs[i*irSize + j][1] = GetLE_ALshort(&data, &datalen) / 32768.0f;
954 }
955 }
956 else if(sampleType == SAMPLETYPE_S24)
957 for(i = 0;i < irCount;i++)
958 {
959 for(j = 0;j < irSize;j++)
960 {
961 coeffs[i*irSize + j][0] = GetLE_ALint24(&data, &datalen) / 8388608.0f;
962 coeffs[i*irSize + j][1] = GetLE_ALint24(&data, &datalen) / 8388608.0f;
963 }
964 }
965
966 for(i = 0;i < irCount;i++)
967 {
968 delays[i][0] = GetLE_ALubyte(&data, &datalen);
969 if(delays[i][0] > MAX_HRIR_DELAY)
970 {
971 ERR("Invalid delays[%d][0]: %d (%d)\n", i, delays[i][0], MAX_HRIR_DELAY);
972 failed = AL_TRUE;
973 }
974 delays[i][1] = GetLE_ALubyte(&data, &datalen);
975 if(delays[i][1] > MAX_HRIR_DELAY)
976 {
977 ERR("Invalid delays[%d][1]: %d (%d)\n", i, delays[i][1], MAX_HRIR_DELAY);
978 failed = AL_TRUE;
979 }
980 }
981 }
982 }
983
984 if(!failed)
985 {
986 if(channelType == CHANTYPE_LEFTONLY)
987 {
988 /* Mirror the left ear responses to the right ear. */
989 for(i = 0;i < evCount;i++)
990 {
991 ALushort evoffset = evOffset[i];
992 ALubyte azcount = azCount[i];
993 for(j = 0;j < azcount;j++)
994 {
995 ALsizei lidx = evoffset + j;
996 ALsizei ridx = evoffset + ((azcount-j) % azcount);
997 ALsizei k;
998
999 for(k = 0;k < irSize;k++)
1000 coeffs[ridx*irSize + k][1] = coeffs[lidx*irSize + k][0];
1001 delays[ridx][1] = delays[lidx][0];
1002 }
1003 }
1004 }
1005
1006 Hrtf = CreateHrtfStore(rate, irSize,
1007 (ALfloat)distance / 1000.0f, evCount, irCount, azCount, evOffset,
1008 coeffs, delays, filename
1009 );
1010 }
1011
1012 free(evOffset);
1013 free(coeffs);
1014 free(delays);
1015 return Hrtf;
1016 }
1017
1018
1019 static void AddFileEntry(vector_EnumeratedHrtf *list, const_al_string filename)
1020 {
1021 EnumeratedHrtf entry = { AL_STRING_INIT_STATIC(), NULL };
1022 struct HrtfEntry *loaded_entry;
1023 const EnumeratedHrtf *iter;
1024 const char *name;
1025 const char *ext;
1026 int i;
1027
1028 /* Check if this file has already been loaded globally. */
1029 loaded_entry = LoadedHrtfs;
1030 while(loaded_entry)
1031 {
1032 if(alstr_cmp_cstr(filename, loaded_entry->filename) == 0)
1033 {
1034 /* Check if this entry has already been added to the list. */
1035 #define MATCH_ENTRY(i) (loaded_entry == (i)->hrtf)
1036 VECTOR_FIND_IF(iter, const EnumeratedHrtf, *list, MATCH_ENTRY);
1037 if(iter != VECTOR_END(*list))
1038 {
1039 TRACE("Skipping duplicate file entry %s\n", alstr_get_cstr(filename));
1040 return;
1041 }
1042 #undef MATCH_FNAME
1043
1044 break;
1045 }
1046 loaded_entry = loaded_entry->next;
1047 }
1048
1049 if(!loaded_entry)
1050 {
1051 TRACE("Got new file \"%s\"\n", alstr_get_cstr(filename));
1052
1053 loaded_entry = al_calloc(DEF_ALIGN,
1054 FAM_SIZE(struct HrtfEntry, filename, alstr_length(filename)+1)
1055 );
1056 loaded_entry->next = LoadedHrtfs;
1057 loaded_entry->handle = NULL;
1058 strcpy(loaded_entry->filename, alstr_get_cstr(filename));
1059 LoadedHrtfs = loaded_entry;
1060 }
1061
1062 /* TODO: Get a human-readable name from the HRTF data (possibly coming in a
1063 * format update). */
1064 name = strrchr(alstr_get_cstr(filename), '/');
1065 if(!name) name = strrchr(alstr_get_cstr(filename), '\\');
1066 if(!name) name = alstr_get_cstr(filename);
1067 else ++name;
1068
1069 ext = strrchr(name, '.');
1070
1071 i = 0;
1072 do {
1073 if(!ext)
1074 alstr_copy_cstr(&entry.name, name);
1075 else
1076 alstr_copy_range(&entry.name, name, ext);
1077 if(i != 0)
1078 {
1079 char str[64];
1080 snprintf(str, sizeof(str), " #%d", i+1);
1081 alstr_append_cstr(&entry.name, str);
1082 }
1083 ++i;
1084
1085 #define MATCH_NAME(i) (alstr_cmp(entry.name, (i)->name) == 0)
1086 VECTOR_FIND_IF(iter, const EnumeratedHrtf, *list, MATCH_NAME);
1087 #undef MATCH_NAME
1088 } while(iter != VECTOR_END(*list));
1089 entry.hrtf = loaded_entry;
1090
1091 TRACE("Adding entry \"%s\" from file \"%s\"\n", alstr_get_cstr(entry.name),
1092 alstr_get_cstr(filename));
1093 VECTOR_PUSH_BACK(*list, entry);
1094 }
1095
1096 /* Unfortunate that we have to duplicate AddFileEntry to take a memory buffer
1097 * for input instead of opening the given filename.
1098 */
1099 static void AddBuiltInEntry(vector_EnumeratedHrtf *list, const_al_string filename, ALuint residx)
1100 {
1101 EnumeratedHrtf entry = { AL_STRING_INIT_STATIC(), NULL };
1102 struct HrtfEntry *loaded_entry;
1103 struct Hrtf *hrtf = NULL;
1104 const EnumeratedHrtf *iter;
1105 const char *name;
1106 const char *ext;
1107 int i;
1108
1109 loaded_entry = LoadedHrtfs;
1110 while(loaded_entry)
1111 {
1112 if(alstr_cmp_cstr(filename, loaded_entry->filename) == 0)
1113 {
1114 #define MATCH_ENTRY(i) (loaded_entry == (i)->hrtf)
1115 VECTOR_FIND_IF(iter, const EnumeratedHrtf, *list, MATCH_ENTRY);
1116 if(iter != VECTOR_END(*list))
1117 {
1118 TRACE("Skipping duplicate file entry %s\n", alstr_get_cstr(filename));
1119 return;
1120 }
1121 #undef MATCH_FNAME
1122
1123 break;
1124 }
1125 loaded_entry = loaded_entry->next;
1126 }
1127
1128 if(!loaded_entry)
1129 {
1130 size_t namelen = alstr_length(filename)+32;
1131
1132 TRACE("Got new file \"%s\"\n", alstr_get_cstr(filename));
1133
1134 loaded_entry = al_calloc(DEF_ALIGN,
1135 FAM_SIZE(struct HrtfEntry, filename, namelen)
1136 );
1137 loaded_entry->next = LoadedHrtfs;
1138 loaded_entry->handle = hrtf;
1139 snprintf(loaded_entry->filename, namelen, "!%u_%s",
1140 residx, alstr_get_cstr(filename));
1141 LoadedHrtfs = loaded_entry;
1142 }
1143
1144 /* TODO: Get a human-readable name from the HRTF data (possibly coming in a
1145 * format update). */
1146 name = strrchr(alstr_get_cstr(filename), '/');
1147 if(!name) name = strrchr(alstr_get_cstr(filename), '\\');
1148 if(!name) name = alstr_get_cstr(filename);
1149 else ++name;
1150
1151 ext = strrchr(name, '.');
1152
1153 i = 0;
1154 do {
1155 if(!ext)
1156 alstr_copy_cstr(&entry.name, name);
1157 else
1158 alstr_copy_range(&entry.name, name, ext);
1159 if(i != 0)
1160 {
1161 char str[64];
1162 snprintf(str, sizeof(str), " #%d", i+1);
1163 alstr_append_cstr(&entry.name, str);
1164 }
1165 ++i;
1166
1167 #define MATCH_NAME(i) (alstr_cmp(entry.name, (i)->name) == 0)
1168 VECTOR_FIND_IF(iter, const EnumeratedHrtf, *list, MATCH_NAME);
1169 #undef MATCH_NAME
1170 } while(iter != VECTOR_END(*list));
1171 entry.hrtf = loaded_entry;
1172
1173 TRACE("Adding built-in entry \"%s\"\n", alstr_get_cstr(entry.name));
1174 VECTOR_PUSH_BACK(*list, entry);
1175 }
1176
1177
1178 #define IDR_DEFAULT_44100_MHR 1
1179 #define IDR_DEFAULT_48000_MHR 2
1180
1181 #ifndef ALSOFT_EMBED_HRTF_DATA
1182
1183 static const ALubyte *GetResource(int UNUSED(name), size_t *size)
1184 {
1185 *size = 0;
1186 return NULL;
1187 }
1188
1189 #else
1190
1191 #include "default-44100.mhr.h"
1192 #include "default-48000.mhr.h"
1193
1194 static const ALubyte *GetResource(int name, size_t *size)
1195 {
1196 if(name == IDR_DEFAULT_44100_MHR)
1197 {
1198 *size = sizeof(hrtf_default_44100);
1199 return hrtf_default_44100;
1200 }
1201 if(name == IDR_DEFAULT_48000_MHR)
1202 {
1203 *size = sizeof(hrtf_default_48000);
1204 return hrtf_default_48000;
1205 }
1206 *size = 0;
1207 return NULL;
1208 }
1209 #endif
1210
1211 vector_EnumeratedHrtf EnumerateHrtf(const_al_string devname)
1212 {
1213 vector_EnumeratedHrtf list = VECTOR_INIT_STATIC();
1214 const char *defaulthrtf = "";
1215 const char *pathlist = "";
1216 bool usedefaults = true;
1217
1218 if(ConfigValueStr(alstr_get_cstr(devname), NULL, "hrtf-paths", &pathlist))
1219 {
1220 al_string pname = AL_STRING_INIT_STATIC();
1221 while(pathlist && *pathlist)
1222 {
1223 const char *next, *end;
1224
1225 while(isspace(*pathlist) || *pathlist == ',')
1226 pathlist++;
1227 if(*pathlist == '\0')
1228 continue;
1229
1230 next = strchr(pathlist, ',');
1231 if(next)
1232 end = next++;
1233 else
1234 {
1235 end = pathlist + strlen(pathlist);
1236 usedefaults = false;
1237 }
1238
1239 while(end != pathlist && isspace(*(end-1)))
1240 --end;
1241 if(end != pathlist)
1242 {
1243 vector_al_string flist;
1244 size_t i;
1245
1246 alstr_copy_range(&pname, pathlist, end);
1247
1248 flist = SearchDataFiles(".mhr", alstr_get_cstr(pname));
1249 for(i = 0;i < VECTOR_SIZE(flist);i++)
1250 AddFileEntry(&list, VECTOR_ELEM(flist, i));
1251 VECTOR_FOR_EACH(al_string, flist, alstr_reset);
1252 VECTOR_DEINIT(flist);
1253 }
1254
1255 pathlist = next;
1256 }
1257
1258 alstr_reset(&pname);
1259 }
1260 else if(ConfigValueExists(alstr_get_cstr(devname), NULL, "hrtf_tables"))
1261 ERR("The hrtf_tables option is deprecated, please use hrtf-paths instead.\n");
1262
1263 if(usedefaults)
1264 {
1265 al_string ename = AL_STRING_INIT_STATIC();
1266 vector_al_string flist;
1267 const ALubyte *rdata;
1268 size_t rsize, i;
1269
1270 flist = SearchDataFiles(".mhr", "openal/hrtf");
1271 for(i = 0;i < VECTOR_SIZE(flist);i++)
1272 AddFileEntry(&list, VECTOR_ELEM(flist, i));
1273 VECTOR_FOR_EACH(al_string, flist, alstr_reset);
1274 VECTOR_DEINIT(flist);
1275
1276 rdata = GetResource(IDR_DEFAULT_44100_MHR, &rsize);
1277 if(rdata != NULL && rsize > 0)
1278 {
1279 alstr_copy_cstr(&ename, "Built-In 44100hz");
1280 AddBuiltInEntry(&list, ename, IDR_DEFAULT_44100_MHR);
1281 }
1282
1283 rdata = GetResource(IDR_DEFAULT_48000_MHR, &rsize);
1284 if(rdata != NULL && rsize > 0)
1285 {
1286 alstr_copy_cstr(&ename, "Built-In 48000hz");
1287 AddBuiltInEntry(&list, ename, IDR_DEFAULT_48000_MHR);
1288 }
1289 alstr_reset(&ename);
1290 }
1291
1292 if(VECTOR_SIZE(list) > 1 && ConfigValueStr(alstr_get_cstr(devname), NULL, "default-hrtf", &defaulthrtf))
1293 {
1294 const EnumeratedHrtf *iter;
1295 /* Find the preferred HRTF and move it to the front of the list. */
1296 #define FIND_ENTRY(i) (alstr_cmp_cstr((i)->name, defaulthrtf) == 0)
1297 VECTOR_FIND_IF(iter, const EnumeratedHrtf, list, FIND_ENTRY);
1298 #undef FIND_ENTRY
1299 if(iter == VECTOR_END(list))
1300 WARN("Failed to find default HRTF \"%s\"\n", defaulthrtf);
1301 else if(iter != VECTOR_BEGIN(list))
1302 {
1303 EnumeratedHrtf entry = *iter;
1304 memmove(&VECTOR_ELEM(list,1), &VECTOR_ELEM(list,0),
1305 (iter-VECTOR_BEGIN(list))*sizeof(EnumeratedHrtf));
1306 VECTOR_ELEM(list,0) = entry;
1307 }
1308 }
1309
1310 return list;
1311 }
1312
1313 void FreeHrtfList(vector_EnumeratedHrtf *list)
1314 {
1315 #define CLEAR_ENTRY(i) alstr_reset(&(i)->name)
1316 VECTOR_FOR_EACH(EnumeratedHrtf, *list, CLEAR_ENTRY);
1317 VECTOR_DEINIT(*list);
1318 #undef CLEAR_ENTRY
1319 }
1320
1321 struct Hrtf *GetLoadedHrtf(struct HrtfEntry *entry)
1322 {
1323 struct Hrtf *hrtf = NULL;
1324 struct FileMapping fmap;
1325 const ALubyte *rdata;
1326 const char *name;
1327 ALuint residx;
1328 size_t rsize;
1329 char ch;
1330
1331 while(ATOMIC_FLAG_TEST_AND_SET(&LoadedHrtfLock, almemory_order_seq_cst))
1332 althrd_yield();
1333
1334 if(entry->handle)
1335 {
1336 hrtf = entry->handle;
1337 Hrtf_IncRef(hrtf);
1338 goto done;
1339 }
1340
1341 fmap.ptr = NULL;
1342 fmap.len = 0;
1343 if(sscanf(entry->filename, "!%u%c", &residx, &ch) == 2 && ch == '_')
1344 {
1345 name = strchr(entry->filename, ch)+1;
1346
1347 TRACE("Loading %s...\n", name);
1348 rdata = GetResource(residx, &rsize);
1349 if(rdata == NULL || rsize == 0)
1350 {
1351 ERR("Could not get resource %u, %s\n", residx, name);
1352 goto done;
1353 }
1354 }
1355 else
1356 {
1357 name = entry->filename;
1358
1359 TRACE("Loading %s...\n", entry->filename);
1360 fmap = MapFileToMem(entry->filename);
1361 if(fmap.ptr == NULL)
1362 {
1363 ERR("Could not open %s\n", entry->filename);
1364 goto done;
1365 }
1366
1367 rdata = fmap.ptr;
1368 rsize = fmap.len;
1369 }
1370
1371 if(rsize < sizeof(magicMarker02))
1372 ERR("%s data is too short ("SZFMT" bytes)\n", name, rsize);
1373 else if(memcmp(rdata, magicMarker02, sizeof(magicMarker02)) == 0)
1374 {
1375 TRACE("Detected data set format v2\n");
1376 hrtf = LoadHrtf02(rdata+sizeof(magicMarker02),
1377 rsize-sizeof(magicMarker02), name
1378 );
1379 }
1380 else if(memcmp(rdata, magicMarker01, sizeof(magicMarker01)) == 0)
1381 {
1382 TRACE("Detected data set format v1\n");
1383 hrtf = LoadHrtf01(rdata+sizeof(magicMarker01),
1384 rsize-sizeof(magicMarker01), name
1385 );
1386 }
1387 else if(memcmp(rdata, magicMarker00, sizeof(magicMarker00)) == 0)
1388 {
1389 TRACE("Detected data set format v0\n");
1390 hrtf = LoadHrtf00(rdata+sizeof(magicMarker00),
1391 rsize-sizeof(magicMarker00), name
1392 );
1393 }
1394 else
1395 ERR("Invalid header in %s: \"%.8s\"\n", name, (const char*)rdata);
1396 if(fmap.ptr)
1397 UnmapFileMem(&fmap);
1398
1399 if(!hrtf)
1400 {
1401 ERR("Failed to load %s\n", name);
1402 goto done;
1403 }
1404 entry->handle = hrtf;
1405 Hrtf_IncRef(hrtf);
1406
1407 TRACE("Loaded HRTF support for format: %s %uhz\n",
1408 DevFmtChannelsString(DevFmtStereo), hrtf->sampleRate);
1409
1410 done:
1411 ATOMIC_FLAG_CLEAR(&LoadedHrtfLock, almemory_order_seq_cst);
1412 return hrtf;
1413 }
1414
1415
1416 void Hrtf_IncRef(struct Hrtf *hrtf)
1417 {
1418 uint ref = IncrementRef(&hrtf->ref);
1419 TRACEREF("%p increasing refcount to %u\n", hrtf, ref);
1420 }
1421
1422 void Hrtf_DecRef(struct Hrtf *hrtf)
1423 {
1424 struct HrtfEntry *Hrtf;
1425 uint ref = DecrementRef(&hrtf->ref);
1426 TRACEREF("%p decreasing refcount to %u\n", hrtf, ref);
1427 if(ref == 0)
1428 {
1429 while(ATOMIC_FLAG_TEST_AND_SET(&LoadedHrtfLock, almemory_order_seq_cst))
1430 althrd_yield();
1431
1432 Hrtf = LoadedHrtfs;
1433 while(Hrtf != NULL)
1434 {
1435 /* Need to double-check that it's still unused, as another device
1436 * could've reacquired this HRTF after its reference went to 0 and
1437 * before the lock was taken.
1438 */
1439 if(hrtf == Hrtf->handle && ReadRef(&hrtf->ref) == 0)
1440 {
1441 al_free(Hrtf->handle);
1442 Hrtf->handle = NULL;
1443 TRACE("Unloaded unused HRTF %s\n", Hrtf->filename);
1444 }
1445 Hrtf = Hrtf->next;
1446 }
1447
1448 ATOMIC_FLAG_CLEAR(&LoadedHrtfLock, almemory_order_seq_cst);
1449 }
1450 }
1451
1452
1453 void FreeHrtfs(void)
1454 {
1455 struct HrtfEntry *Hrtf = LoadedHrtfs;
1456 LoadedHrtfs = NULL;
1457
1458 while(Hrtf != NULL)
1459 {
1460 struct HrtfEntry *next = Hrtf->next;
1461 al_free(Hrtf->handle);
1462 al_free(Hrtf);
1463 Hrtf = next;
1464 }
1465 }
Software OpenAL implementation