src/third_party/WebKit/Source/platform/audio/HRTFDatabase.cpp

   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 "config.h"
  30
  31 #if ENABLE(WEB_AUDIO)
  32
  33 #include "platform/audio/HRTFDatabase.h"
  34
  35
  36 using namespace std;
  37
  38 namespace blink {
  39
  40 const int HRTFDatabase::MinElevation = -45;
  41 const int HRTFDatabase::MaxElevation = 90;
  42 const unsigned HRTFDatabase::RawElevationAngleSpacing = 15;
  43 const unsigned HRTFDatabase::NumberOfRawElevations = 10; // -45 -> +90 (each 15 degrees)
  44 const unsigned HRTFDatabase::InterpolationFactor = 1;
  45 const unsigned HRTFDatabase::NumberOfTotalElevations = NumberOfRawElevations * InterpolationFactor;
  46
  47 PassOwnPtr<HRTFDatabase> HRTFDatabase::create(float sampleRate)
  48 {
  49     OwnPtr<HRTFDatabase> hrtfDatabase = adoptPtr(new HRTFDatabase(sampleRate));
  50     return hrtfDatabase.release();
  51 }
  52
  53 HRTFDatabase::HRTFDatabase(float sampleRate)
  54     : m_elevations(NumberOfTotalElevations)
  55     , m_sampleRate(sampleRate)
  56 {
  57     unsigned elevationIndex = 0;
  58     for (int elevation = MinElevation; elevation <= MaxElevation; elevation += RawElevationAngleSpacing) {
  59         OwnPtr<HRTFElevation> hrtfElevation = HRTFElevation::createForSubject("Composite", elevation, sampleRate);
  60         ASSERT(hrtfElevation.get());
  61         if (!hrtfElevation.get())
  62             return;
  63
  64         m_elevations[elevationIndex] = hrtfElevation.release();
  65         elevationIndex += InterpolationFactor;
  66     }
  67
  68     // Now, go back and interpolate elevations.
  69     if (InterpolationFactor > 1) {
  70         for (unsigned i = 0; i < NumberOfTotalElevations; i += InterpolationFactor) {
  71             unsigned j = (i + InterpolationFactor);
  72             if (j >= NumberOfTotalElevations)
  73                 j = i; // for last elevation interpolate with itself
  74
  75             // Create the interpolated convolution kernels and delays.
  76             for (unsigned jj = 1; jj < InterpolationFactor; ++jj) {
  77                 float x = static_cast<float>(jj) / static_cast<float>(InterpolationFactor);
  78                 m_elevations[i + jj] = HRTFElevation::createByInterpolatingSlices(m_elevations[i].get(), m_elevations[j].get(), x, sampleRate);
  79                 ASSERT(m_elevations[i + jj].get());
  80             }
  81         }
  82     }
  83 }
  84
  85 void HRTFDatabase::getKernelsFromAzimuthElevation(double azimuthBlend, unsigned azimuthIndex, double elevationAngle, HRTFKernel* &kernelL, HRTFKernel* &kernelR,
  86                                                   double& frameDelayL, double& frameDelayR)
  87 {
  88     unsigned elevationIndex = indexFromElevationAngle(elevationAngle);
  89     ASSERT_WITH_SECURITY_IMPLICATION(elevationIndex < m_elevations.size() && m_elevations.size() > 0);
  90
  91     if (!m_elevations.size()) {
  92         kernelL = 0;
  93         kernelR = 0;
  94         return;
  95     }
  96
  97     if (elevationIndex > m_elevations.size() - 1)
  98         elevationIndex = m_elevations.size() - 1;
  99
 100     HRTFElevation* hrtfElevation = m_elevations[elevationIndex].get();
 101     ASSERT(hrtfElevation);
 102     if (!hrtfElevation) {
 103         kernelL = 0;
 104         kernelR = 0;
 105         return;
 106     }
 107
 108     hrtfElevation->getKernelsFromAzimuth(azimuthBlend, azimuthIndex, kernelL, kernelR, frameDelayL, frameDelayR);
 109 }
 110
 111 unsigned HRTFDatabase::indexFromElevationAngle(double elevationAngle)
 112 {
 113     // Clamp to allowed range.
 114     elevationAngle = max(static_cast<double>(MinElevation), elevationAngle);
 115     elevationAngle = min(static_cast<double>(MaxElevation), elevationAngle);
 116
 117     unsigned elevationIndex = static_cast<int>(InterpolationFactor * (elevationAngle - MinElevation) / RawElevationAngleSpacing);
 118     return elevationIndex;
 119 }
 120
 121 } // namespace blink
 122
 123 #endif // ENABLE(WEB_AUDIO)