/*
* maximilian.cpp
* platform independent synthesis library using portaudio or rtaudio
*
* Created by Mick Grierson on 29/12/2009.
* Copyright 2009 Mick Grierson & Strangeloop Limited. All rights reserved.
* Thanks to the Goldsmiths Creative Computing Team.
* Special thanks to Arturo Castro for the PortAudio implementation.
*
* 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.
*
*/
#include "maximilian.h"
#include "math.h"
/* Maximilian can be configured to load ogg vorbis format files using the
* loadOgg() method.
* Uncomment the following to include Sean Barrett's Ogg Vorbis decoder.
* If you're on windows, make sure to add the files std_vorbis.c and std_vorbis.h to your project*/
//#define VORBIS
#ifdef VORBIS
extern "C" {
#include "stb_vorbis.h"
}
#endif
//This used to be important for dealing with multichannel playback
float chandiv= 1;
int maxiSettings::sampleRate = 44100;
int maxiSettings::channels = 2;
int maxiSettings::bufferSize = 1024;
//this is a 514-point sinewave table that has many uses.
double sineBuffer[514]={0,0.012268,0.024536,0.036804,0.049042,0.06131,0.073547,0.085785,0.097992,0.1102,0.12241,0.13455,0.1467,0.15884,0.17093,0.18301,0.19507,0.20709,0.21909,0.23105,0.24295,0.25485,0.26669,0.2785,0.29025,0.30197,0.31366,0.32529,0.33685,0.34839,0.35986,0.37128,0.38266,0.39395,0.40521,0.41641,0.42752,0.4386,0.44958,0.46051,0.47137,0.48215,0.49286,0.50351,0.51407,0.52457,0.53497,0.54529,0.55554,0.5657,0.57578,0.58575,0.59567,0.60547,0.6152,0.62482,0.63437,0.6438,0.65314,0.66238,0.67151,0.68057,0.68951,0.69833,0.70706,0.7157,0.72421,0.7326,0.74091,0.74908,0.75717,0.76514,0.77298,0.7807,0.7883,0.79581,0.80316,0.81042,0.81754,0.82455,0.83142,0.8382,0.84482,0.85132,0.8577,0.86392,0.87006,0.87604,0.88187,0.8876,0.89319,0.89862,0.90396,0.90912,0.91415,0.91907,0.92383,0.92847,0.93295,0.93729,0.9415,0.94556,0.94949,0.95325,0.95691,0.96039,0.96375,0.96692,0.97,0.9729,0.97565,0.97827,0.98074,0.98306,0.98523,0.98724,0.98914,0.99084,0.99243,0.99387,0.99515,0.99628,0.99725,0.99808,0.99875,0.99927,0.99966,0.99988,0.99997,0.99988,0.99966,0.99927,0.99875,0.99808,0.99725,0.99628,0.99515,0.99387,0.99243,0.99084,0.98914,0.98724,0.98523,0.98306,0.98074,0.97827,0.97565,0.9729,0.97,0.96692,0.96375,0.96039,0.95691,0.95325,0.94949,0.94556,0.9415,0.93729,0.93295,0.92847,0.92383,0.91907,0.91415,0.90912,0.90396,0.89862,0.89319,0.8876,0.88187,0.87604,0.87006,0.86392,0.8577,0.85132,0.84482,0.8382,0.83142,0.82455,0.81754,0.81042,0.80316,0.79581,0.7883,0.7807,0.77298,0.76514,0.75717,0.74908,0.74091,0.7326,0.72421,0.7157,0.70706,0.69833,0.68951,0.68057,0.67151,0.66238,0.65314,0.6438,0.63437,0.62482,0.6152,0.60547,0.59567,0.58575,0.57578,0.5657,0.55554,0.54529,0.53497,0.52457,0.51407,0.50351,0.49286,0.48215,0.47137,0.46051,0.44958,0.4386,0.42752,0.41641,0.40521,0.39395,0.38266,0.37128,0.35986,0.34839,0.33685,0.32529,0.31366,0.30197,0.29025,0.2785,0.26669,0.25485,0.24295,0.23105,0.21909,0.20709,0.19507,0.18301,0.17093,0.15884,0.1467,0.13455,0.12241,0.1102,0.097992,0.085785,0.073547,0.06131,0.049042,0.036804,0.024536,0.012268,0,-0.012268,-0.024536,-0.036804,-0.049042,-0.06131,-0.073547,-0.085785,-0.097992,-0.1102,-0.12241,-0.13455,-0.1467,-0.15884,-0.17093,-0.18301,-0.19507,-0.20709,-0.21909,-0.23105,-0.24295,-0.25485,-0.26669,-0.2785,-0.29025,-0.30197,-0.31366,-0.32529,-0.33685,-0.34839,-0.35986,-0.37128,-0.38266,-0.39395,-0.40521,-0.41641,-0.42752,-0.4386,-0.44958,-0.46051,-0.47137,-0.48215,-0.49286,-0.50351,-0.51407,-0.52457,-0.53497,-0.54529,-0.55554,-0.5657,-0.57578,-0.58575,-0.59567,-0.60547,-0.6152,-0.62482,-0.63437,-0.6438,-0.65314,-0.66238,-0.67151,-0.68057,-0.68951,-0.69833,-0.70706,-0.7157,-0.72421,-0.7326,-0.74091,-0.74908,-0.75717,-0.76514,-0.77298,-0.7807,-0.7883,-0.79581,-0.80316,-0.81042,-0.81754,-0.82455,-0.83142,-0.8382,-0.84482,-0.85132,-0.8577,-0.86392,-0.87006,-0.87604,-0.88187,-0.8876,-0.89319,-0.89862,-0.90396,-0.90912,-0.91415,-0.91907,-0.92383,-0.92847,-0.93295,-0.93729,-0.9415,-0.94556,-0.94949,-0.95325,-0.95691,-0.96039,-0.96375,-0.96692,-0.97,-0.9729,-0.97565,-0.97827,-0.98074,-0.98306,-0.98523,-0.98724,-0.98914,-0.99084,-0.99243,-0.99387,-0.99515,-0.99628,-0.99725,-0.99808,-0.99875,-0.99927,-0.99966,-0.99988,-0.99997,-0.99988,-0.99966,-0.99927,-0.99875,-0.99808,-0.99725,-0.99628,-0.99515,-0.99387,-0.99243,-0.99084,-0.98914,-0.98724,-0.98523,-0.98306,-0.98074,-0.97827,-0.97565,-0.9729,-0.97,-0.96692,-0.96375,-0.96039,-0.95691,-0.95325,-0.94949,-0.94556,-0.9415,-0.93729,-0.93295,-0.92847,-0.92383,-0.91907,-0.91415,-0.90912,-0.90396,-0.89862,-0.89319,-0.8876,-0.88187,-0.87604,-0.87006,-0.86392,-0.8577,-0.85132,-0.84482,-0.8382,-0.83142,-0.82455,-0.81754,-0.81042,-0.80316,-0.79581,-0.7883,-0.7807,-0.77298,-0.76514,-0.75717,-0.74908,-0.74091,-0.7326,-0.72421,-0.7157,-0.70706,-0.69833,-0.68951,-0.68057,-0.67151,-0.66238,-0.65314,-0.6438,-0.63437,-0.62482,-0.6152,-0.60547,-0.59567,-0.58575,-0.57578,-0.5657,-0.55554,-0.54529,-0.53497,-0.52457,-0.51407,-0.50351,-0.49286,-0.48215,-0.47137,-0.46051,-0.44958,-0.4386,-0.42752,-0.41641,-0.40521,-0.39395,-0.38266,-0.37128,-0.35986,-0.34839,-0.33685,-0.32529,-0.31366,-0.30197,-0.29025,-0.2785,-0.26669,-0.25485,-0.24295,-0.23105,-0.21909,-0.20709,-0.19507,-0.18301,-0.17093,-0.15884,-0.1467,-0.13455,-0.12241,-0.1102,-0.097992,-0.085785,-0.073547,-0.06131,-0.049042,-0.036804,-0.024536,-0.012268,0,0.012268
};
// This is a transition table that helps with bandlimited oscs.
double transition[1001]={-0.500003,-0.500003,-0.500023,-0.500063,-0.500121,-0.500179,-0.500259,
-0.50036,-0.500476,-0.500591,-0.500732,-0.500893,-0.501066,-0.501239,
-0.50144,-0.501661,-0.501891,-0.502123,-0.502382,-0.502662,-0.502949,
-0.50324,-0.503555,-0.503895,-0.504238,-0.504587,-0.504958,-0.505356,
-0.505754,-0.506162,-0.506589,-0.507042,-0.507495,-0.50796,-0.508444,
-0.508951,-0.509458,-0.509979,-0.510518,-0.511079,-0.511638,-0.512213,
-0.512808,-0.51342,-0.51403,-0.514659,-0.515307,-0.51597,-0.51663,-0.517312,
-0.518012,-0.518724,-0.519433,-0.520166,-0.520916,-0.521675,-0.522432,
-0.523214,-0.524013,-0.524819,-0.525624,-0.526451,-0.527298,-0.528147,
-0.528999,-0.52987,-0.530762,-0.531654,-0.532551,-0.533464,-0.534399,
-0.535332,-0.536271,-0.537226,-0.538202,-0.539172,-0.540152,-0.541148,
-0.542161,-0.543168,-0.544187,-0.54522,-0.546269,-0.54731,-0.548365,
-0.549434,-0.550516,-0.55159,-0.552679,-0.553781,-0.554893,-0.555997,
-0.557118,-0.558252,-0.559391,-0.560524,-0.561674,-0.562836,-0.564001,
-0.565161,-0.566336,-0.567524,-0.568712,-0.569896,-0.571095,-0.572306,
-0.573514,-0.574721,-0.575939,-0.577171,-0.578396,-0.579622,-0.580858,
-0.582108,-0.583348,-0.58459,-0.585842,-0.587106,-0.588358,-0.589614,
-0.590879,-0.592154,-0.593415,-0.594682,-0.595957,-0.59724,-0.598507,
-0.599782,-0.601064,-0.602351,-0.603623,-0.604902,-0.606189,-0.607476,
-0.60875,-0.610032,-0.611319,-0.612605,-0.613877,-0.615157,-0.616443,
-0.617723,-0.618992,-0.620268,-0.621548,-0.62282,-0.624083,-0.62535,
-0.626622,-0.627882,-0.629135,-0.630391,-0.631652,-0.632898,-0.634138,
-0.63538,-0.636626,-0.637854,-0.639078,-0.640304,-0.641531,-0.642739,
-0.643943,-0.645149,-0.646355,-0.647538,-0.64872,-0.649903,-0.651084,
-0.652241,-0.653397,-0.654553,-0.655705,-0.656834,-0.657961,-0.659087,
-0.660206,-0.661304,-0.662399,-0.663492,-0.664575,-0.665639,-0.666699,
-0.667756,-0.6688,-0.669827,-0.670849,-0.671866,-0.672868,-0.673854,
-0.674835,-0.675811,-0.676767,-0.677709,-0.678646,-0.679576,-0.680484,
-0.68138,-0.682269,-0.683151,-0.684008,-0.684854,-0.685693,-0.686524,
-0.687327,-0.688119,-0.688905,-0.689682,-0.690428,-0.691164,-0.691893,
-0.692613,-0.6933,-0.693978,-0.694647,-0.695305,-0.695932,-0.696549,
-0.697156,-0.697748,-0.698313,-0.698865,-0.699407,-0.699932,-0.700431,
-0.700917,-0.701391,-0.701845,-0.702276,-0.702693,-0.703097,-0.703478,
-0.703837,-0.704183,-0.704514,-0.704819,-0.705105,-0.705378,-0.705633,
-0.70586,-0.706069,-0.706265,-0.706444,-0.706591,-0.706721,-0.706837,
-0.706938,-0.707003,-0.707051,-0.707086,-0.707106,-0.707086,-0.707051,
-0.707001,-0.706935,-0.706832,-0.706711,-0.706576,-0.706421,-0.706233,
-0.706025,-0.705802,-0.705557,-0.705282,-0.704984,-0.704671,-0.704334,
-0.703969,-0.703582,-0.703176,-0.702746,-0.702288,-0.70181,-0.701312,
-0.700785,-0.700234,-0.699664,-0.69907,-0.698447,-0.6978,-0.697135,
-0.696446,-0.695725,-0.694981,-0.694219,-0.693435,-0.692613,-0.691771,
-0.690911,-0.69003,-0.689108,-0.688166,-0.687206,-0.686227,-0.685204,
-0.684162,-0.683101,-0.682019,-0.680898,-0.679755,-0.678592,-0.677407,
-0.676187,-0.674941,-0.673676,-0.672386,-0.671066,-0.669718,-0.66835,
-0.666955,-0.665532,-0.664083,-0.662611,-0.661112,-0.659585,-0.658035,
-0.656459,-0.654854,-0.653223,-0.651572,-0.649892,-0.648181,-0.646446,
-0.644691,-0.642909,-0.641093,-0.639253,-0.637393,-0.63551,-0.633588,
-0.631644,-0.62968,-0.627695,-0.625668,-0.623621,-0.621553,-0.619464,
-0.617334,-0.615183,-0.613011,-0.610817,-0.608587,-0.606333,-0.604058,
-0.60176,-0.599429,-0.597072,-0.594695,-0.592293,-0.589862,-0.587404,
-0.584925,-0.58242,-0.579888,-0.577331,-0.574751,-0.572145,-0.569512,
-0.566858,-0.564178,-0.561471,-0.558739,-0.555988,-0.553209,-0.550402,
-0.547572,-0.544723,-0.54185,-0.538944,-0.536018,-0.533072,-0.530105,
-0.527103,-0.524081,-0.52104,-0.51798,-0.514883,-0.511767,-0.508633,
-0.505479,-0.502291,-0.499083,-0.495857,-0.492611,-0.489335,-0.486037,
-0.48272,-0.479384,-0.476021,-0.472634,-0.46923,-0.465805,-0.462356,
-0.458884,-0.455394,-0.451882,-0.448348,-0.444795,-0.44122,-0.437624,
-0.434008,-0.430374,-0.426718,-0.423041,-0.419344,-0.415631,-0.411897,
-0.40814,-0.404365,-0.400575,-0.396766,-0.392933,-0.389082,-0.385217,
-0.381336,-0.377428,-0.373505,-0.369568,-0.365616,-0.361638,-0.357645,
-0.353638,-0.349617,-0.345572,-0.341512,-0.337438,-0.33335,-0.329242,
-0.325118,-0.32098,-0.316829,-0.31266,-0.308474,-0.304276,-0.300063,
-0.295836,-0.291593,-0.287337,-0.283067,-0.278783,-0.274487,-0.270176,
-0.265852,-0.261515,-0.257168,-0.252806,-0.248431,-0.244045,-0.239649,
-0.23524,-0.230817,-0.226385,-0.221943,-0.21749,-0.213024,-0.208548,
-0.204064,-0.199571,-0.195064,-0.190549,-0.186026,-0.181495,-0.176952,
-0.1724,-0.167842,-0.163277,-0.1587,-0.154117,-0.149527,-0.14493,-0.140325,
-0.135712,-0.131094,-0.12647,-0.121839,-0.117201,-0.112559,-0.10791,
-0.103257,-0.0985979,-0.0939343,-0.0892662,-0.0845935,-0.079917,-0.0752362,
-0.0705516,-0.0658635,-0.0611729,-0.0564786,-0.0517814,-0.0470818,-0.0423802,
-0.0376765,-0.0329703,-0.0282629,-0.0235542,-0.0188445,-0.0141335,-0.00942183,
-0.00470983,2.41979e-06,0.00471481,0.00942681,0.0141384,0.0188494,0.023559,
0.028268,0.0329754,0.0376813,0.0423851,0.0470868,0.0517863,0.0564836,
0.0611777,0.0658683,0.0705566,0.0752412,0.0799218,0.0845982,0.0892712,
0.0939393,0.0986028,0.103262,0.107915,0.112563,0.117206,0.121844,0.126475,
0.131099,0.135717,0.14033,0.144935,0.149531,0.154122,0.158705,0.163281,
0.167847,0.172405,0.176956,0.1815,0.18603,0.190553,0.195069,0.199576,
0.204068,0.208552,0.213028,0.217495,0.221947,0.226389,0.230822,0.235245,
0.239653,0.244049,0.248436,0.252811,0.257173,0.26152,0.265857,0.270181,
0.274491,0.278788,0.283071,0.287341,0.291597,0.29584,0.300068,0.30428,
0.308478,0.312664,0.316833,0.320984,0.325122,0.329246,0.333354,0.337442,
0.341516,0.345576,0.34962,0.353642,0.357649,0.361642,0.36562,0.369572,
0.373509,0.377432,0.38134,0.385221,0.389086,0.392936,0.39677,0.400579,
0.404369,0.408143,0.4119,0.415634,0.419347,0.423044,0.426721,0.430377,
0.434011,0.437627,0.441223,0.444798,0.448351,0.451885,0.455397,0.458887,
0.462359,0.465807,0.469232,0.472637,0.476024,0.479386,0.482723,0.486039,
0.489338,0.492613,0.49586,0.499086,0.502294,0.505481,0.508635,0.511769,
0.514885,0.517982,0.521042,0.524083,0.527105,0.530107,0.533074,0.53602,
0.538946,0.541851,0.544725,0.547574,0.550404,0.553211,0.555989,0.55874,
0.561472,0.564179,0.566859,0.569514,0.572146,0.574753,0.577332,0.579889,
0.582421,0.584926,0.587405,0.589863,0.592294,0.594696,0.597073,0.59943,
0.60176,0.604059,0.606333,0.608588,0.610818,0.613012,0.615183,0.617335,
0.619464,0.621553,0.623621,0.625669,0.627696,0.629681,0.631645,0.633588,
0.63551,0.637393,0.639253,0.641093,0.642909,0.644691,0.646446,0.648181,
0.649892,0.651572,0.653223,0.654854,0.656459,0.658035,0.659585,0.661112,
0.662611,0.664083,0.665532,0.666955,0.66835,0.669718,0.671066,0.672386,
0.673676,0.674941,0.676187,0.677407,0.678592,0.679755,0.680898,0.682019,
0.683101,0.684162,0.685204,0.686227,0.687206,0.688166,0.689108,0.69003,
0.690911,0.691771,0.692613,0.693435,0.694219,0.694981,0.695725,0.696447,
0.697135,0.6978,0.698447,0.69907,0.699664,0.700234,0.700786,0.701312,
0.70181,0.702288,0.702746,0.703177,0.703582,0.703969,0.704334,0.704671,
0.704984,0.705282,0.705557,0.705802,0.706025,0.706233,0.706422,0.706576,
0.706712,0.706832,0.706936,0.707002,0.707051,0.707086,0.707106,0.707086,
0.707051,0.707003,0.706939,0.706838,0.706721,0.706592,0.706445,0.706265,
0.70607,0.705861,0.705634,0.705378,0.705105,0.70482,0.704515,0.704184,
0.703837,0.703478,0.703097,0.702694,0.702276,0.701846,0.701392,0.700917,
0.700432,0.699932,0.699408,0.698866,0.698314,0.697749,0.697156,0.696549,
0.695933,0.695305,0.694648,0.693979,0.693301,0.692613,0.691894,0.691165,
0.690428,0.689683,0.688905,0.68812,0.687327,0.686525,0.685693,0.684854,
0.684009,0.683152,0.68227,0.68138,0.680485,0.679577,0.678647,0.67771,
0.676768,0.675811,0.674836,0.673855,0.672869,0.671867,0.670849,0.669827,
0.668801,0.667757,0.6667,0.66564,0.664576,0.663493,0.6624,0.661305,
0.660207,0.659088,0.657962,0.656834,0.655705,0.654553,0.653398,0.652241,
0.651085,0.649903,0.648721,0.647539,0.646356,0.645149,0.643944,0.642739,
0.641532,0.640304,0.639079,0.637855,0.636626,0.635381,0.634139,0.632899,
0.631652,0.630392,0.629136,0.627883,0.626622,0.62535,0.624083,0.62282,
0.621548,0.620268,0.618993,0.617724,0.616443,0.615158,0.613878,0.612605,
0.61132,0.610032,0.608751,0.607477,0.606189,0.604903,0.603623,0.602351,
0.601065,0.599782,0.598508,0.59724,0.595957,0.594682,0.593415,0.592154,
0.59088,0.589615,0.588359,0.587106,0.585843,0.584591,0.583349,0.582108,
0.580859,0.579623,0.578397,0.577172,0.575939,0.574721,0.573515,0.572307,
0.571095,0.569897,0.568713,0.567525,0.566337,0.565161,0.564002,0.562837,
0.561674,0.560525,0.559392,0.558252,0.557119,0.555998,0.554893,0.553782,
0.552679,0.55159,0.550516,0.549434,0.548365,0.54731,0.546269,0.54522,
0.544187,0.543168,0.542161,0.541148,0.540153,0.539173,0.538202,0.537226,
0.536271,0.535332,0.5344,0.533464,0.532551,0.531654,0.530762,0.52987,
0.528999,0.528147,0.527298,0.526451,0.525624,0.524819,0.524014,0.523214,
0.522432,0.521675,0.520916,0.520166,0.519433,0.518724,0.518012,0.517312,
0.51663,0.51597,0.515307,0.51466,0.51403,0.51342,0.512808,0.512213,
0.511638,0.511079,0.510518,0.509979,0.509458,0.508951,0.508444,0.50796,
0.507495,0.507042,0.506589,0.506162,0.505754,0.505356,0.504958,0.504587,
0.504237,0.503895,0.503555,0.50324,0.502949,0.502662,0.502382,0.502123,
0.501891,0.501661,0.50144,0.501239,0.501066,0.500893,0.500732,0.500591,
0.500476,0.50036,0.500259,0.500179,0.500121,0.500063,0.500023,0.500003,0.500003};
//This is a lookup table for converting midi to frequency
double mtofarray[129]={0, 8.661957, 9.177024, 9.722718, 10.3, 10.913383, 11.562325, 12.25, 12.978271, 13.75, 14.567617, 15.433853, 16.351599, 17.323914, 18.354048, 19.445436, 20.601723, 21.826765, 23.124651, 24.5, 25.956543, 27.5, 29.135235, 30.867706, 32.703197, 34.647827, 36.708096, 38.890873, 41.203445, 43.65353, 46.249302, 49., 51.913086, 55., 58.27047, 61.735413, 65.406395, 69.295654, 73.416191, 77.781746, 82.406891, 87.30706, 92.498604, 97.998856, 103.826172, 110., 116.540939, 123.470825, 130.81279, 138.591309, 146.832382, 155.563492, 164.813782, 174.61412, 184.997208, 195.997711, 207.652344, 220., 233.081879, 246.94165, 261.62558, 277.182617,293.664764, 311.126984, 329.627563, 349.228241, 369.994415, 391.995422, 415.304688, 440., 466.163757, 493.883301, 523.25116, 554.365234, 587.329529, 622.253967, 659.255127, 698.456482, 739.988831, 783.990845, 830.609375, 880., 932.327515, 987.766602, 1046.502319, 1108.730469, 1174.659058, 1244.507935, 1318.510254, 1396.912964, 1479.977661, 1567.981689, 1661.21875, 1760., 1864.655029, 1975.533203, 2093.004639, 2217.460938, 2349.318115, 2489.015869, 2637.020508, 2793.825928, 2959.955322, 3135.963379, 3322.4375, 3520., 3729.31, 3951.066406, 4186.009277, 4434.921875, 4698.63623, 4978.031738, 5274.041016, 5587.651855, 5919.910645, 6271.926758, 6644.875, 7040., 7458.620117, 7902.132812, 8372.018555, 8869.84375, 9397.272461, 9956.063477, 10548.082031, 11175.303711, 11839.821289, 12543.853516, 13289.75};
void setup();//use this to do any initialisation if you want.
void play(double *channels);//run dac!
maxiOsc::maxiOsc(){
//When you create an oscillator, the constructor sets the phase of the oscillator to 0.
phase = 0.0;
}
double maxiOsc::noise() {
//White Noise
//always the same unless you seed it.
float r = rand()/(float)RAND_MAX;
output=r*2-1;
return(output);
}
void maxiOsc::phaseReset(double phaseIn) {
//This allows you to set the phase of the oscillator to anything you like.
phase=phaseIn;
}
double maxiOsc::sinewave(double frequency) {
//This is a sinewave oscillator
output=sin (phase*(TWOPI));
if ( phase >= 1.0 ) phase -= 1.0;
phase += (1./(maxiSettings::sampleRate/(frequency)));
return(output);
}
double maxiOsc::sinebuf4(double frequency) {
//This is a sinewave oscillator that uses 4 point interpolation on a 514 point buffer
double remainder;
double a,b,c,d,a1,a2,a3;
phase += 512./(maxiSettings::sampleRate/(frequency));
if ( phase >= 511 ) phase -=512;
remainder = phase - floor(phase);
if (phase==0) {
a=sineBuffer[(long) 512];
b=sineBuffer[(long) phase];
c=sineBuffer[(long) phase+1];
d=sineBuffer[(long) phase+2];
} else {
a=sineBuffer[(long) phase-1];
b=sineBuffer[(long) phase];
c=sineBuffer[(long) phase+1];
d=sineBuffer[(long) phase+2];
}
a1 = 0.5f * (c - a);
a2 = a - 2.5 * b + 2.f * c - 0.5f * d;
a3 = 0.5f * (d - a) + 1.5f * (b - c);
output = double (((a3 * remainder + a2) * remainder + a1) * remainder + b);
return(output);
}
double maxiOsc::sinebuf(double frequency) { //specify the frequency of the oscillator in Hz / cps etc.
//This is a sinewave oscillator that uses linear interpolation on a 514 point buffer
double remainder;
phase += 512./(maxiSettings::sampleRate/(frequency*chandiv));
if ( phase >= 511 ) phase -=512;
remainder = phase - floor(phase);
output = (double) ((1-remainder) * sineBuffer[1+ (long) phase] + remainder * sineBuffer[2+(long) phase]);
return(output);
}
double maxiOsc::coswave(double frequency) {
//This is a cosine oscillator
output=cos (phase*(TWOPI));
if ( phase >= 1.0 ) phase -= 1.0;
phase += (1./(maxiSettings::sampleRate/(frequency)));
return(output);
}
double maxiOsc::phasor(double frequency) {
//This produces a floating point linear ramp between 0 and 1 at the desired frequency
output=phase;
if ( phase >= 1.0 ) phase -= 1.0;
phase += (1./(maxiSettings::sampleRate/(frequency)));
return(output);
}
double maxiOsc::square(double frequency) {
//This is a square wave
if (phase<0.5) output=-1;
if (phase>0.5) output=1;
if ( phase >= 1.0 ) phase -= 1.0;
phase += (1./(maxiSettings::sampleRate/(frequency)));
return(output);
}
double maxiOsc::pulse(double frequency, double duty) {
//This is a pulse generator that creates a signal between -1 and 1.
if (duty<0.) duty=0;
if (duty>1.) duty=1;
if ( phase >= 1.0 ) phase -= 1.0;
phase += (1./(maxiSettings::sampleRate/(frequency)));
if (phase<duty) output=-1.;
if (phase>duty) output=1.;
return(output);
}
double maxiOsc::phasor(double frequency, double startphase, double endphase) {
//This is a phasor that takes a value for the start and end of the ramp.
output=phase;
if (phase<startphase) {
phase=startphase;
}
if ( phase >= endphase ) phase = startphase;
phase += ((endphase-startphase)/(maxiSettings::sampleRate/(frequency)));
return(output);
}
double maxiOsc::saw(double frequency) {
//Sawtooth generator. This is like a phasor but goes between -1 and 1
output=phase;
if ( phase >= 1.0 ) phase -= 2.0;
phase += (1./(maxiSettings::sampleRate/(frequency)));
return(output);
}
double maxiOsc::sawn(double frequency) {
//Bandlimited sawtooth generator. Woohoo.
if ( phase >= 0.5 ) phase -= 1.0;
phase += (1./(maxiSettings::sampleRate/(frequency)));
double temp=(8820.22/frequency)*phase;
if (temp<-0.5) {
temp=-0.5;
}
if (temp>0.5) {
temp=0.5;
}
temp*=1000.0f;
temp+=500.0f;
double remainder = temp - floor(temp);
output = (double) ((1.0f-remainder) * transition[(long)temp] + remainder * transition[1+(long)temp]) - phase;
return(output);
}
double maxiOsc::rect(double frequency, double duty) {
return (output);
}
double maxiOsc::triangle(double frequency) {
//This is a triangle wave.
if ( phase >= 1.0 ) phase -= 1.0;
phase += (1./(maxiSettings::sampleRate/(frequency)));
if (phase <= 0.5 ) {
output =(phase - 0.25) * 4;
} else {
output =((1.0-phase) - 0.25) * 4;
}
return(output);
}
// don't use this nonsense. Use ramps instead.
// ..er... I mean "This method is deprecated"
double maxiEnvelope::line(int numberofsegments,double segments[1000]) {
//This is a basic multi-segment ramp generator that you can use for more or less anything.
//However, it's not that intuitive.
if (isPlaying==1) {//only make a sound once you've been triggered
period=4./(segments[valindex+1]*0.0044);
nextval=segments[valindex+2];
currentval=segments[valindex];
if (currentval-amplitude > 0.0000001 && valindex < numberofsegments) {
amplitude += ((currentval-startVal)/(maxiSettings::sampleRate/period));
} else if (currentval-amplitude < -0.0000001 && valindex < numberofsegments) {
amplitude -= (((currentval-startVal)*(-1))/(maxiSettings::sampleRate/period));
} else if (valindex >numberofsegments-1) {
valindex=numberofsegments-2;
} else {
valindex=valindex+2;
startVal=currentval;
}
output=amplitude;
}
else {
output=0;
}
return(output);
}
//and this is also deprecated
void maxiEnvelope::trigger(int index, double amp) {
isPlaying=1;//ok the envelope is being used now.
valindex=index;
amplitude=amp;
}
void maxiEnvelope::trigger(bool noteOn) {
if (noteOn) trig=1;
if (noteOn==false) trig=0;
}
double maxiEnvelope::ramp(double startVal, double endVal, double duration){
if (trig!=0) {
phase=startVal;
isPlaying=true;
trig=0;
}
if (isPlaying) {
if (startVal<endVal) {
phase += ((endVal-startVal)/(maxiSettings::sampleRate/(1./duration)));
if ( phase >= endVal ) phase = endVal;
}
if (startVal > endVal) {
phase += ((endVal-startVal)/(maxiSettings::sampleRate/(1./duration)));
if ( phase <= endVal ) phase = endVal;
}
return(phase);
} else {
return(0);
}
}
double maxiEnvelope::ramps(std::vector<double> rampsArray){
if (trig!=0) {
valindex=0;
endVal=rampsArray[valindex+1];
isPlaying=true;
trig=0;
}
if (isPlaying) {
if (valindex>0 && rampsArray[valindex-1]==rampsArray[valindex+1]) {
period += (1/(maxiSettings::sampleRate/(1./rampsArray[valindex])));
if (period>=1) {
phase=endVal;
startVal=phase;
if (valindex+2<rampsArray.size()){
valindex+=2;
endVal=rampsArray[valindex+1];
period=0;
}
} output=phase;
}
if (valindex==0 && output==endVal) {
period += (1/(maxiSettings::sampleRate/(1./rampsArray[valindex])));
if (period>=1) {
phase=endVal;
startVal=phase;
if (valindex+2<rampsArray.size()){
valindex+=2;
endVal=rampsArray[valindex+1];
period=0;
}
} output=phase;
}
if (phase<endVal) {
phase += ((endVal-startVal)/(maxiSettings::sampleRate/(1./rampsArray[valindex])));
if ( phase >= endVal ) {
phase=endVal;
startVal=phase;
if (valindex+2<rampsArray.size()){
valindex+=2;
endVal=rampsArray[valindex+1];
}
}
output=phase;
}
if (phase > endVal) {
phase += ((endVal-startVal)/(maxiSettings::sampleRate/(1./rampsArray[valindex])));
if ( phase <= endVal ) {
phase=endVal;
startVal=phase;
if (valindex+2<rampsArray.size()){
valindex+=2;
endVal=rampsArray[valindex+1];
}
} output=phase;
}
return(output);
} else {
return(0);
}
}
double maxiEnvelope::ar(double attack, double release) {
if (trig!=0) {
//phase=0;
releaseMode=false;
trig=0;
}
if (phase<1 && releaseMode==false) {
phase += ((1)/(maxiSettings::sampleRate/(1./attack)));
if ( phase >= 1 ) {
phase=1;
releaseMode=true;
};
}
if (releaseMode==true) {
phase += ((-1)/(maxiSettings::sampleRate/(1./release)));
if ( phase <= 0 ) phase = 0;
}
return phase;
}
double maxiEnvelope::adsr(double attack, double decay, double sustain, double release) {
if (trig!=0 && !attackMode) {
// phase=0.;
releaseMode=false;
decayMode=false;
sustainMode=false;
attackMode=true;
trig=0;
}
if (attackMode) {
phase += ((1)/(maxiSettings::sampleRate/(1./attack)));
if ( phase >= 1 ) {
phase=1;
attackMode=false;
decayMode=true;
};
}
if (decayMode) {
phase += ((-1)/(maxiSettings::sampleRate/(1./decay)));
if ( phase <= sustain ) {
phase=sustain;
decayMode=false;
sustainMode=true;
};
}
if (sustainMode) {
if (noteOn) {
phase=sustain;
}
if (!noteOn) {
sustainMode=false;
releaseMode=true;
}
}
if (releaseMode) {
phase += ((-sustain)/(maxiSettings::sampleRate/(1./release)));
if ( phase <= 0 ) {
phase = 0;
releaseMode=false;
}
}
return phase;
}
//Delay with feedback
maxiDelayline::maxiDelayline() {
memset( memory, 0, 88200*sizeof (double) );
}
double maxiDelayline::dl(double input, int size, double feedback) {
if ( phase >=size ) {
phase = 0;
}
output=memory[phase];
memory[phase]=(memory[phase]*feedback)+(input*feedback)*0.5;
phase+=1;
return(output);
}
double maxiDelayline::dl(double input, int size, double feedback, int position) {
if ( phase >=size ) phase = 0;
if ( position >=size ) position = 0;
output=memory[position];
memory[phase]=(memory[phase]*feedback)+(input*feedback)*chandiv;
phase+=1;
return(output);
}
maxiFractionalDelay::maxiFractionalDelay ( void ) {
memset( memory, 0, delaySize*sizeof (double) );
}
double maxiFractionalDelay::dl ( double sig, double delayTime, double feedback )
{
// Set delay time
delayTime = fmin(fabs(delayTime), delaySize);
int32_t delay = delayTime; // Truncated
double fractAmount = delayTime - delay; // Fractional remainder
double truncAmount = 1.0f - fractAmount; // Inverse fractional remainder
// Update read pointer
readPointer = writePointer - delay;
if (readPointer < 0)
readPointer += delaySize;
int readPointerFractPart = readPointer-1;
if (readPointerFractPart < 0)
readPointerFractPart += delaySize;
// Get interpolated sample
double y = memory[readPointer] * truncAmount + memory[readPointerFractPart] * fractAmount;
// Write new sample
memory[writePointer] = y * feedback + sig;
// Increment write pointer
if (++writePointer >= delaySize)
writePointer -= delaySize;
return y;
}
//I particularly like these. cutoff between 0 and 1
double maxiFilter::lopass(double input, double cutoff) {
output=outputs[0] + cutoff*(input-outputs[0]);
outputs[0]=output;
return(output);
}
//as above
double maxiFilter::hipass(double input, double cutoff) {
output=input-(outputs[0] + cutoff*(input-outputs[0]));
outputs[0]=output;
return(output);
}
//awesome. cuttof is freq in hz. res is between 1 and whatever. Watch out!
double maxiFilter::lores(double input,double cutoff1, double resonance) {
cutoff=cutoff1;
if (cutoff<10) cutoff=10;
if (cutoff>(maxiSettings::sampleRate)) cutoff=(maxiSettings::sampleRate);
if (resonance<1.) resonance = 1.;
z=cos(TWOPI*cutoff/maxiSettings::sampleRate);
c=2-2*z;
double r=(sqrt(2.0)*sqrt(-pow((z-1.0),3.0))+resonance*(z-1))/(resonance*(z-1));
x=x+(input-y)*c;
y=y+x;
x=x*r;
output=y;
return(output);
}
//working hires filter
double maxiFilter::hires(double input,double cutoff1, double resonance) {
cutoff=cutoff1;
if (cutoff<10) cutoff=10;
if (cutoff>(maxiSettings::sampleRate)) cutoff=(maxiSettings::sampleRate);
if (resonance<1.) resonance = 1.;
z=cos(TWOPI*cutoff/maxiSettings::sampleRate);
c=2-2*z;
double r=(sqrt(2.0)*sqrt(-pow((z-1.0),3.0))+resonance*(z-1))/(resonance*(z-1));
x=x+(input-y)*c;
y=y+x;
x=x*r;
output=input-y;
return(output);
}
//This works a bit. Needs attention.
double maxiFilter::bandpass(double input,double cutoff1, double resonance) {
cutoff=cutoff1;
if (cutoff>(maxiSettings::sampleRate*0.5)) cutoff=(maxiSettings::sampleRate*0.5);
if (resonance>=1.) resonance=0.999999;
z=cos(TWOPI*cutoff/maxiSettings::sampleRate);
inputs[0] = (1-resonance)*(sqrt(resonance*(resonance-4.0*pow(z,2.0)+2.0)+1));
inputs[1] = 2*z*resonance;
inputs[2] = pow((resonance*-1),2);
output=inputs[0]*input+inputs[1]*outputs[1]+inputs[2]*outputs[2];
outputs[2]=outputs[1];
outputs[1]=output;
return(output);
}
//stereo bus
double *maxiMix::stereo(double input,double two[2],double x) {
if (x>1) x=1;
if (x<0) x=0;
two[0]=input*sqrt(1.0-x);
two[1]=input*sqrt(x);
return(two);
}
//quad bus
double *maxiMix::quad(double input,double four[4],double x,double y) {
if (x>1) x=1;
if (x<0) x=0;
if (y>1) y=1;
if (y<0) y=0;
four[0]=input*sqrt((1.0-x)*y);
four[1]=input*sqrt((1.0-x)*(1.0-y));
four[2]=input*sqrt(x*y);
four[3]=input*sqrt(x*(1.0-y));
return(four);
}
//ambisonic bus
double *maxiMix::ambisonic(double input,double eight[8],double x,double y,double z) {
if (x>1) x=1;
if (x<0) x=0;
if (y>1) y=1;
if (y<0) y=0;
if (z>1) y=1;
if (z<0) y=0;
eight[0]=input*(sqrt((1.0-x)*y)*1.0-z);
eight[1]=input*(sqrt((1.0-x)*(1.0-y))*1.0-z);
eight[2]=input*(sqrt(x*y)*1.0-z);
eight[3]=input*(sqrt(x*(1.0-y))*1.0-z);
eight[4]=input*(sqrt((1.0-x)*y)*z);
eight[5]=input*(sqrt((1.0-x)*(1.0-y))*z);
eight[6]=input*sqrt((x*y)*z);
eight[7]=input*sqrt((x*(1.0-y))*z);
return(eight);
}
//This is the maxiSample load function. It just calls read.
bool maxiSample::load(string fileName, int channel) {
myPath = fileName;
readChannel=channel;
return read();
}
// This is for OGG loading
bool maxiSample::loadOgg(string fileName, int channel) {
#ifdef VORBIS
bool result;
readChannel=channel;
int channelx;
// cout << fileName << endl;
free(temp);
myDataSize = stb_vorbis_decode_filename(const_cast<char*>(fileName.c_str()), &channelx, &temp);
result = myDataSize > 0;
printf("\nchannels = %d\nlength = %d",channelx,myDataSize);
printf("\n");
myChannels=(short)channelx;
length=myDataSize;
mySampleRate=44100;
if (myChannels>1) {
int position=0;
int channel=readChannel;
for (int i=channel;i<myDataSize*2;i+=myChannels) {
temp[position]=temp[i];
position++;
}
}
return result; // this should probably be something more descriptive
#else
assert(false); // called but VORBIS not defined!
#endif
return 0;
}
//This sets the playback position to the start of a sample
void maxiSample::trigger() {
position = 0;
recordPosition = 0;
}
//This is the main read function.
bool maxiSample::read()
{
bool result;
ifstream inFile( myPath.c_str(), ios::in | ios::binary);
result = inFile.is_open();
if (result) {
bool datafound = false;
inFile.seekg(4, ios::beg);
inFile.read( (char*) &myChunkSize, 4 ); // read the ChunkSize
inFile.seekg(16, ios::beg);
inFile.read( (char*) &mySubChunk1Size, 4 ); // read the SubChunk1Size
//inFile.seekg(20, ios::beg);
inFile.read( (char*) &myFormat, sizeof(short) ); // read the file format. This should be 1 for PCM
//inFile.seekg(22, ios::beg);
inFile.read( (char*) &myChannels, sizeof(short) ); // read the # of channels (1 or 2)
//inFile.seekg(24, ios::beg);
inFile.read( (char*) &mySampleRate, sizeof(int) ); // read the samplerate
//inFile.seekg(28, ios::beg);
inFile.read( (char*) &myByteRate, sizeof(int) ); // read the byterate
//inFile.seekg(32, ios::beg);
inFile.read( (char*) &myBlockAlign, sizeof(short) ); // read the blockalign
//inFile.seekg(34, ios::beg);
inFile.read( (char*) &myBitsPerSample, sizeof(short) ); // read the bitspersample
//ignore any extra chunks
char chunkID[5]="";
chunkID[4] = 0;
int filePos = 20 + mySubChunk1Size;
while(!datafound && !inFile.eof()) {
inFile.seekg(filePos, ios::beg);
inFile.read((char*) &chunkID, sizeof(char) * 4);
inFile.seekg(filePos + 4, ios::beg);
inFile.read( (char*) &myDataSize, sizeof(int) ); // read the size of the data
filePos += 8;
if (strcmp(chunkID,"data") == 0) {
datafound = true;
}else{
filePos += myDataSize;
}
}
// read the data chunk
char * myData = (char*) malloc(myDataSize * sizeof(char));
inFile.seekg(filePos, ios::beg);
inFile.read(myData, myDataSize);
length=myDataSize*(0.5/myChannels);
inFile.close(); // close the input file
cout << "Ch: " << myChannels << ", len: " << length << endl;
if (myChannels>1) {
int position=0;
int channel=readChannel*2;
for (int i=channel;i<myDataSize;i+=(myChannels*2)) {
myData[position]=myData[i];
myData[position+1]=myData[i+1];
position+=2;
}
}
free(temp);
temp = (short*) malloc(myDataSize * sizeof(char));
memcpy(temp, myData, myDataSize * sizeof(char));
free(myData);
}else {
// cout << "ERROR: Could not load sample: " <<myPath << endl; //This line seems to be hated by windows
printf("ERROR: Could not load sample.");
}
return result; // this should probably be something more descriptive
}
//This plays back at the correct speed. Always loops.
double maxiSample::play() {
position++;
if ((long) position >= length) position=0;
output = (double) temp[(long)position]/32767.0;
return output;
}
void maxiSample::setPosition(double newPos) {
position = maxiMap::clamp<double>(newPos, 0.0, 1.0) * length;
}
//start end and points are between 0 and 1
double maxiSample::playLoop(double start, double end) {
position++;
if (position < length * start) position = length * start;
if ((long) position >= length * end) position = length * start;
output = (double) temp[(long)position]/32767.0;
return output;
}
double maxiSample::playUntil(double end) {
position++;
if ((long) position<length * end)
output = (double) temp[(long)position]/32767.0;
else {
output=0;
}
return output;
}
//This plays back at the correct speed. Only plays once. To retrigger, you have to manually reset the position
double maxiSample::playOnce() {
position++;
if ((long) position<length)
output = (double) temp[(long)position]/32767.0;
else {
output=0;
}
return output;
}
//Same as above but takes a speed value specified as a ratio, with 1.0 as original speed
double maxiSample::playOnce(double speed) {
position=position+((speed*chandiv)/(maxiSettings::sampleRate/mySampleRate));
double remainder = position - (long) position;
if ((long) position<length)
output = (double) ((1-remainder) * temp[1+ (long) position] + remainder * temp[2+(long) position])/32767;//linear interpolation
else
output=0;
return(output);
}
//As above but looping
double maxiSample::play(double speed) {
double remainder;
long a,b;
position=position+((speed*chandiv)/(maxiSettings::sampleRate/mySampleRate));
if (speed >=0) {
if ((long) position>=length-1) position=1;
remainder = position - floor(position);
if (position+1<length) {
a=position+1;
}
else {
a=length-1;
}
if (position+2<length)
{
b=position+2;
}
else {
b=length-1;
}
output = (double) ((1-remainder) * temp[a] + remainder * temp[b])/32767;//linear interpolation
} else {
if ((long) position<0) position=length;
remainder = position - floor(position);
if (position-1>=0) {
a=position-1;
}
else {
a=0;
}
if (position-2>=0) {
b=position-2;
}
else {
b=0;
}
output = (double) ((-1-remainder) * temp[a] + remainder * temp[b])/32767;//linear interpolation
}
return(output);
}
//placeholder
double maxiSample::play(double frequency, double start, double end) {
return play(frequency, start, end, position);
}
//This allows you to say how often a second you want a specific chunk of audio to play
double maxiSample::play(double frequency, double start, double end, double &pos) {
double remainder;
if (end>=length) end=length-1;
long a,b;
if (frequency >0.) {
if (pos<start) {
pos=start;
}
if ( pos >= end ) pos = start;
pos += ((end-start)/((maxiSettings::sampleRate)/(frequency*chandiv)));
remainder = pos - floor(pos);
long posl = floor(pos);
if (posl+1<length) {
a=posl+1;
}
else {
a=posl-1;
}
if (posl+2<length) {
b=posl+2;
}
else {
b=length-1;
}
output = (double) ((1-remainder) * temp[a] +
remainder * temp[b])/32767;//linear interpolation
} else {
frequency*=-1.;
if ( pos <= start ) pos = end;
pos -= ((end-start)/(maxiSettings::sampleRate/(frequency*chandiv)));
remainder = pos - floor(pos);
long posl = floor(pos);
if (posl-1>=0) {
a=posl-1;
}
else {
a=0;
}
if (posl-2>=0) {
b=posl-2;
}
else {
b=0;
}
output = (double) ((-1-remainder) * temp[a] +
remainder * temp[b])/32767;//linear interpolation
}
return(output);
}
//Same as above. better cubic inerpolation. Cobbled together from various (pd externals, yehar, other places).
double maxiSample::play4(double frequency, double start, double end) {
double remainder;
double a,b,c,d,a1,a2,a3;
if (frequency >0.) {
if (position<start) {
position=start;
}
if ( position >= end ) position = start;
position += ((end-start)/(maxiSettings::sampleRate/(frequency*chandiv)));
remainder = position - floor(position);
if (position>0) {
a=temp[(int)(floor(position))-1];
} else {
a=temp[0];
}
b=temp[(long) position];
if (position<end-2) {
c=temp[(long) position+1];
} else {
c=temp[0];
}
if (position<end-3) {
d=temp[(long) position+2];
} else {
d=temp[0];
}
a1 = 0.5f * (c - a);
a2 = a - 2.5 * b + 2.f * c - 0.5f * d;
a3 = 0.5f * (d - a) + 1.5f * (b - c);
output = (double) (((a3 * remainder + a2) * remainder + a1) * remainder + b) / 32767;
} else {
frequency*=-1.;
if ( position <= start ) position = end;
position -= ((end-start)/(maxiSettings::sampleRate/(frequency*chandiv)));
remainder = position - floor(position);
if (position>start && position < end-1) {
a=temp[(long) position+1];
} else {
a=temp[0];
}
b=temp[(long) position];
if (position>start) {
c=temp[(long) position-1];
} else {
c=temp[0];
}
if (position>start+1) {
d=temp[(long) position-2];
} else {
d=temp[0];
}
a1 = 0.5f * (c - a);
a2 = a - 2.5 * b + 2.f * c - 0.5f * d;
a3 = 0.5f * (d - a) + 1.5f * (b - c);
output = (double) (((a3 * remainder + a2) * -remainder + a1) * -remainder + b) / 32767;
}
return(output);
}
//You don't need to worry about this stuff.
double maxiSample::bufferPlay(unsigned char &bufferin,long length) {
double remainder;
short* buffer = (short *)&bufferin;
position=(position+1);
remainder = position - (long) position;
if ((long) position>length) position=0;
output = (double) ((1-remainder) * buffer[1+ (long) position] + remainder * buffer[2+(long) position])/32767;//linear interpolation
return(output);
}
double maxiSample::bufferPlay(unsigned char &bufferin,double speed,long length) {
double remainder;
long a,b;
short* buffer = (short *)&bufferin;
position=position+((speed*chandiv)/(maxiSettings::sampleRate/mySampleRate));
if (speed >=0) {
if ((long) position>=length-1) position=1;
remainder = position - floor(position);
if (position+1<length) {
a=position+1;
}
else {
a=length-1;
}
if (position+2<length)
{
b=position+2;
}
else {
b=length-1;
}
output = (double) ((1-remainder) * buffer[a] + remainder * buffer[b])/32767;//linear interpolation
} else {
if ((long) position<0) position=length;
remainder = position - floor(position);
if (position-1>=0) {
a=position-1;
}
else {
a=0;
}
if (position-2>=0) {
b=position-2;
}
else {
b=0;
}
output = (double) ((-1-remainder) * buffer[a] + remainder * buffer[b])/32767;//linear interpolation
}
return(output);
}
double maxiSample::bufferPlay(unsigned char &bufferin,double frequency, double start, double end) {
double remainder;
length=end;
long a,b;
short* buffer = (short *)&bufferin;
if (frequency >0.) {
if (position<start) {
position=start;
}
if ( position >= end ) position = start;
position += ((end-start)/(maxiSettings::sampleRate/(frequency*chandiv)));
remainder = position - floor(position);
long pos = floor(position);
if (pos+1<length) {
a=pos+1;
}
else {
a=pos-1;
}
if (pos+2<length) {
b=pos+2;
}
else {
b=length-1;
}
output = (double) ((1-remainder) * buffer[a] +
remainder * buffer[b])/32767;//linear interpolation
} else {
frequency*=-1.;
if ( position <= start ) position = end;
position -= ((end-start)/(maxiSettings::sampleRate/(frequency*chandiv)));
remainder = position - floor(position);
long pos = floor(position);
if (pos-1>=0) {
a=pos-1;
}
else {
a=0;
}
if (pos-2>=0) {
b=pos-2;
}
else {
b=0;
}
output = (double) ((-1-remainder) * buffer[a] +
remainder * buffer[b])/32767;//linear interpolation
}
return(output);
}
//better cubic inerpolation. Cobbled together from various (pd externals, yehar, other places).
double maxiSample::bufferPlay4(unsigned char &bufferin,double frequency, double start, double end) {
double remainder;
double a,b,c,d,a1,a2,a3;
short* buffer = (short*)&bufferin;
if (frequency >0.) {
if (position<start) {
position=start;
}
if ( position >= end ) position = start;
position += ((end-start)/(maxiSettings::sampleRate/(frequency*chandiv)));
remainder = position - floor(position);
if (position>0) {
a=buffer[(int)(floor(position))-1];
} else {
a=buffer[0];
}
b=buffer[(long) position];
if (position<end-2) {
c=buffer[(long) position+1];
} else {
c=buffer[0];
}
if (position<end-3) {
d=buffer[(long) position+2];
} else {
d=buffer[0];
}
a1 = 0.5f * (c - a);
a2 = a - 2.5 * b + 2.f * c - 0.5f * d;
a3 = 0.5f * (d - a) + 1.5f * (b - c);
output = (double) (((a3 * remainder + a2) * remainder + a1) * remainder + b) / 32767;
} else {
frequency*=-1.;
if ( position <= start ) position = end;
position -= ((end-start)/(maxiSettings::sampleRate/(frequency*chandiv)));
remainder = position - floor(position);
if (position>start && position < end-1) {
a=buffer[(long) position+1];
} else {
a=buffer[0];
}
b=buffer[(long) position];
if (position>start) {
c=buffer[(long) position-1];
} else {
c=buffer[0];
}
if (position>start+1) {
d=buffer[(long) position-2];
} else {
d=buffer[0];
}
a1 = 0.5f * (c - a);
a2 = a - 2.5 * b + 2.f * c - 0.5f * d;
a3 = 0.5f * (d - a) + 1.5f * (b - c);
output = (double) (((a3 * remainder + a2) * -remainder + a1) * -remainder + b) / 32767;
}
return(output);
}
long maxiSample::getLength() {
return (length=myDataSize*0.5);
}
void maxiSample::setLength(unsigned long numSamples) {
cout << "Length: " << numSamples << endl;
short *newData = (short*) malloc(sizeof(short) * numSamples);
if (NULL!=temp) {
unsigned long copyLength = min((unsigned long)length, numSamples);
memcpy(newData, temp, sizeof(short) * copyLength);
}
temp = newData;
myDataSize = int(numSamples * 2);
length=numSamples;
position=0;
recordPosition=0;
}
void maxiSample::clear() {
memset(temp, 0, myDataSize);
}
void maxiSample::reset() {
position=0;
}
void maxiSample::normalise(float maxLevel) {
short maxValue = 0;
for(int i=0; i < length; i++) {
if (abs(temp[i]) > maxValue) {
maxValue = abs(temp[i]);
}
}
float scale = 32767.0 * maxLevel / (float) maxValue;
for(int i=0; i < length; i++) {
temp[i] = round(scale * (float) temp[i]);
}
}
void maxiSample::autoTrim(float alpha, float threshold, bool trimStart, bool trimEnd) {
int startMarker=0;
if(trimStart) {
maxiLagExp<float> startLag(alpha, 0);
while(startMarker < length) {
startLag.addSample(abs(temp[startMarker]));
if (startLag.value() > threshold) {
break;
}
startMarker++;
}
}
int endMarker = int(length-1);
if(trimEnd) {
maxiLagExp<float> endLag(alpha, 0);
while(endMarker > 0) {
endLag.addSample(abs(temp[endMarker]));
if (endLag.value() > threshold) {
break;
}
endMarker--;
}
}
cout << "Autotrim: start: " << startMarker << ", end: " << endMarker << endl;
int newLength = endMarker - startMarker;
if (newLength > 0) {
short *newData = (short*) malloc(sizeof(short) * newLength);
for(int i=0; i < newLength; i++) {
newData[i] = temp[i+startMarker];
}
free(temp);
temp = newData;
myDataSize = newLength * 2;
length=newLength;
position=0;
recordPosition=0;
//envelope the start
int fadeSize=int(min((long)100, length));
for(int i=0; i < fadeSize; i++) {
float factor = i / (float) fadeSize;
temp[i] = round(temp[i] * factor);
temp[length - 1 - i] = round(temp[length - 1 - i] * factor);
}
}
}
/* OK this compressor and gate are now ready to use. The envelopes, like all the envelopes in this recent update, use stupid algorithms for
incrementing - consequently a long attack is something like 0.0001 and a long release is like 0.9999.
Annoyingly, a short attack is 0.1, and a short release is 0.99. I'll sort this out laters */
double maxiDyn::gate(double input, double threshold, long holdtime, double attack, double release) {
if (fabs(input)>threshold && attackphase!=1){
holdcount=0;
releasephase=0;
attackphase=1;
if(amplitude==0) amplitude=0.01;
}
if (attackphase==1 && amplitude<1) {
amplitude*=(1+attack);
output=input*amplitude;
}
if (amplitude>=1) {
attackphase=0;
holdphase=1;
}
if (holdcount<holdtime && holdphase==1) {
output=input;
holdcount++;
}
if (holdcount==holdtime) {
holdphase=0;
releasephase=1;
}
if (releasephase==1 && amplitude>0.) {
output=input*(amplitude*=release);
}
return output;
}
double maxiDyn::compressor(double input, double ratio, double threshold, double attack, double release) {
if (fabs(input)>threshold && attackphase!=1){
holdcount=0;
releasephase=0;
attackphase=1;
if(currentRatio==0) currentRatio=ratio;
}
if (attackphase==1 && currentRatio<ratio-1) {
currentRatio*=(1+attack);
}
if (currentRatio>=ratio-1) {
attackphase=0;
releasephase=1;
}
if (releasephase==1 && currentRatio>0.) {
currentRatio*=release;
}
if (input>0.) {
output = input/(1.+currentRatio);
} else {
output = input/(1.+currentRatio);
}
return output*(1+log(ratio));
}
double maxiDyn::compress(double input) {
if (fabs(input)>threshold && attackphase!=1){
holdcount=0;
releasephase=0;
attackphase=1;
if(currentRatio==0) currentRatio=ratio;
}
if (attackphase==1 && currentRatio<ratio-1) {
currentRatio*=(1+attack);
}
if (currentRatio>=ratio-1) {
attackphase=0;
releasephase=1;
}
if (releasephase==1 && currentRatio>0.) {
currentRatio*=release;
}
if (input>0.) {
output = input/(1.+currentRatio);
} else {
output = input/(1.+currentRatio);
}
return output*(1+log(ratio));
}
/* Lots of people struggle with the envelope generators so here's a new easy one.
It takes mental numbers for attack and release tho. Basically, they're exponentials.
I'll map them out later so that it's a bit more intuitive */
double maxiEnv::ar(double input, double attack, double release, long holdtime, int trigger) {
if (trigger==1 && attackphase!=1 && holdphase!=1){
holdcount=0;
releasephase=0;
attackphase=1;
}
if (attackphase==1) {
amplitude+=(1*attack);
output=input*amplitude;
}
if (amplitude>=1) {
amplitude=1;
attackphase=0;
holdphase=1;
}
if (holdcount<holdtime && holdphase==1) {
output=input;
holdcount++;
}
if (holdcount==holdtime && trigger==1) {
output=input;
}
if (holdcount==holdtime && trigger!=1) {
holdphase=0;
releasephase=1;
}
if (releasephase==1 && amplitude>0.) {
output=input*(amplitude*=release);
}
return output;
}
/* adsr. It's not bad, very simple to use*/
double maxiEnv::adsr(double input, double attack, double decay, double sustain, double release, long holdtime, int trigger) {
if (trigger==1 && attackphase!=1 && holdphase!=1 && decayphase!=1){
holdcount=0;
decayphase=0;
sustainphase=0;
releasephase=0;
attackphase=1;
}
if (attackphase==1) {
releasephase=0;
amplitude+=(1*attack);
output=input*amplitude;
if (amplitude>=1) {
amplitude=1;
attackphase=0;
decayphase=1;
}
}
if (decayphase==1) {
output=input*(amplitude*=decay);
if (amplitude<=sustain) {
decayphase=0;
holdphase=1;
}
}
if (holdcount<holdtime && holdphase==1) {
output=input*amplitude;
holdcount++;
}
if (holdcount>=holdtime && trigger==1) {
output=input*amplitude;
}
if (holdcount>=holdtime && trigger!=1) {
holdphase=0;
releasephase=1;
}
if (releasephase==1 && amplitude>0.) {
output=input*(amplitude*=release);
}
return output;
}
double maxiEnv::adsr(double input, int trigger) {
if (trigger==1 && attackphase!=1 && holdphase!=1 && decayphase!=1){
holdcount=0;
decayphase=0;
sustainphase=0;
releasephase=0;
attackphase=1;
}
if (attackphase==1) {
releasephase=0;
amplitude+=(1*attack);
output=input*amplitude;
if (amplitude>=1) {
amplitude=1;
attackphase=0;
decayphase=1;
}
}
if (decayphase==1) {
output=input*(amplitude*=decay);
if (amplitude<=sustain) {
decayphase=0;
holdphase=1;
}
}
if (holdcount<holdtime && holdphase==1) {
output=input*amplitude;
holdcount++;
}
if (holdcount>=holdtime && trigger==1) {
output=input*amplitude;
}
if (holdcount>=holdtime && trigger!=1) {
holdphase=0;
releasephase=1;
}
if (releasephase==1 && amplitude>0.) {
output=input*(amplitude*=release);
}
return output;
}
void maxiEnv::setAttack(double attackMS) {
attack = 1-pow( 0.01, 1.0 / ( attackMS * maxiSettings::sampleRate * 0.001 ) );
}
void maxiEnv::setRelease(double releaseMS) {
release = pow( 0.01, 1.0 / ( releaseMS * maxiSettings::sampleRate * 0.001 ) );
}
void maxiEnv::setSustain(double sustainL) {
sustain = sustainL;
}
void maxiEnv::setDecay(double decayMS) {
decay = pow( 0.01, 1.0 / ( decayMS * maxiSettings::sampleRate * 0.001 ) );
}
void maxiDyn::setAttack(double attackMS) {
attack = pow( 0.01, 1.0 / ( attackMS * maxiSettings::sampleRate * 0.001 ) );
}
void maxiDyn::setRelease(double releaseMS) {
release = pow( 0.01, 1.0 / ( releaseMS * maxiSettings::sampleRate * 0.001 ) );
}
void maxiDyn::setThreshold(double thresholdI) {
threshold = thresholdI;
}
void maxiDyn::setRatio(double ratioF) {
ratio = ratioF;
}
double convert::mtof(int midinote) {
return mtofarray[midinote];
}
int convert::ftom(double frequency) {
double baseFrequency = 440;
return round(12 * log2(frequency/baseFrequency)) + 69;
}
double convert::atodb(double amplitude) {
return 20 * log10(amplitude);
}
double convert::dbtoa(double decibels) {
return pow(10, (decibels * 0.5));
}
template<> void maxiEnvelopeFollower::setAttack(double attackMS) {
attack = pow( 0.01, 1.0 / ( attackMS * maxiSettings::sampleRate * 0.001 ) );
}
template<> void maxiEnvelopeFollower::setRelease(double releaseMS) {
release = pow( 0.01, 1.0 / ( releaseMS * maxiSettings::sampleRate * 0.001 ) );
}
double pitchRatios[256]= {0.0006517771980725,0.0006905338959768,0.0007315951515920, 0.0007750981021672, 0.0008211878011934, 0.0008700182079338, 0.0009217521874234, 0.0009765623835847, 0.0010346318595111, 0.0010961542138830, 0.0011613349197432, 0.0012303915573284, 0.0013035543961450, 0.0013810677919537, 0.0014631903031841, 0.0015501962043345, 0.0016423756023869, 0.0017400364158675, 0.0018435043748468, 0.0019531247671694, 0.0020692637190223, 0.0021923084277660, 0.0023226698394865, 0.0024607831146568, 0.0026071087922901, 0.0027621355839074, 0.0029263808391988, 0.0031003924086690, 0.0032847514376044, 0.0034800728317350, 0.0036870087496936, 0.0039062497671694, 0.0041385274380445, 0.0043846168555319, 0.0046453396789730, 0.0049215662293136, 0.0052142175845802, 0.0055242711678147, 0.0058527616783977, 0.0062007848173380, 0.0065695028752089, 0.0069601456634700, 0.0073740174993873, 0.0078124995343387, 0.0082770548760891, 0.0087692337110639, 0.0092906802892685, 0.0098431324586272, 0.0104284351691604, 0.0110485423356295, 0.0117055233567953, 0.0124015696346760, 0.0131390057504177, 0.0139202913269401, 0.0147480349987745, 0.0156249990686774, 0.0165541097521782, 0.0175384692847729, 0.0185813605785370, 0.0196862649172544, 0.0208568722009659, 0.0220970865339041, 0.0234110467135906, 0.0248031392693520, 0.0262780115008354, 0.0278405826538801, 0.0294960699975491, 0.0312499981373549, 0.0331082195043564, 0.0350769385695457, 0.0371627211570740, 0.0393725298345089, 0.0417137444019318, 0.0441941730678082, 0.0468220934271812, 0.0496062822639942, 0.0525560230016708, 0.0556811690330505, 0.0589921437203884, 0.0624999962747097, 0.0662164390087128, 0.0701538771390915, 0.0743254423141479, 0.0787450596690178, 0.0834274888038635, 0.0883883461356163, 0.0936441868543625, 0.0992125645279884, 0.1051120460033417, 0.1113623380661011, 0.1179842874407768, 0.1249999925494194, 0.1324328780174255, 0.1403077542781830, 0.1486508846282959, 0.1574901193380356, 0.1668549776077271, 0.1767766922712326, 0.1872883737087250, 0.1984251290559769, 0.2102240920066833, 0.2227246761322021, 0.2359685748815536, 0.2500000000000000, 0.2648657560348511, 0.2806155085563660, 0.2973017692565918, 0.3149802684783936, 0.3337099552154541, 0.3535533845424652, 0.3745767772197723, 0.3968502581119537, 0.4204482138156891, 0.4454493522644043, 0.4719371497631073, 0.5000000000000000, 0.5297315716743469, 0.5612310171127319, 0.5946035385131836, 0.6299605369567871, 0.6674199104309082, 0.7071067690849304, 0.7491535544395447, 0.7937005162239075, 0.8408964276313782, 0.8908987045288086, 0.9438742995262146, 1.0000000000000000, 1.0594631433486938, 1.1224620342254639, 1.1892070770263672, 1.2599210739135742, 1.3348398208618164, 1.4142135381698608, 1.4983071088790894, 1.5874010324478149, 1.6817928552627563, 1.7817974090576172, 1.8877485990524292, 2.0000000000000000, 2.1189262866973877, 2.2449240684509277, 2.3784141540527344, 2.5198421478271484, 2.6696796417236328, 2.8284270763397217, 2.9966142177581787, 3.1748020648956299, 3.3635857105255127, 3.5635950565338135, 3.7754974365234375, 4.0000000000000000, 4.2378525733947754, 4.4898481369018555, 4.7568287849426270, 5.0396842956542969, 5.3393597602844238, 5.6568546295166016, 5.9932284355163574, 6.3496046066284180, 6.7271714210510254, 7.1271901130676270, 7.5509948730468750, 8.0000000000000000, 8.4757051467895508, 8.9796962738037109, 9.5136575698852539, 10.0793685913085938, 10.6787195205688477, 11.3137092590332031, 11.9864568710327148, 12.6992092132568359, 13.4543428421020508, 14.2543802261352539, 15.1019897460937500, 16.0000000000000000, 16.9514102935791016, 17.9593944549560547, 19.0273151397705078, 20.1587371826171875, 21.3574390411376953, 22.6274185180664062, 23.9729137420654297, 25.3984184265136719, 26.9086875915527344, 28.5087604522705078, 30.2039794921875000, 32.0000000000000000, 33.9028205871582031, 35.9187889099121094, 38.0546302795410156, 40.3174743652343750, 42.7148780822753906, 45.2548370361328125, 47.9458274841308594, 50.7968368530273438, 53.8173751831054688, 57.0175209045410156, 60.4079589843750000, 64.0000076293945312, 67.8056411743164062, 71.8375778198242188, 76.1092605590820312, 80.6349563598632812, 85.4297561645507812, 90.5096740722656250, 95.8916625976562500, 101.5936737060546875, 107.6347503662109375, 114.0350418090820312, 120.8159179687500000, 128.0000152587890625, 135.6112823486328125, 143.6751556396484375, 152.2185211181640625, 161.2699127197265625, 170.8595123291015625, 181.0193481445312500, 191.7833251953125000, 203.1873474121093750, 215.2695007324218750, 228.0700836181640625, 241.6318511962890625, 256.0000305175781250, 271.2225646972656250, 287.3503112792968750, 304.4370422363281250, 322.5398254394531250, 341.7190246582031250, 362.0386962890625000, 383.5666503906250000, 406.3746948242187500, 430.5390014648437500, 456.1401977539062500, 483.2637023925781250, 512.0000610351562500, 542.4451293945312500, 574.7006225585937500, 608.8740844726562500, 645.0796508789062500, 683.4380493164062500, 724.0773925781250000, 767.1333007812500000, 812.7494506835937500, 861.0780029296875000, 912.2803955078125000, 966.5274047851562500, 1024.0001220703125000, 1084.8903808593750000, 1149.4012451171875000, 1217.7481689453125000, 1290.1593017578125000, 1366.8762207031250000, 1448.1549072265625000, 1534.2666015625000000, 1625.4989013671875000};
maxiKick::maxiKick(){
maxiKick::envelope.setAttack(0);
maxiKick::setPitch(200);
maxiKick::envelope.setDecay(1);
maxiKick::envelope.setSustain(1);
maxiKick::envelope.setRelease(500);
maxiKick::envelope.holdtime=1;
maxiKick::envelope.trigger=0;
};
double maxiKick::play(){
envOut=envelope.adsr(1.,envelope.trigger);
if (inverse) {
envOut=fabs(1-envOut);
}
output=kick.sinewave(pitch*envOut)*envOut;
if (envelope.trigger==1) {
envelope.trigger=0;
}
if (useDistortion) {
output=distort.fastAtanDist(output, distortion);
}
if (useFilter) {
output=filter.lores(output, cutoff, resonance);
}
if (useLimiter) {
if (output*gain > 1) {
return 1.;
} else if (output*gain < -1) {
return -1.;
} else {
return output*gain;
}
} else {
return output*gain;
}
};
void maxiKick::setRelease(double release) {
envelope.setRelease(release);
}
void maxiKick::setPitch(double newPitch) {
pitch=newPitch;
}
void maxiKick::trigger() {
envelope.trigger=1;
}
maxiSnare::maxiSnare(){
maxiSnare::envelope.setAttack(0);
maxiSnare::setPitch(800);
maxiSnare::envelope.setDecay(20);
maxiSnare::envelope.setSustain(0.05);
maxiSnare::envelope.setRelease(300);
maxiSnare::envelope.holdtime=1;
maxiSnare::envelope.trigger=0;
};
double maxiSnare::play(){
envOut=envelope.adsr(1.,envelope.trigger);
if (inverse) {
envOut=fabs(1-envOut);
}
output=(tone.triangle(pitch*(0.1+(envOut*0.85)))+noise.noise())*envOut;
if (envelope.trigger==1) {
envelope.trigger=0;
}
if (useDistortion) {
output=distort.fastAtanDist(output, distortion);
}
if (useFilter) {
output=filter.lores(output, cutoff, resonance);
}
if (useLimiter) {
if (output*gain > 1) {
return 1.;
} else if (output*gain < -1) {
return -1.;
} else {
return output*gain;
}
} else {
return output*gain;
}
};
void maxiSnare::setRelease(double release) {
envelope.setRelease(release);
}
void maxiSnare::setPitch(double newPitch) {
pitch=newPitch;
}
void maxiSnare::trigger() {
envelope.trigger=1;
}
maxiHats::maxiHats(){
maxiHats::envelope.setAttack(0);
maxiHats::setPitch(12000);
maxiHats::envelope.setDecay(20);
maxiHats::envelope.setSustain(0.1);
maxiHats::envelope.setRelease(300);
maxiHats::envelope.holdtime=1;
maxiHats::envelope.trigger=0;
maxiHats::filter.setCutoff(8000);
maxiHats::filter.setResonance(1);
};
double maxiHats::play(){
envOut=envelope.adsr(1.,envelope.trigger);
if (inverse) {
envOut=fabs(1-envOut);
}
output=(tone.sinebuf(pitch)+noise.noise())*envOut;
if (envelope.trigger==1) {
envelope.trigger=0;
}
if (useDistortion) {
output=distort.fastAtanDist(output, distortion);
}
if (useFilter) {
output=filter.play(output, 0., 0., 1., 0.);
}
if (useLimiter) {
if (output*gain > 1) {
return 1.;
} else if (output*gain < -1) {
return -1.;
} else {
return output*gain;
}
} else {
return output*gain;
}
};
void maxiHats::setRelease(double release) {
envelope.setRelease(release);
}
void maxiHats::setPitch(double newPitch) {
pitch=newPitch;
}
void maxiHats::trigger() {
envelope.trigger=1;
}
maxiClock::maxiClock() {
playHead=0;
currentCount=0;
lastCount=0;
bpm=120;
ticks=1;
maxiClock::setTempo(bpm);
}
void maxiClock::ticker() {
tick=false;
currentCount=floor(timer.phasor(bps));//this sets up a metronome that ticks n times a second
if (lastCount!=currentCount) {//if we have a new timer int this sample,
tick=true;
playHead++;//iterate the playhead
}
}
void maxiClock::setTempo(double bpmIn) {
bpm=bpmIn;
bps=(bpm/60.)*ticks;
}
void maxiClock::setTicksPerBeat(int ticksPerBeat) {
ticks=ticksPerBeat;
maxiClock::setTempo(bpm);
}
maxiSampler::maxiSampler() {
maxiSampler::voices=32;
maxiSampler::currentVoice=0;
for (int i=0;i<voices;i++) {
maxiSampler::envelopes[i].setAttack(0);
maxiSampler::envelopes[i].setDecay(1);
maxiSampler::envelopes[i].setSustain(1.);
maxiSampler::envelopes[i].setRelease(2000);
maxiSampler::envelopes[i].holdtime=1;
maxiSampler::envelopes[i].trigger=0;
maxiSampler::envOut[i]=0;
maxiSampler::pitch[i]=0;
maxiSampler::outputs[i]=0;
}
}
void maxiSampler::setNumVoices(int numVoices) {
voices=numVoices;
}
double maxiSampler::play() {
output=0;
for (int i=0;i<voices;i++) {
envOut[i]=envelopes[i].adsr(envOutGain[i],envelopes[i].trigger);
if (envOut[i]>0.) {
outputs[i]=samples[i].play(pitchRatios[(int)pitch[i]+originalPitch]*((1./samples[i].length)*maxiSettings::sampleRate),0,samples[i].length)*envOut[i];
output+=outputs[i]/voices;
if (envelopes[i].trigger==1 && !sustain) {
envelopes[i].trigger=0;
}
}
} return output;
}
void maxiSampler::load(string inFile, bool setall) {
if (setall) {
for (int i=0;i<voices;i++) {
samples[i].load(inFile);
}
} else {
samples[currentVoice].load(inFile);
}
}
void maxiSampler::setPitch(double pitchIn, bool setall) {
if (setall) {
for (int i=0;i<voices;i++) {
pitch[i]=pitchIn;
}
} else {
pitch[currentVoice]=pitchIn;
}
}
void maxiSampler::midiNoteOn(double pitchIn, double velocity, bool setall) {
if (setall) {
for (int i=0;i<voices;i++) {
pitch[i]=pitchIn;
}
} else {
pitch[currentVoice]=pitchIn;
envOutGain[currentVoice]=velocity/128;
}
}
void maxiSampler::midiNoteOff(double pitchIn, double velocity, bool setall) {
for (int i=0;i<voices;i++){
if (pitch[i]==pitchIn) {
envelopes[i].trigger=0;
}
}
}
void maxiSampler::setAttack(double attackD, bool setall) {
if (setall) {
for (int i=0;i<voices;i++) {
envelopes[i].setAttack(attackD);
}
} else {
envelopes[currentVoice].setAttack(attackD);
}
}
void maxiSampler::setDecay(double decayD, bool setall) {
if (setall) {
for (int i=0;i<voices;i++) {
envelopes[i].setDecay(decayD);
}
} else {
envelopes[currentVoice].setDecay(decayD);
}
}
void maxiSampler::setSustain(double sustainD, bool setall) {
if (setall) {
for (int i=0;i<voices;i++) {
envelopes[i].setSustain(sustainD);
}
} else {
envelopes[currentVoice].setSustain(sustainD);
}
}
void maxiSampler::setRelease(double releaseD, bool setall) {
if (setall) {
for (int i=0;i<voices;i++) {
envelopes[i].setRelease(releaseD);
}
} else {
envelopes[currentVoice].setRelease(releaseD);
}
}
void maxiSampler::setPosition(double positionD, bool setall){
if (setall) {
for (int i=0;i<voices;i++) {
samples[i].setPosition(positionD);
}
} else {
samples[currentVoice].setPosition(positionD);
}
}
void maxiSampler::trigger() {
envelopes[currentVoice].trigger=1;
samples[currentVoice].trigger();
currentVoice++;
currentVoice=currentVoice%voices;
}
///*************************************************************
///
/// Init most variables
///
///*************************************************************
maxiRecorder::maxiRecorder() :
bufferSize(maxiSettings::sampleRate * 2),
bufferQueueSize(3),
bufferIndex(0),
recordedAmountFrames(0),
threadRunning(false)
{
}
///*************************************************************
///
/// Free resources in RAII manner and save to wav when this
/// object lifetime is up!
///
///*************************************************************
maxiRecorder::~maxiRecorder()
{
if (isRecording()) saveToWav();
freeResources();
}
///*************************************************************
///
/// Free resources only if they exist; may be called multiple
/// times in it's lifetime becasuse of destructor.
///
///*************************************************************
void maxiRecorder::freeResources()
{
if (savedBuffers.size() > 0)
{
for (int i = 0; i < savedBuffers.size(); ++i)
{
delete[] savedBuffers.front();
savedBuffers.pop();
}
}
}
///*************************************************************
///
/// Simply return if we are recording or not
///
///*************************************************************
bool maxiRecorder::isRecording() const
{
return doRecord;
}
///*************************************************************
///
/// Set the filename (and path)
///
///*************************************************************
void maxiRecorder::setup(std::string _filename)
{
filename = _filename;
}
///*************************************************************
///
/// You don't need to worry calling about this, this runs in a
/// detached thread after calling startRecording() and means
/// we can allocate new memory without interupting the real
/// time audio callback. Rather than have multiple ifdefs to
/// split up the same function to use different OS's threading
/// services I just redefined the thing because it looked
/// cleaner.
///
/// A pointer to the instanciated object from this class is
/// passed as '_context'. We do this because we are calling this
/// function from a static helper func and need to access object
/// specific variables.
///
///*************************************************************
#if defined(OS_IS_UNIX)
void* maxiRecorder::update(void* _context)
{
pthread_setcanceltype(PTHREAD_CANCEL_ASYNCHRONOUS, NULL);
maxiRecorder* _this = static_cast<maxiRecorder*>(_context);
_this->threadRunning = true;
while (_this->doRecord)
{
usleep((useconds_t)(10000. / bufferSize / maxiSettings::sampleRate));
while (_this->bufferQueueSize > _this->bufferQueue.size())
{
_this->enqueueBuffer();
}
}
std::cout << "\nFinishing Recording..." << std::endl;
return 0;
}
#elif defined(OS_IS_WIN)
void* maxiRecorder::update(void* _context)
{
maxiRecorder* _this = static_cast<maxiRecorder*>(_context);
_this->threadRunning = true;
long time = (long)(1000. / bufferSize / maxiSettings::sampleRate);
while (_this->doRecord)
{
Sleep(time);
while (_this->bufferQueueSize > _this->bufferQueue.size())
{
_this->enqueueBuffer();
}
}
std::cout << "\nFinishing Recording..." << std::endl;
return 0;
}
#endif
///*************************************************************
///
/// This whacks the update method into a detached thread which
/// will ensure there are enough buffers available in the queue.
///
///*************************************************************
void maxiRecorder::startRecording()
{
doRecord = true;
for (int i = 0; i < bufferQueueSize; ++i)
{
enqueueBuffer();
}
#if defined(OS_IS_UNIX)
if (pthread_create(&daemon, NULL, &maxiRecorder::update_pthread_helper, this) != 0)
std::cout << "pthread created incorrectly" << std::endl;
if (pthread_detach(daemon) != 0)
std::cout << "pthread detached incorrectly" << std::endl;
#elif defined(OS_IS_WIN)
unsigned threadID;
daemonHandle = (HANDLE) _beginthreadex(NULL, 0, &maxiRecorder::update_pthread_helper, (void *) this, 0, &threadID);
#endif
}
///*************************************************************
///
/// stops the recording thread by allowing the stack to unwind
/// in update() and stops the Windows thread if necessary
///
///*************************************************************
void maxiRecorder::stopRecording()
{
#if defined(OS_IS_WIN)
CloseHandle(daemonHandle);
#endif
doRecord = false;
}
///*************************************************************
///
/// Shamelessly hacked together largely using this link
/// http://stackoverflow.com/questions/22226872/two-problems-
/// when-writing-to-wav-c
///
/// I haven't bothered to change the syntax much here and the
/// templated write function was too convenient during my
/// struggle to get my head around wav headers so I've kept it
/// very similar to the attached code.
///
/// A clear optimisation would be to make getProcessedData()
/// return a vector of shorts to prevent another loop
///
///*************************************************************
template <typename T>
void maxiRecorder::write(std::ofstream& _stream, const T& _t) {
_stream.write((const char*)&_t, sizeof(T));
}
void maxiRecorder::saveToWav()
{
if (isRecording()) stopRecording();
if (threadRunning)
{
#if defined(OS_IS_UNIX)
pthread_cancel(daemon);
#endif
threadRunning = false;
}
std::vector<double> pcmData = getProcessedData();
std::vector<short> pcmDataInt;
pcmDataInt.resize(pcmData.size());
for (int i = 0; i < pcmData.size(); ++i)
pcmDataInt[i] = (short) (pcmData[i] * 3276.7);
int sampleRate = maxiSettings::sampleRate;
short channels = maxiSettings::channels;
int buffSize = int(pcmDataInt.size()) * 2;
std::ofstream stream(filename.c_str(), std::ios::binary);
if (stream.is_open())
{
/* Header */
stream.write("RIFF", 4);
write<int>(stream, 36 + buffSize);
stream.write("WAVE", 4);
/* Format Chunk */
stream.write("fmt ", 4);
write<int>(stream, 16);
write<short>(stream, 1);
write<short>(stream, channels);
write<int>(stream, sampleRate);
write<int>(stream, sampleRate * channels * sizeof(short));
write<short>(stream, channels * sizeof(short));
write<short>(stream, 8 * sizeof(short));
/* Data Chunk */
stream.write("data", 4);
stream.write((const char*) &buffSize, 4);
stream.write((const char*) pcmDataInt.data(), buffSize);
stream.close();
std::cout << "Wrote " << buffSize
<< " bytes to " << filename.c_str()
<< std::endl;
}
else
{
std::cerr << "Failed to open file: " << strerror(errno) << std::endl;
}
}
///*************************************************************
///
/// This takes the queue of user generated arrays and whacks
/// it into a vector of doubles which is easier to transverse
/// and work with.
///
///*************************************************************
std::vector<double> maxiRecorder::getProcessedData()
{
std::vector<double> userData;
int dataSize = int(savedBuffers.size()) * bufferSize;
userData.resize(dataSize);
int savedIndex = 0;
for (int i = 0; i < dataSize; ++i)
{
if (i > 0 && i % bufferSize == 0)
{
savedIndex = 0;
savedBuffers.pop();
}
userData[i] = savedBuffers.front()[savedIndex];
++savedIndex;
}
double lastFrame = userData[userData.size() - 1];
while (lastFrame == 0)
{
userData.pop_back();
lastFrame = userData[userData.size() - 1];
}
return userData;
}
///*************************************************************
///
/// Pass the buffer of audio to this method and it will write
/// it to the right place. This method is real-time safe and
/// is overriden for ofx's / port's floats array or maximilian's
/// / rtaudio's double
///
///*************************************************************
void maxiRecorder::passData(double* _in, int _inBufferSize)
{
for (int i = 0; i < _inBufferSize; ++i)
{
if (bufferIndex >= bufferSize)
{
bufferQueue.pop();
bufferIndex = 0;
}
bufferQueue.front()[bufferIndex] = _in[i];
++bufferIndex;
++recordedAmountFrames;
}
}
void maxiRecorder::passData(float* _in, int _inBufferSize)
{
for (int i = 0; i < _inBufferSize; ++i)
{
if (bufferIndex >= bufferSize)
{
bufferQueue.pop();
bufferIndex = 0;
}
bufferQueue.front()[bufferIndex] = (double) _in[i];
++bufferIndex;
++recordedAmountFrames;
}
}
///*************************************************************
///
/// The method to instanciate new memory and ensure the correct
/// structures have the correct addresses.
///
///*************************************************************
void maxiRecorder::enqueueBuffer()
{
double* b = new double[bufferSize];
bufferQueue.push(b);
savedBuffers.push(b);
}