From c6ab02a626e15b4a12fc09ecd844eb8b95380c3c Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?=C3=89anna=20=C3=93=20Cath=C3=A1in?= Date: Wed, 7 Apr 2021 14:35:25 +0100 Subject: MLECO-1252 ASR sample application using the public ArmNN C++ API. MIME-Version: 1.0 Content-Type: text/plain; charset=UTF-8 Content-Transfer-Encoding: 8bit Change-Id: I98cd505b8772a8c8fa88308121bc94135bb45068 Signed-off-by: Éanna Ó Catháin --- samples/SpeechRecognition/test/PreprocessTest.cpp | 136 ++++++++++++++++++++++ 1 file changed, 136 insertions(+) create mode 100644 samples/SpeechRecognition/test/PreprocessTest.cpp (limited to 'samples/SpeechRecognition/test/PreprocessTest.cpp') diff --git a/samples/SpeechRecognition/test/PreprocessTest.cpp b/samples/SpeechRecognition/test/PreprocessTest.cpp new file mode 100644 index 0000000000..2b98831fda --- /dev/null +++ b/samples/SpeechRecognition/test/PreprocessTest.cpp @@ -0,0 +1,136 @@ +// +// Copyright © 2020 Arm Ltd and Contributors. All rights reserved. +// SPDX-License-Identifier: MIT +// + +#include +#include + +#include "Preprocess.hpp" +#include "DataStructures.hpp" + +void PopulateTestWavVector(std::vector& vec) +{ + constexpr int int16max = std::numeric_limits::max(); + int val = 0; + for (size_t i = 0; i < vec.size(); ++i, ++val) + { + + /* We want a differential filter response from both - order 1 + * and 2 => Don't have a linear signal here - we use a signal + * using squares for example. Alternate sign flips might work + * just as well and will be computationally less work! */ + int valsq = val * val; + if (valsq > int16max) + { + val = 0; + valsq = 0; + } + vec[i] = valsq; + } +} + +TEST_CASE("Preprocessing calculation INT8") +{ + /*Test Constants: */ + const uint32_t windowLen = 512; + const uint32_t windowStride = 160; + const float quantScale = 0.1410219967365265; + const int quantOffset = -11; + int numMfccVectors = 10; + const int sampFreq = 16000; + const int frameLenMs = 32; + const int frameLenSamples = sampFreq * frameLenMs * 0.001; + const int numMfccFeats = 13; + const int audioDataToPreProcess = 512 + ((numMfccVectors -1) * windowStride); + int outputBufferSize = numMfccVectors * numMfccFeats * 3; + + /* Test wav memory */ + std::vector testWav1((windowStride * numMfccVectors) + + (windowLen - windowStride)); + /* Populate with dummy input */ + PopulateTestWavVector(testWav1); + + MfccParams mfccParams(sampFreq, 128, 0, 8000, numMfccFeats, frameLenSamples, false, numMfccVectors); + + MFCC mfccInst = MFCC(mfccParams); + + std::vector fullAudioData; + + for(int i = 0; i < 4; ++i) + { + for (auto f : testWav1) + { + fullAudioData.emplace_back(static_cast(f) / (1<<15)); + } + } + + Preprocess prep(frameLenSamples, windowStride, mfccInst); + + std::vector outputBuffer(outputBufferSize); + + prep.Invoke(fullAudioData.data(), audioDataToPreProcess, outputBuffer, quantOffset, quantScale); + + int8_t expectedResult[numMfccVectors][numMfccFeats*3] = + { + /* Feature vec 0 */ + -32, 4, -9, -8, -10, -10, -11, -11, -11, -11, -12, -11, -11, /* MFCCs */ + -11, -11, -11, -11, -11, -11, -11, -11, -11, -11, -11, -11, -11, /* Delta 1 */ + -10, -10, -10, -10, -10, -10, -10, -10, -10, -10, -10, -10, -10, /* Delta 2 */ + + /* Feature vec 1 */ + -31, 4, -9, -8, -10, -10, -11, -11, -11, -11, -12, -11, -11, + -11, -11, -11, -11, -11, -11, -11, -11, -11, -11, -11, -11, -11, + -10, -10, -10, -10, -10, -10, -10, -10, -10, -10, -10, -10, -10, + + /* Feature vec 2 */ + -31, 4, -9, -9, -10, -10, -11, -11, -11, -11, -12, -12, -12, + -11, -11, -11, -11, -11, -11, -11, -11, -11, -11, -11, -11, -11, + -10, -10, -10, -10, -10, -10, -10, -10, -10, -10, -10, -10, -10, + + /* Feature vec 3 */ + -31, 4, -9, -9, -10, -10, -11, -11, -11, -11, -11, -12, -12, + -11, -11, -11, -11, -11, -11, -11, -11, -11, -11, -11, -11, -11, + -10, -10, -10, -10, -10, -10, -10, -10, -10, -10, -10, -10, -10, + + /* Feature vec 4 : this should have valid delta 1 and delta 2 */ + -31, 4, -9, -9, -10, -10, -11, -11, -11, -11, -11, -12, -12, + -38, -29, -9, 1, -2, -7, -8, -8, -12, -16, -14, -5, 5, + -68, -50, -13, 5, 0, -9, -9, -8, -13, -20, -19, -3, 15, + + /* Feature vec 5 : this should have valid delta 1 and delta 2 */ + -31, 4, -9, -8, -10, -10, -11, -11, -11, -11, -11, -12, -12, + -62, -45, -11, 5, 0, -8, -9, -8, -12, -19, -17, -3, 13, + -27, -22, -13, -9, -11, -12, -12, -11, -11, -13, -13, -10, -6, + + /* Feature vec 6 */ + -31, 4, -9, -8, -10, -10, -11, -11, -11, -11, -12, -11, -11, + -11, -11, -11, -11, -11, -11, -11, -11, -11, -11, -11, -11, -11, + -10, -10, -10, -10, -10, -10, -10, -10, -10, -10, -10, -10, -10, + + /* Feature vec 7 */ + -32, 4, -9, -8, -10, -10, -11, -11, -11, -12, -12, -11, -11, + -11, -11, -11, -11, -11, -11, -11, -11, -11, -11, -11, -11, -11, + -10, -10, -10, -10, -10, -10, -10, -10, -10, -10, -10, -10, -10, + + /* Feature vec 8 */ + -32, 4, -9, -8, -10, -10, -11, -11, -11, -12, -12, -11, -11, + -11, -11, -11, -11, -11, -11, -11, -11, -11, -11, -11, -11, -11, + -10, -10, -10, -10, -10, -10, -10, -10, -10, -10, -10, -10, -10, + + /* Feature vec 9 */ + -31, 4, -9, -8, -10, -10, -11, -11, -11, -11, -12, -11, -11, + -11, -11, -11, -11, -11, -11, -11, -11, -11, -11, -11, -11, -11, + -10, -10, -10, -10, -10, -10, -10, -10, -10, -10, -10, -10, -10 + }; + + /* Check that the elements have been calculated correctly */ + for (uint32_t j = 0; j < numMfccVectors; ++j) + { + for (uint32_t i = 0; i < numMfccFeats * 3; ++i) + { + size_t tensorIdx = (j * numMfccFeats * 3) + i; + CHECK(static_cast(outputBuffer.at(tensorIdx) == static_cast(expectedResult[j][i]))); + } + } +} -- cgit v1.2.1