| 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 |  |