首頁 探索 說明
登入
TorqueEngine
/
Torque3D-clone
2
0
複刻 0
檔案 問題管理 0 合併請求 0 Wiki
分支: development
分支列表 標籤列表
Torque3D-clone
/
Engine
/
lib
/
openal-soft
/
core
/
ambidefs.cpp

ambidefs.cpp 26 KB
永久連結 文件歷史 原始文件

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210211212213214215216217218219220221222223224225226227228229230231232233234235236237238239240241242243244245246247248249250251252253254255256257258259260261262263264265266267268269270271272273274275276277278279280281282283284285286287288289290291292293294295296297298299300301302303304305306307308 
  1. 

  2. #include "config.h"
    3. 

  3. 

  4. #include "ambidefs.h"
    5. 

  5. 

  6. #include "alnumbers.h"
    7. 

  7. 

  8. 

  9. namespace {
    10. 

  10. 

  11. using AmbiChannelFloatArray = std::array<float,MaxAmbiChannels>;
    12. 

  12. 

  13. constexpr auto inv_sqrt2f = static_cast<float>(1.0/al::numbers::sqrt2);
    14. 

  14. constexpr auto inv_sqrt3f = static_cast<float>(1.0/al::numbers::sqrt3);
    15. 

  15. 

  16. 

  17. /* These HF gains are derived from the same 32-point speaker array. The scale
    18. 

  18.  * factor between orders represents the same scale factors for any (regular)
    19. 

  19.  * speaker array decoder. e.g. Given a first-order source and second-order
    20. 

  20.  * output, applying an HF scale of HFScales[1][0] / HFScales[2][0] to channel 0
    21. 

  21.  * will result in that channel being subsequently decoded for second-order as
    22. 

  22.  * if it was a first-order decoder for that same speaker array.
    23. 

  23.  */
    24. 

  24. constexpr std::array HFScales{
    25. 

  25.     std::array{4.000000000e+00f, 2.309401077e+00f, 1.192569588e+00f, 7.189495850e-01f},
    26. 

  26.     std::array{4.000000000e+00f, 2.309401077e+00f, 1.192569588e+00f, 7.189495850e-01f},
    27. 

  27.     std::array{2.981423970e+00f, 2.309401077e+00f, 1.192569588e+00f, 7.189495850e-01f},
    28. 

  28.     std::array{2.359168820e+00f, 2.031565936e+00f, 1.444598386e+00f, 7.189495850e-01f},
    29. 

  29.     /*std::array{1.947005434e+00f, 1.764337084e+00f, 1.424707344e+00f, 9.755104127e-01f, 4.784482742e-01f}, */
    30. 

  30. };
    31. 

  31. 

  32. /* Same as above, but using a 10-point horizontal-only speaker array. Should
    33. 

  33.  * only be used when the device is mixing in 2D B-Format for horizontal-only
    34. 

  34.  * output.
    35. 

  35.  */
    36. 

  36. constexpr std::array HFScales2D{
    37. 

  37.     std::array{2.236067977e+00f, 1.581138830e+00f, 9.128709292e-01f, 6.050756345e-01f},
    38. 

  38.     std::array{2.236067977e+00f, 1.581138830e+00f, 9.128709292e-01f, 6.050756345e-01f},
    39. 

  39.     std::array{1.825741858e+00f, 1.581138830e+00f, 9.128709292e-01f, 6.050756345e-01f},
    40. 

  40.     std::array{1.581138830e+00f, 1.460781803e+00f, 1.118033989e+00f, 6.050756345e-01f},
    41. 

  41.     /*std::array{1.414213562e+00f, 1.344997024e+00f, 1.144122806e+00f, 8.312538756e-01f, 4.370160244e-01f}, */
    42. 

  42. };
    43. 

  43. 

  44. 

  45. /* This calculates a first-order "upsampler" matrix. It combines a first-order
    46. 

  46.  * decoder matrix with a max-order encoder matrix, creating a matrix that
    47. 

  47.  * behaves as if the B-Format input signal is first decoded to a speaker array
    48. 

  48.  * at first-order, then those speaker feeds are encoded to a higher-order
    49. 

  49.  * signal. While not perfect, this should accurately encode a lower-order
    50. 

  50.  * signal into a higher-order signal.
    51. 

  51.  */
    52. 

  52. constexpr std::array FirstOrderDecoder{
    53. 

  53.     std::array{1.250000000e-01f,  1.250000000e-01f,  1.250000000e-01f,  1.250000000e-01f},
    54. 

  54.     std::array{1.250000000e-01f,  1.250000000e-01f,  1.250000000e-01f, -1.250000000e-01f},
    55. 

  55.     std::array{1.250000000e-01f, -1.250000000e-01f,  1.250000000e-01f,  1.250000000e-01f},
    56. 

  56.     std::array{1.250000000e-01f, -1.250000000e-01f,  1.250000000e-01f, -1.250000000e-01f},
    57. 

  57.     std::array{1.250000000e-01f,  1.250000000e-01f, -1.250000000e-01f,  1.250000000e-01f},
    58. 

  58.     std::array{1.250000000e-01f,  1.250000000e-01f, -1.250000000e-01f, -1.250000000e-01f},
    59. 

  59.     std::array{1.250000000e-01f, -1.250000000e-01f, -1.250000000e-01f,  1.250000000e-01f},
    60. 

  60.     std::array{1.250000000e-01f, -1.250000000e-01f, -1.250000000e-01f, -1.250000000e-01f},
    61. 

  61. };
    62. 

  62. constexpr std::array FirstOrderEncoder{
    63. 

  63.     CalcAmbiCoeffs( inv_sqrt3f,  inv_sqrt3f,  inv_sqrt3f),
    64. 

  64.     CalcAmbiCoeffs( inv_sqrt3f,  inv_sqrt3f, -inv_sqrt3f),
    65. 

  65.     CalcAmbiCoeffs(-inv_sqrt3f,  inv_sqrt3f,  inv_sqrt3f),
    66. 

  66.     CalcAmbiCoeffs(-inv_sqrt3f,  inv_sqrt3f, -inv_sqrt3f),
    67. 

  67.     CalcAmbiCoeffs( inv_sqrt3f, -inv_sqrt3f,  inv_sqrt3f),
    68. 

  68.     CalcAmbiCoeffs( inv_sqrt3f, -inv_sqrt3f, -inv_sqrt3f),
    69. 

  69.     CalcAmbiCoeffs(-inv_sqrt3f, -inv_sqrt3f,  inv_sqrt3f),
    70. 

  70.     CalcAmbiCoeffs(-inv_sqrt3f, -inv_sqrt3f, -inv_sqrt3f),
    71. 

  71. };
    72. 

  72. static_assert(FirstOrderDecoder.size() == FirstOrderEncoder.size(), "First-order mismatch");
    73. 

  73. 

  74. /* This calculates a 2D first-order "upsampler" matrix. Same as the first-order
    75. 

  75.  * matrix, just using a more optimized speaker array for horizontal-only
    76. 

  76.  * content.
    77. 

  77.  */
    78. 

  78. constexpr std::array FirstOrder2DDecoder{
    79. 

  79.     std::array{2.500000000e-01f,  2.041241452e-01f, 0.0f,  2.041241452e-01f},
    80. 

  80.     std::array{2.500000000e-01f,  2.041241452e-01f, 0.0f, -2.041241452e-01f},
    81. 

  81.     std::array{2.500000000e-01f, -2.041241452e-01f, 0.0f,  2.041241452e-01f},
    82. 

  82.     std::array{2.500000000e-01f, -2.041241452e-01f, 0.0f, -2.041241452e-01f},
    83. 

  83. };
    84. 

  84. constexpr std::array FirstOrder2DEncoder{
    85. 

  85.     CalcAmbiCoeffs( inv_sqrt2f, 0.0f,  inv_sqrt2f),
    86. 

  86.     CalcAmbiCoeffs( inv_sqrt2f, 0.0f, -inv_sqrt2f),
    87. 

  87.     CalcAmbiCoeffs(-inv_sqrt2f, 0.0f,  inv_sqrt2f),
    88. 

  88.     CalcAmbiCoeffs(-inv_sqrt2f, 0.0f, -inv_sqrt2f),
    89. 

  89. };
    90. 

  90. static_assert(FirstOrder2DDecoder.size() == FirstOrder2DEncoder.size(), "First-order 2D mismatch");
    91. 

  91. 

  92. 

  93. /* This calculates a second-order "upsampler" matrix. Same as the first-order
    94. 

  94.  * matrix, just using a slightly more dense speaker array suitable for second-
    95. 

  95.  * order content.
    96. 

  96.  */
    97. 

  97. constexpr std::array SecondOrderDecoder{
    98. 

  98.     std::array{8.333333333e-02f,  0.000000000e+00f, -7.588274978e-02f,  1.227808683e-01f,  0.000000000e+00f,  0.000000000e+00f, -1.591525047e-02f, -1.443375673e-01f,  1.167715449e-01f},
    99. 

  99.     std::array{8.333333333e-02f, -1.227808683e-01f,  0.000000000e+00f,  7.588274978e-02f, -1.443375673e-01f,  0.000000000e+00f, -9.316949906e-02f,  0.000000000e+00f, -7.216878365e-02f},
    100. 

  100.     std::array{8.333333333e-02f, -7.588274978e-02f,  1.227808683e-01f,  0.000000000e+00f,  0.000000000e+00f, -1.443375673e-01f,  1.090847495e-01f,  0.000000000e+00f, -4.460276122e-02f},
    101. 

  101.     std::array{8.333333333e-02f,  0.000000000e+00f,  7.588274978e-02f,  1.227808683e-01f,  0.000000000e+00f,  0.000000000e+00f, -1.591525047e-02f,  1.443375673e-01f,  1.167715449e-01f},
    102. 

  102.     std::array{8.333333333e-02f, -1.227808683e-01f,  0.000000000e+00f, -7.588274978e-02f,  1.443375673e-01f,  0.000000000e+00f, -9.316949906e-02f,  0.000000000e+00f, -7.216878365e-02f},
    103. 

  103.     std::array{8.333333333e-02f,  7.588274978e-02f, -1.227808683e-01f,  0.000000000e+00f,  0.000000000e+00f, -1.443375673e-01f,  1.090847495e-01f,  0.000000000e+00f, -4.460276122e-02f},
    104. 

  104.     std::array{8.333333333e-02f,  0.000000000e+00f, -7.588274978e-02f, -1.227808683e-01f,  0.000000000e+00f,  0.000000000e+00f, -1.591525047e-02f,  1.443375673e-01f,  1.167715449e-01f},
    105. 

  105.     std::array{8.333333333e-02f,  1.227808683e-01f,  0.000000000e+00f, -7.588274978e-02f, -1.443375673e-01f,  0.000000000e+00f, -9.316949906e-02f,  0.000000000e+00f, -7.216878365e-02f},
    106. 

  106.     std::array{8.333333333e-02f,  7.588274978e-02f,  1.227808683e-01f,  0.000000000e+00f,  0.000000000e+00f,  1.443375673e-01f,  1.090847495e-01f,  0.000000000e+00f, -4.460276122e-02f},
    107. 

  107.     std::array{8.333333333e-02f,  0.000000000e+00f,  7.588274978e-02f, -1.227808683e-01f,  0.000000000e+00f,  0.000000000e+00f, -1.591525047e-02f, -1.443375673e-01f,  1.167715449e-01f},
    108. 

  108.     std::array{8.333333333e-02f,  1.227808683e-01f,  0.000000000e+00f,  7.588274978e-02f,  1.443375673e-01f,  0.000000000e+00f, -9.316949906e-02f,  0.000000000e+00f, -7.216878365e-02f},
    109. 

  109.     std::array{8.333333333e-02f, -7.588274978e-02f, -1.227808683e-01f,  0.000000000e+00f,  0.000000000e+00f,  1.443375673e-01f,  1.090847495e-01f,  0.000000000e+00f, -4.460276122e-02f},
    110. 

  110. };
    111. 

  111. constexpr std::array SecondOrderEncoder{
    112. 

  112.     CalcAmbiCoeffs( 0.000000000e+00f, -5.257311121e-01f,  8.506508084e-01f),
    113. 

  113.     CalcAmbiCoeffs(-8.506508084e-01f,  0.000000000e+00f,  5.257311121e-01f),
    114. 

  114.     CalcAmbiCoeffs(-5.257311121e-01f,  8.506508084e-01f,  0.000000000e+00f),
    115. 

  115.     CalcAmbiCoeffs( 0.000000000e+00f,  5.257311121e-01f,  8.506508084e-01f),
    116. 

  116.     CalcAmbiCoeffs(-8.506508084e-01f,  0.000000000e+00f, -5.257311121e-01f),
    117. 

  117.     CalcAmbiCoeffs( 5.257311121e-01f, -8.506508084e-01f,  0.000000000e+00f),
    118. 

  118.     CalcAmbiCoeffs( 0.000000000e+00f, -5.257311121e-01f, -8.506508084e-01f),
    119. 

  119.     CalcAmbiCoeffs( 8.506508084e-01f,  0.000000000e+00f, -5.257311121e-01f),
    120. 

  120.     CalcAmbiCoeffs( 5.257311121e-01f,  8.506508084e-01f,  0.000000000e+00f),
    121. 

  121.     CalcAmbiCoeffs( 0.000000000e+00f,  5.257311121e-01f, -8.506508084e-01f),
    122. 

  122.     CalcAmbiCoeffs( 8.506508084e-01f,  0.000000000e+00f,  5.257311121e-01f),
    123. 

  123.     CalcAmbiCoeffs(-5.257311121e-01f, -8.506508084e-01f,  0.000000000e+00f),
    124. 

  124. };
    125. 

  125. static_assert(SecondOrderDecoder.size() == SecondOrderEncoder.size(), "Second-order mismatch");
    126. 

  126. 

  127. /* This calculates a 2D second-order "upsampler" matrix. Same as the second-
    128. 

  128.  * order matrix, just using a more optimized speaker array for horizontal-only
    129. 

  129.  * content.
    130. 

  130.  */
    131. 

  131. constexpr std::array SecondOrder2DDecoder{
    132. 

  132.     std::array{1.666666667e-01f, -9.622504486e-02f, 0.0f,  1.666666667e-01f, -1.490711985e-01f, 0.0f, 0.0f, 0.0f,  8.606629658e-02f},
    133. 

  133.     std::array{1.666666667e-01f, -1.924500897e-01f, 0.0f,  0.000000000e+00f,  0.000000000e+00f, 0.0f, 0.0f, 0.0f, -1.721325932e-01f},
    134. 

  134.     std::array{1.666666667e-01f, -9.622504486e-02f, 0.0f, -1.666666667e-01f,  1.490711985e-01f, 0.0f, 0.0f, 0.0f,  8.606629658e-02f},
    135. 

  135.     std::array{1.666666667e-01f,  9.622504486e-02f, 0.0f, -1.666666667e-01f, -1.490711985e-01f, 0.0f, 0.0f, 0.0f,  8.606629658e-02f},
    136. 

  136.     std::array{1.666666667e-01f,  1.924500897e-01f, 0.0f,  0.000000000e+00f,  0.000000000e+00f, 0.0f, 0.0f, 0.0f, -1.721325932e-01f},
    137. 

  137.     std::array{1.666666667e-01f,  9.622504486e-02f, 0.0f,  1.666666667e-01f,  1.490711985e-01f, 0.0f, 0.0f, 0.0f,  8.606629658e-02f},
    138. 

  138. };
    139. 

  139. constexpr std::array SecondOrder2DEncoder{
    140. 

  140.     CalcAmbiCoeffs(-0.50000000000f, 0.0f,  0.86602540379f),
    141. 

  141.     CalcAmbiCoeffs(-1.00000000000f, 0.0f,  0.00000000000f),
    142. 

  142.     CalcAmbiCoeffs(-0.50000000000f, 0.0f, -0.86602540379f),
    143. 

  143.     CalcAmbiCoeffs( 0.50000000000f, 0.0f, -0.86602540379f),
    144. 

  144.     CalcAmbiCoeffs( 1.00000000000f, 0.0f,  0.00000000000f),
    145. 

  145.     CalcAmbiCoeffs( 0.50000000000f, 0.0f,  0.86602540379f),
    146. 

  146. };
    147. 

  147. static_assert(SecondOrder2DDecoder.size() == SecondOrder2DEncoder.size(),
    148. 

  148.     "Second-order 2D mismatch");
    149. 

  149. 

  150. 

  151. /* This calculates a third-order "upsampler" matrix. Same as the first-order
    152. 

  152.  * matrix, just using a more dense speaker array suitable for third-order
    153. 

  153.  * content.
    154. 

  154.  */
    155. 

  155. constexpr std::array ThirdOrderDecoder{
    156. 

  156.     std::array{5.000000000e-02f,  3.090169944e-02f,  8.090169944e-02f,  0.000000000e+00f,  0.000000000e+00f,  6.454972244e-02f,  9.045084972e-02f,  0.000000000e+00f, -1.232790000e-02f, -1.256118221e-01f,  0.000000000e+00f,  1.126112056e-01f,  7.944389175e-02f,  0.000000000e+00f,  2.421151497e-02f,  0.000000000e+00f},
    157. 

  157.     std::array{5.000000000e-02f, -3.090169944e-02f,  8.090169944e-02f,  0.000000000e+00f,  0.000000000e+00f, -6.454972244e-02f,  9.045084972e-02f,  0.000000000e+00f, -1.232790000e-02f,  1.256118221e-01f,  0.000000000e+00f, -1.126112056e-01f,  7.944389175e-02f,  0.000000000e+00f,  2.421151497e-02f,  0.000000000e+00f},
    158. 

  158.     std::array{5.000000000e-02f,  3.090169944e-02f, -8.090169944e-02f,  0.000000000e+00f,  0.000000000e+00f, -6.454972244e-02f,  9.045084972e-02f,  0.000000000e+00f, -1.232790000e-02f, -1.256118221e-01f,  0.000000000e+00f,  1.126112056e-01f, -7.944389175e-02f,  0.000000000e+00f, -2.421151497e-02f,  0.000000000e+00f},
    159. 

  159.     std::array{5.000000000e-02f, -3.090169944e-02f, -8.090169944e-02f,  0.000000000e+00f,  0.000000000e+00f,  6.454972244e-02f,  9.045084972e-02f,  0.000000000e+00f, -1.232790000e-02f,  1.256118221e-01f,  0.000000000e+00f, -1.126112056e-01f, -7.944389175e-02f,  0.000000000e+00f, -2.421151497e-02f,  0.000000000e+00f},
    160. 

  160.     std::array{5.000000000e-02f,  8.090169944e-02f,  0.000000000e+00f,  3.090169944e-02f,  6.454972244e-02f,  0.000000000e+00f, -5.590169944e-02f,  0.000000000e+00f, -7.216878365e-02f, -7.763237543e-02f,  0.000000000e+00f, -2.950836627e-02f,  0.000000000e+00f, -1.497759251e-01f,  0.000000000e+00f, -7.763237543e-02f},
    161. 

  161.     std::array{5.000000000e-02f,  8.090169944e-02f,  0.000000000e+00f, -3.090169944e-02f, -6.454972244e-02f,  0.000000000e+00f, -5.590169944e-02f,  0.000000000e+00f, -7.216878365e-02f, -7.763237543e-02f,  0.000000000e+00f, -2.950836627e-02f,  0.000000000e+00f,  1.497759251e-01f,  0.000000000e+00f,  7.763237543e-02f},
    162. 

  162.     std::array{5.000000000e-02f, -8.090169944e-02f,  0.000000000e+00f,  3.090169944e-02f, -6.454972244e-02f,  0.000000000e+00f, -5.590169944e-02f,  0.000000000e+00f, -7.216878365e-02f,  7.763237543e-02f,  0.000000000e+00f,  2.950836627e-02f,  0.000000000e+00f, -1.497759251e-01f,  0.000000000e+00f, -7.763237543e-02f},
    163. 

  163.     std::array{5.000000000e-02f, -8.090169944e-02f,  0.000000000e+00f, -3.090169944e-02f,  6.454972244e-02f,  0.000000000e+00f, -5.590169944e-02f,  0.000000000e+00f, -7.216878365e-02f,  7.763237543e-02f,  0.000000000e+00f,  2.950836627e-02f,  0.000000000e+00f,  1.497759251e-01f,  0.000000000e+00f,  7.763237543e-02f},
    164. 

  164.     std::array{5.000000000e-02f,  0.000000000e+00f,  3.090169944e-02f,  8.090169944e-02f,  0.000000000e+00f,  0.000000000e+00f, -3.454915028e-02f,  6.454972244e-02f,  8.449668365e-02f,  0.000000000e+00f,  0.000000000e+00f,  0.000000000e+00f,  3.034486645e-02f, -6.779013272e-02f,  1.659481923e-01f,  4.797944664e-02f},
    165. 

  165.     std::array{5.000000000e-02f,  0.000000000e+00f,  3.090169944e-02f, -8.090169944e-02f,  0.000000000e+00f,  0.000000000e+00f, -3.454915028e-02f, -6.454972244e-02f,  8.449668365e-02f,  0.000000000e+00f,  0.000000000e+00f,  0.000000000e+00f,  3.034486645e-02f,  6.779013272e-02f,  1.659481923e-01f, -4.797944664e-02f},
    166. 

  166.     std::array{5.000000000e-02f,  0.000000000e+00f, -3.090169944e-02f,  8.090169944e-02f,  0.000000000e+00f,  0.000000000e+00f, -3.454915028e-02f, -6.454972244e-02f,  8.449668365e-02f,  0.000000000e+00f,  0.000000000e+00f,  0.000000000e+00f, -3.034486645e-02f, -6.779013272e-02f, -1.659481923e-01f,  4.797944664e-02f},
    167. 

  167.     std::array{5.000000000e-02f,  0.000000000e+00f, -3.090169944e-02f, -8.090169944e-02f,  0.000000000e+00f,  0.000000000e+00f, -3.454915028e-02f,  6.454972244e-02f,  8.449668365e-02f,  0.000000000e+00f,  0.000000000e+00f,  0.000000000e+00f, -3.034486645e-02f,  6.779013272e-02f, -1.659481923e-01f, -4.797944664e-02f},
    168. 

  168.     std::array{5.000000000e-02f,  5.000000000e-02f,  5.000000000e-02f,  5.000000000e-02f,  6.454972244e-02f,  6.454972244e-02f,  0.000000000e+00f,  6.454972244e-02f,  0.000000000e+00f,  1.016220987e-01f,  6.338656910e-02f, -1.092600649e-02f, -7.364853795e-02f,  1.011266756e-01f, -7.086833869e-02f, -1.482646439e-02f},
    169. 

  169.     std::array{5.000000000e-02f,  5.000000000e-02f,  5.000000000e-02f, -5.000000000e-02f, -6.454972244e-02f,  6.454972244e-02f,  0.000000000e+00f, -6.454972244e-02f,  0.000000000e+00f,  1.016220987e-01f, -6.338656910e-02f, -1.092600649e-02f, -7.364853795e-02f, -1.011266756e-01f, -7.086833869e-02f,  1.482646439e-02f},
    170. 

  170.     std::array{5.000000000e-02f, -5.000000000e-02f,  5.000000000e-02f,  5.000000000e-02f, -6.454972244e-02f, -6.454972244e-02f,  0.000000000e+00f,  6.454972244e-02f,  0.000000000e+00f, -1.016220987e-01f, -6.338656910e-02f,  1.092600649e-02f, -7.364853795e-02f,  1.011266756e-01f, -7.086833869e-02f, -1.482646439e-02f},
    171. 

  171.     std::array{5.000000000e-02f, -5.000000000e-02f,  5.000000000e-02f, -5.000000000e-02f,  6.454972244e-02f, -6.454972244e-02f,  0.000000000e+00f, -6.454972244e-02f,  0.000000000e+00f, -1.016220987e-01f,  6.338656910e-02f,  1.092600649e-02f, -7.364853795e-02f, -1.011266756e-01f, -7.086833869e-02f,  1.482646439e-02f},
    172. 

  172.     std::array{5.000000000e-02f,  5.000000000e-02f, -5.000000000e-02f,  5.000000000e-02f,  6.454972244e-02f, -6.454972244e-02f,  0.000000000e+00f, -6.454972244e-02f,  0.000000000e+00f,  1.016220987e-01f, -6.338656910e-02f, -1.092600649e-02f,  7.364853795e-02f,  1.011266756e-01f,  7.086833869e-02f, -1.482646439e-02f},
    173. 

  173.     std::array{5.000000000e-02f,  5.000000000e-02f, -5.000000000e-02f, -5.000000000e-02f, -6.454972244e-02f, -6.454972244e-02f,  0.000000000e+00f,  6.454972244e-02f,  0.000000000e+00f,  1.016220987e-01f,  6.338656910e-02f, -1.092600649e-02f,  7.364853795e-02f, -1.011266756e-01f,  7.086833869e-02f,  1.482646439e-02f},
    174. 

  174.     std::array{5.000000000e-02f, -5.000000000e-02f, -5.000000000e-02f,  5.000000000e-02f, -6.454972244e-02f,  6.454972244e-02f,  0.000000000e+00f, -6.454972244e-02f,  0.000000000e+00f, -1.016220987e-01f,  6.338656910e-02f,  1.092600649e-02f,  7.364853795e-02f,  1.011266756e-01f,  7.086833869e-02f, -1.482646439e-02f},
    175. 

  175.     std::array{5.000000000e-02f, -5.000000000e-02f, -5.000000000e-02f, -5.000000000e-02f,  6.454972244e-02f,  6.454972244e-02f,  0.000000000e+00f,  6.454972244e-02f,  0.000000000e+00f, -1.016220987e-01f, -6.338656910e-02f,  1.092600649e-02f,  7.364853795e-02f, -1.011266756e-01f,  7.086833869e-02f,  1.482646439e-02f},
    176. 

  176. };
    177. 

  177. constexpr std::array ThirdOrderEncoder{
    178. 

  178.     CalcAmbiCoeffs( 0.35682208976f,  0.93417235897f,  0.00000000000f),
    179. 

  179.     CalcAmbiCoeffs(-0.35682208976f,  0.93417235897f,  0.00000000000f),
    180. 

  180.     CalcAmbiCoeffs( 0.35682208976f, -0.93417235897f,  0.00000000000f),
    181. 

  181.     CalcAmbiCoeffs(-0.35682208976f, -0.93417235897f,  0.00000000000f),
    182. 

  182.     CalcAmbiCoeffs( 0.93417235897f,  0.00000000000f,  0.35682208976f),
    183. 

  183.     CalcAmbiCoeffs( 0.93417235897f,  0.00000000000f, -0.35682208976f),
    184. 

  184.     CalcAmbiCoeffs(-0.93417235897f,  0.00000000000f,  0.35682208976f),
    185. 

  185.     CalcAmbiCoeffs(-0.93417235897f,  0.00000000000f, -0.35682208976f),
    186. 

  186.     CalcAmbiCoeffs( 0.00000000000f,  0.35682208976f,  0.93417235897f),
    187. 

  187.     CalcAmbiCoeffs( 0.00000000000f,  0.35682208976f, -0.93417235897f),
    188. 

  188.     CalcAmbiCoeffs( 0.00000000000f, -0.35682208976f,  0.93417235897f),
    189. 

  189.     CalcAmbiCoeffs( 0.00000000000f, -0.35682208976f, -0.93417235897f),
    190. 

  190.     CalcAmbiCoeffs(     inv_sqrt3f,      inv_sqrt3f,      inv_sqrt3f),
    191. 

  191.     CalcAmbiCoeffs(     inv_sqrt3f,      inv_sqrt3f,     -inv_sqrt3f),
    192. 

  192.     CalcAmbiCoeffs(    -inv_sqrt3f,      inv_sqrt3f,      inv_sqrt3f),
    193. 

  193.     CalcAmbiCoeffs(    -inv_sqrt3f,      inv_sqrt3f,     -inv_sqrt3f),
    194. 

  194.     CalcAmbiCoeffs(     inv_sqrt3f,     -inv_sqrt3f,      inv_sqrt3f),
    195. 

  195.     CalcAmbiCoeffs(     inv_sqrt3f,     -inv_sqrt3f,     -inv_sqrt3f),
    196. 

  196.     CalcAmbiCoeffs(    -inv_sqrt3f,     -inv_sqrt3f,      inv_sqrt3f),
    197. 

  197.     CalcAmbiCoeffs(    -inv_sqrt3f,     -inv_sqrt3f,     -inv_sqrt3f),
    198. 

  198. };
    199. 

  199. static_assert(ThirdOrderDecoder.size() == ThirdOrderEncoder.size(), "Third-order mismatch");
    200. 

  200. 

  201. /* This calculates a 2D third-order "upsampler" matrix. Same as the third-order
    202. 

  202.  * matrix, just using a more optimized speaker array for horizontal-only
    203. 

  203.  * content.
    204. 

  204.  */
    205. 

  205. constexpr std::array ThirdOrder2DDecoder{
    206. 

  206.     std::array{1.250000000e-01f, -5.523559567e-02f, 0.0f,  1.333505242e-01f, -9.128709292e-02f, 0.0f, 0.0f, 0.0f,  9.128709292e-02f, -1.104247249e-01f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,  4.573941867e-02f},
    207. 

  207.     std::array{1.250000000e-01f, -1.333505242e-01f, 0.0f,  5.523559567e-02f, -9.128709292e-02f, 0.0f, 0.0f, 0.0f, -9.128709292e-02f,  4.573941867e-02f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, -1.104247249e-01f},
    208. 

  208.     std::array{1.250000000e-01f, -1.333505242e-01f, 0.0f, -5.523559567e-02f,  9.128709292e-02f, 0.0f, 0.0f, 0.0f, -9.128709292e-02f,  4.573941867e-02f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,  1.104247249e-01f},
    209. 

  209.     std::array{1.250000000e-01f, -5.523559567e-02f, 0.0f, -1.333505242e-01f,  9.128709292e-02f, 0.0f, 0.0f, 0.0f,  9.128709292e-02f, -1.104247249e-01f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, -4.573941867e-02f},
    210. 

  210.     std::array{1.250000000e-01f,  5.523559567e-02f, 0.0f, -1.333505242e-01f, -9.128709292e-02f, 0.0f, 0.0f, 0.0f,  9.128709292e-02f,  1.104247249e-01f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, -4.573941867e-02f},
    211. 

  211.     std::array{1.250000000e-01f,  1.333505242e-01f, 0.0f, -5.523559567e-02f, -9.128709292e-02f, 0.0f, 0.0f, 0.0f, -9.128709292e-02f, -4.573941867e-02f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,  1.104247249e-01f},
    212. 

  212.     std::array{1.250000000e-01f,  1.333505242e-01f, 0.0f,  5.523559567e-02f,  9.128709292e-02f, 0.0f, 0.0f, 0.0f, -9.128709292e-02f, -4.573941867e-02f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, -1.104247249e-01f},
    213. 

  213.     std::array{1.250000000e-01f,  5.523559567e-02f, 0.0f,  1.333505242e-01f,  9.128709292e-02f, 0.0f, 0.0f, 0.0f,  9.128709292e-02f,  1.104247249e-01f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,  4.573941867e-02f},
    214. 

  214. };
    215. 

  215. constexpr std::array ThirdOrder2DEncoder{
    216. 

  216.     CalcAmbiCoeffs(-0.38268343237f, 0.0f,  0.92387953251f),
    217. 

  217.     CalcAmbiCoeffs(-0.92387953251f, 0.0f,  0.38268343237f),
    218. 

  218.     CalcAmbiCoeffs(-0.92387953251f, 0.0f, -0.38268343237f),
    219. 

  219.     CalcAmbiCoeffs(-0.38268343237f, 0.0f, -0.92387953251f),
    220. 

  220.     CalcAmbiCoeffs( 0.38268343237f, 0.0f, -0.92387953251f),
    221. 

  221.     CalcAmbiCoeffs( 0.92387953251f, 0.0f, -0.38268343237f),
    222. 

  222.     CalcAmbiCoeffs( 0.92387953251f, 0.0f,  0.38268343237f),
    223. 

  223.     CalcAmbiCoeffs( 0.38268343237f, 0.0f,  0.92387953251f),
    224. 

  224. };
    225. 

  225. static_assert(ThirdOrder2DDecoder.size() == ThirdOrder2DEncoder.size(), "Third-order 2D mismatch");
    226. 

  226. 

  227. 

  228. /* This calculates a 2D fourth-order "upsampler" matrix. There is no 3D fourth-
    229. 

  229.  * order upsampler since fourth-order is the max order we'll be supporting for
    230. 

  230.  * the foreseeable future. This is only necessary for mixing horizontal-only
    231. 

  231.  * fourth-order content to 3D.
    232. 

  232.  */
    233. 

  233. constexpr std::array FourthOrder2DDecoder{
    234. 

  234.     std::array{1.000000000e-01f,  3.568220898e-02f, 0.0f,  1.098185471e-01f,  6.070619982e-02f, 0.0f, 0.0f, 0.0f,  8.355491589e-02f,  7.735682057e-02f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,  5.620301997e-02f,  8.573754253e-02f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,  2.785781628e-02f},
    235. 

  235.     std::array{1.000000000e-01f,  9.341723590e-02f, 0.0f,  6.787159473e-02f,  9.822469464e-02f, 0.0f, 0.0f, 0.0f, -3.191513794e-02f,  2.954767620e-02f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, -9.093839659e-02f, -5.298871540e-02f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, -7.293270986e-02f},
    236. 

  236.     std::array{1.000000000e-01f,  1.154700538e-01f, 0.0f,  0.000000000e+00f,  0.000000000e+00f, 0.0f, 0.0f, 0.0f, -1.032795559e-01f, -9.561828875e-02f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,  0.000000000e+00f,  0.000000000e+00f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,  9.014978717e-02f},
    237. 

  237.     std::array{1.000000000e-01f,  9.341723590e-02f, 0.0f, -6.787159473e-02f, -9.822469464e-02f, 0.0f, 0.0f, 0.0f, -3.191513794e-02f,  2.954767620e-02f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,  9.093839659e-02f,  5.298871540e-02f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, -7.293270986e-02f},
    238. 

  238.     std::array{1.000000000e-01f,  3.568220898e-02f, 0.0f, -1.098185471e-01f, -6.070619982e-02f, 0.0f, 0.0f, 0.0f,  8.355491589e-02f,  7.735682057e-02f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, -5.620301997e-02f, -8.573754253e-02f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,  2.785781628e-02f},
    239. 

  239.     std::array{1.000000000e-01f, -3.568220898e-02f, 0.0f, -1.098185471e-01f,  6.070619982e-02f, 0.0f, 0.0f, 0.0f,  8.355491589e-02f, -7.735682057e-02f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, -5.620301997e-02f,  8.573754253e-02f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,  2.785781628e-02f},
    240. 

  240.     std::array{1.000000000e-01f, -9.341723590e-02f, 0.0f, -6.787159473e-02f,  9.822469464e-02f, 0.0f, 0.0f, 0.0f, -3.191513794e-02f, -2.954767620e-02f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,  9.093839659e-02f, -5.298871540e-02f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, -7.293270986e-02f},
    241. 

  241.     std::array{1.000000000e-01f, -1.154700538e-01f, 0.0f,  0.000000000e+00f,  0.000000000e+00f, 0.0f, 0.0f, 0.0f, -1.032795559e-01f,  9.561828875e-02f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,  0.000000000e+00f,  0.000000000e+00f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,  9.014978717e-02f},
    242. 

  242.     std::array{1.000000000e-01f, -9.341723590e-02f, 0.0f,  6.787159473e-02f, -9.822469464e-02f, 0.0f, 0.0f, 0.0f, -3.191513794e-02f, -2.954767620e-02f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, -9.093839659e-02f,  5.298871540e-02f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, -7.293270986e-02f},
    243. 

  243.     std::array{1.000000000e-01f, -3.568220898e-02f, 0.0f,  1.098185471e-01f, -6.070619982e-02f, 0.0f, 0.0f, 0.0f,  8.355491589e-02f, -7.735682057e-02f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,  5.620301997e-02f, -8.573754253e-02f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f,  2.785781628e-02f},
    244. 

  244. };
    245. 

  245. constexpr std::array FourthOrder2DEncoder{
    246. 

  246.     CalcAmbiCoeffs( 3.090169944e-01f,  0.000000000e+00f,  9.510565163e-01f),
    247. 

  247.     CalcAmbiCoeffs( 8.090169944e-01f,  0.000000000e+00f,  5.877852523e-01f),
    248. 

  248.     CalcAmbiCoeffs( 1.000000000e+00f,  0.000000000e+00f,  0.000000000e+00f),
    249. 

  249.     CalcAmbiCoeffs( 8.090169944e-01f,  0.000000000e+00f, -5.877852523e-01f),
    250. 

  250.     CalcAmbiCoeffs( 3.090169944e-01f,  0.000000000e+00f, -9.510565163e-01f),
    251. 

  251.     CalcAmbiCoeffs(-3.090169944e-01f,  0.000000000e+00f, -9.510565163e-01f),
    252. 

  252.     CalcAmbiCoeffs(-8.090169944e-01f,  0.000000000e+00f, -5.877852523e-01f),
    253. 

  253.     CalcAmbiCoeffs(-1.000000000e+00f,  0.000000000e+00f,  0.000000000e+00f),
    254. 

  254.     CalcAmbiCoeffs(-8.090169944e-01f,  0.000000000e+00f,  5.877852523e-01f),
    255. 

  255.     CalcAmbiCoeffs(-3.090169944e-01f,  0.000000000e+00f,  9.510565163e-01f),
    256. 

  256. };
    257. 

  257. static_assert(FourthOrder2DDecoder.size() == FourthOrder2DEncoder.size(), "Fourth-order 2D mismatch");
    258. 

  258. 

  259. 

  260. template<size_t N, size_t M>
    261. 

  261. constexpr auto CalcAmbiUpsampler(const std::array<std::array<float,N>,M> &decoder,
    262. 

  262.     const std::array<AmbiChannelFloatArray,M> &encoder)
    263. 

  263. {
    264. 

  264.     std::array<AmbiChannelFloatArray,N> res{};
    265. 

  265. 

  266.     for(size_t i{0};i < decoder[0].size();++i)
    267. 

  267.     {
    268. 

  268.         for(size_t j{0};j < encoder[0].size();++j)
    269. 

  269.         {
    270. 

  270.             double sum{0.0};
    271. 

  271.             for(size_t k{0};k < decoder.size();++k)
    272. 

  272.                 sum += double{decoder[k][i]} * encoder[k][j];
    273. 

  273.             res[i][j] = static_cast<float>(sum);
    274. 

  274.         }
    275. 

  275.     }
    276. 

  276. 

  277.     return res;
    278. 

  278. }
    279. 

  279. 

  280. } // namespace
    281. 

  281. 

  282. const std::array<std::array<float,MaxAmbiChannels>,4> AmbiScale::FirstOrderUp{CalcAmbiUpsampler(FirstOrderDecoder, FirstOrderEncoder)};
    283. 

  283. const std::array<std::array<float,MaxAmbiChannels>,4> AmbiScale::FirstOrder2DUp{CalcAmbiUpsampler(FirstOrder2DDecoder, FirstOrder2DEncoder)};
    284. 

  284. const std::array<std::array<float,MaxAmbiChannels>,9> AmbiScale::SecondOrderUp{CalcAmbiUpsampler(SecondOrderDecoder, SecondOrderEncoder)};
    285. 

  285. const std::array<std::array<float,MaxAmbiChannels>,9> AmbiScale::SecondOrder2DUp{CalcAmbiUpsampler(SecondOrder2DDecoder, SecondOrder2DEncoder)};
    286. 

  286. const std::array<std::array<float,MaxAmbiChannels>,16> AmbiScale::ThirdOrderUp{CalcAmbiUpsampler(ThirdOrderDecoder, ThirdOrderEncoder)};
    287. 

  287. const std::array<std::array<float,MaxAmbiChannels>,16> AmbiScale::ThirdOrder2DUp{CalcAmbiUpsampler(ThirdOrder2DDecoder, ThirdOrder2DEncoder)};
    288. 

  288. const std::array<std::array<float,MaxAmbiChannels>,25> AmbiScale::FourthOrder2DUp{CalcAmbiUpsampler(FourthOrder2DDecoder, FourthOrder2DEncoder)};
    289. 

  289. 

  290. 

  291. std::array<float,MaxAmbiOrder+1> AmbiScale::GetHFOrderScales(const uint src_order,
    292. 

  292.     const uint dev_order, const bool horizontalOnly) noexcept
    293. 

  293. {
    294. 

  294.     std::array<float,MaxAmbiOrder+1> res{};
    295. 

  295. 

  296.     if(!horizontalOnly)
    297. 

  297.     {
    298. 

  298.         for(size_t i{0};i < MaxAmbiOrder+1;++i)
    299. 

  299.             res[i] = HFScales[src_order][i] / HFScales[dev_order][i];
    300. 

  300.     }
    301. 

  301.     else
    302. 

  302.     {
    303. 

  303.         for(size_t i{0};i < MaxAmbiOrder+1;++i)
    304. 

  304.             res[i] = HFScales2D[src_order][i] / HFScales2D[dev_order][i];
    305. 

  305.     }
    306. 

  306. 

  307.     return res;
    308. 

  308. }
    309.