From a4dd0ad63c00f4dee3b86dfd3075d1d61b2b3180 Mon Sep 17 00:00:00 2001 From: sanine Date: Sat, 27 Aug 2022 23:52:56 -0500 Subject: add plibsys --- portaudio/examples/paex_ocean_shore.c | 533 ---------------------------------- 1 file changed, 533 deletions(-) delete mode 100644 portaudio/examples/paex_ocean_shore.c (limited to 'portaudio/examples/paex_ocean_shore.c') diff --git a/portaudio/examples/paex_ocean_shore.c b/portaudio/examples/paex_ocean_shore.c deleted file mode 100644 index 9424e8b..0000000 --- a/portaudio/examples/paex_ocean_shore.c +++ /dev/null @@ -1,533 +0,0 @@ -/** @file paex_ocean_shore.c - @ingroup examples_src - @brief Generate Pink Noise using Gardner method, and make "waves". Provides an example of how to - post stuff to/from the audio callback using lock-free FIFOs implemented by the PA ringbuffer. - - Optimization suggested by James McCartney uses a tree - to select which random value to replace. -
-    x x x x x x x x x x x x x x x x
-    x   x   x   x   x   x   x   x
-    x       x       x       x
-     x               x
-       x
-
- Tree is generated by counting trailing zeros in an increasing index. - When the index is zero, no random number is selected. - - @author Phil Burk http://www.softsynth.com - Robert Bielik -*/ -/* - * $Id$ - * - * This program uses the PortAudio Portable Audio Library. - * For more information see: http://www.portaudio.com - * Copyright (c) 1999-2000 Ross Bencina and Phil Burk - * - * Permission is hereby granted, free of charge, to any person obtaining - * a copy of this software and associated documentation files - * (the "Software"), to deal in the Software without restriction, - * including without limitation the rights to use, copy, modify, merge, - * publish, distribute, sublicense, and/or sell copies of the Software, - * and to permit persons to whom the Software is furnished to do so, - * subject to the following conditions: - * - * The above copyright notice and this permission notice shall be - * included in all copies or substantial portions of the Software. - * - * THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, - * EXPRESS OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF - * MERCHANTABILITY, FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. - * IN NO EVENT SHALL THE AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR - * ANY CLAIM, DAMAGES OR OTHER LIABILITY, WHETHER IN AN ACTION OF - * CONTRACT, TORT OR OTHERWISE, ARISING FROM, OUT OF OR IN CONNECTION - * WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE SOFTWARE. - */ - -/* - * The text above constitutes the entire PortAudio license; however, - * the PortAudio community also makes the following non-binding requests: - * - * Any person wishing to distribute modifications to the Software is - * requested to send the modifications to the original developer so that - * they can be incorporated into the canonical version. It is also - * requested that these non-binding requests be included along with the - * license above. - */ - -#include -#include -#include -#include -#include - -#include "portaudio.h" -#include "pa_ringbuffer.h" -#include "pa_util.h" - -#define PINK_MAX_RANDOM_ROWS (30) -#define PINK_RANDOM_BITS (24) -#define PINK_RANDOM_SHIFT ((sizeof(long)*8)-PINK_RANDOM_BITS) - -typedef struct -{ - long pink_Rows[PINK_MAX_RANDOM_ROWS]; - long pink_RunningSum; /* Used to optimize summing of generators. */ - int pink_Index; /* Incremented each sample. */ - int pink_IndexMask; /* Index wrapped by ANDing with this mask. */ - float pink_Scalar; /* Used to scale within range of -1.0 to +1.0 */ -} -PinkNoise; - -typedef struct -{ - float bq_b0; - float bq_b1; - float bq_b2; - float bq_a1; - float bq_a2; -} BiQuad; - -typedef enum -{ - State_kAttack, - State_kPreDecay, - State_kDecay, - State_kCnt, -} EnvState; - -typedef struct -{ - PinkNoise wave_left; - PinkNoise wave_right; - - BiQuad wave_bq_coeffs; - float wave_bq_left[2]; - float wave_bq_right[2]; - - EnvState wave_envelope_state; - float wave_envelope_level; - float wave_envelope_max_level; - float wave_pan_left; - float wave_pan_right; - float wave_attack_incr; - float wave_decay_incr; - -} OceanWave; - -/* Prototypes */ -static unsigned long GenerateRandomNumber( void ); -void InitializePinkNoise( PinkNoise *pink, int numRows ); -float GeneratePinkNoise( PinkNoise *pink ); -unsigned GenerateWave( OceanWave* wave, float* output, unsigned noOfFrames); - -/************************************************************/ -/* Calculate pseudo-random 32 bit number based on linear congruential method. */ -static unsigned long GenerateRandomNumber( void ) -{ - /* Change this seed for different random sequences. */ - static unsigned long randSeed = 22222; - randSeed = (randSeed * 196314165) + 907633515; - return randSeed; -} - -/************************************************************/ -/* Setup PinkNoise structure for N rows of generators. */ -void InitializePinkNoise( PinkNoise *pink, int numRows ) -{ - int i; - long pmax; - pink->pink_Index = 0; - pink->pink_IndexMask = (1<pink_Scalar = 1.0f / pmax; - /* Initialize rows. */ - for( i=0; ipink_Rows[i] = 0; - pink->pink_RunningSum = 0; -} - -/* Generate Pink noise values between -1.0 and +1.0 */ -float GeneratePinkNoise( PinkNoise *pink ) -{ - long newRandom; - long sum; - float output; - /* Increment and mask index. */ - pink->pink_Index = (pink->pink_Index + 1) & pink->pink_IndexMask; - /* If index is zero, don't update any random values. */ - if( pink->pink_Index != 0 ) - { - /* Determine how many trailing zeros in PinkIndex. */ - /* This algorithm will hang if n==0 so test first. */ - int numZeros = 0; - int n = pink->pink_Index; - while( (n & 1) == 0 ) - { - n = n >> 1; - numZeros++; - } - /* Replace the indexed ROWS random value. - * Subtract and add back to RunningSum instead of adding all the random - * values together. Only one changes each time. - */ - pink->pink_RunningSum -= pink->pink_Rows[numZeros]; - newRandom = ((long)GenerateRandomNumber()) >> PINK_RANDOM_SHIFT; - pink->pink_RunningSum += newRandom; - pink->pink_Rows[numZeros] = newRandom; - } - - /* Add extra white noise value. */ - newRandom = ((long)GenerateRandomNumber()) >> PINK_RANDOM_SHIFT; - sum = pink->pink_RunningSum + newRandom; - /* Scale to range of -1.0 to 0.9999. */ - output = pink->pink_Scalar * sum; - return output; -} - -float ProcessBiquad(const BiQuad* coeffs, float* memory, float input) -{ - float w = input - coeffs->bq_a1 * memory[0] - coeffs->bq_a2 * memory[1]; - float out = coeffs->bq_b1 * memory[0] + coeffs->bq_b2 * memory[1] + coeffs->bq_b0 * w; - memory[1] = memory[0]; - memory[0] = w; - return out; -} - -static const float one_over_2Q_LP = 0.3f; -static const float one_over_2Q_HP = 1.0f; - -unsigned GenerateWave( OceanWave* wave, float* output, unsigned noOfFrames ) -{ - unsigned retval=0,i; - float targetLevel, levelIncr, currentLevel; - switch (wave->wave_envelope_state) - { - case State_kAttack: - targetLevel = noOfFrames * wave->wave_attack_incr + wave->wave_envelope_level; - if (targetLevel >= wave->wave_envelope_max_level) - { - /* Go to decay state */ - wave->wave_envelope_state = State_kPreDecay; - targetLevel = wave->wave_envelope_max_level; - } - /* Calculate lowpass biquad coeffs - - alpha = sin(w0)/(2*Q) - - b0 = (1 - cos(w0))/2 - b1 = 1 - cos(w0) - b2 = (1 - cos(w0))/2 - a0 = 1 + alpha - a1 = -2*cos(w0) - a2 = 1 - alpha - - w0 = [0 - pi[ - */ - { - const float w0 = 3.141592654f * targetLevel / wave->wave_envelope_max_level; - const float alpha = sinf(w0) * one_over_2Q_LP; - const float cosw0 = cosf(w0); - const float a0_fact = 1.0f / (1.0f + alpha); - wave->wave_bq_coeffs.bq_b1 = (1.0f - cosw0) * a0_fact; - wave->wave_bq_coeffs.bq_b0 = wave->wave_bq_coeffs.bq_b1 * 0.5f; - wave->wave_bq_coeffs.bq_b2 = wave->wave_bq_coeffs.bq_b0; - wave->wave_bq_coeffs.bq_a2 = (1.0f - alpha) * a0_fact; - wave->wave_bq_coeffs.bq_a1 = -2.0f * cosw0 * a0_fact; - } - break; - - case State_kPreDecay: - /* Reset biquad state */ - memset(wave->wave_bq_left, 0, 2 * sizeof(float)); - memset(wave->wave_bq_right, 0, 2 * sizeof(float)); - wave->wave_envelope_state = State_kDecay; - - /* Deliberate fall-through */ - - case State_kDecay: - targetLevel = noOfFrames * wave->wave_decay_incr + wave->wave_envelope_level; - if (targetLevel < 0.001f) - { - /* < -60 dB, we're done */ - wave->wave_envelope_state = 3; - retval = 1; - } - /* Calculate highpass biquad coeffs - - alpha = sin(w0)/(2*Q) - - b0 = (1 + cos(w0))/2 - b1 = -(1 + cos(w0)) - b2 = (1 + cos(w0))/2 - a0 = 1 + alpha - a1 = -2*cos(w0) - a2 = 1 - alpha - - w0 = [0 - pi/2[ - */ - { - const float v = targetLevel / wave->wave_envelope_max_level; - const float w0 = 1.5707963f * (1.0f - (v*v)); - const float alpha = sinf(w0) * one_over_2Q_HP; - const float cosw0 = cosf(w0); - const float a0_fact = 1.0f / (1.0f + alpha); - wave->wave_bq_coeffs.bq_b1 = (float)(- (1 + cosw0) * a0_fact); - wave->wave_bq_coeffs.bq_b0 = -wave->wave_bq_coeffs.bq_b1 * 0.5f; - wave->wave_bq_coeffs.bq_b2 = wave->wave_bq_coeffs.bq_b0; - wave->wave_bq_coeffs.bq_a2 = (float)((1.0 - alpha) * a0_fact); - wave->wave_bq_coeffs.bq_a1 = (float)(-2.0 * cosw0 * a0_fact); - } - break; - - default: - break; - } - - currentLevel = wave->wave_envelope_level; - wave->wave_envelope_level = targetLevel; - levelIncr = (targetLevel - currentLevel) / noOfFrames; - - for (i = 0; i < noOfFrames; ++i, currentLevel += levelIncr) - { - (*output++) += ProcessBiquad(&wave->wave_bq_coeffs, wave->wave_bq_left, (GeneratePinkNoise(&wave->wave_left))) * currentLevel * wave->wave_pan_left; - (*output++) += ProcessBiquad(&wave->wave_bq_coeffs, wave->wave_bq_right, (GeneratePinkNoise(&wave->wave_right))) * currentLevel * wave->wave_pan_right; - } - - return retval; -} - - -/*******************************************************************/ - -/* Context for callback routine. */ -typedef struct -{ - OceanWave* waves[16]; /* Maximum 16 waves */ - unsigned noOfActiveWaves; - - /* Ring buffer (FIFO) for "communicating" towards audio callback */ - PaUtilRingBuffer rBufToRT; - void* rBufToRTData; - - /* Ring buffer (FIFO) for "communicating" from audio callback */ - PaUtilRingBuffer rBufFromRT; - void* rBufFromRTData; -} -paTestData; - -/* This routine will be called by the PortAudio engine when audio is needed. -** It may called at interrupt level on some machines so don't do anything -** that could mess up the system like calling malloc() or free(). -*/ -static int patestCallback(const void* inputBuffer, - void* outputBuffer, - unsigned long framesPerBuffer, - const PaStreamCallbackTimeInfo* timeInfo, - PaStreamCallbackFlags statusFlags, - void* userData) -{ - int i; - paTestData *data = (paTestData*)userData; - float *out = (float*)outputBuffer; - (void) inputBuffer; /* Prevent "unused variable" warnings. */ - - /* Reset output data first */ - memset(out, 0, framesPerBuffer * 2 * sizeof(float)); - - for (i = 0; i < 16; ++i) - { - /* Consume the input queue */ - if (data->waves[i] == 0 && PaUtil_GetRingBufferReadAvailable(&data->rBufToRT)) - { - OceanWave* ptr = 0; - PaUtil_ReadRingBuffer(&data->rBufToRT, &ptr, 1); - data->waves[i] = ptr; - } - - if (data->waves[i] != 0) - { - if (GenerateWave(data->waves[i], out, framesPerBuffer)) - { - /* If wave is "done", post it back to the main thread for deletion */ - PaUtil_WriteRingBuffer(&data->rBufFromRT, &data->waves[i], 1); - data->waves[i] = 0; - } - } - } - return paContinue; -} - -#define NEW_ROW_SIZE (12 + (8*rand())/RAND_MAX) - -OceanWave* InitializeWave(double SR, float attackInSeconds, float maxLevel, float positionLeftRight) -{ - OceanWave* wave = NULL; - static unsigned lastNoOfRows = 12; - unsigned newNoOfRows; - - wave = (OceanWave*)PaUtil_AllocateMemory(sizeof(OceanWave)); - if (wave != NULL) - { - InitializePinkNoise(&wave->wave_left, lastNoOfRows); - while ((newNoOfRows = NEW_ROW_SIZE) == lastNoOfRows); - InitializePinkNoise(&wave->wave_right, newNoOfRows); - lastNoOfRows = newNoOfRows; - - wave->wave_envelope_state = State_kAttack; - wave->wave_envelope_level = 0.f; - wave->wave_envelope_max_level = maxLevel; - wave->wave_attack_incr = wave->wave_envelope_max_level / (attackInSeconds * (float)SR); - wave->wave_decay_incr = - wave->wave_envelope_max_level / (attackInSeconds * 4 * (float)SR); - - wave->wave_pan_left = sqrtf(1.0f - positionLeftRight); - wave->wave_pan_right = sqrtf(positionLeftRight); - } - return wave; -} - -static float GenerateFloatRandom(float minValue, float maxValue) -{ - return minValue + ((maxValue - minValue) * rand()) / RAND_MAX; -} - -/*******************************************************************/ -int main(void); -int main(void) -{ - PaStream* stream; - PaError err; - paTestData data = {0}; - PaStreamParameters outputParameters; - double tstamp; - double tstart; - double tdelta = 0; - static const double SR = 44100.0; - static const int FPB = 128; /* Frames per buffer: 2.9 ms buffers. */ - - /* Initialize communication buffers (queues) */ - data.rBufToRTData = PaUtil_AllocateMemory(sizeof(OceanWave*) * 256); - if (data.rBufToRTData == NULL) - { - return 1; - } - PaUtil_InitializeRingBuffer(&data.rBufToRT, sizeof(OceanWave*), 256, data.rBufToRTData); - - data.rBufFromRTData = PaUtil_AllocateMemory(sizeof(OceanWave*) * 256); - if (data.rBufFromRTData == NULL) - { - return 1; - } - PaUtil_InitializeRingBuffer(&data.rBufFromRT, sizeof(OceanWave*), 256, data.rBufFromRTData); - - err = Pa_Initialize(); - if( err != paNoError ) goto error; - - /* Open a stereo PortAudio stream so we can hear the result. */ - outputParameters.device = Pa_GetDefaultOutputDevice(); /* Take the default output device. */ - if (outputParameters.device == paNoDevice) { - fprintf(stderr,"Error: No default output device.\n"); - goto error; - } - outputParameters.channelCount = 2; /* Stereo output, most likely supported. */ - outputParameters.hostApiSpecificStreamInfo = NULL; - outputParameters.sampleFormat = paFloat32; /* 32 bit floating point output. */ - outputParameters.suggestedLatency = Pa_GetDeviceInfo(outputParameters.device)->defaultLowOutputLatency; - err = Pa_OpenStream(&stream, - NULL, /* No input. */ - &outputParameters, - SR, /* Sample rate. */ - FPB, /* Frames per buffer. */ - paDitherOff, /* Clip but don't dither */ - patestCallback, - &data); - if( err != paNoError ) goto error; - - err = Pa_StartStream( stream ); - if( err != paNoError ) goto error; - - printf("Stereo \"ocean waves\" for one minute...\n"); - - tstart = PaUtil_GetTime(); - tstamp = tstart; - srand( (unsigned)time(NULL) ); - - while( ( err = Pa_IsStreamActive( stream ) ) == 1 ) - { - const double tcurrent = PaUtil_GetTime(); - - /* Delete "waves" that the callback is finished with */ - while (PaUtil_GetRingBufferReadAvailable(&data.rBufFromRT) > 0) - { - OceanWave* ptr = 0; - PaUtil_ReadRingBuffer(&data.rBufFromRT, &ptr, 1); - if (ptr != 0) - { - printf("Wave is deleted...\n"); - PaUtil_FreeMemory(ptr); - --data.noOfActiveWaves; - } - } - - if (tcurrent - tstart < 60.0) /* Only start new "waves" during one minute */ - { - if (tcurrent >= tstamp) - { - double tdelta = GenerateFloatRandom(1.0f, 4.0f); - tstamp += tdelta; - - if (data.noOfActiveWaves<16) - { - const float attackTime = GenerateFloatRandom(2.0f, 6.0f); - const float level = GenerateFloatRandom(0.1f, 1.0f); - const float pos = GenerateFloatRandom(0.0f, 1.0f); - OceanWave* p = InitializeWave(SR, attackTime, level, pos); - if (p != NULL) - { - /* Post wave to audio callback */ - PaUtil_WriteRingBuffer(&data.rBufToRT, &p, 1); - ++data.noOfActiveWaves; - - printf("Starting wave at level = %.2f, attack = %.2lf, pos = %.2lf\n", level, attackTime, pos); - } - } - } - } - else - { - if (data.noOfActiveWaves == 0) - { - printf("All waves finished!\n"); - break; - } - } - - Pa_Sleep(100); - } - if( err < 0 ) goto error; - - err = Pa_CloseStream( stream ); - if( err != paNoError ) goto error; - - if (data.rBufToRTData) - { - PaUtil_FreeMemory(data.rBufToRTData); - } - if (data.rBufFromRTData) - { - PaUtil_FreeMemory(data.rBufFromRTData); - } - - Pa_Sleep(1000); - - Pa_Terminate(); - return 0; - -error: - Pa_Terminate(); - fprintf( stderr, "An error occurred while using the portaudio stream\n" ); - fprintf( stderr, "Error number: %d\n", err ); - fprintf( stderr, "Error message: %s\n", Pa_GetErrorText( err ) ); - return 0; -} -- cgit v1.2.1