diff options
Diffstat (limited to 'util/adt_generate_qt.m')
| -rw-r--r-- | util/adt_generate_qt.m | 232 |
1 files changed, 232 insertions, 0 deletions
diff --git a/util/adt_generate_qt.m b/util/adt_generate_qt.m new file mode 100644 index 000000000..45b3dcd03 --- /dev/null +++ b/util/adt_generate_qt.m @@ -0,0 +1,232 @@ +%{ +**************************************************************************** +** +** Copyright (C) 2021 The Qt Company Ltd. +** Contact: https://www.qt.io/licensing/ +** +** This file is the build configuration utility of the Qt Toolkit. +** +** $QT_BEGIN_LICENSE:GPL-EXCEPT$ +** Commercial License Usage +** Licensees holding valid commercial Qt licenses may use this file in +** accordance with the commercial license agreement provided with the +** Software or, alternatively, in accordance with the terms contained in +** a written agreement between you and The Qt Company. For licensing terms +** and conditions see https://www.qt.io/terms-conditions. For further +** information use the contact form at https://www.qt.io/contact-us. +** +** GNU General Public License Usage +** Alternatively, this file may be used under the terms of the GNU +** General Public License version 3 as published by the Free Software +** Foundation with exceptions as appearing in the file LICENSE.GPL3-EXCEPT +** included in the packaging of this file. Please review the following +** information to ensure the GNU General Public License requirements will +** be met: https://www.gnu.org/licenses/gpl-3.0.html. +** +** $QT_END_LICENSE$ +** +**************************************************************************** +%} + +% Create mapping matrixes to convert ambisonic to different speaker layouts supported in Qt Multimedia +% +% This uses the ADT octave framework from https://bitbucket.org/ambidecodertoolbox/adt/src/master/ +% to generate conversion matrixes between ambisonic formats and various speaker configurations +% +% We're generating two band conversion matrices with a cutoff frequency of 380 Herz. Filtering low +% and high frequencies differently is important to get a decent spatial reproduction. For details +% see the "Is my decoder Ambisonic?" paper (https://ambisonics.dreamhosters.com/BLaH3.pdf) +% +function adt_generate_qt() + % assume speakers are 2 meters from listener. This should avoid near-field + % effects and should work for most room setups + radius = 2; + + [outfile,msg] = fopen("qambisonicdecoderdata_p.h",'w'); + fprintf(outfile, "/****************************************************************************\n"); + fprintf(outfile, "**\n"); + fprintf(outfile, "** Copyright (C) 2016 The Qt Company Ltd.\n"); + fprintf(outfile, "** Contact: https://www.qt.io/licensing/\n"); + fprintf(outfile, "**\n"); + fprintf(outfile, "** This file is part of the Qt Toolkit.\n"); + fprintf(outfile, "**\n"); + fprintf(outfile, "** $QT_BEGIN_LICENSE:LGPL$\n"); + fprintf(outfile, "** Commercial License Usage\n"); + fprintf(outfile, "** Licensees holding valid commercial Qt licenses may use this file in\n"); + fprintf(outfile, "** accordance with the commercial license agreement provided with the\n"); + fprintf(outfile, "** Software or, alternatively, in accordance with the terms contained in\n"); + fprintf(outfile, "** a written agreement between you and The Qt Company. For licensing terms\n"); + fprintf(outfile, "** and conditions see https://www.qt.io/terms-conditions. For further\n"); + fprintf(outfile, "** information use the contact form at https://www.qt.io/contact-us.\n"); + fprintf(outfile, "**\n"); + fprintf(outfile, "** GNU Lesser General Public License Usage\n"); + fprintf(outfile, "** Alternatively, this file may be used under the terms of the GNU Lesser\n"); + fprintf(outfile, "** General Public License version 3 as published by the Free Software\n"); + fprintf(outfile, "** Foundation and appearing in the file LICENSE.LGPL3 included in the\n"); + fprintf(outfile, "** packaging of this file. Please review the following information to\n"); + fprintf(outfile, "** ensure the GNU Lesser General Public License version 3 requirements\n"); + fprintf(outfile, "** will be met: https://www.gnu.org/licenses/lgpl-3.0.html.\n"); + fprintf(outfile, "**\n"); + fprintf(outfile, "** GNU General Public License Usage\n"); + fprintf(outfile, "** Alternatively, this file may be used under the terms of the GNU\n"); + fprintf(outfile, "** General Public License version 2.0 or (at your option) the GNU General\n"); + fprintf(outfile, "** Public license version 3 or any later version approved by the KDE Free\n"); + fprintf(outfile, "** Qt Foundation. The licenses are as published by the Free Software\n"); + fprintf(outfile, "** Foundation and appearing in the file LICENSE.GPL2 and LICENSE.GPL3\n"); + fprintf(outfile, "** included in the packaging of this file. Please review the following\n"); + fprintf(outfile, "** information to ensure the GNU General Public License requirements will\n"); + fprintf(outfile, "** be met: https://www.gnu.org/licenses/gpl-2.0.html and\n"); + fprintf(outfile, "** https://www.gnu.org/licenses/gpl-3.0.html.\n"); + fprintf(outfile, "**\n"); + fprintf(outfile, "** $QT_END_LICENSE$\n"); + fprintf(outfile, "**\n"); + fprintf(outfile, "****************************************************************************/\n"); + fprintf(outfile, "#ifndef QAMBISONICDECODERDATA_P_H\n"); + fprintf(outfile, "#define QAMBISONICDECODERDATA_P_H\n\n"); + fprintf(outfile, "#include <qtmultimediaglobal_p.h>\n\n"); + fprintf(outfile, "// W A R N I N G\n"); + fprintf(outfile, "// -------------\n"); + fprintf(outfile, "//\n"); + fprintf(outfile, "// This file is not part of the Qt API. It exists purely as an\n"); + fprintf(outfile, "// implementation detail. This header file may change from version to\n"); + fprintf(outfile, "// version without notice, or even be removed.\n"); + fprintf(outfile, "//\n"); + fprintf(outfile, "// We mean it.\n"); + fprintf(outfile, "//\n\n"); + fprintf(outfile, "// This file is generated by the matlab/octave file adt_generate_qt.m\n"); + fprintf(outfile, "// using the Ambisonic Decoder Toolbox (https://bitbucket.org/ambidecodertoolbox/adt/src/master/)\n"); + fprintf(outfile, "\n\n"); + fprintf(outfile, "QT_BEGIN_NAMESPACE\n\n"); + + % cover top/bottom and back for mono and stereo + imag_speakers = [0,0,radius; 0,0,-radius; 0,radius,0; 0,-radius,0; -radius,0,0]; + + % Mono, one speaker up front + S = ambi_spkr_array(... + ... % array name + 'mono', ... + ... % coordinate codes, unit codes + ... % Azimuth, Elevation, Radius; Degrees, Degrees, Meters + 'AER', 'DDM', ... + ... % speaker name, [azimuth, elevation, radius] + 'C', [ 0, 0, radius] ... + ); + createDecoders(S, imag_speakers, outfile); + + % Stereo, assume -30 and 30 degree speakers + S = ambi_spkr_array(... + 'stereo', ... + 'AER', 'DDM', ... + 'L', [ 30, 0, radius], ... + 'R', [ -30, 0, radius] ... + ); + createDecoders(S, imag_speakers, outfile); + + S.lfeRow = 3; + S.name = "2dot1"; + createDecoders(S, imag_speakers, outfile); + + % cover top/bottom for surround + imag_speakers = [0,0,radius; 0,0,-radius]; + + % https://www.dolby.com/us/en/guide/surround-sound-speaker-setup/5-1-setup.html + % Dolby 5.1 -- F: 22-30, S: 110-120, C in same plane as F + S = ambi_spkr_array(... + '5dot0', ... + 'AER', 'DDM', ... + 'L', [ 30, 0, radius], ... + 'R', [ -30, 0, radius], ... + 'C', [ 0, 0, radius], ... + 'Ls', [ 110, 0, radius], ... + 'Rs', [-110, 0, radius] ... + ); + createDecoders(S, imag_speakers, outfile); + + S.lfeRow = 4; + S.name = "5dot1"; + createDecoders(S, imag_speakers, outfile); + + % https://www.dolby.com/us/en/guide/surround-sound-speaker-setup/7-1-setup.html + % Dolby 7.1 F: 22-30, S: 90-110, B: 135-150 + S = ambi_spkr_array(... + ... % array name + '7dot0', ... + ... % coordinate codes, unit codes + ... % Azimuth, Elevation, Radius; Degrees, Degrees, Meters + 'AER', 'DDM', ... + ... % speaker name, [azimuth, elevation, radius] + 'L', [ 30, 0, radius], ... + 'R', [ -30, 0, radius], ... + 'C', [ 0, 0, radius], ... + 'Ls', [ 90, 0, radius], ... + 'Rs', [ -90, 0, radius], ... + 'Lb', [ 150, 0, radius], ... + 'Rb', [-150, 0, radius] ... + ); + createDecoders(S, imag_speakers, outfile); + + S.lfeRow = 4; + S.name = "7dot1"; + createDecoders(S, imag_speakers, outfile); + + fprintf(outfile, "QT_END_NAMESPACE\n\n"); + fprintf(outfile, "#endif\n\n"); + fclose(outfile); +end + +% remove rounding errors +function m = trimMatrix(m) + for i = 1:rows(m) + for j = 1:columns(m) + if (abs(m(i, j)) < 1e-4) + m(i, j) = 0; + endif + endfor + endfor +end + +function writeLFERow(outfile, m, suffix) + if (strcmp(suffix, 'hf')) + fprintf(outfile, "0.0f, "); % no need for High frequency data on the LFE channel + else + fprintf(outfile, "0.5f, "); + endif + for i = 2:columns(m) + fprintf(outfile, "0.0f, "); + endfor + fprintf(outfile, "// LFE\n"); +end + +function writeMatrix(outfile, level, S, M, suffix) + m = trimMatrix(M); + hasLFE = isfield(S, "lfeRow"); + r = rows(m); + c = columns(m); + if (hasLFE) + r = r + 1; + endif + fprintf(outfile, "// Decoder matrix for %s, ambisonic level %d\n", S.name, level); + fprintf(outfile, "static constexpr float decoderMatrix_%s_%d_%s[%d*%d] = {\n", S.name, level, suffix, r, c); + for i = 1:rows(S.id) + fprintf(outfile, "%ff, ", m(i, :)); + fprintf(outfile, "// %s\n", S.id(i, 1){1}); + if (hasLFE && S.lfeRow == i + 1) + writeLFERow(outfile, m, suffix); + endif + endfor + fprintf(outfile, "};\n\n"); +end + +function createOneDecoder(S, imag_speakers, outfile, level) + [D,S,M,C] = ambi_run_allrad(S, level, imag_speakers, [S.name '_' int2str(level)], false, "amb", 1, 3); + writeMatrix(outfile, level, S, M.lf, "lf"); + m = ambi_apply_gamma(M.hf, D.hf_gains, C); + writeMatrix(outfile, level, S, m, "hf"); +end + +function createDecoders(S, imag_speakers, outfile) + for level = [1:3] + createOneDecoder(S, imag_speakers, outfile, level) + endfor +end + |