|
1 /* |
|
2 * Copyright (C) 2010 Google Inc. All rights reserved. |
|
3 * |
|
4 * Redistribution and use in source and binary forms, with or without |
|
5 * modification, are permitted provided that the following conditions |
|
6 * are met: |
|
7 * |
|
8 * 1. Redistributions of source code must retain the above copyright |
|
9 * notice, this list of conditions and the following disclaimer. |
|
10 * 2. Redistributions in binary form must reproduce the above copyright |
|
11 * notice, this list of conditions and the following disclaimer in the |
|
12 * documentation and/or other materials provided with the distribution. |
|
13 * 3. Neither the name of Apple Computer, Inc. ("Apple") nor the names of |
|
14 * its contributors may be used to endorse or promote products derived |
|
15 * from this software without specific prior written permission. |
|
16 * |
|
17 * THIS SOFTWARE IS PROVIDED BY APPLE AND ITS CONTRIBUTORS "AS IS" AND ANY |
|
18 * EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED |
|
19 * WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE ARE |
|
20 * DISCLAIMED. IN NO EVENT SHALL APPLE OR ITS CONTRIBUTORS BE LIABLE FOR ANY |
|
21 * DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES |
|
22 * (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; |
|
23 * LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND |
|
24 * ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT |
|
25 * (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF |
|
26 * THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE. |
|
27 */ |
|
28 |
|
29 #include "HRTFDatabase.h" |
|
30 |
|
31 #include "HRTFElevation.h" |
|
32 |
|
33 using namespace std; |
|
34 |
|
35 namespace WebCore { |
|
36 |
|
37 const int HRTFDatabase::MinElevation = -45; |
|
38 const int HRTFDatabase::MaxElevation = 90; |
|
39 const unsigned HRTFDatabase::RawElevationAngleSpacing = 15; |
|
40 const unsigned HRTFDatabase::NumberOfRawElevations = 10; // -45 -> +90 (each 15 degrees) |
|
41 const unsigned HRTFDatabase::InterpolationFactor = 1; |
|
42 const unsigned HRTFDatabase::NumberOfTotalElevations = NumberOfRawElevations * InterpolationFactor; |
|
43 |
|
44 nsReturnRef<HRTFDatabase> HRTFDatabase::create(float sampleRate) |
|
45 { |
|
46 return nsReturnRef<HRTFDatabase>(new HRTFDatabase(sampleRate)); |
|
47 } |
|
48 |
|
49 HRTFDatabase::HRTFDatabase(float sampleRate) |
|
50 : m_sampleRate(sampleRate) |
|
51 { |
|
52 m_elevations.SetLength(NumberOfTotalElevations); |
|
53 |
|
54 unsigned elevationIndex = 0; |
|
55 for (int elevation = MinElevation; elevation <= MaxElevation; elevation += RawElevationAngleSpacing) { |
|
56 nsAutoRef<HRTFElevation> hrtfElevation(HRTFElevation::createBuiltin(elevation, sampleRate)); |
|
57 MOZ_ASSERT(hrtfElevation.get()); |
|
58 if (!hrtfElevation.get()) |
|
59 return; |
|
60 |
|
61 m_elevations[elevationIndex] = hrtfElevation.out(); |
|
62 elevationIndex += InterpolationFactor; |
|
63 } |
|
64 |
|
65 // Now, go back and interpolate elevations. |
|
66 if (InterpolationFactor > 1) { |
|
67 for (unsigned i = 0; i < NumberOfTotalElevations; i += InterpolationFactor) { |
|
68 unsigned j = (i + InterpolationFactor); |
|
69 if (j >= NumberOfTotalElevations) |
|
70 j = i; // for last elevation interpolate with itself |
|
71 |
|
72 // Create the interpolated convolution kernels and delays. |
|
73 for (unsigned jj = 1; jj < InterpolationFactor; ++jj) { |
|
74 float x = static_cast<float>(jj) / static_cast<float>(InterpolationFactor); |
|
75 m_elevations[i + jj] = HRTFElevation::createByInterpolatingSlices(m_elevations[i].get(), m_elevations[j].get(), x, sampleRate); |
|
76 MOZ_ASSERT(m_elevations[i + jj].get()); |
|
77 } |
|
78 } |
|
79 } |
|
80 } |
|
81 |
|
82 size_t HRTFDatabase::sizeOfIncludingThis(mozilla::MallocSizeOf aMallocSizeOf) const |
|
83 { |
|
84 size_t amount = aMallocSizeOf(this); |
|
85 amount += m_elevations.SizeOfExcludingThis(aMallocSizeOf); |
|
86 for (size_t i = 0; i < m_elevations.Length(); i++) { |
|
87 amount += m_elevations[i]->sizeOfIncludingThis(aMallocSizeOf); |
|
88 } |
|
89 |
|
90 return amount; |
|
91 } |
|
92 |
|
93 void HRTFDatabase::getKernelsFromAzimuthElevation(double azimuthBlend, unsigned azimuthIndex, double elevationAngle, HRTFKernel* &kernelL, HRTFKernel* &kernelR, |
|
94 double& frameDelayL, double& frameDelayR) |
|
95 { |
|
96 unsigned elevationIndex = indexFromElevationAngle(elevationAngle); |
|
97 MOZ_ASSERT(elevationIndex < m_elevations.Length() && m_elevations.Length() > 0); |
|
98 |
|
99 if (!m_elevations.Length()) { |
|
100 kernelL = 0; |
|
101 kernelR = 0; |
|
102 return; |
|
103 } |
|
104 |
|
105 if (elevationIndex > m_elevations.Length() - 1) |
|
106 elevationIndex = m_elevations.Length() - 1; |
|
107 |
|
108 HRTFElevation* hrtfElevation = m_elevations[elevationIndex].get(); |
|
109 MOZ_ASSERT(hrtfElevation); |
|
110 if (!hrtfElevation) { |
|
111 kernelL = 0; |
|
112 kernelR = 0; |
|
113 return; |
|
114 } |
|
115 |
|
116 hrtfElevation->getKernelsFromAzimuth(azimuthBlend, azimuthIndex, kernelL, kernelR, frameDelayL, frameDelayR); |
|
117 } |
|
118 |
|
119 unsigned HRTFDatabase::indexFromElevationAngle(double elevationAngle) |
|
120 { |
|
121 // Clamp to allowed range. |
|
122 elevationAngle = mozilla::clamped(elevationAngle, |
|
123 static_cast<double>(MinElevation), |
|
124 static_cast<double>(MaxElevation)); |
|
125 |
|
126 unsigned elevationIndex = static_cast<int>(InterpolationFactor * (elevationAngle - MinElevation) / RawElevationAngleSpacing); |
|
127 return elevationIndex; |
|
128 } |
|
129 |
|
130 } // namespace WebCore |