/* * Compressed floating-point Twiddle Factor. * * This file implements the non-inline functions declared in * fpr.h, as well as the constants for FFT / iFFT. * * ============================================================================= * Copyright (c) 2023 by Cryptographic Engineering Research Group (CERG) * ECE Department, George Mason University * Fairfax, VA, U.S.A. * Author: Duc Tri Nguyen * Licensed under the Apache License, Version 2.0 (the "License"); * you may not use this file except in compliance with the License. * You may obtain a copy of the License at * http://www.apache.org/licenses/LICENSE-2.0 * Unless required by applicable law or agreed to in writing, software * distributed under the License is distributed on an "AS IS" BASIS, * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. * See the License for the specific language governing permissions and * limitations under the License. * ============================================================================= * @author Duc Tri Nguyen , */ #include "inner.h" const fpr fpr_p2_tab[] = { 2.00000000000, 1.00000000000, 0.50000000000, 0.25000000000, 0.12500000000, 0.06250000000, 0.03125000000, 0.01562500000, 0.00781250000, 0.00390625000, 0.00195312500 }; const fpr fpr_tab_log2[] = { 0.707106781186547524400844362, 0.707106781186547524400844362, // 4, 5 }; const fpr fpr_tab_log3[] = { 0.923879532511286756128183189, 0.382683432365089771728459984, // 8, 9 -0.382683432365089771728459984, 0.923879532511286756128183189, }; const fpr fpr_tab_log4[] = { 0.980785280403230449126182236, 0.195090322016128267848284868, // 16 0.555570233019602224742830814, 0.831469612302545237078788378, // 20 }; const fpr fpr_tab_log5[] = { 0.995184726672196886244836953, 0.098017140329560601994195564, // 32 0.634393284163645498215171613, 0.773010453362736960810906610, // 36 0.881921264348355029712756864, 0.471396736825997648556387626, // 40 0.290284677254462367636192376, 0.956940335732208864935797887, // 44 }; const fpr fpr_tab_log6[] = { 0.998795456205172392714771605, 0.049067674327418014254954977, // 64 0.671558954847018400625376850, 0.740951125354959091175616897, // 68 0.903989293123443331586200297, 0.427555093430282094320966857, // 72 0.336889853392220050689253213, 0.941544065183020778412509403, // 76 0.970031253194543992603984207, 0.242980179903263889948274162, // 80 0.514102744193221726593693839, 0.857728610000272069902269984, // 84 0.803207531480644909806676513, 0.595699304492433343467036529, // 88 0.146730474455361751658850130, 0.989176509964780973451673738, // 92 }; const fpr fpr_tab_log7[] = { 0.999698818696204220115765650, 0.024541228522912288031734529, // 128 0.689540544737066924616730630, 0.724247082951466920941069243, // 132 0.914209755703530654635014829, 0.405241314004989870908481306, // 136 0.359895036534988148775104572, 0.932992798834738887711660256, // 140 0.975702130038528544460395766, 0.219101240156869797227737547, // 144 0.534997619887097210663076905, 0.844853565249707073259571205, // 148 0.817584813151583696504920884, 0.575808191417845300745972454, // 152 0.170961888760301226363642357, 0.985277642388941244774018433, // 156 0.992479534598709998156767252, 0.122410675199216198498704474, // 160 0.615231590580626845484913563, 0.788346427626606262009164705, // 164 0.870086991108711418652292404, 0.492898192229784036873026689, // 168 0.266712757474898386325286515, 0.963776065795439866686464356, // 172 0.949528180593036667195936074, 0.313681740398891476656478846, // 176 0.449611329654606600046294579, 0.893224301195515320342416447, // 180 0.757208846506484547575464054, 0.653172842953776764084203014, // 184 0.073564563599667423529465622, 0.997290456678690216135597140, // 188 }; const fpr fpr_tab_log8[] = { 0.999924701839144540921646491, 0.012271538285719926079408262, // 256 0.698376249408972853554813503, 0.715730825283818654125532623, // 260 0.919113851690057743908477789, 0.393992040061048108596188661, // 264 0.371317193951837543411934967, 0.928506080473215565937167396, // 268 0.978317370719627633106240097, 0.207111376192218549708116020, // 272 0.545324988422046422313987347, 0.838224705554838043186996856, // 276 0.824589302785025264474803737, 0.565731810783613197389765011, // 280 0.183039887955140958516532578, 0.983105487431216327180301155, // 284 0.993906970002356041546922813, 0.110222207293883058807899140, // 288 0.624859488142386377084072816, 0.780737228572094478301588484, // 292 0.876070094195406607095844268, 0.482183772079122748517344481, // 296 0.278519689385053105207848526, 0.960430519415565811199035138, // 300 0.953306040354193836916740383, 0.302005949319228067003463232, // 304 0.460538710958240023633181487, 0.887639620402853947760181617, // 308 0.765167265622458925888815999, 0.643831542889791465068086063, // 312 0.085797312344439890461556332, 0.996312612182778012627226190, // 316 0.998118112900149207125155861, 0.061320736302208577782614593, // 320 0.662415777590171761113069817, 0.749136394523459325469203257, // 324 0.898674465693953843041976744, 0.438616238538527637647025738, // 328 0.325310292162262934135954708, 0.945607325380521325730945387, // 332 0.966976471044852109087220226, 0.254865659604514571553980779, // 336 0.503538383725717558691867071, 0.863972856121586737918147054, // 340 0.795836904608883536262791915, 0.605511041404325513920626941, // 344 0.134580708507126186316358409, 0.990902635427780025108237011, // 348 0.987301418157858382399815802, 0.158858143333861441684385360, // 352 0.585797857456438860328080838, 0.810457198252594791726703434, // 356 0.851355193105265142261290312, 0.524589682678468906215098464, // 360 0.231058108280671119643236018, 0.972939952205560145467720114, // 364 0.937339011912574923201899593, 0.348418680249434568419308588, // 368 0.416429560097637182562598911, 0.909167983090522376563884788, // 372 0.732654271672412834615546649, 0.680600997795453050594430464, // 376 0.036807222941358832324332691, 0.999322384588349500896221011, // 380 }; const fpr fpr_tab_log9[] = { 0.999981175282601142656990438, 0.006135884649154475359640235, // 512 0.702754744457225302452914421, 0.711432195745216441522130290, // 516 0.921514039342041943465396332, 0.388345046698826291624993541, // 520 0.377007410216418256726567823, 0.926210242138311341974793388, // 524 0.979569765685440534439326110, 0.201104634842091911558443546, // 528 0.550457972936604802977289893, 0.834862874986380056304401383, // 532 0.828045045257755752067527592, 0.560661576197336023839710223, // 536 0.189068664149806212754997837, 0.981963869109555264072848154, // 540 0.994564570734255452119106243, 0.104121633872054579120943880, // 544 0.629638238914927025372981341, 0.776888465673232450040827983, // 548 0.879012226428633477831323711, 0.476799230063322133342158117, // 552 0.284407537211271843618310615, 0.958703474895871555374645792, // 556 0.955141168305770721498157712, 0.296150888243623824121786128, // 560 0.465976495767966177902756065, 0.884797098430937780104007041, // 564 0.769103337645579639346626069, 0.639124444863775743801488193, // 568 0.091908956497132728624990979, 0.995767414467659793982495643, // 572 0.998475580573294752208559038, 0.055195244349689939809447526, // 576 0.666999922303637506650154222, 0.745057785441465962407907310, // 580 0.901348847046022014570746093, 0.433093818853151968484222638, // 584 0.331106305759876401737190737, 0.943593458161960361495301445, // 588 0.968522094274417316221088329, 0.248927605745720168110682816, // 592 0.508830142543107036931749324, 0.860866938637767279344583877, // 596 0.799537269107905033500246232, 0.600616479383868926653875896, // 600 0.140658239332849230714788846, 0.990058210262297105505906464, // 604 0.988257567730749491404792538, 0.152797185258443427720336613, // 608 0.590759701858874228423887908, 0.806847553543799272206514313, // 612 0.854557988365400520767862276, 0.519355990165589587361829932, // 616 0.237023605994367206867735915, 0.971503890986251775537099622, // 620 0.939459223602189911962669246, 0.342660717311994397592781983, // 624 0.422000270799799685941287941, 0.906595704514915365332960588, // 628 0.736816568877369875090132520, 0.676092703575315960360419228, // 632 0.042938256934940823077124540, 0.999077727752645382888781997, // 636 0.999529417501093163079703322, 0.030674803176636625934021028, // 640 0.685083667772700381362052545, 0.728464390448225196492035438, // 644 0.911706032005429851404397325, 0.410843171057903942183466675, // 648 0.354163525420490382357395796, 0.935183509938947577642207480, // 652 0.974339382785575860518721668, 0.225083911359792835991642120, // 656 0.529803624686294668216054671, 0.848120344803297251279133563, // 660 0.814036329705948361654516690, 0.580813958095764545075595272, // 664 0.164913120489969921418189113, 0.986308097244598647863297524, // 668 0.991709753669099522860049931, 0.128498110793793172624415589, // 672 0.610382806276309452716352152, 0.792106577300212351782342879, // 676 0.867046245515692651480195629, 0.498227666972781852410983869, // 680 0.260794117915275518280186509, 0.965394441697689374550843858, // 684 0.947585591017741134653387321, 0.319502030816015677901518272, // 688 0.444122144570429231642069418, 0.895966249756185155914560282, // 692 0.753186799043612482483430486, 0.657806693297078656931182264, // 696 0.067443919563664057897972422, 0.997723066644191609848546728, // 700 0.996820299291165714972629398, 0.079682437971430121147120656, // 704 0.648514401022112445084560551, 0.761202385484261814029709836, // 708 0.890448723244757889952150560, 0.455083587126343823535869268, // 712 0.307849640041534893682063646, 0.951435020969008369549175569, // 716 0.962121404269041595429604316, 0.272621355449948984493347477, // 720 0.487550160148435954641485027, 0.873094978418290098636085973, // 724 0.784556597155575233023892575, 0.620057211763289178646268191, // 728 0.116318630911904767252544319, 0.993211949234794533104601012, // 732 0.984210092386929073193874387, 0.177004220412148756196839844, // 736 0.570780745886967280232652864, 0.821102514991104679060430820, // 740 0.841554977436898409603499520, 0.540171472729892881297845480, // 744 0.213110319916091373967757518, 0.977028142657754351485866211, // 748 0.930766961078983731944872340, 0.365612997804773870011745909, // 752 0.399624199845646828544117031, 0.916679059921042663116457013, // 756 0.720002507961381629076682999, 0.693971460889654009003734389, // 760 0.018406729905804820927366313, 0.999830581795823422015722275, // 764 }; const fpr fpr_tab_log10[] = { 0.999995293809576171511580126, 0.003067956762965976270145365, // 1024 0.704934080375904908852523758, 0.709272826438865651316533772, // 1028 0.922701128333878570437264227, 0.385516053843918864075607949, // 1032 0.379847208924051170576281147, 0.925049240782677590302371869, // 1036 0.980182135968117392690210009, 0.198098410717953586179324918, // 1040 0.553016705580027531764226988, 0.833170164701913186439915922, // 1044 0.829761233794523042469023765, 0.558118531220556115693702964, // 1048 0.192080397049892441679288205, 0.981379193313754574318224190, // 1052 0.994879330794805620591166107, 0.101069862754827824987887585, // 1056 0.632018735939809021909403706, 0.774953106594873878359129282, // 1060 0.880470889052160770806542929, 0.474100214650550014398580015, // 1064 0.287347459544729526477331841, 0.957826413027532890321037029, // 1068 0.956045251349996443270479823, 0.293219162694258650606608599, // 1072 0.468688822035827933697617870, 0.883363338665731594736308015, // 1076 0.771060524261813773200605759, 0.636761861236284230413943435, // 1080 0.094963495329638998938034312, 0.995480755491926941769171600, // 1084 0.998640218180265222418199049, 0.052131704680283321236358216, // 1088 0.669282588346636065720696366, 0.743007952135121693517362293, // 1092 0.902673318237258806751502391, 0.430326481340082633908199031, // 1096 0.333999651442009404650865481, 0.942573197601446879280758735, // 1100 0.969281235356548486048290738, 0.245955050335794611599924709, // 1104 0.511468850437970399504391001, 0.859301818357008404783582139, // 1108 0.801376171723140219430247777, 0.598160706996342311724958652, // 1112 0.143695033150294454819773349, 0.989622017463200834623694454, // 1116 0.988721691960323767604516485, 0.149764534677321517229695737, // 1120 0.593232295039799808047809426, 0.805031331142963597922659282, // 1124 0.856147328375194481019630732, 0.516731799017649881508753876, // 1128 0.240003022448741486568922365, 0.970772140728950302138169611, // 1132 0.940506070593268323787291309, 0.339776884406826857828825803, // 1136 0.424779681209108833357226189, 0.905296759318118774354048329, // 1140 0.738887324460615147933116508, 0.673829000378756060917568372, // 1144 0.046003182130914628814301788, 0.998941293186856850633930266, // 1148 0.999618822495178597116830637, 0.027608145778965741612354872, // 1152 0.687315340891759108199186948, 0.726359155084345976817494315, // 1156 0.912962190428398164628018233, 0.408044162864978680820747499, // 1160 0.357030961233430032614954036, 0.934092550404258914729877883, // 1164 0.975025345066994146844913468, 0.222093620973203534094094721, // 1168 0.532403127877197971442805218, 0.846490938774052078300544488, // 1172 0.815814410806733789010772660, 0.578313796411655563342245019, // 1176 0.167938294974731178054745536, 0.985797509167567424700995000, // 1180 0.992099313142191757112085445, 0.125454983411546238542336453, // 1184 0.612810082429409703935211936, 0.790230221437310055030217152, // 1188 0.868570705971340895340449876, 0.495565261825772531150266670, // 1192 0.263754678974831383611349322, 0.964589793289812723836432159, // 1196 0.948561349915730288158494826, 0.316593375556165867243047035, // 1200 0.446868840162374195353044389, 0.894599485631382678433072126, // 1204 0.755201376896536527598710756, 0.655492852999615385312679701, // 1208 0.070504573389613863027351471, 0.997511456140303459699448390, // 1212 0.997060070339482978987989949, 0.076623861392031492278332463, // 1216 0.650846684996380915068975573, 0.759209188978388033485525443, // 1220 0.891840709392342727796478697, 0.452349587233770874133026703, // 1224 0.310767152749611495835997250, 0.950486073949481721759926101, // 1228 0.962953266873683886347921481, 0.269668325572915106525464462, // 1232 0.490226483288291154229598449, 0.871595086655951034842481435, // 1236 0.786455213599085757522319464, 0.617647307937803932403979402, // 1240 0.119365214810991364593637790, 0.992850414459865090793563344, // 1244 0.984748501801904218556553176, 0.173983873387463827950700807, // 1248 0.573297166698042212820171239, 0.819347520076796960824689637, // 1252 0.843208239641845437161743865, 0.537587076295645482502214932, // 1256 0.216106797076219509948385131, 0.976369731330021149312732194, // 1260 0.931884265581668106718557199, 0.362755724367397216204854462, // 1264 0.402434650859418441082533934, 0.915448716088267819566431292, // 1268 0.722128193929215321243607198, 0.691759258364157774906734132, // 1272 0.021474080275469507418374898, 0.999769405351215321657617036, // 1276 0.999882347454212525633049627, 0.015339206284988101044151868, // 1280 0.696177131491462944788582591, 0.717870045055731736211325329, // 1284 0.917900775621390457642276297, 0.396809987416710328595290911, // 1288 0.368466829953372331712746222, 0.929640895843181265457918066, // 1292 0.977677357824509979943404762, 0.210111836880469621717489972, // 1296 0.542750784864515906586768661, 0.839893794195999504583383987, // 1300 0.822849781375826332046780034, 0.568258952670131549790548489, // 1304 0.180022901405699522679906590, 0.983662419211730274396237776, // 1308 0.993564135520595333782021697, 0.113270952177564349018228733, // 1312 0.622461279374149972519166721, 0.782650596166575738458949301, // 1316 0.874586652278176112634431897, 0.484869248000791101822951699, // 1320 0.275571819310958163076425168, 0.961280485811320641748659653, // 1324 0.952375012719765858529893608, 0.304929229735402406490728633, // 1328 0.457813303598877221904961155, 0.889048355854664562540777729, // 1332 0.763188417263381271704838297, 0.646176012983316364832802220, // 1336 0.082740264549375693111987083, 0.996571145790554847093566910, // 1340 0.997925286198596012623025462, 0.064382630929857460819324537, // 1344 0.660114342067420478559490747, 0.751165131909686411205819422, // 1348 0.897324580705418281231391836, 0.441371268731716692879988968, // 1352 0.322407678801069848384807478, 0.946600913083283570044599823, // 1356 0.966190003445412555433832961, 0.257831102162159005614471295, // 1360 0.500885382611240786241285004, 0.865513624090569082825488358, // 1364 0.793975477554337164895083757, 0.607949784967773667243642671, // 1368 0.131540028702883111103387493, 0.991310859846115418957349799, // 1372 0.986809401814185476970235952, 0.161886393780111837641387995, // 1376 0.583308652937698294392830961, 0.812250586585203913049744181, // 1380 0.849741768000852489471268395, 0.527199134781901348464274575, // 1384 0.228072083170885739254457379, 0.973644249650811925318383912, // 1388 0.936265667170278246576310996, 0.351292756085567125601307623, // 1392 0.413638312238434547471944324, 0.910441292258067196934095369, // 1396 0.730562769227827561177758850, 0.682845546385248068164596123, // 1400 0.033741171851377584833716112, 0.999430604555461772019008327, // 1404 0.999204758618363895492950001, 0.039872927587739811128578738, // 1408 0.678350043129861486873655042, 0.734738878095963464563223604, // 1412 0.907886116487666212038681480, 0.419216888363223956433010020, // 1416 0.345541324963989065539191723, 0.938403534063108112192420774, // 1420 0.972226497078936305708321144, 0.234041958583543423191242045, // 1424 0.521975292937154342694258318, 0.852960604930363657746588082, // 1428 0.808656181588174991946968128, 0.588281548222645304786439813, // 1432 0.155828397654265235743101486, 0.987784141644572154230969032, // 1436 0.990485084256457037998682243, 0.137620121586486044948441663, // 1440 0.603066598540348201693430617, 0.797690840943391108362662755, // 1444 0.862423956111040538690933878, 0.506186645345155291048942344, // 1448 0.251897818154216950498106628, 0.967753837093475465243391912, // 1452 0.944604837261480265659265493, 0.328209843579092526107916817, // 1456 0.435857079922255491032544080, 0.900015892016160228714535267, // 1460 0.747100605980180144323078847, 0.664710978203344868130324985, // 1464 0.058258264500435759613979782, 0.998301544933892840738782163, // 1468 0.996044700901251989887944810, 0.088853552582524596561586535, // 1472 0.641481012808583151988739898, 0.767138911935820381181694573, // 1476 0.886222530148880631647990821, 0.463259783551860197390719637, // 1480 0.299079826308040476750336973, 0.954228095109105629780430732, // 1484 0.959571513081984528335528181, 0.281464937925757984095231007, // 1488 0.479493757660153026679839798, 0.877545290207261291668470750, // 1492 0.778816512381475953374724325, 0.627251815495144113509622565, // 1496 0.107172424956808849175529148, 0.994240449453187946358413442, // 1500 0.982539302287441255907040396, 0.186055151663446648105438304, // 1504 0.563199344013834115007363772, 0.826321062845663480311195452, // 1508 0.836547727223511984524285790, 0.547894059173100165608820571, // 1512 0.204108966092816874181696950, 0.978948175319062194715480124, // 1516 0.927362525650401087274536959, 0.374164062971457997104393020, // 1520 0.391170384302253888687512949, 0.920318276709110566440076541, // 1524 0.713584868780793592903125099, 0.700568793943248366792866380, // 1528 0.009203754782059819315102378, 0.999957644551963866333120920, // 1532 }; const fpr *fpr_table[] = { NULL, NULL, fpr_tab_log2, fpr_tab_log3, fpr_tab_log4, fpr_tab_log5, fpr_tab_log6, fpr_tab_log7, fpr_tab_log8, fpr_tab_log9, fpr_tab_log10, };