| 1 | /**************************************************************************** |
| 2 | ** |
| 3 | ** Copyright (C) 2017 The Qt Company Ltd. |
| 4 | ** Contact: https://www.qt.io/licensing/ |
| 5 | ** |
| 6 | ** This file is part of the examples of the Qt Toolkit. |
| 7 | ** |
| 8 | ** $QT_BEGIN_LICENSE:BSD$ |
| 9 | ** Commercial License Usage |
| 10 | ** Licensees holding valid commercial Qt licenses may use this file in |
| 11 | ** accordance with the commercial license agreement provided with the |
| 12 | ** Software or, alternatively, in accordance with the terms contained in |
| 13 | ** a written agreement between you and The Qt Company. For licensing terms |
| 14 | ** and conditions see https://www.qt.io/terms-conditions. For further |
| 15 | ** information use the contact form at https://www.qt.io/contact-us. |
| 16 | ** |
| 17 | ** BSD License Usage |
| 18 | ** Alternatively, you may use this file under the terms of the BSD license |
| 19 | ** as follows: |
| 20 | ** |
| 21 | ** "Redistribution and use in source and binary forms, with or without |
| 22 | ** modification, are permitted provided that the following conditions are |
| 23 | ** met: |
| 24 | ** * Redistributions of source code must retain the above copyright |
| 25 | ** notice, this list of conditions and the following disclaimer. |
| 26 | ** * Redistributions in binary form must reproduce the above copyright |
| 27 | ** notice, this list of conditions and the following disclaimer in |
| 28 | ** the documentation and/or other materials provided with the |
| 29 | ** distribution. |
| 30 | ** * Neither the name of The Qt Company Ltd nor the names of its |
| 31 | ** contributors may be used to endorse or promote products derived |
| 32 | ** from this software without specific prior written permission. |
| 33 | ** |
| 34 | ** |
| 35 | ** THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS |
| 36 | ** "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT |
| 37 | ** LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR |
| 38 | ** A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT |
| 39 | ** OWNER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, |
| 40 | ** SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT |
| 41 | ** LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, |
| 42 | ** DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY |
| 43 | ** THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT |
| 44 | ** (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE |
| 45 | ** OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE." |
| 46 | ** |
| 47 | ** $QT_END_LICENSE$ |
| 48 | ** |
| 49 | ****************************************************************************/ |
| 50 | |
| 51 | #include "spectrum.h" |
| 52 | #include "utils.h" |
| 53 | #include <QByteArray> |
| 54 | #include <QAudioFormat> |
| 55 | #include <qmath.h> |
| 56 | #include <qendian.h> |
| 57 | |
| 58 | void generateTone(const SweptTone &tone, const QAudioFormat &format, QByteArray &buffer) |
| 59 | { |
| 60 | Q_ASSERT(isPCMS16LE(format)); |
| 61 | |
| 62 | const int channelBytes = format.sampleSize() / 8; |
| 63 | const int sampleBytes = format.channelCount() * channelBytes; |
| 64 | int length = buffer.size(); |
| 65 | const int numSamples = buffer.size() / sampleBytes; |
| 66 | |
| 67 | Q_ASSERT(length % sampleBytes == 0); |
| 68 | Q_UNUSED(sampleBytes) // suppress warning in release builds |
| 69 | |
| 70 | unsigned char *ptr = reinterpret_cast<unsigned char *>(buffer.data()); |
| 71 | |
| 72 | qreal phase = 0.0; |
| 73 | |
| 74 | const qreal d = 2 * M_PI / format.sampleRate(); |
| 75 | |
| 76 | // We can't generate a zero-frequency sine wave |
| 77 | const qreal startFreq = tone.startFreq ? tone.startFreq : 1.0; |
| 78 | |
| 79 | // Amount by which phase increases on each sample |
| 80 | qreal phaseStep = d * startFreq; |
| 81 | |
| 82 | // Amount by which phaseStep increases on each sample |
| 83 | // If this is non-zero, the output is a frequency-swept tone |
| 84 | const qreal phaseStepStep = d * (tone.endFreq - startFreq) / numSamples; |
| 85 | |
| 86 | while (length) { |
| 87 | const qreal x = tone.amplitude * qSin(v: phase); |
| 88 | const qint16 value = realToPcm(real: x); |
| 89 | for (int i=0; i<format.channelCount(); ++i) { |
| 90 | qToLittleEndian<qint16>(src: value, dest: ptr); |
| 91 | ptr += channelBytes; |
| 92 | length -= channelBytes; |
| 93 | } |
| 94 | |
| 95 | phase += phaseStep; |
| 96 | while (phase > 2 * M_PI) |
| 97 | phase -= 2 * M_PI; |
| 98 | phaseStep += phaseStepStep; |
| 99 | } |
| 100 | } |
| 101 | |