Mercurial > hg > beaglert
diff examples/04-Audio/filter-FIR/render.cpp @ 464:8fcfbfb32aa0 prerelease
Examples reorder with subdirectories. Added header to each project. Moved Doxygen to bottom of render.cpp.
author | Robert Jack <robert.h.jack@gmail.com> |
---|---|
date | Mon, 20 Jun 2016 16:20:38 +0100 |
parents | |
children | b935f890e512 |
line wrap: on
line diff
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/examples/04-Audio/filter-FIR/render.cpp Mon Jun 20 16:20:38 2016 +0100 @@ -0,0 +1,172 @@ +/* + ____ _____ _ _ +| __ )| ____| | / \ +| _ \| _| | | / _ \ +| |_) | |___| |___ / ___ \ +|____/|_____|_____/_/ \_\ + +The platform for ultra-low latency audio and sensor processing + +http://bela.io + +A project of the Augmented Instruments Laboratory within the +Centre for Digital Music at Queen Mary University of London. +http://www.eecs.qmul.ac.uk/~andrewm + +(c) 2016 Augmented Instruments Laboratory: Andrew McPherson, + Astrid Bin, Liam Donovan, Christian Heinrichs, Robert Jack, + Giulio Moro, Laurel Pardue, Victor Zappi. All rights reserved. + +The Bela software is distributed under the GNU Lesser General Public License +(LGPL 3.0), available here: https://www.gnu.org/licenses/lgpl-3.0.txt +*/ + + +#define ENABLE_NE10_FIR_FLOAT_NEON // Define needed for Ne10 library + +#include <Bela.h> +#include <cmath> +#include <ne10/NE10.h> // neon library +#include "SampleData.h" +#include "FIRfilter.h" + +SampleData gSampleData; // User defined structure to get complex data from main +int gReadPtr; // Position of last read sample from file + +// filter vars +ne10_fir_instance_f32_t gFIRfilter; +ne10_float32_t *gFIRfilterIn; +ne10_float32_t *gFIRfilterOut; +ne10_uint32_t blockSize; +ne10_float32_t *gFIRfilterState; + +void initialise_filter(BelaContext *context); + +// Task for handling the update of the frequencies using the matrix +AuxiliaryTask gTriggerSamplesTask; + +bool initialise_trigger(); +void trigger_samples(); + +bool setup(BelaContext *context, void *userData) +{ + + // Retrieve a parameter passed in from the initAudio() call + gSampleData = *(SampleData *)userData; + + gReadPtr = -1; + + initialise_filter(context); + + // Initialise auxiliary tasks + if(!initialise_trigger()) + return false; + + return true; +} + +void render(BelaContext *context, void *userData) +{ + for(unsigned int n = 0; n < context->audioFrames; n++) { + float in = 0; + + // If triggered... + if(gReadPtr != -1) + in += gSampleData.samples[gReadPtr++]; // ...read each sample... + + if(gReadPtr >= gSampleData.sampleLen) + gReadPtr = -1; + + gFIRfilterIn[n] = in; + } + + ne10_fir_float_neon(&gFIRfilter, gFIRfilterIn, gFIRfilterOut, blockSize); + + for(unsigned int n = 0; n < context->audioFrames; n++) { + for(unsigned int channel = 0; channel < context->audioChannels; channel++) + context->audioOut[n * context->audioChannels + channel] = gFIRfilterOut[n]; // ...and put it in both left and right channel + } + + + // Request that the lower-priority task run at next opportunity + Bela_scheduleAuxiliaryTask(gTriggerSamplesTask); +} + +// Initialise NE10 data structures to define FIR filter + +void initialise_filter(BelaContext *context) +{ + blockSize = context->audioFrames; + gFIRfilterState = (ne10_float32_t *) NE10_MALLOC ((FILTER_TAP_NUM+blockSize-1) * sizeof (ne10_float32_t)); + gFIRfilterIn = (ne10_float32_t *) NE10_MALLOC (blockSize * sizeof (ne10_float32_t)); + gFIRfilterOut = (ne10_float32_t *) NE10_MALLOC (blockSize * sizeof (ne10_float32_t)); + ne10_fir_init_float(&gFIRfilter, FILTER_TAP_NUM, filterTaps, gFIRfilterState, blockSize); +} + + +// Initialise the auxiliary task +// and print info + +bool initialise_trigger() +{ + if((gTriggerSamplesTask = Bela_createAuxiliaryTask(&trigger_samples, 50, "bela-trigger-samples")) == 0) + return false; + + rt_printf("Press 'a' to trigger sample, 's' to stop\n"); + rt_printf("Press 'q' to quit\n"); + + return true; +} + +// This is a lower-priority call to periodically read keyboard input +// and trigger samples. By placing it at a lower priority, +// it has minimal effect on the audio performance but it will take longer to +// complete if the system is under heavy audio load. + +void trigger_samples() +{ + // This is not a real-time task! + // Cos getchar is a system call, not handled by Xenomai. + // This task will be automatically down graded. + + char keyStroke = '.'; + + keyStroke = getchar(); + while(getchar()!='\n'); // to read the first stroke + + switch (keyStroke) + { + case 'a': + gReadPtr = 0; + break; + case 's': + gReadPtr = -1; + break; + case 'q': + gShouldStop = true; + break; + default: + break; + } +} + + +void cleanup(BelaContext *context, void *userData) +{ + delete[] gSampleData.samples; + + NE10_FREE(gFIRfilterState); + NE10_FREE(gFIRfilterIn); + NE10_FREE(gFIRfilterOut); +} + +/* ------------ Project Explantation ------------ */ + +/** +\example 04-filter-FIR + +Finite Impulse Response Filter +------------------------------ + +This is an example of a finite impulse response filter implementation. +*/