/* Amplitude scaling factors */ #define M17_0DB 0.1414213562 #define M6_99DB 0.4472135955 #define M4_77DB 0.5773502692 #define M3_01DB 0.7071067812 #define M1_76DB 0.8164965809 #define DELAYBUFLEN 1024 /* Length of the delay buffer */ #define HRTFFILTLEN 64 /* HRTF filter length */ #define IRTHRESH 0.001 /* Impulse response pruning thresh. */ #define AMPLNORM M6_99DB /* Overall amplitude renormalization */ #define BASSFILTFREQ 180 /* Bass compensation filter cut (Hz) */ #define BASSFILTLEN 193 /* Bass compensation filter length */ #define BASSGAIN M4_77DB /* Bass compensation gain */ #define BASSCROSS 0.35 /* Bass cross talk */ #define CFECHODELAY 360 /* Center front echo delay (samples) */ #define CFECHOAMPL M17_0DB /* Center front echo amplitude */ #define STEXPAND2 0.07 /* Stereo expansion / 2 */ /* Head related impulse response (HRIR) derived from KEMAR measurement data by Bill Gardner and Keith Martin URL: http://sound.media.mit.edu/KEMAR.html Distributed under GPL with authors' permission */ /* EQUALIZED KEMAR HRIR How to generate these data: 1. You need the MIT Media Lab's KEMAR data, read it into an software capable of signal/time series analysis (like Mathematica, Matlab, ...) 2. Construct an equalizing FIR inverse filter by calculating the transfer function of the front, same side impulse response, then take 1 over the absolute magnitude. 3. Cut the poles in the inverse filter's transfer function specification by limiting the amplification to 2.5 (note, this number assumes that you have correct signal processing normalization of the Fourier transform). 4. Design the FIR inverse filter by calculating the inverse Fourier transform, then chopping the coefficients down to a reasonable number (N = 15 is used here), apply the Kaiser-Bessel window (alpha = 2 is used here). Note the objective is remove the color bias only (as if you are using an equalizer), _not_ to do a full inverse filtering. (Note: beta = pi*alpha in other notation.) For N = 15, alpha = 2, you should get the following impulse response: 0.001001558668605168, 0.00698457265741865, 0.040453643039829436, 0.012230541722147855, -0.11939760844854072, 0.16468099899755967, -0.30297563073747436, 1.3140211791355982, -0.30297563073747436, 0.16468099899755967, -0.11939760844854072, 0.012230541722147855, 0.040453643039829436, 0.00698457265741865, 0.001001558668605168 5. Linearly convolve all KEMAR HRIR with this inverse filter. 6. Resample from 44.1 kHz sampling frequency to 48 kHz. */ /* Center front (-5 degree) - not 0 degree in order to create a clear front image from a finite distance */ float cf_filt[128] = { -0.00008638082319075036, 0.0003198059946385229, -0.0005010631339162132, 0.0011424741331126876, -0.001584220794688753, 0.001742715363246275, -0.0011080796626780694, 0.0001651829990860167, 0.005235028337314985, 0.0035223828473357776, 0.010057681388431303, -0.033469432129545514, 0.013391253316233523, 0.004858462839827063, 0.08172161220103627, 0.26158596134500023, -0.12420314583323326, -0.14298458356097565, 0.14421897280453896, -0.1635792507629016, -0.02187136722480014, 0.2426863044711817, 0.07229814207917194, 0.0942742651913879, 0.29856830878076834, 0.2944146162057754, -0.12122157003421209, -0.19640092165631157, 0.11623836502034968, -0.05794027397995521, -0.34313138341973776, -0.19487516249168105, 0.010118993953802401, -0.09460218797710966, -0.16761521117359582, 0.004154461610153861, 0.052768641758969316, -0.00041823982226147407, 0.021634960445143514, 0.07562793486871108, 0.08895407129506479, 0.039857755093416214, 0.044257936180422945, 0.061557584906101664, 0.015547268541895703, -0.023908191934932484, 0.009498030443468223, 0.03816269003221337, -0.009820500607303615, -0.042003975527908084, -0.03335447117311547, -0.029294510859746596, -0.05212623136198511, -0.073427547153205, -0.061190797824958836, -0.04620925059966413, -0.04204619420702159, -0.02331915902615157, 0.00211481411477094, 0.00852563995740107, 0.008766809731743399, 0.008666632180812078, 0.018279202191625352, 0.02924332751289675, 0.022293148257836494, 0.012362146008584188, 0.008572582458807008, 0.006491370763597344, -0.0019366944997535774, -0.006318669309634434, -0.006457921690218669, -0.015050265524669974, -0.02110660282616213, -0.017027809096377904, -0.01651052305334348, -0.022770064150046673, -0.01999875754219472, -0.012294792027337775, -0.011506057031057188, -0.011448970577312903, -0.004823572302580925, 0.0022451134042777883, 0.004145473526859826, 0.005629030064546135, 0.008588029213398785, 0.010092048834844231, 0.007182013245552008, 0.0014600979508720656, -0.0038314646272511756, -0.003443901997881347, -0.0029483418254804047, -0.007609357112679647, -0.006518368948030822, -0.004495803701497202, -0.007109113004849672, -0.008346237278084265, -0.005560847336252453, -0.002993453167553188, -0.005122897816824269, -0.004389782626604215, -0.0010912633695218108, -0.0019712029474458835, -0.005870162265802235, -0.005626159534954128, -0.00027254977910844407, 0.0013794425431633785, -0.0005919083190430062, -0.0007861203545416682, -0.0007049560240893946, -0.0032720188494468868, -0.004460645567968504, -0.0032018528193571696, -0.0030579229375062105, -0.003593998902656612, -0.0038032977997776445, -0.004613776010454773, -0.0059796549143736845, -0.00420126194319768, -0.0012374419948287222, 0.0008572699213050608, 0.0021490971020081094, 0.00047295283198381995, -0.0009670277915884887, -0.001354440866080231, -0.002962902746547851, -0.00533935813338141, -0.005469203016468759, -0.004355784273189485 }; /* Front, same side (30 degree) */ float af_filt[128] = { -0.004140580614755493, 0.005790934614385445, 0.003318916682081112, 0.014257145544366063, 0.007328442487127339, -0.06550381777876194, 0.03502225818161845, -0.013090579770708259, 0.2517776798694195, 0.420770489950659, -0.3418854608834852, -0.24371032493696737, 0.04901356150030018, -0.1778083521632833, 0.26448004245714163, 0.23245199964546834, 0.033053145803936305, 0.46811222821062415, 0.5359265986255611, -0.011912195468533787, -0.39763432601411647, -0.034482864386898314, 0.029445398240649626, -0.3850940407863262, -0.3272328478175581, -0.14701421403616477, -0.08522137400169517, -0.14936851633336035, -0.09432605283433723, 0.0991200405937827, 0.011075012089917331, -0.0051036489980274685, 0.0933903289749412, 0.1344189369609565, 0.10507466913017807, 0.04240159324684365, 0.06559270110638656, 0.026984119875617524, -0.03359846103238096, -0.018000197099174275, 0.042031818548178244, 0.03849039666888434, -0.02450829674011345, -0.03407882403088576, -0.029230189282961977, -0.046964865291761734, -0.09458258700116245, -0.1527349330901158, -0.15411577687826097, -0.08761679790956928, -0.033623549089171874, -0.007204768531481949, 0.008273425020444852, 0.021368717994908505, 0.04366608267875025, 0.05660907333076205, 0.06775726495503939, 0.05621881735271431, 0.03576231950669927, 0.02500825834889175, 0.015423811076054601, 0.007903258334503761, -0.0053873014137761945, -0.006987955469434698, -0.012027972007598602, -0.025228281243816594, -0.026225091797257318, -0.023809293997344882, -0.03250172017712295, -0.03195696301067249, -0.01784813952189948, -0.01663267233760342, -0.016184530450468065, -0.011659883749357463, -0.0035378511240219163, -0.0005485800790443406, 0.0018432660108168625, 0.011634844139907534, 0.018333603402051105, 0.020447379185133056, 0.00850783664147828, 0.0004694148911037838, -0.0017047130409786676, -0.0022409152834483997, -0.000860472174892845, -0.004111075059198666, -0.003527843382056666, -0.009640160874903018, -0.01750044574231376, -0.015613389403672443, -0.010790028120953001, -0.0095313499955768, -0.007469721416726809, -0.0019186578145117315, -0.00014977322572890802, -0.0029803838028179728, -0.006520567233727221, 0.000035015132033882596, 0.009245098100543752, 0.009896930052308656, 0.008316744929565786, 0.004575207140193997, -0.0000647420103997081, -0.004502916832871627, -0.004225962213251224, -0.002886014126381486, -0.006416834142585976, -0.007156609995423569, -0.008840274447579472, -0.01441763751386817, -0.015435817484659574, -0.00924487254924743, -0.0021571721940235205, 0.0028540722992305453, 0.00273577475088536, -0.000788412365513805, -0.0032650029728365907, -0.003880217646231338, -0.0035302087299613778, -0.0038531436176586246, -0.0011921632190514074, -0.0020722967099011938, -0.004985351145629344, -0.0042375588844648735, -0.003030360463006021, -0.0014161075428041471, -0.0005083025643192044, 0.00035096963769606926 }; /* Front, opposite (-30 degree) */ float of_filt[128] = { -0.000013472538374193126, -0.00008048061877079751, 0.000043927265781258155, -0.000017931700794858892, -0.000034774602476112886, -0.00009576223008735474, 0.0001557797638630691, -0.00018742885883751094, 0.00026512448626705716, -0.0001451040203319678, -0.00008263233117758043, 0.0006486245853639179, -0.0010631408451846698, 0.0026571994100746143, 0.0014179177997092787, 0.0062326502956616256, -0.008194149324545333, -0.006568029415878379, 0.009538759710818582, 0.012309193558632693, 0.12336638055838955, 0.046164307101829005, -0.10228706407884815, 0.04047687260345798, -0.00296595313977046, -0.07057949208414134, 0.08172114840714612, 0.08736490764127891, 0.05105250431333021, 0.11627179512747428, 0.20717888490340705, 0.09375052213570291, -0.09784374168330194, -0.010493571845901443, 0.053131894303891716, -0.10157443971694806, -0.16264032634244974, -0.05402369511361273, -0.0274403608654217, -0.09860277022495063, -0.06841875821090282, -0.004434574400066223, -0.0005222661652743502, -0.006231881259827263, 0.014410397820340159, 0.04885649512730243, 0.04361962569042684, 0.03399214029009391, 0.04961073933475931, 0.04067325604132289, 0.007850647519227257, 0.004564440466905299, 0.02257107958021618, 0.008183791928884486, -0.014913479343180557, -0.018685938460856224, -0.01745737397226911, -0.02327177054233603, -0.03723048632685227, -0.044739390162299685, -0.042651220125613766, -0.03730017561004743, -0.029039465434276192, -0.01885087458914294, -0.01207127752277769, -0.006779800724164512, -0.001930416967444157, 0.000029454577995528385, 0.0013822760965755472, 0.0014799128583230202, 0.0002068200609199832, 0.0022254295286201083, 0.005143858159434566, 0.0018580542060917013, -0.0019426046325146259, -0.0014464042108543495, -0.0034430083560735582, -0.009692758426099499, -0.011840035292593485, -0.010716508855893968, -0.012939889036853034, -0.0121846427935653, -0.006198503315630782, -0.0023186723099380305, -0.002679872498314837, -0.003086020446226295, -0.0015709623347698936, -0.0008147490468332398, -0.0012384575726770983, -0.0005212877089109362, 0.0017707578744906142, 0.001324932723905786, -0.0017023653780617696, -0.0045108927752919425, -0.005422155613096912, -0.0039489323837623835, -0.005295995750506547, -0.00629706566356189, -0.004685732198036754, -0.0048076735568143, -0.005978864279217503, -0.005928999306332966, -0.004187703549017582, -0.003213999896976475, -0.0028068699816073414, -0.0010889703907593833, 0.0003276714243495386, -0.0013015007040186994, -0.003208050402434782, -0.0025115319088208545, -0.0013787553006401076, -0.0018279087370218635, -0.0025904836507747754, -0.002071221947222004, -0.0026424212922485594, -0.0039837031817577, -0.0041635566057380965, -0.004355223489150822, -0.004350395332709937, -0.0036693292471930263, -0.003386384394185026, -0.003972568655001128, -0.004332336840023821, -0.002648767912111827, -0.001384410080218114, -0.0011353792711849466, -0.0013726264946164232, -0.0020075119315034313 }; /* Rear, same side (135 degree) */ float ar_filt[128] = { 0.004573315040810066, 0.0013592578059426913, 0.01553271930902704, -0.0002356117224941317, -0.05746098219774702, 0.03430688963370445, 0.00808371687447385, 0.21893535841158596, 0.2984357591724814, -0.3302799746504719, -0.3194029149806245, 0.21633225051331056, 0.24371260938097083, -0.08843705549751085, 0.03939684701343366, 0.45386926431114494, 0.07599118140753386, -0.18114706160474578, 0.285640624686038, 0.4049515236666218, -0.05347890222071792, -0.31464359045319074, -0.1033502246468194, -0.04553593949283157, -0.1880747731157464, -0.13629090230626037, -0.10435789106123239, -0.19818232801888755, -0.16701805476330397, -0.022793111199284, 0.058049696762683685, 0.007048321372693906, -0.002966419183225961, 0.10140569697797763, 0.11648999956673124, 0.05218347182779882, 0.028427001212735392, 0.04151900310166159, -0.0006960604221423734, -0.05898623212226975, -0.03801934531173312, -0.029306970535287986, -0.04549125782835908, -0.0599222718506552, -0.058299618975430116, -0.03765579129720727, -0.03559302657499581, -0.020647901025903054, -0.005720957338744348, -0.0041915732688915545, -0.0011470880098346143, 0.008737404798553, 0.023444168098121512, 0.024204226042172663, 0.01894897166475026, 0.020807655257479588, 0.021570431128040954, 0.006800556178576289, -0.009000089216921362, -0.010969824547067934, -0.0033653428332822374, -0.012676936164668659, -0.026739938673413587, -0.023427869194287573, -0.023302007105117244, -0.023647155590533712, -0.021289317515613106, -0.009120487305069884, 0.0009251551667728967, 0.00004285344125653763, -0.00009042365479456271, 0.00022573242339446494, 0.00720168491586098, 0.007111875505402431, 0.003186514817683482, 0.00810087718334745, 0.012619557025922575, 0.007854726400013397, -0.0024013592881066267, -0.001452457473161119, -0.0025535188366093945, -0.012428911627809337, -0.013729251536694145, -0.0070099675146427344, -0.007165284278706593, -0.01639289295622301, -0.015831795079778305, -0.007305768485523729, -0.003608863157004021, -0.0032640528878698084, 0.0030901263998481944, 0.00749497566124848, 0.002515185532327241, 0.00004840875738621367, 0.0017596043486043966, 0.0046229941553338144, 0.0034259167322926096, 0.003707347634186093, 0.0035584806528586328, -0.0019078936035275198, -0.006057343815214898, -0.0069262470468817, -0.004345020728618624, -0.004177623355574794, -0.005373506556122508, -0.006624933928893836, -0.008679541408588839, -0.010718719681595322, -0.011392246979594496, -0.007893917064389902, -0.0027572935365832536, -0.00006064707149834412, -0.0012512537319656323, -0.0024501501002409786, -0.0022106788572895998, -0.00234124933370301, -0.0008953445167066823, 0.0005393670625637734, -0.00033175600142209297, -0.004023994309351289, -0.008655472335784443, -0.009899957354849682, -0.008664952919996412, -0.00553483124503576, -0.003735336089277662, -0.002754824348643885, -0.0026884314856593368, -0.004084181815125924 }; /* Rear, opposite (-135 degree) */ float or_filt[128] = { 0.0001220944028243897, -0.000021785381808441314, 5.823057988603169e-6, -0.00001217768176447613, -0.00006123604397345513, 5.574117262531134e-6, -0.00004935331914366778, 1.1771577934768211e-6, -0.000059236211621095756, 9.503536190497286e-6, -0.0001494445696103564, 0.00012248858284145305, -0.0000963975321456313, 6.017905197665205e-6, 0.00003353395360402643, -0.0001931511015359506, 0.0005113536523931485, -0.0005676652619386114, 0.0012057159755477467, 0.0009370492250339692, 0.004596472288877981, -0.0018831773384237068, -0.008208535225621212, 0.0038178646400751056, 0.008726517739105965, 0.06664363898418262, 0.06788684221502142, -0.04492315162807267, -0.04019906311255255, 0.026203059677375153, 0.013678129114847544, -0.014334962223993527, 0.010141709596167392, 0.11559131576945537, 0.1586402064538425, 0.059975334707967023, 0.004671725963777715, 0.031498678282775874, 0.014338626006524587, -0.014749719448472231, -0.02055508237941379, -0.05407690143992048, -0.07767559836886143, -0.05029091786216801, -0.030808335706574427, -0.03401958135442541, -0.030520368430288967, -0.014168302104259355, 0.011907621259989802, 0.014286081013069, 0.006377467879613449, 0.018546823568277478, 0.028297012771618273, 0.025222339408338186, 0.021931611353415138, 0.019708894333646355, 0.01729258494072014, 0.017468204169564034, 0.009729094845051928, -0.002976992018531901, -0.00956986166277019, -0.016125733548332074, -0.02934094241442545, -0.04133767871051455, -0.043536981145416466, -0.0385966307108608, -0.02784453599342459, -0.018995135307247116, -0.012849534096536747, -0.004437491064613308, 0.00028385411598204655, 0.003776874988516643, 0.008069432041547833, 0.008764754183751848, 0.008020908861878062, 0.006830351461360802, 0.002218330884267235, -0.0020478725582339444, -0.003997428121462543, -0.007066287373515421, -0.00940847412544698, -0.010938998446237963, -0.011775483016151306, -0.011391103919484287, -0.010586061195163017, -0.009842793078929053, -0.007753202010139829, -0.00569213732353025, -0.006506783349722073, -0.005346134281903736, -0.003913089814898934, -0.0036091443854759727, -0.0020328564301266287, 0.00017932870773467683, 0.0032779786679056357, 0.003969695813293966, 0.0020339334412434987, -0.00011345940675415259, -0.0018344103399567666, -0.003556764701666365, -0.004263523639408391, -0.002940568582022133, -0.0034341188272627556, -0.006023399920020824, -0.0077456903203677865, -0.007912219312377842, -0.00625202770436523, -0.00530785086116117, -0.005569722659634311, -0.004664448462594344, -0.0037747773914077747, -0.004175649656985592, -0.004659601521384289, -0.005008602967819641, -0.004730625079902729, -0.0034039554356604146, -0.0017110333873406587, -0.0006091938771510242, -0.0016051679050678297, -0.003312864664007262, -0.004505512715977288, -0.004152222189861692, -0.003218596419678823, -0.0027277806209877343, -0.001715005444317267, -0.0012589960071233749, -0.001852908777923165, -0.002540339553144362 }; /* Center rear (180 degree) */ float cr_filt[128] = { -0.00005989110716536726, -0.00022790291829128702, 0.0002659166098971966, -0.0003774772716776257, 0.0004540309551867803, -0.000420238187386368, 0.00025518536450885686, 0.00028285526288953955, -0.001016391007574093, 0.0028634984299063795, 0.0021574799687976045, 0.01035121276682072, -0.010481720917298163, -0.013197198495899292, 0.0031928225328717195, 0.02903137843618603, 0.1632429772511569, 0.1047487989875262, -0.10464685060623742, -0.09260196288035998, -0.007514241993554443, 0.013596249226741712, -0.019876166508450258, 0.1255626123599804, 0.3648170359521724, 0.19458249753223478, -0.04434070930031298, 0.046582528121935265, 0.09484427347230277, -0.03137795311969644, -0.10297437925363695, -0.09351091015917065, -0.1129521092162064, -0.14925322995658827, -0.1231466295584665, -0.06356719756705227, -0.05442277895126282, -0.07568433015661316, -0.023314932828602003, 0.04260950721318558, 0.02249026315598923, 0.02048195669571197, 0.05651342117268278, 0.05885038917623213, 0.03797102097397795, 0.011767394419953451, 0.00560502503429481, 0.005051125343961189, -0.012925933188033823, -0.023918884651306566, -0.013251659441678816, -0.010694772488866284, -0.03080486448617846, -0.03661278237783158, -0.0379227303416262, -0.042189005718490775, -0.026595666344596286, -0.009759025956801257, -0.002064986663513004, -0.002420117028098389, -0.006629991977552491, 0.004619970897631026, 0.019450642967537877, 0.0173521119057514, 0.017641425439988062, 0.02270029598048491, 0.018976431925275348, 0.009299852902290885, -0.001695039371619912, -0.00675162574265618, -0.009380968871003034, -0.011208396125485165, -0.01308640049201482, -0.0165636375633249, -0.022004099870933345, -0.025173458684139286, -0.016918759559175375, -0.00865150653575917, -0.006999929082792643, -0.005454830010518988, -0.0021129521131095317, 0.00018717090054046307, -0.0002864344788569993, 0.0017615225381095569, 0.006985907557802283, 0.010577308310476465, 0.006466104789306027, -0.0014988738575948326, -0.0039669755229277195, -0.0065156971200080235, -0.009343206924192169, -0.0076430644693577495, -0.004395214976600924, -0.003052735340422483, -0.007019103043066595, -0.00974109267696527, -0.007968015032797376, -0.007801513845528344, -0.007535748903681969, -0.003543341967287925, 0.0015083125553729722, 0.0023345972556147025, -0.0010043623069557037, -0.0025295765105203746, -0.0023701840891643634, -0.0005908186035024362, 0.0029826252289082847, 0.004829048542117764, 0.004488360022902081, 0.00002643748103005408, -0.0042100779212597295, -0.006170600558114495, -0.007267149164680168, -0.006825522903494639, -0.006899834372739123, -0.0073493916110062675, -0.009554351265163382, -0.011790297433830197, -0.010645796603734424, -0.0064661575394022106, -0.002026743466524137, -0.0004337034584909932, -0.0011172647031654614, -0.0017947816283674731, -0.00255615052036616, -0.0017721562881944813, -0.0002379619297227554, 0.0007130120121089036 };