libm-test.inc 158 KB

12345678910111213141516171819202122232425262728293031323334353637383940414243444546474849505152535455565758596061626364656667686970717273747576777879808182838485868788899091929394959697989910010110210310410510610710810911011111211311411511611711811912012112212312412512612712812913013113213313413513613713813914014114214314414514614714814915015115215315415515615715815916016116216316416516616716816917017117217317417517617717817918018118218318418518618718818919019119219319419519619719819920020120220320420520620720820921021121221321421521621721821922022122222322422522622722822923023123223323423523623723823924024124224324424524624724824925025125225325425525625725825926026126226326426526626726826927027127227327427527627727827928028128228328428528628728828929029129229329429529629729829930030130230330430530630730830931031131231331431531631731831932032132232332432532632732832933033133233333433533633733833934034134234334434534634734834935035135235335435535635735835936036136236336436536636736836937037137237337437537637737837938038138238338438538638738838939039139239339439539639739839940040140240340440540640740840941041141241341441541641741841942042142242342442542642742842943043143243343443543643743843944044144244344444544644744844945045145245345445545645745845946046146246346446546646746846947047147247347447547647747847948048148248348448548648748848949049149249349449549649749849950050150250350450550650750850951051151251351451551651751851952052152252352452552652752852953053153253353453553653753853954054154254354454554654754854955055155255355455555655755855956056156256356456556656756856957057157257357457557657757857958058158258358458558658758858959059159259359459559659759859960060160260360460560660760860961061161261361461561661761861962062162262362462562662762862963063163263363463563663763863964064164264364464564664764864965065165265365465565665765865966066166266366466566666766866967067167267367467567667767867968068168268368468568668768868969069169269369469569669769869970070170270370470570670770870971071171271371471571671771871972072172272372472572672772872973073173273373473573673773873974074174274374474574674774874975075175275375475575675775875976076176276376476576676776876977077177277377477577677777877978078178278378478578678778878979079179279379479579679779879980080180280380480580680780880981081181281381481581681781881982082182282382482582682782882983083183283383483583683783883984084184284384484584684784884985085185285385485585685785885986086186286386486586686786886987087187287387487587687787887988088188288388488588688788888989089189289389489589689789889990090190290390490590690790890991091191291391491591691791891992092192292392492592692792892993093193293393493593693793893994094194294394494594694794894995095195295395495595695795895996096196296396496596696796896997097197297397497597697797897998098198298398498598698798898999099199299399499599699799899910001001100210031004100510061007100810091010101110121013101410151016101710181019102010211022102310241025102610271028102910301031103210331034103510361037103810391040104110421043104410451046104710481049105010511052105310541055105610571058105910601061106210631064106510661067106810691070107110721073107410751076107710781079108010811082108310841085108610871088108910901091109210931094109510961097109810991100110111021103110411051106110711081109111011111112111311141115111611171118111911201121112211231124112511261127112811291130113111321133113411351136113711381139114011411142114311441145114611471148114911501151115211531154115511561157115811591160116111621163116411651166116711681169117011711172117311741175117611771178117911801181118211831184118511861187118811891190119111921193119411951196119711981199120012011202120312041205120612071208120912101211121212131214121512161217121812191220122112221223122412251226122712281229123012311232123312341235123612371238123912401241124212431244124512461247124812491250125112521253125412551256125712581259126012611262126312641265126612671268126912701271127212731274127512761277127812791280128112821283128412851286128712881289129012911292129312941295129612971298129913001301130213031304130513061307130813091310131113121313131413151316131713181319132013211322132313241325132613271328132913301331133213331334133513361337133813391340134113421343134413451346134713481349135013511352135313541355135613571358135913601361136213631364136513661367136813691370137113721373137413751376137713781379138013811382138313841385138613871388138913901391139213931394139513961397139813991400140114021403140414051406140714081409141014111412141314141415141614171418141914201421142214231424142514261427142814291430143114321433143414351436143714381439144014411442144314441445144614471448144914501451145214531454145514561457145814591460146114621463146414651466146714681469147014711472147314741475147614771478147914801481148214831484148514861487148814891490149114921493149414951496149714981499150015011502150315041505150615071508150915101511151215131514151515161517151815191520152115221523152415251526152715281529153015311532153315341535153615371538153915401541154215431544154515461547154815491550155115521553155415551556155715581559156015611562156315641565156615671568156915701571157215731574157515761577157815791580158115821583158415851586158715881589159015911592159315941595159615971598159916001601160216031604160516061607160816091610161116121613161416151616161716181619162016211622162316241625162616271628162916301631163216331634163516361637163816391640164116421643164416451646164716481649165016511652165316541655165616571658165916601661166216631664166516661667166816691670167116721673167416751676167716781679168016811682168316841685168616871688168916901691169216931694169516961697169816991700170117021703170417051706170717081709171017111712171317141715171617171718171917201721172217231724172517261727172817291730173117321733173417351736173717381739174017411742174317441745174617471748174917501751175217531754175517561757175817591760176117621763176417651766176717681769177017711772177317741775177617771778177917801781178217831784178517861787178817891790179117921793179417951796179717981799180018011802180318041805180618071808180918101811181218131814181518161817181818191820182118221823182418251826182718281829183018311832183318341835183618371838183918401841184218431844184518461847184818491850185118521853185418551856185718581859186018611862186318641865186618671868186918701871187218731874187518761877187818791880188118821883188418851886188718881889189018911892189318941895189618971898189919001901190219031904190519061907190819091910191119121913191419151916191719181919192019211922192319241925192619271928192919301931193219331934193519361937193819391940194119421943194419451946194719481949195019511952195319541955195619571958195919601961196219631964196519661967196819691970197119721973197419751976197719781979198019811982198319841985198619871988198919901991199219931994199519961997199819992000200120022003200420052006200720082009201020112012201320142015201620172018201920202021202220232024202520262027202820292030203120322033203420352036203720382039204020412042204320442045204620472048204920502051205220532054205520562057205820592060206120622063206420652066206720682069207020712072207320742075207620772078207920802081208220832084208520862087208820892090209120922093209420952096209720982099210021012102210321042105210621072108210921102111211221132114211521162117211821192120212121222123212421252126212721282129213021312132213321342135213621372138213921402141214221432144214521462147214821492150215121522153215421552156215721582159216021612162216321642165216621672168216921702171217221732174217521762177217821792180218121822183218421852186218721882189219021912192219321942195219621972198219922002201220222032204220522062207220822092210221122122213221422152216221722182219222022212222222322242225222622272228222922302231223222332234223522362237223822392240224122422243224422452246224722482249225022512252225322542255225622572258225922602261226222632264226522662267226822692270227122722273227422752276227722782279228022812282228322842285228622872288228922902291229222932294229522962297229822992300230123022303230423052306230723082309231023112312231323142315231623172318231923202321232223232324232523262327232823292330233123322333233423352336233723382339234023412342234323442345234623472348234923502351235223532354235523562357235823592360236123622363236423652366236723682369237023712372237323742375237623772378237923802381238223832384238523862387238823892390239123922393239423952396239723982399240024012402240324042405240624072408240924102411241224132414241524162417241824192420242124222423242424252426242724282429243024312432243324342435243624372438243924402441244224432444244524462447244824492450245124522453245424552456245724582459246024612462246324642465246624672468246924702471247224732474247524762477247824792480248124822483248424852486248724882489249024912492249324942495249624972498249925002501250225032504250525062507250825092510251125122513251425152516251725182519252025212522252325242525252625272528252925302531253225332534253525362537253825392540254125422543254425452546254725482549255025512552255325542555255625572558255925602561256225632564256525662567256825692570257125722573257425752576257725782579258025812582258325842585258625872588258925902591259225932594259525962597259825992600260126022603260426052606260726082609261026112612261326142615261626172618261926202621262226232624262526262627262826292630263126322633263426352636263726382639264026412642264326442645264626472648264926502651265226532654265526562657265826592660266126622663266426652666266726682669267026712672267326742675267626772678267926802681268226832684268526862687268826892690269126922693269426952696269726982699270027012702270327042705270627072708270927102711271227132714271527162717271827192720272127222723272427252726272727282729273027312732273327342735273627372738273927402741274227432744274527462747274827492750275127522753275427552756275727582759276027612762276327642765276627672768276927702771277227732774277527762777277827792780278127822783278427852786278727882789279027912792279327942795279627972798279928002801280228032804280528062807280828092810281128122813281428152816281728182819282028212822282328242825282628272828282928302831283228332834283528362837283828392840284128422843284428452846284728482849285028512852285328542855285628572858285928602861286228632864286528662867286828692870287128722873287428752876287728782879288028812882288328842885288628872888288928902891289228932894289528962897289828992900290129022903290429052906290729082909291029112912291329142915291629172918291929202921292229232924292529262927292829292930293129322933293429352936293729382939294029412942294329442945294629472948294929502951295229532954295529562957295829592960296129622963296429652966296729682969297029712972297329742975297629772978297929802981298229832984298529862987298829892990299129922993299429952996299729982999300030013002300330043005300630073008300930103011301230133014301530163017301830193020302130223023302430253026302730283029303030313032303330343035303630373038303930403041304230433044304530463047304830493050305130523053305430553056305730583059306030613062306330643065306630673068306930703071307230733074307530763077307830793080308130823083308430853086308730883089309030913092309330943095309630973098309931003101310231033104310531063107310831093110311131123113311431153116311731183119312031213122312331243125312631273128312931303131313231333134313531363137313831393140314131423143314431453146314731483149315031513152315331543155315631573158315931603161316231633164316531663167316831693170317131723173317431753176317731783179318031813182318331843185318631873188318931903191319231933194319531963197319831993200320132023203320432053206320732083209321032113212321332143215321632173218321932203221322232233224322532263227322832293230323132323233323432353236323732383239324032413242324332443245324632473248324932503251325232533254325532563257325832593260326132623263326432653266326732683269327032713272327332743275327632773278327932803281328232833284328532863287328832893290329132923293329432953296329732983299330033013302330333043305330633073308330933103311331233133314331533163317331833193320332133223323332433253326332733283329333033313332333333343335333633373338333933403341334233433344334533463347334833493350335133523353335433553356335733583359336033613362336333643365336633673368336933703371337233733374337533763377337833793380338133823383338433853386338733883389339033913392339333943395339633973398339934003401340234033404340534063407340834093410341134123413341434153416341734183419342034213422342334243425342634273428342934303431343234333434343534363437343834393440344134423443344434453446344734483449345034513452345334543455345634573458345934603461346234633464346534663467346834693470347134723473347434753476347734783479348034813482348334843485348634873488348934903491349234933494349534963497349834993500350135023503350435053506350735083509351035113512351335143515351635173518351935203521352235233524352535263527352835293530353135323533353435353536353735383539354035413542354335443545354635473548354935503551355235533554355535563557355835593560356135623563356435653566356735683569357035713572357335743575357635773578357935803581358235833584358535863587358835893590359135923593359435953596359735983599360036013602360336043605360636073608360936103611361236133614361536163617361836193620362136223623362436253626362736283629363036313632363336343635363636373638363936403641364236433644364536463647364836493650365136523653365436553656365736583659366036613662366336643665366636673668366936703671367236733674367536763677367836793680368136823683368436853686368736883689369036913692369336943695369636973698369937003701370237033704370537063707370837093710371137123713371437153716371737183719372037213722372337243725372637273728372937303731373237333734373537363737373837393740374137423743374437453746374737483749375037513752375337543755375637573758375937603761376237633764376537663767376837693770377137723773377437753776377737783779378037813782378337843785378637873788378937903791379237933794379537963797379837993800380138023803380438053806380738083809381038113812381338143815381638173818381938203821382238233824382538263827382838293830383138323833383438353836383738383839384038413842384338443845384638473848384938503851385238533854385538563857385838593860386138623863386438653866386738683869387038713872387338743875387638773878387938803881388238833884388538863887388838893890389138923893389438953896389738983899390039013902390339043905390639073908390939103911391239133914391539163917391839193920392139223923392439253926392739283929393039313932393339343935393639373938393939403941394239433944394539463947394839493950395139523953395439553956395739583959396039613962396339643965396639673968396939703971397239733974397539763977397839793980398139823983398439853986398739883989399039913992399339943995399639973998399940004001400240034004400540064007400840094010401140124013401440154016401740184019402040214022402340244025402640274028402940304031403240334034403540364037403840394040404140424043404440454046404740484049405040514052405340544055405640574058405940604061406240634064406540664067406840694070407140724073407440754076407740784079408040814082408340844085408640874088408940904091409240934094409540964097409840994100410141024103410441054106410741084109411041114112411341144115411641174118411941204121412241234124412541264127412841294130413141324133413441354136413741384139414041414142414341444145414641474148414941504151415241534154415541564157415841594160416141624163416441654166416741684169417041714172417341744175417641774178417941804181418241834184418541864187418841894190419141924193419441954196419741984199420042014202420342044205420642074208420942104211421242134214421542164217421842194220422142224223422442254226422742284229423042314232423342344235423642374238423942404241424242434244424542464247424842494250425142524253425442554256425742584259426042614262426342644265426642674268426942704271427242734274427542764277427842794280428142824283428442854286428742884289429042914292429342944295429642974298429943004301430243034304430543064307430843094310431143124313431443154316431743184319432043214322432343244325432643274328432943304331433243334334433543364337433843394340434143424343434443454346434743484349435043514352435343544355435643574358435943604361436243634364436543664367436843694370437143724373437443754376437743784379438043814382438343844385438643874388438943904391439243934394439543964397439843994400440144024403440444054406440744084409441044114412441344144415441644174418441944204421442244234424442544264427442844294430443144324433443444354436443744384439444044414442444344444445444644474448444944504451445244534454445544564457445844594460446144624463446444654466446744684469447044714472447344744475447644774478447944804481448244834484448544864487448844894490449144924493449444954496449744984499450045014502450345044505450645074508450945104511451245134514451545164517451845194520452145224523452445254526452745284529453045314532453345344535453645374538453945404541454245434544454545464547454845494550455145524553455445554556455745584559456045614562456345644565456645674568456945704571457245734574457545764577457845794580458145824583458445854586458745884589459045914592459345944595459645974598459946004601460246034604460546064607460846094610461146124613461446154616461746184619462046214622462346244625462646274628462946304631463246334634463546364637463846394640464146424643464446454646464746484649465046514652465346544655465646574658465946604661466246634664466546664667466846694670467146724673467446754676467746784679468046814682468346844685468646874688468946904691469246934694469546964697469846994700470147024703470447054706470747084709471047114712471347144715471647174718471947204721472247234724472547264727472847294730473147324733473447354736473747384739474047414742474347444745474647474748474947504751475247534754475547564757475847594760476147624763476447654766476747684769477047714772477347744775477647774778477947804781478247834784478547864787478847894790479147924793479447954796479747984799480048014802480348044805480648074808480948104811481248134814481548164817481848194820482148224823482448254826482748284829483048314832483348344835483648374838483948404841484248434844484548464847484848494850485148524853485448554856485748584859486048614862486348644865486648674868486948704871487248734874487548764877487848794880488148824883488448854886488748884889489048914892489348944895489648974898489949004901490249034904490549064907490849094910491149124913491449154916491749184919492049214922492349244925492649274928492949304931493249334934493549364937493849394940494149424943494449454946494749484949495049514952495349544955495649574958495949604961496249634964496549664967496849694970497149724973497449754976497749784979498049814982498349844985498649874988498949904991499249934994499549964997499849995000500150025003500450055006500750085009501050115012501350145015501650175018501950205021502250235024502550265027502850295030503150325033503450355036503750385039504050415042504350445045504650475048504950505051505250535054505550565057505850595060506150625063506450655066506750685069507050715072507350745075507650775078507950805081508250835084508550865087508850895090509150925093509450955096509750985099510051015102510351045105
  1. /* Copyright (C) 1997-2002, 2003, 2004, 2005, 2006 Free Software Foundation, Inc.
  2. This file is part of the GNU C Library.
  3. Contributed by Andreas Jaeger <aj@suse.de>, 1997.
  4. The GNU C Library is free software; you can redistribute it and/or
  5. modify it under the terms of the GNU Lesser General Public
  6. License as published by the Free Software Foundation; either
  7. version 2.1 of the License, or (at your option) any later version.
  8. The GNU C Library is distributed in the hope that it will be useful,
  9. but WITHOUT ANY WARRANTY; without even the implied warranty of
  10. MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
  11. Lesser General Public License for more details.
  12. You should have received a copy of the GNU Lesser General Public
  13. License along with the GNU C Library; if not, write to the Free
  14. Software Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA
  15. 02111-1307 USA. */
  16. /* Part of testsuite for libm.
  17. This file is processed by a perl script. The resulting file has to
  18. be included by a master file that defines:
  19. Macros:
  20. FUNC(function): converts general function name (like cos) to
  21. name with correct suffix (e.g. cosl or cosf)
  22. MATHCONST(x): like FUNC but for constants (e.g convert 0.0 to 0.0L)
  23. FLOAT: floating point type to test
  24. - TEST_MSG: informal message to be displayed
  25. CHOOSE(Clongdouble,Cdouble,Cfloat,Cinlinelongdouble,Cinlinedouble,Cinlinefloat):
  26. chooses one of the parameters as delta for testing
  27. equality
  28. PRINTF_EXPR Floating point conversion specification to print a variable
  29. of type FLOAT with printf. PRINTF_EXPR just contains
  30. the specifier, not the percent and width arguments,
  31. e.g. "f".
  32. PRINTF_XEXPR Like PRINTF_EXPR, but print in hexadecimal format.
  33. PRINTF_NEXPR Like PRINTF_EXPR, but print nice. */
  34. /* This testsuite has currently tests for:
  35. acos, acosh, asin, asinh, atan, atan2, atanh,
  36. cbrt, ceil, copysign, cos, cosh, erf, erfc, exp, exp10, exp2, expm1,
  37. fabs, fdim, floor, fma, fmax, fmin, fmod, fpclassify,
  38. frexp, gamma, hypot,
  39. ilogb, isfinite, isinf, isnan, isnormal,
  40. isless, islessequal, isgreater, isgreaterequal, islessgreater, isunordered,
  41. j0, j1, jn,
  42. ldexp, lgamma, log, log10, log1p, log2, logb,
  43. modf, nearbyint, nextafter,
  44. pow, remainder, remquo, rint, lrint, llrint,
  45. round, lround, llround,
  46. scalb, scalbn, scalbln, signbit, sin, sincos, sinh, sqrt, tan, tanh, tgamma, trunc,
  47. y0, y1, yn, significand
  48. and for the following complex math functions:
  49. cabs, cacos, cacosh, carg, casin, casinh, catan, catanh,
  50. ccos, ccosh, cexp, clog, cpow, cproj, csin, csinh, csqrt, ctan, ctanh.
  51. At the moment the following functions aren't tested:
  52. drem, nan
  53. Parameter handling is primitive in the moment:
  54. --verbose=[0..3] for different levels of output:
  55. 0: only error count
  56. 1: basic report on failed tests (default)
  57. 2: full report on all tests
  58. -v for full output (equals --verbose=3)
  59. -u for generation of an ULPs file
  60. */
  61. /* "Philosophy":
  62. This suite tests some aspects of the correct implementation of
  63. mathematical functions in libm. Some simple, specific parameters
  64. are tested for correctness but there's no exhaustive
  65. testing. Handling of specific inputs (e.g. infinity, not-a-number)
  66. is also tested. Correct handling of exceptions is checked
  67. against. These implemented tests should check all cases that are
  68. specified in ISO C99.
  69. Exception testing: At the moment only divide-by-zero and invalid
  70. exceptions are tested. Overflow/underflow and inexact exceptions
  71. aren't checked at the moment.
  72. NaN values: There exist signalling and quiet NaNs. This implementation
  73. only uses quiet NaN as parameter but does not differenciate
  74. between the two kinds of NaNs as result.
  75. Inline functions: Inlining functions should give an improvement in
  76. speed - but not in precission. The inlined functions return
  77. reasonable values for a reasonable range of input values. The
  78. result is not necessarily correct for all values and exceptions are
  79. not correctly raised in all cases. Problematic input and return
  80. values are infinity, not-a-number and minus zero. This suite
  81. therefore does not check these specific inputs and the exception
  82. handling for inlined mathematical functions - just the "reasonable"
  83. values are checked.
  84. Beware: The tests might fail for any of the following reasons:
  85. - Tests are wrong
  86. - Functions are wrong
  87. - Floating Point Unit not working properly
  88. - Compiler has errors
  89. With e.g. gcc 2.7.2.2 the test for cexp fails because of a compiler error.
  90. To Do: All parameter should be numbers that can be represented as
  91. exact floating point values. Currently some values cannot be
  92. represented exactly and therefore the result is not the expected
  93. result. For this we will use 36 digits so that numbers can be
  94. represented exactly. */
  95. #ifndef _GNU_SOURCE
  96. # define _GNU_SOURCE
  97. #endif
  98. #include "libm-test-ulps.h"
  99. #include <complex.h>
  100. #include <math.h>
  101. #include <float.h>
  102. #include "fenv.h"
  103. #include <limits.h>
  104. #include <errno.h>
  105. #include <stdlib.h>
  106. #include <stdio.h>
  107. #include <string.h>
  108. #include <getopt.h>
  109. /* Possible exceptions */
  110. #define NO_EXCEPTION 0x0
  111. #define INVALID_EXCEPTION 0x1
  112. #define DIVIDE_BY_ZERO_EXCEPTION 0x2
  113. /* The next flags signals that those exceptions are allowed but not required. */
  114. #define INVALID_EXCEPTION_OK 0x4
  115. #define DIVIDE_BY_ZERO_EXCEPTION_OK 0x8
  116. #define EXCEPTIONS_OK INVALID_EXCEPTION_OK+DIVIDE_BY_ZERO_EXCEPTION_OK
  117. /* Some special test flags, passed togther with exceptions. */
  118. #define IGNORE_ZERO_INF_SIGN 0x10
  119. /* Various constants (we must supply them precalculated for accuracy). */
  120. #define M_PI_6l .52359877559829887307710723054658383L
  121. #define M_E2l 7.389056098930650227230427460575008L
  122. #define M_E3l 20.085536923187667740928529654581719L
  123. #define M_2_SQRT_PIl 3.5449077018110320545963349666822903L /* 2 sqrt (M_PIl) */
  124. #define M_SQRT_PIl 1.7724538509055160272981674833411451L /* sqrt (M_PIl) */
  125. #define M_LOG_SQRT_PIl 0.57236494292470008707171367567652933L /* log(sqrt(M_PIl)) */
  126. #define M_LOG_2_SQRT_PIl 1.265512123484645396488945797134706L /* log(2*sqrt(M_PIl)) */
  127. #define M_PI_34l (M_PIl - M_PI_4l) /* 3*pi/4 */
  128. #define M_PI_34_LOG10El (M_PIl - M_PI_4l) * M_LOG10El
  129. #define M_PI2_LOG10El M_PI_2l * M_LOG10El
  130. #define M_PI4_LOG10El M_PI_4l * M_LOG10El
  131. #define M_PI_LOG10El M_PIl * M_LOG10El
  132. #define M_SQRT_2_2 0.70710678118654752440084436210484903L /* sqrt (2) / 2 */
  133. static FILE *ulps_file; /* File to document difference. */
  134. static int output_ulps; /* Should ulps printed? */
  135. static int noErrors; /* number of errors */
  136. static int noTests; /* number of tests (without testing exceptions) */
  137. static int noExcTests; /* number of tests for exception flags */
  138. static int noXFails; /* number of expected failures. */
  139. static int noXPasses; /* number of unexpected passes. */
  140. static int verbose;
  141. static int output_max_error; /* Should the maximal errors printed? */
  142. static int output_points; /* Should the single function results printed? */
  143. static int ignore_max_ulp; /* Should we ignore max_ulp? */
  144. static FLOAT minus_zero, plus_zero;
  145. static FLOAT plus_infty, minus_infty, nan_value, max_value, min_value;
  146. static FLOAT max_error, real_max_error, imag_max_error;
  147. #define BUILD_COMPLEX(real, imag) \
  148. ({ __complex__ FLOAT __retval; \
  149. __real__ __retval = (real); \
  150. __imag__ __retval = (imag); \
  151. __retval; })
  152. #define BUILD_COMPLEX_INT(real, imag) \
  153. ({ __complex__ int __retval; \
  154. __real__ __retval = (real); \
  155. __imag__ __retval = (imag); \
  156. __retval; })
  157. #define MANT_DIG CHOOSE ((LDBL_MANT_DIG-1), (DBL_MANT_DIG-1), (FLT_MANT_DIG-1), \
  158. (LDBL_MANT_DIG-1), (DBL_MANT_DIG-1), (FLT_MANT_DIG-1))
  159. static void
  160. init_max_error (void)
  161. {
  162. max_error = 0;
  163. real_max_error = 0;
  164. imag_max_error = 0;
  165. feclearexcept (FE_ALL_EXCEPT);
  166. }
  167. static void
  168. set_max_error (FLOAT current, FLOAT *curr_max_error)
  169. {
  170. if (current > *curr_max_error)
  171. *curr_max_error = current;
  172. }
  173. /* Should the message print to screen? This depends on the verbose flag,
  174. and the test status. */
  175. static int
  176. print_screen (int ok, int xfail)
  177. {
  178. if (output_points
  179. && (verbose > 1
  180. || (verbose == 1 && ok == xfail)))
  181. return 1;
  182. return 0;
  183. }
  184. /* Should the message print to screen? This depends on the verbose flag,
  185. and the test status. */
  186. static int
  187. print_screen_max_error (int ok, int xfail)
  188. {
  189. if (output_max_error
  190. && (verbose > 1
  191. || ((verbose == 1) && (ok == xfail))))
  192. return 1;
  193. return 0;
  194. }
  195. /* Update statistic counters. */
  196. static void
  197. update_stats (int ok, int xfail)
  198. {
  199. ++noTests;
  200. if (ok && xfail)
  201. ++noXPasses;
  202. else if (!ok && xfail)
  203. ++noXFails;
  204. else if (!ok && !xfail)
  205. ++noErrors;
  206. }
  207. static void
  208. print_ulps (const char *test_name, FLOAT ulp)
  209. {
  210. if (output_ulps)
  211. {
  212. fprintf (ulps_file, "Test \"%s\":\n", test_name);
  213. fprintf (ulps_file, "%s: %.0" PRINTF_NEXPR "\n",
  214. CHOOSE("ldouble", "double", "float",
  215. "ildouble", "idouble", "ifloat"),
  216. FUNC(ceil) (ulp));
  217. }
  218. }
  219. static void
  220. print_function_ulps (const char *function_name, FLOAT ulp)
  221. {
  222. if (output_ulps)
  223. {
  224. fprintf (ulps_file, "Function: \"%s\":\n", function_name);
  225. fprintf (ulps_file, "%s: %.0" PRINTF_NEXPR "\n",
  226. CHOOSE("ldouble", "double", "float",
  227. "ildouble", "idouble", "ifloat"),
  228. FUNC(ceil) (ulp));
  229. }
  230. }
  231. static void
  232. print_complex_function_ulps (const char *function_name, FLOAT real_ulp,
  233. FLOAT imag_ulp)
  234. {
  235. if (output_ulps)
  236. {
  237. if (real_ulp != 0.0)
  238. {
  239. fprintf (ulps_file, "Function: Real part of \"%s\":\n", function_name);
  240. fprintf (ulps_file, "%s: %.0" PRINTF_NEXPR "\n",
  241. CHOOSE("ldouble", "double", "float",
  242. "ildouble", "idouble", "ifloat"),
  243. FUNC(ceil) (real_ulp));
  244. }
  245. if (imag_ulp != 0.0)
  246. {
  247. fprintf (ulps_file, "Function: Imaginary part of \"%s\":\n", function_name);
  248. fprintf (ulps_file, "%s: %.0" PRINTF_NEXPR "\n",
  249. CHOOSE("ldouble", "double", "float",
  250. "ildouble", "idouble", "ifloat"),
  251. FUNC(ceil) (imag_ulp));
  252. }
  253. }
  254. }
  255. /* Test if Floating-Point stack hasn't changed */
  256. static void
  257. fpstack_test (const char *test_name)
  258. {
  259. #ifdef i386
  260. static int old_stack;
  261. int sw;
  262. __asm__ ("fnstsw" : "=a" (sw));
  263. sw >>= 11;
  264. sw &= 7;
  265. if (sw != old_stack)
  266. {
  267. printf ("FP-Stack wrong after test %s (%d, should be %d)\n",
  268. test_name, sw, old_stack);
  269. ++noErrors;
  270. old_stack = sw;
  271. }
  272. #endif
  273. }
  274. static void
  275. print_max_error (const char *func_name, FLOAT allowed, int xfail)
  276. {
  277. int ok = 0;
  278. if (max_error == 0.0 || (max_error <= allowed && !ignore_max_ulp))
  279. {
  280. ok = 1;
  281. }
  282. if (!ok)
  283. print_function_ulps (func_name, max_error);
  284. if (print_screen_max_error (ok, xfail))
  285. {
  286. printf ("Maximal error of `%s'\n", func_name);
  287. printf (" is : %.0" PRINTF_NEXPR " ulp\n", FUNC(ceil) (max_error));
  288. printf (" accepted: %.0" PRINTF_NEXPR " ulp\n", FUNC(ceil) (allowed));
  289. }
  290. update_stats (ok, xfail);
  291. }
  292. static void
  293. print_complex_max_error (const char *func_name, __complex__ FLOAT allowed,
  294. __complex__ int xfail)
  295. {
  296. int ok = 0;
  297. if ((real_max_error == 0 && imag_max_error == 0)
  298. || (real_max_error <= __real__ allowed
  299. && imag_max_error <= __imag__ allowed
  300. && !ignore_max_ulp))
  301. {
  302. ok = 1;
  303. }
  304. if (!ok)
  305. print_complex_function_ulps (func_name, real_max_error, imag_max_error);
  306. if (print_screen_max_error (ok, xfail))
  307. {
  308. printf ("Maximal error of real part of: %s\n", func_name);
  309. printf (" is : %.0" PRINTF_NEXPR " ulp\n",
  310. FUNC(ceil) (real_max_error));
  311. printf (" accepted: %.0" PRINTF_NEXPR " ulp\n",
  312. FUNC(ceil) (__real__ allowed));
  313. printf ("Maximal error of imaginary part of: %s\n", func_name);
  314. printf (" is : %.0" PRINTF_NEXPR " ulp\n",
  315. FUNC(ceil) (imag_max_error));
  316. printf (" accepted: %.0" PRINTF_NEXPR " ulp\n",
  317. FUNC(ceil) (__imag__ allowed));
  318. }
  319. update_stats (ok, xfail);
  320. }
  321. /* Test whether a given exception was raised. */
  322. static void
  323. test_single_exception (const char *test_name,
  324. int exception,
  325. int exc_flag,
  326. int fe_flag,
  327. const char *flag_name)
  328. {
  329. #ifndef TEST_INLINE
  330. int ok = 1;
  331. if (exception & exc_flag)
  332. {
  333. if (fetestexcept (fe_flag))
  334. {
  335. if (print_screen (1, 0))
  336. printf ("Pass: %s: Exception \"%s\" set\n", test_name, flag_name);
  337. }
  338. else
  339. {
  340. ok = 0;
  341. if (print_screen (0, 0))
  342. printf ("Failure: %s: Exception \"%s\" not set\n",
  343. test_name, flag_name);
  344. }
  345. }
  346. else
  347. {
  348. if (fetestexcept (fe_flag))
  349. {
  350. ok = 0;
  351. if (print_screen (0, 0))
  352. printf ("Failure: %s: Exception \"%s\" set\n",
  353. test_name, flag_name);
  354. }
  355. else
  356. {
  357. if (print_screen (1, 0))
  358. printf ("%s: Exception \"%s\" not set\n", test_name,
  359. flag_name);
  360. }
  361. }
  362. if (!ok)
  363. ++noErrors;
  364. #endif
  365. }
  366. /* Test whether exceptions given by EXCEPTION are raised. Ignore thereby
  367. allowed but not required exceptions.
  368. */
  369. static void
  370. test_exceptions (const char *test_name, int exception)
  371. {
  372. ++noExcTests;
  373. #ifdef FE_DIVBYZERO
  374. if ((exception & DIVIDE_BY_ZERO_EXCEPTION_OK) == 0)
  375. test_single_exception (test_name, exception,
  376. DIVIDE_BY_ZERO_EXCEPTION, FE_DIVBYZERO,
  377. "Divide by zero");
  378. #endif
  379. #ifdef FE_INVALID
  380. if ((exception & INVALID_EXCEPTION_OK) == 0)
  381. test_single_exception (test_name, exception, INVALID_EXCEPTION, FE_INVALID,
  382. "Invalid operation");
  383. #endif
  384. feclearexcept (FE_ALL_EXCEPT);
  385. }
  386. static void
  387. check_float_internal (const char *test_name, FLOAT computed, FLOAT expected,
  388. FLOAT max_ulp, int xfail, int exceptions,
  389. FLOAT *curr_max_error)
  390. {
  391. int ok = 0;
  392. int print_diff = 0;
  393. FLOAT diff = 0;
  394. FLOAT ulp = 0;
  395. test_exceptions (test_name, exceptions);
  396. if (isnan (computed) && isnan (expected))
  397. ok = 1;
  398. else if (isinf (computed) && isinf (expected))
  399. {
  400. /* Test for sign of infinities. */
  401. if ((exceptions & IGNORE_ZERO_INF_SIGN) == 0
  402. && signbit (computed) != signbit (expected))
  403. {
  404. ok = 0;
  405. printf ("infinity has wrong sign.\n");
  406. }
  407. else
  408. ok = 1;
  409. }
  410. /* Don't calc ulp for NaNs or infinities. */
  411. else if (isinf (computed) || isnan (computed) || isinf (expected) || isnan (expected))
  412. ok = 0;
  413. else
  414. {
  415. diff = FUNC(fabs) (computed - expected);
  416. /* ilogb (0) isn't allowed. */
  417. if (expected == 0.0)
  418. ulp = diff / FUNC(ldexp) (1.0, - MANT_DIG);
  419. else
  420. ulp = diff / FUNC(ldexp) (1.0, FUNC(ilogb) (expected) - MANT_DIG);
  421. set_max_error (ulp, curr_max_error);
  422. print_diff = 1;
  423. if ((exceptions & IGNORE_ZERO_INF_SIGN) == 0
  424. && computed == 0.0 && expected == 0.0
  425. && signbit(computed) != signbit (expected))
  426. ok = 0;
  427. else if (ulp <= 0.5 || (ulp <= max_ulp && !ignore_max_ulp))
  428. ok = 1;
  429. else
  430. {
  431. ok = 0;
  432. print_ulps (test_name, ulp);
  433. }
  434. }
  435. if (print_screen (ok, xfail))
  436. {
  437. if (!ok)
  438. printf ("Failure: ");
  439. printf ("Test: %s\n", test_name);
  440. printf ("Result:\n");
  441. printf (" is: % .20" PRINTF_EXPR " % .20" PRINTF_XEXPR "\n",
  442. computed, computed);
  443. printf (" should be: % .20" PRINTF_EXPR " % .20" PRINTF_XEXPR "\n",
  444. expected, expected);
  445. if (print_diff)
  446. {
  447. printf (" difference: % .20" PRINTF_EXPR " % .20" PRINTF_XEXPR
  448. "\n", diff, diff);
  449. printf (" ulp : % .4" PRINTF_NEXPR "\n", ulp);
  450. printf (" max.ulp : % .4" PRINTF_NEXPR "\n", max_ulp);
  451. }
  452. }
  453. update_stats (ok, xfail);
  454. fpstack_test (test_name);
  455. }
  456. static void
  457. check_float (const char *test_name, FLOAT computed, FLOAT expected,
  458. FLOAT max_ulp, int xfail, int exceptions)
  459. {
  460. check_float_internal (test_name, computed, expected, max_ulp, xfail,
  461. exceptions, &max_error);
  462. }
  463. static void
  464. check_complex (const char *test_name, __complex__ FLOAT computed,
  465. __complex__ FLOAT expected,
  466. __complex__ FLOAT max_ulp, __complex__ int xfail,
  467. int exception)
  468. {
  469. FLOAT part_comp, part_exp, part_max_ulp;
  470. int part_xfail;
  471. char str[200];
  472. sprintf (str, "Real part of: %s", test_name);
  473. part_comp = __real__ computed;
  474. part_exp = __real__ expected;
  475. part_max_ulp = __real__ max_ulp;
  476. part_xfail = __real__ xfail;
  477. check_float_internal (str, part_comp, part_exp, part_max_ulp, part_xfail,
  478. exception, &real_max_error);
  479. sprintf (str, "Imaginary part of: %s", test_name);
  480. part_comp = __imag__ computed;
  481. part_exp = __imag__ expected;
  482. part_max_ulp = __imag__ max_ulp;
  483. part_xfail = __imag__ xfail;
  484. /* Don't check again for exceptions, just pass through the
  485. zero/inf sign test. */
  486. check_float_internal (str, part_comp, part_exp, part_max_ulp, part_xfail,
  487. exception & IGNORE_ZERO_INF_SIGN,
  488. &imag_max_error);
  489. }
  490. /* Check that computed and expected values are equal (int values). */
  491. static void
  492. check_int (const char *test_name, int computed, int expected, int max_ulp,
  493. int xfail, int exceptions)
  494. {
  495. int diff = computed - expected;
  496. int ok = 0;
  497. test_exceptions (test_name, exceptions);
  498. noTests++;
  499. if (abs (diff) <= max_ulp)
  500. ok = 1;
  501. if (!ok)
  502. print_ulps (test_name, diff);
  503. if (print_screen (ok, xfail))
  504. {
  505. if (!ok)
  506. printf ("Failure: ");
  507. printf ("Test: %s\n", test_name);
  508. printf ("Result:\n");
  509. printf (" is: %d\n", computed);
  510. printf (" should be: %d\n", expected);
  511. }
  512. update_stats (ok, xfail);
  513. fpstack_test (test_name);
  514. }
  515. /* Check that computed and expected values are equal (long int values). */
  516. static void
  517. check_long (const char *test_name, long int computed, long int expected,
  518. long int max_ulp, int xfail, int exceptions)
  519. {
  520. long int diff = computed - expected;
  521. int ok = 0;
  522. test_exceptions (test_name, exceptions);
  523. noTests++;
  524. if (labs (diff) <= max_ulp)
  525. ok = 1;
  526. if (!ok)
  527. print_ulps (test_name, diff);
  528. if (print_screen (ok, xfail))
  529. {
  530. if (!ok)
  531. printf ("Failure: ");
  532. printf ("Test: %s\n", test_name);
  533. printf ("Result:\n");
  534. printf (" is: %ld\n", computed);
  535. printf (" should be: %ld\n", expected);
  536. }
  537. update_stats (ok, xfail);
  538. fpstack_test (test_name);
  539. }
  540. /* Check that computed value is true/false. */
  541. static void
  542. check_bool (const char *test_name, int computed, int expected,
  543. long int max_ulp, int xfail, int exceptions)
  544. {
  545. int ok = 0;
  546. test_exceptions (test_name, exceptions);
  547. noTests++;
  548. if ((computed == 0) == (expected == 0))
  549. ok = 1;
  550. if (print_screen (ok, xfail))
  551. {
  552. if (!ok)
  553. printf ("Failure: ");
  554. printf ("Test: %s\n", test_name);
  555. printf ("Result:\n");
  556. printf (" is: %d\n", computed);
  557. printf (" should be: %d\n", expected);
  558. }
  559. update_stats (ok, xfail);
  560. fpstack_test (test_name);
  561. }
  562. /* check that computed and expected values are equal (long int values) */
  563. static void
  564. check_longlong (const char *test_name, long long int computed,
  565. long long int expected,
  566. long long int max_ulp, int xfail,
  567. int exceptions)
  568. {
  569. long long int diff = computed - expected;
  570. int ok = 0;
  571. test_exceptions (test_name, exceptions);
  572. noTests++;
  573. if (llabs (diff) <= max_ulp)
  574. ok = 1;
  575. if (!ok)
  576. print_ulps (test_name, diff);
  577. if (print_screen (ok, xfail))
  578. {
  579. if (!ok)
  580. printf ("Failure:");
  581. printf ("Test: %s\n", test_name);
  582. printf ("Result:\n");
  583. printf (" is: %lld\n", computed);
  584. printf (" should be: %lld\n", expected);
  585. }
  586. update_stats (ok, xfail);
  587. fpstack_test (test_name);
  588. }
  589. /* This is to prevent messages from the SVID libm emulation. */
  590. int
  591. matherr (struct exception *x __attribute__ ((unused)))
  592. {
  593. return 1;
  594. }
  595. /****************************************************************************
  596. Tests for single functions of libm.
  597. Please keep them alphabetically sorted!
  598. ****************************************************************************/
  599. static void
  600. acos_test (void)
  601. {
  602. errno = 0;
  603. FUNC(acos) (0);
  604. if (errno == ENOSYS)
  605. /* Function not implemented. */
  606. return;
  607. START (acos);
  608. TEST_f_f (acos, plus_infty, nan_value, INVALID_EXCEPTION);
  609. TEST_f_f (acos, minus_infty, nan_value, INVALID_EXCEPTION);
  610. TEST_f_f (acos, nan_value, nan_value);
  611. /* |x| > 1: */
  612. TEST_f_f (acos, 1.125L, nan_value, INVALID_EXCEPTION);
  613. TEST_f_f (acos, -1.125L, nan_value, INVALID_EXCEPTION);
  614. TEST_f_f (acos, 0, M_PI_2l);
  615. TEST_f_f (acos, minus_zero, M_PI_2l);
  616. TEST_f_f (acos, 1, 0);
  617. TEST_f_f (acos, -1, M_PIl);
  618. TEST_f_f (acos, 0.5, M_PI_6l*2.0);
  619. TEST_f_f (acos, -0.5, M_PI_6l*4.0);
  620. TEST_f_f (acos, 0.75L, 0.722734247813415611178377352641333362L);
  621. TEST_f_f (acos, 2e-17L, 1.57079632679489659923132169163975144L);
  622. TEST_f_f (acos, 0.0625L, 1.50825556499840522843072005474337068L);
  623. END (acos);
  624. }
  625. static void
  626. acosh_test (void)
  627. {
  628. errno = 0;
  629. FUNC(acosh) (7);
  630. if (errno == ENOSYS)
  631. /* Function not implemented. */
  632. return;
  633. START (acosh);
  634. TEST_f_f (acosh, plus_infty, plus_infty);
  635. TEST_f_f (acosh, minus_infty, nan_value, INVALID_EXCEPTION);
  636. /* x < 1: */
  637. TEST_f_f (acosh, -1.125L, nan_value, INVALID_EXCEPTION);
  638. TEST_f_f (acosh, 1, 0);
  639. TEST_f_f (acosh, 7, 2.63391579384963341725009269461593689L);
  640. END (acosh);
  641. }
  642. static void
  643. asin_test (void)
  644. {
  645. errno = 0;
  646. FUNC(asin) (0);
  647. if (errno == ENOSYS)
  648. /* Function not implemented. */
  649. return;
  650. START (asin);
  651. TEST_f_f (asin, plus_infty, nan_value, INVALID_EXCEPTION);
  652. TEST_f_f (asin, minus_infty, nan_value, INVALID_EXCEPTION);
  653. TEST_f_f (asin, nan_value, nan_value);
  654. /* asin x == NaN plus invalid exception for |x| > 1. */
  655. TEST_f_f (asin, 1.125L, nan_value, INVALID_EXCEPTION);
  656. TEST_f_f (asin, -1.125L, nan_value, INVALID_EXCEPTION);
  657. TEST_f_f (asin, 0, 0);
  658. TEST_f_f (asin, minus_zero, minus_zero);
  659. TEST_f_f (asin, 0.5, M_PI_6l);
  660. TEST_f_f (asin, -0.5, -M_PI_6l);
  661. TEST_f_f (asin, 1.0, M_PI_2l);
  662. TEST_f_f (asin, -1.0, -M_PI_2l);
  663. TEST_f_f (asin, 0.75L, 0.848062078981481008052944338998418080L);
  664. END (asin);
  665. }
  666. static void
  667. asinh_test (void)
  668. {
  669. errno = 0;
  670. FUNC(asinh) (0.7L);
  671. if (errno == ENOSYS)
  672. /* Function not implemented. */
  673. return;
  674. START (asinh);
  675. TEST_f_f (asinh, 0, 0);
  676. TEST_f_f (asinh, minus_zero, minus_zero);
  677. #ifndef TEST_INLINE
  678. TEST_f_f (asinh, plus_infty, plus_infty);
  679. TEST_f_f (asinh, minus_infty, minus_infty);
  680. #endif
  681. TEST_f_f (asinh, nan_value, nan_value);
  682. TEST_f_f (asinh, 0.75L, 0.693147180559945309417232121458176568L);
  683. END (asinh);
  684. }
  685. static void
  686. atan_test (void)
  687. {
  688. errno = 0;
  689. FUNC(atan) (0);
  690. if (errno == ENOSYS)
  691. /* Function not implemented. */
  692. return;
  693. START (atan);
  694. TEST_f_f (atan, 0, 0);
  695. TEST_f_f (atan, minus_zero, minus_zero);
  696. TEST_f_f (atan, plus_infty, M_PI_2l);
  697. TEST_f_f (atan, minus_infty, -M_PI_2l);
  698. TEST_f_f (atan, nan_value, nan_value);
  699. TEST_f_f (atan, 1, M_PI_4l);
  700. TEST_f_f (atan, -1, -M_PI_4l);
  701. TEST_f_f (atan, 0.75L, 0.643501108793284386802809228717322638L);
  702. END (atan);
  703. }
  704. static void
  705. atanh_test (void)
  706. {
  707. errno = 0;
  708. FUNC(atanh) (0.7L);
  709. if (errno == ENOSYS)
  710. /* Function not implemented. */
  711. return;
  712. START (atanh);
  713. TEST_f_f (atanh, 0, 0);
  714. TEST_f_f (atanh, minus_zero, minus_zero);
  715. TEST_f_f (atanh, 1, plus_infty, DIVIDE_BY_ZERO_EXCEPTION);
  716. TEST_f_f (atanh, -1, minus_infty, DIVIDE_BY_ZERO_EXCEPTION);
  717. TEST_f_f (atanh, nan_value, nan_value);
  718. /* atanh (x) == NaN plus invalid exception if |x| > 1. */
  719. TEST_f_f (atanh, 1.125L, nan_value, INVALID_EXCEPTION);
  720. TEST_f_f (atanh, -1.125L, nan_value, INVALID_EXCEPTION);
  721. TEST_f_f (atanh, 0.75L, 0.972955074527656652552676371721589865L);
  722. END (atanh);
  723. }
  724. static void
  725. atan2_test (void)
  726. {
  727. errno = 0;
  728. FUNC(atan2) (-0, 1);
  729. if (errno == ENOSYS)
  730. /* Function not implemented. */
  731. return;
  732. START (atan2);
  733. /* atan2 (0,x) == 0 for x > 0. */
  734. TEST_ff_f (atan2, 0, 1, 0);
  735. /* atan2 (-0,x) == -0 for x > 0. */
  736. TEST_ff_f (atan2, minus_zero, 1, minus_zero);
  737. TEST_ff_f (atan2, 0, 0, 0);
  738. TEST_ff_f (atan2, minus_zero, 0, minus_zero);
  739. /* atan2 (+0,x) == +pi for x < 0. */
  740. TEST_ff_f (atan2, 0, -1, M_PIl);
  741. /* atan2 (-0,x) == -pi for x < 0. */
  742. TEST_ff_f (atan2, minus_zero, -1, -M_PIl);
  743. TEST_ff_f (atan2, 0, minus_zero, M_PIl);
  744. TEST_ff_f (atan2, minus_zero, minus_zero, -M_PIl);
  745. /* atan2 (y,+0) == pi/2 for y > 0. */
  746. TEST_ff_f (atan2, 1, 0, M_PI_2l);
  747. /* atan2 (y,-0) == pi/2 for y > 0. */
  748. TEST_ff_f (atan2, 1, minus_zero, M_PI_2l);
  749. /* atan2 (y,+0) == -pi/2 for y < 0. */
  750. TEST_ff_f (atan2, -1, 0, -M_PI_2l);
  751. /* atan2 (y,-0) == -pi/2 for y < 0. */
  752. TEST_ff_f (atan2, -1, minus_zero, -M_PI_2l);
  753. /* atan2 (y,inf) == +0 for finite y > 0. */
  754. TEST_ff_f (atan2, 1, plus_infty, 0);
  755. /* atan2 (y,inf) == -0 for finite y < 0. */
  756. TEST_ff_f (atan2, -1, plus_infty, minus_zero);
  757. /* atan2(+inf, x) == pi/2 for finite x. */
  758. TEST_ff_f (atan2, plus_infty, -1, M_PI_2l);
  759. /* atan2(-inf, x) == -pi/2 for finite x. */
  760. TEST_ff_f (atan2, minus_infty, 1, -M_PI_2l);
  761. /* atan2 (y,-inf) == +pi for finite y > 0. */
  762. TEST_ff_f (atan2, 1, minus_infty, M_PIl);
  763. /* atan2 (y,-inf) == -pi for finite y < 0. */
  764. TEST_ff_f (atan2, -1, minus_infty, -M_PIl);
  765. TEST_ff_f (atan2, plus_infty, plus_infty, M_PI_4l);
  766. TEST_ff_f (atan2, minus_infty, plus_infty, -M_PI_4l);
  767. TEST_ff_f (atan2, plus_infty, minus_infty, M_PI_34l);
  768. TEST_ff_f (atan2, minus_infty, minus_infty, -M_PI_34l);
  769. TEST_ff_f (atan2, nan_value, nan_value, nan_value);
  770. TEST_ff_f (atan2, 0.75L, 1, 0.643501108793284386802809228717322638L);
  771. TEST_ff_f (atan2, -0.75L, 1.0L, -0.643501108793284386802809228717322638L);
  772. TEST_ff_f (atan2, 0.75L, -1.0L, 2.49809154479650885165983415456218025L);
  773. TEST_ff_f (atan2, -0.75L, -1.0L, -2.49809154479650885165983415456218025L);
  774. TEST_ff_f (atan2, 0.390625L, .00029L, 1.57005392693128974780151246612928941L);
  775. TEST_ff_f (atan2, 1.390625L, 0.9296875L, 0.981498387184244311516296577615519772L);
  776. TEST_ff_f (atan2, -0.00756827042671106339L, -.001792735857538728036L, -1.80338464113663849327153994379639112L);
  777. END (atan2);
  778. }
  779. static void
  780. cabs_test (void)
  781. {
  782. errno = 0;
  783. FUNC(cabs) (BUILD_COMPLEX (0.7L, 12.4L));
  784. if (errno == ENOSYS)
  785. /* Function not implemented. */
  786. return;
  787. START (cabs);
  788. /* cabs (x + iy) is specified as hypot (x,y) */
  789. /* cabs (+inf + i x) == +inf. */
  790. TEST_c_f (cabs, plus_infty, 1.0, plus_infty);
  791. /* cabs (-inf + i x) == +inf. */
  792. TEST_c_f (cabs, minus_infty, 1.0, plus_infty);
  793. TEST_c_f (cabs, minus_infty, nan_value, plus_infty);
  794. TEST_c_f (cabs, minus_infty, nan_value, plus_infty);
  795. TEST_c_f (cabs, nan_value, nan_value, nan_value);
  796. /* cabs (x,y) == cabs (y,x). */
  797. TEST_c_f (cabs, 0.75L, 12.390625L, 12.4133028598606664302388810868156657L);
  798. /* cabs (x,y) == cabs (-x,y). */
  799. TEST_c_f (cabs, -12.390625L, 0.75L, 12.4133028598606664302388810868156657L);
  800. /* cabs (x,y) == cabs (-y,x). */
  801. TEST_c_f (cabs, -0.75L, 12.390625L, 12.4133028598606664302388810868156657L);
  802. /* cabs (x,y) == cabs (-x,-y). */
  803. TEST_c_f (cabs, -12.390625L, -0.75L, 12.4133028598606664302388810868156657L);
  804. /* cabs (x,y) == cabs (-y,-x). */
  805. TEST_c_f (cabs, -0.75L, -12.390625L, 12.4133028598606664302388810868156657L);
  806. /* cabs (x,0) == fabs (x). */
  807. TEST_c_f (cabs, -0.75L, 0, 0.75L);
  808. TEST_c_f (cabs, 0.75L, 0, 0.75L);
  809. TEST_c_f (cabs, -1.0L, 0, 1.0L);
  810. TEST_c_f (cabs, 1.0L, 0, 1.0L);
  811. TEST_c_f (cabs, -5.7e7L, 0, 5.7e7L);
  812. TEST_c_f (cabs, 5.7e7L, 0, 5.7e7L);
  813. TEST_c_f (cabs, 0.75L, 1.25L, 1.45773797371132511771853821938639577L);
  814. END (cabs);
  815. }
  816. #if 0
  817. static void
  818. cacos_test (void)
  819. {
  820. errno = 0;
  821. FUNC(cacos) (BUILD_COMPLEX (0.7L, 1.2L));
  822. if (errno == ENOSYS)
  823. /* Function not implemented. */
  824. return;
  825. START (cacos);
  826. TEST_c_c (cacos, 0, 0, M_PI_2l, minus_zero);
  827. TEST_c_c (cacos, minus_zero, 0, M_PI_2l, minus_zero);
  828. TEST_c_c (cacos, minus_zero, minus_zero, M_PI_2l, 0.0);
  829. TEST_c_c (cacos, 0, minus_zero, M_PI_2l, 0.0);
  830. TEST_c_c (cacos, minus_infty, plus_infty, M_PI_34l, minus_infty);
  831. TEST_c_c (cacos, minus_infty, minus_infty, M_PI_34l, plus_infty);
  832. TEST_c_c (cacos, plus_infty, plus_infty, M_PI_4l, minus_infty);
  833. TEST_c_c (cacos, plus_infty, minus_infty, M_PI_4l, plus_infty);
  834. TEST_c_c (cacos, -10.0, plus_infty, M_PI_2l, minus_infty);
  835. TEST_c_c (cacos, -10.0, minus_infty, M_PI_2l, plus_infty);
  836. TEST_c_c (cacos, 0, plus_infty, M_PI_2l, minus_infty);
  837. TEST_c_c (cacos, 0, minus_infty, M_PI_2l, plus_infty);
  838. TEST_c_c (cacos, 0.1L, plus_infty, M_PI_2l, minus_infty);
  839. TEST_c_c (cacos, 0.1L, minus_infty, M_PI_2l, plus_infty);
  840. TEST_c_c (cacos, minus_infty, 0, M_PIl, minus_infty);
  841. TEST_c_c (cacos, minus_infty, minus_zero, M_PIl, plus_infty);
  842. TEST_c_c (cacos, minus_infty, 100, M_PIl, minus_infty);
  843. TEST_c_c (cacos, minus_infty, -100, M_PIl, plus_infty);
  844. TEST_c_c (cacos, plus_infty, 0, 0.0, minus_infty);
  845. TEST_c_c (cacos, plus_infty, minus_zero, 0.0, plus_infty);
  846. TEST_c_c (cacos, plus_infty, 0.5, 0.0, minus_infty);
  847. TEST_c_c (cacos, plus_infty, -0.5, 0.0, plus_infty);
  848. TEST_c_c (cacos, plus_infty, nan_value, nan_value, plus_infty, IGNORE_ZERO_INF_SIGN);
  849. TEST_c_c (cacos, minus_infty, nan_value, nan_value, plus_infty, IGNORE_ZERO_INF_SIGN);
  850. TEST_c_c (cacos, 0, nan_value, M_PI_2l, nan_value);
  851. TEST_c_c (cacos, minus_zero, nan_value, M_PI_2l, nan_value);
  852. TEST_c_c (cacos, nan_value, plus_infty, nan_value, minus_infty);
  853. TEST_c_c (cacos, nan_value, minus_infty, nan_value, plus_infty);
  854. TEST_c_c (cacos, 10.5, nan_value, nan_value, nan_value, INVALID_EXCEPTION_OK);
  855. TEST_c_c (cacos, -10.5, nan_value, nan_value, nan_value, INVALID_EXCEPTION_OK);
  856. TEST_c_c (cacos, nan_value, 0.75, nan_value, nan_value, INVALID_EXCEPTION_OK);
  857. TEST_c_c (cacos, nan_value, -0.75, nan_value, nan_value, INVALID_EXCEPTION_OK);
  858. TEST_c_c (cacos, nan_value, nan_value, nan_value, nan_value);
  859. TEST_c_c (cacos, 0.75L, 1.25L, 1.11752014915610270578240049553777969L, -1.13239363160530819522266333696834467L);
  860. TEST_c_c (cacos, -2, -3, 2.1414491111159960199416055713254211L, 1.9833870299165354323470769028940395L);
  861. END (cacos, complex);
  862. }
  863. static void
  864. cacosh_test (void)
  865. {
  866. errno = 0;
  867. FUNC(cacosh) (BUILD_COMPLEX (0.7L, 1.2L));
  868. if (errno == ENOSYS)
  869. /* Function not implemented. */
  870. return;
  871. START (cacosh);
  872. TEST_c_c (cacosh, 0, 0, 0.0, M_PI_2l);
  873. TEST_c_c (cacosh, minus_zero, 0, 0.0, M_PI_2l);
  874. TEST_c_c (cacosh, 0, minus_zero, 0.0, -M_PI_2l);
  875. TEST_c_c (cacosh, minus_zero, minus_zero, 0.0, -M_PI_2l);
  876. TEST_c_c (cacosh, minus_infty, plus_infty, plus_infty, M_PI_34l);
  877. TEST_c_c (cacosh, minus_infty, minus_infty, plus_infty, -M_PI_34l);
  878. TEST_c_c (cacosh, plus_infty, plus_infty, plus_infty, M_PI_4l);
  879. TEST_c_c (cacosh, plus_infty, minus_infty, plus_infty, -M_PI_4l);
  880. TEST_c_c (cacosh, -10.0, plus_infty, plus_infty, M_PI_2l);
  881. TEST_c_c (cacosh, -10.0, minus_infty, plus_infty, -M_PI_2l);
  882. TEST_c_c (cacosh, 0, plus_infty, plus_infty, M_PI_2l);
  883. TEST_c_c (cacosh, 0, minus_infty, plus_infty, -M_PI_2l);
  884. TEST_c_c (cacosh, 0.1L, plus_infty, plus_infty, M_PI_2l);
  885. TEST_c_c (cacosh, 0.1L, minus_infty, plus_infty, -M_PI_2l);
  886. TEST_c_c (cacosh, minus_infty, 0, plus_infty, M_PIl);
  887. TEST_c_c (cacosh, minus_infty, minus_zero, plus_infty, -M_PIl);
  888. TEST_c_c (cacosh, minus_infty, 100, plus_infty, M_PIl);
  889. TEST_c_c (cacosh, minus_infty, -100, plus_infty, -M_PIl);
  890. TEST_c_c (cacosh, plus_infty, 0, plus_infty, 0.0);
  891. TEST_c_c (cacosh, plus_infty, minus_zero, plus_infty, minus_zero);
  892. TEST_c_c (cacosh, plus_infty, 0.5, plus_infty, 0.0);
  893. TEST_c_c (cacosh, plus_infty, -0.5, plus_infty, minus_zero);
  894. TEST_c_c (cacosh, plus_infty, nan_value, plus_infty, nan_value);
  895. TEST_c_c (cacosh, minus_infty, nan_value, plus_infty, nan_value);
  896. TEST_c_c (cacosh, 0, nan_value, nan_value, nan_value);
  897. TEST_c_c (cacosh, minus_zero, nan_value, nan_value, nan_value);
  898. TEST_c_c (cacosh, nan_value, plus_infty, plus_infty, nan_value);
  899. TEST_c_c (cacosh, nan_value, minus_infty, plus_infty, nan_value);
  900. TEST_c_c (cacosh, 10.5, nan_value, nan_value, nan_value, INVALID_EXCEPTION_OK);
  901. TEST_c_c (cacosh, -10.5, nan_value, nan_value, nan_value, INVALID_EXCEPTION_OK);
  902. TEST_c_c (cacosh, nan_value, 0.75, nan_value, nan_value, INVALID_EXCEPTION_OK);
  903. TEST_c_c (cacosh, nan_value, -0.75, nan_value, nan_value, INVALID_EXCEPTION_OK);
  904. TEST_c_c (cacosh, nan_value, nan_value, nan_value, nan_value);
  905. TEST_c_c (cacosh, 0.75L, 1.25L, 1.13239363160530819522266333696834467L, 1.11752014915610270578240049553777969L);
  906. TEST_c_c (cacosh, -2, -3, 1.9833870299165354323470769028940395L, -2.1414491111159960199416055713254211L);
  907. END (cacosh, complex);
  908. }
  909. static void
  910. carg_test (void)
  911. {
  912. START (carg);
  913. /* carg (x + iy) is specified as atan2 (y, x) */
  914. /* carg (x + i 0) == 0 for x > 0. */
  915. TEST_c_f (carg, 2.0, 0, 0);
  916. /* carg (x - i 0) == -0 for x > 0. */
  917. TEST_c_f (carg, 2.0, minus_zero, minus_zero);
  918. TEST_c_f (carg, 0, 0, 0);
  919. TEST_c_f (carg, 0, minus_zero, minus_zero);
  920. /* carg (x + i 0) == +pi for x < 0. */
  921. TEST_c_f (carg, -2.0, 0, M_PIl);
  922. /* carg (x - i 0) == -pi for x < 0. */
  923. TEST_c_f (carg, -2.0, minus_zero, -M_PIl);
  924. TEST_c_f (carg, minus_zero, 0, M_PIl);
  925. TEST_c_f (carg, minus_zero, minus_zero, -M_PIl);
  926. /* carg (+0 + i y) == pi/2 for y > 0. */
  927. TEST_c_f (carg, 0, 2.0, M_PI_2l);
  928. /* carg (-0 + i y) == pi/2 for y > 0. */
  929. TEST_c_f (carg, minus_zero, 2.0, M_PI_2l);
  930. /* carg (+0 + i y) == -pi/2 for y < 0. */
  931. TEST_c_f (carg, 0, -2.0, -M_PI_2l);
  932. /* carg (-0 + i y) == -pi/2 for y < 0. */
  933. TEST_c_f (carg, minus_zero, -2.0, -M_PI_2l);
  934. /* carg (inf + i y) == +0 for finite y > 0. */
  935. TEST_c_f (carg, plus_infty, 2.0, 0);
  936. /* carg (inf + i y) == -0 for finite y < 0. */
  937. TEST_c_f (carg, plus_infty, -2.0, minus_zero);
  938. /* carg(x + i inf) == pi/2 for finite x. */
  939. TEST_c_f (carg, 10.0, plus_infty, M_PI_2l);
  940. /* carg(x - i inf) == -pi/2 for finite x. */
  941. TEST_c_f (carg, 10.0, minus_infty, -M_PI_2l);
  942. /* carg (-inf + i y) == +pi for finite y > 0. */
  943. TEST_c_f (carg, minus_infty, 10.0, M_PIl);
  944. /* carg (-inf + i y) == -pi for finite y < 0. */
  945. TEST_c_f (carg, minus_infty, -10.0, -M_PIl);
  946. TEST_c_f (carg, plus_infty, plus_infty, M_PI_4l);
  947. TEST_c_f (carg, plus_infty, minus_infty, -M_PI_4l);
  948. TEST_c_f (carg, minus_infty, plus_infty, 3 * M_PI_4l);
  949. TEST_c_f (carg, minus_infty, minus_infty, -3 * M_PI_4l);
  950. TEST_c_f (carg, nan_value, nan_value, nan_value);
  951. END (carg);
  952. }
  953. static void
  954. casin_test (void)
  955. {
  956. errno = 0;
  957. FUNC(casin) (BUILD_COMPLEX (0.7L, 1.2L));
  958. if (errno == ENOSYS)
  959. /* Function not implemented. */
  960. return;
  961. START (casin);
  962. TEST_c_c (casin, 0, 0, 0.0, 0.0);
  963. TEST_c_c (casin, minus_zero, 0, minus_zero, 0.0);
  964. TEST_c_c (casin, 0, minus_zero, 0.0, minus_zero);
  965. TEST_c_c (casin, minus_zero, minus_zero, minus_zero, minus_zero);
  966. TEST_c_c (casin, plus_infty, plus_infty, M_PI_4l, plus_infty);
  967. TEST_c_c (casin, plus_infty, minus_infty, M_PI_4l, minus_infty);
  968. TEST_c_c (casin, minus_infty, plus_infty, -M_PI_4l, plus_infty);
  969. TEST_c_c (casin, minus_infty, minus_infty, -M_PI_4l, minus_infty);
  970. TEST_c_c (casin, -10.0, plus_infty, minus_zero, plus_infty);
  971. TEST_c_c (casin, -10.0, minus_infty, minus_zero, minus_infty);
  972. TEST_c_c (casin, 0, plus_infty, 0.0, plus_infty);
  973. TEST_c_c (casin, 0, minus_infty, 0.0, minus_infty);
  974. TEST_c_c (casin, minus_zero, plus_infty, minus_zero, plus_infty);
  975. TEST_c_c (casin, minus_zero, minus_infty, minus_zero, minus_infty);
  976. TEST_c_c (casin, 0.1L, plus_infty, 0.0, plus_infty);
  977. TEST_c_c (casin, 0.1L, minus_infty, 0.0, minus_infty);
  978. TEST_c_c (casin, minus_infty, 0, -M_PI_2l, plus_infty);
  979. TEST_c_c (casin, minus_infty, minus_zero, -M_PI_2l, minus_infty);
  980. TEST_c_c (casin, minus_infty, 100, -M_PI_2l, plus_infty);
  981. TEST_c_c (casin, minus_infty, -100, -M_PI_2l, minus_infty);
  982. TEST_c_c (casin, plus_infty, 0, M_PI_2l, plus_infty);
  983. TEST_c_c (casin, plus_infty, minus_zero, M_PI_2l, minus_infty);
  984. TEST_c_c (casin, plus_infty, 0.5, M_PI_2l, plus_infty);
  985. TEST_c_c (casin, plus_infty, -0.5, M_PI_2l, minus_infty);
  986. TEST_c_c (casin, nan_value, plus_infty, nan_value, plus_infty);
  987. TEST_c_c (casin, nan_value, minus_infty, nan_value, minus_infty);
  988. TEST_c_c (casin, 0.0, nan_value, 0.0, nan_value);
  989. TEST_c_c (casin, minus_zero, nan_value, minus_zero, nan_value);
  990. TEST_c_c (casin, plus_infty, nan_value, nan_value, plus_infty, IGNORE_ZERO_INF_SIGN);
  991. TEST_c_c (casin, minus_infty, nan_value, nan_value, plus_infty, IGNORE_ZERO_INF_SIGN);
  992. TEST_c_c (casin, nan_value, 10.5, nan_value, nan_value, INVALID_EXCEPTION_OK);
  993. TEST_c_c (casin, nan_value, -10.5, nan_value, nan_value, INVALID_EXCEPTION_OK);
  994. TEST_c_c (casin, 0.75, nan_value, nan_value, nan_value, INVALID_EXCEPTION_OK);
  995. TEST_c_c (casin, -0.75, nan_value, nan_value, nan_value, INVALID_EXCEPTION_OK);
  996. TEST_c_c (casin, nan_value, nan_value, nan_value, nan_value);
  997. TEST_c_c (casin, 0.75L, 1.25L, 0.453276177638793913448921196101971749L, 1.13239363160530819522266333696834467L);
  998. TEST_c_c (casin, -2, -3, -0.57065278432109940071028387968566963L, -1.9833870299165354323470769028940395L);
  999. END (casin, complex);
  1000. }
  1001. static void
  1002. casinh_test (void)
  1003. {
  1004. errno = 0;
  1005. FUNC(casinh) (BUILD_COMPLEX (0.7L, 1.2L));
  1006. if (errno == ENOSYS)
  1007. /* Function not implemented. */
  1008. return;
  1009. START (casinh);
  1010. TEST_c_c (casinh, 0, 0, 0.0, 0.0);
  1011. TEST_c_c (casinh, minus_zero, 0, minus_zero, 0);
  1012. TEST_c_c (casinh, 0, minus_zero, 0.0, minus_zero);
  1013. TEST_c_c (casinh, minus_zero, minus_zero, minus_zero, minus_zero);
  1014. TEST_c_c (casinh, plus_infty, plus_infty, plus_infty, M_PI_4l);
  1015. TEST_c_c (casinh, plus_infty, minus_infty, plus_infty, -M_PI_4l);
  1016. TEST_c_c (casinh, minus_infty, plus_infty, minus_infty, M_PI_4l);
  1017. TEST_c_c (casinh, minus_infty, minus_infty, minus_infty, -M_PI_4l);
  1018. TEST_c_c (casinh, -10.0, plus_infty, minus_infty, M_PI_2l);
  1019. TEST_c_c (casinh, -10.0, minus_infty, minus_infty, -M_PI_2l);
  1020. TEST_c_c (casinh, 0, plus_infty, plus_infty, M_PI_2l);
  1021. TEST_c_c (casinh, 0, minus_infty, plus_infty, -M_PI_2l);
  1022. TEST_c_c (casinh, minus_zero, plus_infty, minus_infty, M_PI_2l);
  1023. TEST_c_c (casinh, minus_zero, minus_infty, minus_infty, -M_PI_2l);
  1024. TEST_c_c (casinh, 0.1L, plus_infty, plus_infty, M_PI_2l);
  1025. TEST_c_c (casinh, 0.1L, minus_infty, plus_infty, -M_PI_2l);
  1026. TEST_c_c (casinh, minus_infty, 0, minus_infty, 0.0);
  1027. TEST_c_c (casinh, minus_infty, minus_zero, minus_infty, minus_zero);
  1028. TEST_c_c (casinh, minus_infty, 100, minus_infty, 0.0);
  1029. TEST_c_c (casinh, minus_infty, -100, minus_infty, minus_zero);
  1030. TEST_c_c (casinh, plus_infty, 0, plus_infty, 0.0);
  1031. TEST_c_c (casinh, plus_infty, minus_zero, plus_infty, minus_zero);
  1032. TEST_c_c (casinh, plus_infty, 0.5, plus_infty, 0.0);
  1033. TEST_c_c (casinh, plus_infty, -0.5, plus_infty, minus_zero);
  1034. TEST_c_c (casinh, plus_infty, nan_value, plus_infty, nan_value);
  1035. TEST_c_c (casinh, minus_infty, nan_value, minus_infty, nan_value);
  1036. TEST_c_c (casinh, nan_value, 0, nan_value, 0.0);
  1037. TEST_c_c (casinh, nan_value, minus_zero, nan_value, minus_zero);
  1038. TEST_c_c (casinh, nan_value, plus_infty, plus_infty, nan_value, IGNORE_ZERO_INF_SIGN);
  1039. TEST_c_c (casinh, nan_value, minus_infty, plus_infty, nan_value, IGNORE_ZERO_INF_SIGN);
  1040. TEST_c_c (casinh, 10.5, nan_value, nan_value, nan_value, INVALID_EXCEPTION_OK);
  1041. TEST_c_c (casinh, -10.5, nan_value, nan_value, nan_value, INVALID_EXCEPTION_OK);
  1042. TEST_c_c (casinh, nan_value, 0.75, nan_value, nan_value, INVALID_EXCEPTION_OK);
  1043. TEST_c_c (casinh, -0.75, nan_value, nan_value, nan_value, INVALID_EXCEPTION_OK);
  1044. TEST_c_c (casinh, nan_value, nan_value, nan_value, nan_value);
  1045. TEST_c_c (casinh, 0.75L, 1.25L, 1.03171853444778027336364058631006594L, 0.911738290968487636358489564316731207L);
  1046. TEST_c_c (casinh, -2, -3, -1.9686379257930962917886650952454982L, -0.96465850440760279204541105949953237L);
  1047. END (casinh, complex);
  1048. }
  1049. static void
  1050. catan_test (void)
  1051. {
  1052. errno = 0;
  1053. FUNC(catan) (BUILD_COMPLEX (0.7L, 1.2L));
  1054. if (errno == ENOSYS)
  1055. /* Function not implemented. */
  1056. return;
  1057. START (catan);
  1058. TEST_c_c (catan, 0, 0, 0, 0);
  1059. TEST_c_c (catan, minus_zero, 0, minus_zero, 0);
  1060. TEST_c_c (catan, 0, minus_zero, 0, minus_zero);
  1061. TEST_c_c (catan, minus_zero, minus_zero, minus_zero, minus_zero);
  1062. TEST_c_c (catan, plus_infty, plus_infty, M_PI_2l, 0);
  1063. TEST_c_c (catan, plus_infty, minus_infty, M_PI_2l, minus_zero);
  1064. TEST_c_c (catan, minus_infty, plus_infty, -M_PI_2l, 0);
  1065. TEST_c_c (catan, minus_infty, minus_infty, -M_PI_2l, minus_zero);
  1066. TEST_c_c (catan, plus_infty, -10.0, M_PI_2l, minus_zero);
  1067. TEST_c_c (catan, minus_infty, -10.0, -M_PI_2l, minus_zero);
  1068. TEST_c_c (catan, plus_infty, minus_zero, M_PI_2l, minus_zero);
  1069. TEST_c_c (catan, minus_infty, minus_zero, -M_PI_2l, minus_zero);
  1070. TEST_c_c (catan, plus_infty, 0.0, M_PI_2l, 0);
  1071. TEST_c_c (catan, minus_infty, 0.0, -M_PI_2l, 0);
  1072. TEST_c_c (catan, plus_infty, 0.1L, M_PI_2l, 0);
  1073. TEST_c_c (catan, minus_infty, 0.1L, -M_PI_2l, 0);
  1074. TEST_c_c (catan, 0.0, minus_infty, M_PI_2l, minus_zero);
  1075. TEST_c_c (catan, minus_zero, minus_infty, -M_PI_2l, minus_zero);
  1076. TEST_c_c (catan, 100.0, minus_infty, M_PI_2l, minus_zero);
  1077. TEST_c_c (catan, -100.0, minus_infty, -M_PI_2l, minus_zero);
  1078. TEST_c_c (catan, 0.0, plus_infty, M_PI_2l, 0);
  1079. TEST_c_c (catan, minus_zero, plus_infty, -M_PI_2l, 0);
  1080. TEST_c_c (catan, 0.5, plus_infty, M_PI_2l, 0);
  1081. TEST_c_c (catan, -0.5, plus_infty, -M_PI_2l, 0);
  1082. TEST_c_c (catan, nan_value, 0.0, nan_value, 0);
  1083. TEST_c_c (catan, nan_value, minus_zero, nan_value, minus_zero);
  1084. TEST_c_c (catan, nan_value, plus_infty, nan_value, 0);
  1085. TEST_c_c (catan, nan_value, minus_infty, nan_value, minus_zero);
  1086. TEST_c_c (catan, 0.0, nan_value, nan_value, nan_value);
  1087. TEST_c_c (catan, minus_zero, nan_value, nan_value, nan_value);
  1088. TEST_c_c (catan, plus_infty, nan_value, M_PI_2l, 0, IGNORE_ZERO_INF_SIGN);
  1089. TEST_c_c (catan, minus_infty, nan_value, -M_PI_2l, 0, IGNORE_ZERO_INF_SIGN);
  1090. TEST_c_c (catan, nan_value, 10.5, nan_value, nan_value, INVALID_EXCEPTION_OK);
  1091. TEST_c_c (catan, nan_value, -10.5, nan_value, nan_value, INVALID_EXCEPTION_OK);
  1092. TEST_c_c (catan, 0.75, nan_value, nan_value, nan_value, INVALID_EXCEPTION_OK);
  1093. TEST_c_c (catan, -0.75, nan_value, nan_value, nan_value, INVALID_EXCEPTION_OK);
  1094. TEST_c_c (catan, nan_value, nan_value, nan_value, nan_value);
  1095. TEST_c_c (catan, 0.75L, 1.25L, 1.10714871779409050301706546017853704L, 0.549306144334054845697622618461262852L);
  1096. TEST_c_c (catan, -2, -3, -1.4099210495965755225306193844604208L, -0.22907268296853876629588180294200276L);
  1097. END (catan, complex);
  1098. }
  1099. static void
  1100. catanh_test (void)
  1101. {
  1102. errno = 0;
  1103. FUNC(catanh) (BUILD_COMPLEX (0.7L, 1.2L));
  1104. if (errno == ENOSYS)
  1105. /* Function not implemented. */
  1106. return;
  1107. START (catanh);
  1108. TEST_c_c (catanh, 0, 0, 0.0, 0.0);
  1109. TEST_c_c (catanh, minus_zero, 0, minus_zero, 0.0);
  1110. TEST_c_c (catanh, 0, minus_zero, 0.0, minus_zero);
  1111. TEST_c_c (catanh, minus_zero, minus_zero, minus_zero, minus_zero);
  1112. TEST_c_c (catanh, plus_infty, plus_infty, 0.0, M_PI_2l);
  1113. TEST_c_c (catanh, plus_infty, minus_infty, 0.0, -M_PI_2l);
  1114. TEST_c_c (catanh, minus_infty, plus_infty, minus_zero, M_PI_2l);
  1115. TEST_c_c (catanh, minus_infty, minus_infty, minus_zero, -M_PI_2l);
  1116. TEST_c_c (catanh, -10.0, plus_infty, minus_zero, M_PI_2l);
  1117. TEST_c_c (catanh, -10.0, minus_infty, minus_zero, -M_PI_2l);
  1118. TEST_c_c (catanh, minus_zero, plus_infty, minus_zero, M_PI_2l);
  1119. TEST_c_c (catanh, minus_zero, minus_infty, minus_zero, -M_PI_2l);
  1120. TEST_c_c (catanh, 0, plus_infty, 0.0, M_PI_2l);
  1121. TEST_c_c (catanh, 0, minus_infty, 0.0, -M_PI_2l);
  1122. TEST_c_c (catanh, 0.1L, plus_infty, 0.0, M_PI_2l);
  1123. TEST_c_c (catanh, 0.1L, minus_infty, 0.0, -M_PI_2l);
  1124. TEST_c_c (catanh, minus_infty, 0, minus_zero, M_PI_2l);
  1125. TEST_c_c (catanh, minus_infty, minus_zero, minus_zero, -M_PI_2l);
  1126. TEST_c_c (catanh, minus_infty, 100, minus_zero, M_PI_2l);
  1127. TEST_c_c (catanh, minus_infty, -100, minus_zero, -M_PI_2l);
  1128. TEST_c_c (catanh, plus_infty, 0, 0.0, M_PI_2l);
  1129. TEST_c_c (catanh, plus_infty, minus_zero, 0.0, -M_PI_2l);
  1130. TEST_c_c (catanh, plus_infty, 0.5, 0.0, M_PI_2l);
  1131. TEST_c_c (catanh, plus_infty, -0.5, 0.0, -M_PI_2l);
  1132. TEST_c_c (catanh, 0, nan_value, 0.0, nan_value);
  1133. TEST_c_c (catanh, minus_zero, nan_value, minus_zero, nan_value);
  1134. TEST_c_c (catanh, plus_infty, nan_value, 0.0, nan_value);
  1135. TEST_c_c (catanh, minus_infty, nan_value, minus_zero, nan_value);
  1136. TEST_c_c (catanh, nan_value, 0, nan_value, nan_value);
  1137. TEST_c_c (catanh, nan_value, minus_zero, nan_value, nan_value);
  1138. TEST_c_c (catanh, nan_value, plus_infty, 0.0, M_PI_2l, IGNORE_ZERO_INF_SIGN);
  1139. TEST_c_c (catanh, nan_value, minus_infty, 0.0, -M_PI_2l, IGNORE_ZERO_INF_SIGN);
  1140. TEST_c_c (catanh, 10.5, nan_value, nan_value, nan_value, INVALID_EXCEPTION_OK);
  1141. TEST_c_c (catanh, -10.5, nan_value, nan_value, nan_value, INVALID_EXCEPTION_OK);
  1142. TEST_c_c (catanh, nan_value, 0.75, nan_value, nan_value, INVALID_EXCEPTION_OK);
  1143. TEST_c_c (catanh, nan_value, -0.75, nan_value, nan_value, INVALID_EXCEPTION_OK);
  1144. TEST_c_c (catanh, nan_value, nan_value, nan_value, nan_value);
  1145. TEST_c_c (catanh, 0.75L, 1.25L, 0.261492138795671927078652057366532140L, 0.996825126463918666098902241310446708L);
  1146. TEST_c_c (catanh, -2, -3, -0.14694666622552975204743278515471595L, -1.3389725222944935611241935759091443L);
  1147. END (catanh, complex);
  1148. }
  1149. #endif
  1150. static void
  1151. cbrt_test (void)
  1152. {
  1153. errno = 0;
  1154. FUNC(cbrt) (8);
  1155. if (errno == ENOSYS)
  1156. /* Function not implemented. */
  1157. return;
  1158. START (cbrt);
  1159. TEST_f_f (cbrt, 0.0, 0.0);
  1160. TEST_f_f (cbrt, minus_zero, minus_zero);
  1161. TEST_f_f (cbrt, plus_infty, plus_infty);
  1162. TEST_f_f (cbrt, minus_infty, minus_infty);
  1163. TEST_f_f (cbrt, nan_value, nan_value);
  1164. TEST_f_f (cbrt, -0.001L, -0.1L);
  1165. TEST_f_f (cbrt, 8, 2);
  1166. TEST_f_f (cbrt, -27.0, -3.0);
  1167. TEST_f_f (cbrt, 0.9921875L, 0.997389022060725270579075195353955217L);
  1168. TEST_f_f (cbrt, 0.75L, 0.908560296416069829445605878163630251L);
  1169. END (cbrt);
  1170. }
  1171. #if 0
  1172. static void
  1173. ccos_test (void)
  1174. {
  1175. errno = 0;
  1176. FUNC(ccos) (BUILD_COMPLEX (0, 0));
  1177. if (errno == ENOSYS)
  1178. /* Function not implemented. */
  1179. return;
  1180. START (ccos);
  1181. TEST_c_c (ccos, 0.0, 0.0, 1.0, minus_zero);
  1182. TEST_c_c (ccos, minus_zero, 0.0, 1.0, 0.0);
  1183. TEST_c_c (ccos, 0.0, minus_zero, 1.0, 0.0);
  1184. TEST_c_c (ccos, minus_zero, minus_zero, 1.0, minus_zero);
  1185. TEST_c_c (ccos, plus_infty, 0.0, nan_value, 0.0, INVALID_EXCEPTION|IGNORE_ZERO_INF_SIGN);
  1186. TEST_c_c (ccos, plus_infty, minus_zero, nan_value, 0.0, INVALID_EXCEPTION|IGNORE_ZERO_INF_SIGN);
  1187. TEST_c_c (ccos, minus_infty, 0.0, nan_value, 0.0, INVALID_EXCEPTION|IGNORE_ZERO_INF_SIGN);
  1188. TEST_c_c (ccos, minus_infty, minus_zero, nan_value, 0.0, INVALID_EXCEPTION|IGNORE_ZERO_INF_SIGN);
  1189. TEST_c_c (ccos, 0.0, plus_infty, plus_infty, minus_zero);
  1190. TEST_c_c (ccos, 0.0, minus_infty, plus_infty, 0.0);
  1191. TEST_c_c (ccos, minus_zero, plus_infty, plus_infty, 0.0);
  1192. TEST_c_c (ccos, minus_zero, minus_infty, plus_infty, minus_zero);
  1193. TEST_c_c (ccos, plus_infty, plus_infty, plus_infty, nan_value, INVALID_EXCEPTION);
  1194. TEST_c_c (ccos, minus_infty, plus_infty, plus_infty, nan_value, INVALID_EXCEPTION);
  1195. TEST_c_c (ccos, plus_infty, minus_infty, plus_infty, nan_value, INVALID_EXCEPTION);
  1196. TEST_c_c (ccos, minus_infty, minus_infty, plus_infty, nan_value, INVALID_EXCEPTION);
  1197. TEST_c_c (ccos, 4.625, plus_infty, minus_infty, plus_infty);
  1198. TEST_c_c (ccos, 4.625, minus_infty, minus_infty, minus_infty);
  1199. TEST_c_c (ccos, -4.625, plus_infty, minus_infty, minus_infty);
  1200. TEST_c_c (ccos, -4.625, minus_infty, minus_infty, plus_infty);
  1201. TEST_c_c (ccos, plus_infty, 6.75, nan_value, nan_value, INVALID_EXCEPTION);
  1202. TEST_c_c (ccos, plus_infty, -6.75, nan_value, nan_value, INVALID_EXCEPTION);
  1203. TEST_c_c (ccos, minus_infty, 6.75, nan_value, nan_value, INVALID_EXCEPTION);
  1204. TEST_c_c (ccos, minus_infty, -6.75, nan_value, nan_value, INVALID_EXCEPTION);
  1205. TEST_c_c (ccos, nan_value, 0.0, nan_value, 0.0, IGNORE_ZERO_INF_SIGN);
  1206. TEST_c_c (ccos, nan_value, minus_zero, nan_value, 0.0, IGNORE_ZERO_INF_SIGN);
  1207. TEST_c_c (ccos, nan_value, plus_infty, plus_infty, nan_value);
  1208. TEST_c_c (ccos, nan_value, minus_infty, plus_infty, nan_value);
  1209. TEST_c_c (ccos, nan_value, 9.0, nan_value, nan_value, INVALID_EXCEPTION_OK);
  1210. TEST_c_c (ccos, nan_value, -9.0, nan_value, nan_value, INVALID_EXCEPTION_OK);
  1211. TEST_c_c (ccos, 0.0, nan_value, nan_value, 0.0, IGNORE_ZERO_INF_SIGN);
  1212. TEST_c_c (ccos, minus_zero, nan_value, nan_value, 0.0, IGNORE_ZERO_INF_SIGN);
  1213. TEST_c_c (ccos, 10.0, nan_value, nan_value, nan_value, INVALID_EXCEPTION_OK);
  1214. TEST_c_c (ccos, -10.0, nan_value, nan_value, nan_value, INVALID_EXCEPTION_OK);
  1215. TEST_c_c (ccos, plus_infty, nan_value, nan_value, nan_value, INVALID_EXCEPTION_OK);
  1216. TEST_c_c (ccos, minus_infty, nan_value, nan_value, nan_value, INVALID_EXCEPTION_OK);
  1217. TEST_c_c (ccos, nan_value, nan_value, nan_value, nan_value);
  1218. TEST_c_c (ccos, 0.75L, 1.25L, 1.38173873063425888530729933139078645L, -1.09193013555397466170919531722024128L);
  1219. TEST_c_c (ccos, -2, -3, -4.18962569096880723013255501961597373L, -9.10922789375533659797919726277886212L);
  1220. END (ccos, complex);
  1221. }
  1222. static void
  1223. ccosh_test (void)
  1224. {
  1225. errno = 0;
  1226. FUNC(ccosh) (BUILD_COMPLEX (0.7L, 1.2L));
  1227. if (errno == ENOSYS)
  1228. /* Function not implemented. */
  1229. return;
  1230. START (ccosh);
  1231. TEST_c_c (ccosh, 0.0, 0.0, 1.0, 0.0);
  1232. TEST_c_c (ccosh, minus_zero, 0.0, 1.0, minus_zero);
  1233. TEST_c_c (ccosh, 0.0, minus_zero, 1.0, minus_zero);
  1234. TEST_c_c (ccosh, minus_zero, minus_zero, 1.0, 0.0);
  1235. TEST_c_c (ccosh, 0.0, plus_infty, nan_value, 0.0, INVALID_EXCEPTION|IGNORE_ZERO_INF_SIGN);
  1236. TEST_c_c (ccosh, minus_zero, plus_infty, nan_value, 0.0, INVALID_EXCEPTION|IGNORE_ZERO_INF_SIGN);
  1237. TEST_c_c (ccosh, 0.0, minus_infty, nan_value, 0.0, INVALID_EXCEPTION|IGNORE_ZERO_INF_SIGN);
  1238. TEST_c_c (ccosh, minus_zero, minus_infty, nan_value, 0.0, INVALID_EXCEPTION|IGNORE_ZERO_INF_SIGN);
  1239. TEST_c_c (ccosh, plus_infty, 0.0, plus_infty, 0.0);
  1240. TEST_c_c (ccosh, minus_infty, 0.0, plus_infty, minus_zero);
  1241. TEST_c_c (ccosh, plus_infty, minus_zero, plus_infty, minus_zero);
  1242. TEST_c_c (ccosh, minus_infty, minus_zero, plus_infty, 0.0);
  1243. TEST_c_c (ccosh, plus_infty, plus_infty, plus_infty, nan_value, INVALID_EXCEPTION);
  1244. TEST_c_c (ccosh, minus_infty, plus_infty, plus_infty, nan_value, INVALID_EXCEPTION);
  1245. TEST_c_c (ccosh, plus_infty, minus_infty, plus_infty, nan_value, INVALID_EXCEPTION);
  1246. TEST_c_c (ccosh, minus_infty, minus_infty, plus_infty, nan_value, INVALID_EXCEPTION);
  1247. TEST_c_c (ccosh, plus_infty, 4.625, minus_infty, minus_infty);
  1248. TEST_c_c (ccosh, minus_infty, 4.625, minus_infty, plus_infty);
  1249. TEST_c_c (ccosh, plus_infty, -4.625, minus_infty, plus_infty);
  1250. TEST_c_c (ccosh, minus_infty, -4.625, minus_infty, minus_infty);
  1251. TEST_c_c (ccosh, 6.75, plus_infty, nan_value, nan_value, INVALID_EXCEPTION);
  1252. TEST_c_c (ccosh, -6.75, plus_infty, nan_value, nan_value, INVALID_EXCEPTION);
  1253. TEST_c_c (ccosh, 6.75, minus_infty, nan_value, nan_value, INVALID_EXCEPTION);
  1254. TEST_c_c (ccosh, -6.75, minus_infty, nan_value, nan_value, INVALID_EXCEPTION);
  1255. TEST_c_c (ccosh, 0.0, nan_value, nan_value, 0.0, IGNORE_ZERO_INF_SIGN);
  1256. TEST_c_c (ccosh, minus_zero, nan_value, nan_value, 0.0, IGNORE_ZERO_INF_SIGN);
  1257. TEST_c_c (ccosh, plus_infty, nan_value, plus_infty, nan_value);
  1258. TEST_c_c (ccosh, minus_infty, nan_value, plus_infty, nan_value);
  1259. TEST_c_c (ccosh, 9.0, nan_value, nan_value, nan_value, INVALID_EXCEPTION_OK);
  1260. TEST_c_c (ccosh, -9.0, nan_value, nan_value, nan_value, INVALID_EXCEPTION_OK);
  1261. TEST_c_c (ccosh, nan_value, 0.0, nan_value, 0.0, IGNORE_ZERO_INF_SIGN);
  1262. TEST_c_c (ccosh, nan_value, minus_zero, nan_value, 0.0, IGNORE_ZERO_INF_SIGN);
  1263. TEST_c_c (ccosh, nan_value, 10.0, nan_value, nan_value, INVALID_EXCEPTION_OK);
  1264. TEST_c_c (ccosh, nan_value, -10.0, nan_value, nan_value, INVALID_EXCEPTION_OK);
  1265. TEST_c_c (ccosh, nan_value, plus_infty, nan_value, nan_value, INVALID_EXCEPTION_OK);
  1266. TEST_c_c (ccosh, nan_value, minus_infty, nan_value, nan_value, INVALID_EXCEPTION_OK);
  1267. TEST_c_c (ccosh, nan_value, nan_value, nan_value, nan_value);
  1268. TEST_c_c (ccosh, 0.75L, 1.25L, 0.408242591877968807788852146397499084L, 0.780365930845853240391326216300863152L);
  1269. TEST_c_c (ccosh, -2, -3, -3.72454550491532256547397070325597253L, 0.511822569987384608834463849801875634L);
  1270. END (ccosh, complex);
  1271. }
  1272. #endif
  1273. static void
  1274. ceil_test (void)
  1275. {
  1276. START (ceil);
  1277. TEST_f_f (ceil, 0.0, 0.0);
  1278. TEST_f_f (ceil, minus_zero, minus_zero);
  1279. TEST_f_f (ceil, plus_infty, plus_infty);
  1280. TEST_f_f (ceil, minus_infty, minus_infty);
  1281. TEST_f_f (ceil, nan_value, nan_value);
  1282. TEST_f_f (ceil, M_PIl, 4.0);
  1283. TEST_f_f (ceil, -M_PIl, -3.0);
  1284. TEST_f_f (ceil, 0.25, 1.0);
  1285. TEST_f_f (ceil, -0.25, minus_zero);
  1286. #ifdef TEST_LDOUBLE
  1287. /* The result can only be represented in long double. */
  1288. TEST_f_f (ceil, 4503599627370495.5L, 4503599627370496.0L);
  1289. TEST_f_f (ceil, 4503599627370496.25L, 4503599627370497.0L);
  1290. TEST_f_f (ceil, 4503599627370496.5L, 4503599627370497.0L);
  1291. TEST_f_f (ceil, 4503599627370496.75L, 4503599627370497.0L);
  1292. TEST_f_f (ceil, 4503599627370497.5L, 4503599627370498.0L);
  1293. TEST_f_f (ceil, -4503599627370495.5L, -4503599627370495.0L);
  1294. TEST_f_f (ceil, -4503599627370496.25L, -4503599627370496.0L);
  1295. TEST_f_f (ceil, -4503599627370496.5L, -4503599627370496.0L);
  1296. TEST_f_f (ceil, -4503599627370496.75L, -4503599627370496.0L);
  1297. TEST_f_f (ceil, -4503599627370497.5L, -4503599627370497.0L);
  1298. TEST_f_f (ceil, 9007199254740991.5L, 9007199254740992.0L);
  1299. TEST_f_f (ceil, 9007199254740992.25L, 9007199254740993.0L);
  1300. TEST_f_f (ceil, 9007199254740992.5L, 9007199254740993.0L);
  1301. TEST_f_f (ceil, 9007199254740992.75L, 9007199254740993.0L);
  1302. TEST_f_f (ceil, 9007199254740993.5L, 9007199254740994.0L);
  1303. TEST_f_f (ceil, -9007199254740991.5L, -9007199254740991.0L);
  1304. TEST_f_f (ceil, -9007199254740992.25L, -9007199254740992.0L);
  1305. TEST_f_f (ceil, -9007199254740992.5L, -9007199254740992.0L);
  1306. TEST_f_f (ceil, -9007199254740992.75L, -9007199254740992.0L);
  1307. TEST_f_f (ceil, -9007199254740993.5L, -9007199254740993.0L);
  1308. TEST_f_f (ceil, 72057594037927935.5L, 72057594037927936.0L);
  1309. TEST_f_f (ceil, 72057594037927936.25L, 72057594037927937.0L);
  1310. TEST_f_f (ceil, 72057594037927936.5L, 72057594037927937.0L);
  1311. TEST_f_f (ceil, 72057594037927936.75L, 72057594037927937.0L);
  1312. TEST_f_f (ceil, 72057594037927937.5L, 72057594037927938.0L);
  1313. TEST_f_f (ceil, -72057594037927935.5L, -72057594037927935.0L);
  1314. TEST_f_f (ceil, -72057594037927936.25L, -72057594037927936.0L);
  1315. TEST_f_f (ceil, -72057594037927936.5L, -72057594037927936.0L);
  1316. TEST_f_f (ceil, -72057594037927936.75L, -72057594037927936.0L);
  1317. TEST_f_f (ceil, -72057594037927937.5L, -72057594037927937.0L);
  1318. TEST_f_f (ceil, 10141204801825835211973625643007.5L, 10141204801825835211973625643008.0L);
  1319. TEST_f_f (ceil, 10141204801825835211973625643008.25L, 10141204801825835211973625643009.0L);
  1320. TEST_f_f (ceil, 10141204801825835211973625643008.5L, 10141204801825835211973625643009.0L);
  1321. TEST_f_f (ceil, 10141204801825835211973625643008.75L, 10141204801825835211973625643009.0L);
  1322. TEST_f_f (ceil, 10141204801825835211973625643009.5L, 10141204801825835211973625643010.0L);
  1323. #endif
  1324. END (ceil);
  1325. }
  1326. #if 0
  1327. static void
  1328. cexp_test (void)
  1329. {
  1330. errno = 0;
  1331. FUNC(cexp) (BUILD_COMPLEX (0, 0));
  1332. if (errno == ENOSYS)
  1333. /* Function not implemented. */
  1334. return;
  1335. START (cexp);
  1336. TEST_c_c (cexp, plus_zero, plus_zero, 1, 0.0);
  1337. TEST_c_c (cexp, minus_zero, plus_zero, 1, 0.0);
  1338. TEST_c_c (cexp, plus_zero, minus_zero, 1, minus_zero);
  1339. TEST_c_c (cexp, minus_zero, minus_zero, 1, minus_zero);
  1340. TEST_c_c (cexp, plus_infty, plus_zero, plus_infty, 0.0);
  1341. TEST_c_c (cexp, plus_infty, minus_zero, plus_infty, minus_zero);
  1342. TEST_c_c (cexp, minus_infty, plus_zero, 0.0, 0.0);
  1343. TEST_c_c (cexp, minus_infty, minus_zero, 0.0, minus_zero);
  1344. TEST_c_c (cexp, 0.0, plus_infty, nan_value, nan_value, INVALID_EXCEPTION);
  1345. TEST_c_c (cexp, minus_zero, plus_infty, nan_value, nan_value, INVALID_EXCEPTION);
  1346. TEST_c_c (cexp, 0.0, minus_infty, nan_value, nan_value, INVALID_EXCEPTION);
  1347. TEST_c_c (cexp, minus_zero, minus_infty, nan_value, nan_value, INVALID_EXCEPTION);
  1348. TEST_c_c (cexp, 100.0, plus_infty, nan_value, nan_value, INVALID_EXCEPTION);
  1349. TEST_c_c (cexp, -100.0, plus_infty, nan_value, nan_value, INVALID_EXCEPTION);
  1350. TEST_c_c (cexp, 100.0, minus_infty, nan_value, nan_value, INVALID_EXCEPTION);
  1351. TEST_c_c (cexp, -100.0, minus_infty, nan_value, nan_value, INVALID_EXCEPTION);
  1352. TEST_c_c (cexp, minus_infty, 2.0, minus_zero, 0.0);
  1353. TEST_c_c (cexp, minus_infty, 4.0, minus_zero, minus_zero);
  1354. TEST_c_c (cexp, plus_infty, 2.0, minus_infty, plus_infty);
  1355. TEST_c_c (cexp, plus_infty, 4.0, minus_infty, minus_infty);
  1356. TEST_c_c (cexp, plus_infty, plus_infty, plus_infty, nan_value, INVALID_EXCEPTION|IGNORE_ZERO_INF_SIGN);
  1357. TEST_c_c (cexp, plus_infty, minus_infty, plus_infty, nan_value, INVALID_EXCEPTION|IGNORE_ZERO_INF_SIGN);
  1358. TEST_c_c (cexp, minus_infty, plus_infty, 0.0, 0.0, IGNORE_ZERO_INF_SIGN);
  1359. TEST_c_c (cexp, minus_infty, minus_infty, 0.0, minus_zero, IGNORE_ZERO_INF_SIGN);
  1360. TEST_c_c (cexp, minus_infty, nan_value, 0, 0, IGNORE_ZERO_INF_SIGN);
  1361. TEST_c_c (cexp, plus_infty, nan_value, plus_infty, nan_value);
  1362. TEST_c_c (cexp, nan_value, 0.0, nan_value, nan_value, INVALID_EXCEPTION_OK);
  1363. TEST_c_c (cexp, nan_value, 1.0, nan_value, nan_value, INVALID_EXCEPTION_OK);
  1364. TEST_c_c (cexp, nan_value, plus_infty, nan_value, nan_value, INVALID_EXCEPTION_OK);
  1365. TEST_c_c (cexp, 0, nan_value, nan_value, nan_value, INVALID_EXCEPTION_OK);
  1366. TEST_c_c (cexp, 1, nan_value, nan_value, nan_value, INVALID_EXCEPTION_OK);
  1367. TEST_c_c (cexp, nan_value, nan_value, nan_value, nan_value);
  1368. TEST_c_c (cexp, 0.75L, 1.25L, 0.667537446429131586942201977015932112L, 2.00900045494094876258347228145863909L);
  1369. TEST_c_c (cexp, -2.0, -3.0, -0.13398091492954261346140525546115575L, -0.019098516261135196432576240858800925L);
  1370. END (cexp, complex);
  1371. }
  1372. static void
  1373. cimag_test (void)
  1374. {
  1375. START (cimag);
  1376. TEST_c_f (cimag, 1.0, 0.0, 0.0);
  1377. TEST_c_f (cimag, 1.0, minus_zero, minus_zero);
  1378. TEST_c_f (cimag, 1.0, nan_value, nan_value);
  1379. TEST_c_f (cimag, nan_value, nan_value, nan_value);
  1380. TEST_c_f (cimag, 1.0, plus_infty, plus_infty);
  1381. TEST_c_f (cimag, 1.0, minus_infty, minus_infty);
  1382. TEST_c_f (cimag, 2.0, 3.0, 3.0);
  1383. END (cimag);
  1384. }
  1385. static void
  1386. clog_test (void)
  1387. {
  1388. errno = 0;
  1389. FUNC(clog) (BUILD_COMPLEX (-2, -3));
  1390. if (errno == ENOSYS)
  1391. /* Function not implemented. */
  1392. return;
  1393. START (clog);
  1394. TEST_c_c (clog, minus_zero, 0, minus_infty, M_PIl, DIVIDE_BY_ZERO_EXCEPTION);
  1395. TEST_c_c (clog, minus_zero, minus_zero, minus_infty, -M_PIl, DIVIDE_BY_ZERO_EXCEPTION);
  1396. TEST_c_c (clog, 0, 0, minus_infty, 0.0, DIVIDE_BY_ZERO_EXCEPTION);
  1397. TEST_c_c (clog, 0, minus_zero, minus_infty, minus_zero, DIVIDE_BY_ZERO_EXCEPTION);
  1398. TEST_c_c (clog, minus_infty, plus_infty, plus_infty, M_PI_34l);
  1399. TEST_c_c (clog, minus_infty, minus_infty, plus_infty, -M_PI_34l);
  1400. TEST_c_c (clog, plus_infty, plus_infty, plus_infty, M_PI_4l);
  1401. TEST_c_c (clog, plus_infty, minus_infty, plus_infty, -M_PI_4l);
  1402. TEST_c_c (clog, 0, plus_infty, plus_infty, M_PI_2l);
  1403. TEST_c_c (clog, 3, plus_infty, plus_infty, M_PI_2l);
  1404. TEST_c_c (clog, minus_zero, plus_infty, plus_infty, M_PI_2l);
  1405. TEST_c_c (clog, -3, plus_infty, plus_infty, M_PI_2l);
  1406. TEST_c_c (clog, 0, minus_infty, plus_infty, -M_PI_2l);
  1407. TEST_c_c (clog, 3, minus_infty, plus_infty, -M_PI_2l);
  1408. TEST_c_c (clog, minus_zero, minus_infty, plus_infty, -M_PI_2l);
  1409. TEST_c_c (clog, -3, minus_infty, plus_infty, -M_PI_2l);
  1410. TEST_c_c (clog, minus_infty, 0, plus_infty, M_PIl);
  1411. TEST_c_c (clog, minus_infty, 1, plus_infty, M_PIl);
  1412. TEST_c_c (clog, minus_infty, minus_zero, plus_infty, -M_PIl);
  1413. TEST_c_c (clog, minus_infty, -1, plus_infty, -M_PIl);
  1414. TEST_c_c (clog, plus_infty, 0, plus_infty, 0.0);
  1415. TEST_c_c (clog, plus_infty, 1, plus_infty, 0.0);
  1416. TEST_c_c (clog, plus_infty, minus_zero, plus_infty, minus_zero);
  1417. TEST_c_c (clog, plus_infty, -1, plus_infty, minus_zero);
  1418. TEST_c_c (clog, plus_infty, nan_value, plus_infty, nan_value);
  1419. TEST_c_c (clog, minus_infty, nan_value, plus_infty, nan_value);
  1420. TEST_c_c (clog, nan_value, plus_infty, plus_infty, nan_value);
  1421. TEST_c_c (clog, nan_value, minus_infty, plus_infty, nan_value);
  1422. TEST_c_c (clog, 0, nan_value, nan_value, nan_value, INVALID_EXCEPTION_OK);
  1423. TEST_c_c (clog, 3, nan_value, nan_value, nan_value, INVALID_EXCEPTION_OK);
  1424. TEST_c_c (clog, minus_zero, nan_value, nan_value, nan_value, INVALID_EXCEPTION_OK);
  1425. TEST_c_c (clog, -3, nan_value, nan_value, nan_value, INVALID_EXCEPTION_OK);
  1426. TEST_c_c (clog, nan_value, 0, nan_value, nan_value, INVALID_EXCEPTION_OK);
  1427. TEST_c_c (clog, nan_value, 5, nan_value, nan_value, INVALID_EXCEPTION_OK);
  1428. TEST_c_c (clog, nan_value, minus_zero, nan_value, nan_value, INVALID_EXCEPTION_OK);
  1429. TEST_c_c (clog, nan_value, -5, nan_value, nan_value, INVALID_EXCEPTION_OK);
  1430. TEST_c_c (clog, nan_value, nan_value, nan_value, nan_value);
  1431. TEST_c_c (clog, 0.75L, 1.25L, 0.376885901188190075998919126749298416L, 1.03037682652431246378774332703115153L);
  1432. TEST_c_c (clog, -2, -3, 1.2824746787307683680267437207826593L, -2.1587989303424641704769327722648368L);
  1433. END (clog, complex);
  1434. }
  1435. static void
  1436. clog10_test (void)
  1437. {
  1438. errno = 0;
  1439. FUNC(clog10) (BUILD_COMPLEX (0.7L, 1.2L));
  1440. if (errno == ENOSYS)
  1441. /* Function not implemented. */
  1442. return;
  1443. START (clog10);
  1444. TEST_c_c (clog10, minus_zero, 0, minus_infty, M_PIl, DIVIDE_BY_ZERO_EXCEPTION);
  1445. TEST_c_c (clog10, minus_zero, minus_zero, minus_infty, -M_PIl, DIVIDE_BY_ZERO_EXCEPTION);
  1446. TEST_c_c (clog10, 0, 0, minus_infty, 0.0, DIVIDE_BY_ZERO_EXCEPTION);
  1447. TEST_c_c (clog10, 0, minus_zero, minus_infty, minus_zero, DIVIDE_BY_ZERO_EXCEPTION);
  1448. TEST_c_c (clog10, minus_infty, plus_infty, plus_infty, M_PI_34_LOG10El);
  1449. TEST_c_c (clog10, plus_infty, plus_infty, plus_infty, M_PI4_LOG10El);
  1450. TEST_c_c (clog10, plus_infty, minus_infty, plus_infty, -M_PI4_LOG10El);
  1451. TEST_c_c (clog10, 0, plus_infty, plus_infty, M_PI2_LOG10El);
  1452. TEST_c_c (clog10, 3, plus_infty, plus_infty, M_PI2_LOG10El);
  1453. TEST_c_c (clog10, minus_zero, plus_infty, plus_infty, M_PI2_LOG10El);
  1454. TEST_c_c (clog10, -3, plus_infty, plus_infty, M_PI2_LOG10El);
  1455. TEST_c_c (clog10, 0, minus_infty, plus_infty, -M_PI2_LOG10El);
  1456. TEST_c_c (clog10, 3, minus_infty, plus_infty, -M_PI2_LOG10El);
  1457. TEST_c_c (clog10, minus_zero, minus_infty, plus_infty, -M_PI2_LOG10El);
  1458. TEST_c_c (clog10, -3, minus_infty, plus_infty, -M_PI2_LOG10El);
  1459. TEST_c_c (clog10, minus_infty, 0, plus_infty, M_PI_LOG10El);
  1460. TEST_c_c (clog10, minus_infty, 1, plus_infty, M_PI_LOG10El);
  1461. TEST_c_c (clog10, minus_infty, minus_zero, plus_infty, -M_PI_LOG10El);
  1462. TEST_c_c (clog10, minus_infty, -1, plus_infty, -M_PI_LOG10El);
  1463. TEST_c_c (clog10, plus_infty, 0, plus_infty, 0.0);
  1464. TEST_c_c (clog10, plus_infty, 1, plus_infty, 0.0);
  1465. TEST_c_c (clog10, plus_infty, minus_zero, plus_infty, minus_zero);
  1466. TEST_c_c (clog10, plus_infty, -1, plus_infty, minus_zero);
  1467. TEST_c_c (clog10, plus_infty, nan_value, plus_infty, nan_value);
  1468. TEST_c_c (clog10, minus_infty, nan_value, plus_infty, nan_value);
  1469. TEST_c_c (clog10, nan_value, plus_infty, plus_infty, nan_value);
  1470. TEST_c_c (clog10, nan_value, minus_infty, plus_infty, nan_value);
  1471. TEST_c_c (clog10, 0, nan_value, nan_value, nan_value, INVALID_EXCEPTION_OK);
  1472. TEST_c_c (clog10, 3, nan_value, nan_value, nan_value, INVALID_EXCEPTION_OK);
  1473. TEST_c_c (clog10, minus_zero, nan_value, nan_value, nan_value, INVALID_EXCEPTION_OK);
  1474. TEST_c_c (clog10, -3, nan_value, nan_value, nan_value, INVALID_EXCEPTION_OK);
  1475. TEST_c_c (clog10, nan_value, 0, nan_value, nan_value, INVALID_EXCEPTION_OK);
  1476. TEST_c_c (clog10, nan_value, 5, nan_value, nan_value, INVALID_EXCEPTION_OK);
  1477. TEST_c_c (clog10, nan_value, minus_zero, nan_value, nan_value, INVALID_EXCEPTION_OK);
  1478. TEST_c_c (clog10, nan_value, -5, nan_value, nan_value, INVALID_EXCEPTION_OK);
  1479. TEST_c_c (clog10, nan_value, nan_value, nan_value, nan_value);
  1480. TEST_c_c (clog10, 0.75L, 1.25L, 0.163679467193165171449476605077428975L, 0.447486970040493067069984724340855636L);
  1481. TEST_c_c (clog10, -2, -3, 0.556971676153418384603252578971164214L, -0.937554462986374708541507952140189646L);
  1482. END (clog10, complex);
  1483. }
  1484. #endif
  1485. #if 0
  1486. static void
  1487. conj_test (void)
  1488. {
  1489. START (conj);
  1490. TEST_c_c (conj, 0.0, 0.0, 0.0, minus_zero);
  1491. TEST_c_c (conj, 0.0, minus_zero, 0.0, 0.0);
  1492. TEST_c_c (conj, nan_value, nan_value, nan_value, nan_value);
  1493. TEST_c_c (conj, plus_infty, minus_infty, plus_infty, plus_infty);
  1494. TEST_c_c (conj, plus_infty, plus_infty, plus_infty, minus_infty);
  1495. TEST_c_c (conj, 1.0, 2.0, 1.0, -2.0);
  1496. TEST_c_c (conj, 3.0, -4.0, 3.0, 4.0);
  1497. END (conj, complex);
  1498. }
  1499. #endif
  1500. static void
  1501. copysign_test (void)
  1502. {
  1503. START (copysign);
  1504. TEST_ff_f (copysign, 0, 4, 0);
  1505. TEST_ff_f (copysign, 0, -4, minus_zero);
  1506. TEST_ff_f (copysign, minus_zero, 4, 0);
  1507. TEST_ff_f (copysign, minus_zero, -4, minus_zero);
  1508. TEST_ff_f (copysign, plus_infty, 0, plus_infty);
  1509. TEST_ff_f (copysign, plus_infty, minus_zero, minus_infty);
  1510. TEST_ff_f (copysign, minus_infty, 0, plus_infty);
  1511. TEST_ff_f (copysign, minus_infty, minus_zero, minus_infty);
  1512. TEST_ff_f (copysign, 0, plus_infty, 0);
  1513. TEST_ff_f (copysign, 0, minus_zero, minus_zero);
  1514. TEST_ff_f (copysign, minus_zero, plus_infty, 0);
  1515. TEST_ff_f (copysign, minus_zero, minus_zero, minus_zero);
  1516. /* XXX More correctly we would have to check the sign of the NaN. */
  1517. TEST_ff_f (copysign, nan_value, 0, nan_value);
  1518. TEST_ff_f (copysign, nan_value, minus_zero, nan_value);
  1519. TEST_ff_f (copysign, -nan_value, 0, nan_value);
  1520. TEST_ff_f (copysign, -nan_value, minus_zero, nan_value);
  1521. END (copysign);
  1522. }
  1523. static void
  1524. cos_test (void)
  1525. {
  1526. errno = 0;
  1527. FUNC(cos) (0);
  1528. if (errno == ENOSYS)
  1529. /* Function not implemented. */
  1530. return;
  1531. START (cos);
  1532. TEST_f_f (cos, 0, 1);
  1533. TEST_f_f (cos, minus_zero, 1);
  1534. TEST_f_f (cos, plus_infty, nan_value, INVALID_EXCEPTION);
  1535. TEST_f_f (cos, minus_infty, nan_value, INVALID_EXCEPTION);
  1536. TEST_f_f (cos, nan_value, nan_value);
  1537. TEST_f_f (cos, M_PI_6l * 2.0, 0.5);
  1538. TEST_f_f (cos, M_PI_6l * 4.0, -0.5);
  1539. TEST_f_f (cos, M_PI_2l, 0);
  1540. TEST_f_f (cos, 0.75L, 0.731688868873820886311838753000084544L);
  1541. #ifdef TEST_DOUBLE
  1542. TEST_f_f (cos, 0.80190127184058835, 0.69534156199418473);
  1543. #endif
  1544. END (cos);
  1545. }
  1546. static void
  1547. cosh_test (void)
  1548. {
  1549. errno = 0;
  1550. FUNC(cosh) (0.7L);
  1551. if (errno == ENOSYS)
  1552. /* Function not implemented. */
  1553. return;
  1554. START (cosh);
  1555. TEST_f_f (cosh, 0, 1);
  1556. TEST_f_f (cosh, minus_zero, 1);
  1557. #ifndef TEST_INLINE
  1558. TEST_f_f (cosh, plus_infty, plus_infty);
  1559. TEST_f_f (cosh, minus_infty, plus_infty);
  1560. #endif
  1561. TEST_f_f (cosh, nan_value, nan_value);
  1562. TEST_f_f (cosh, 0.75L, 1.29468328467684468784170818539018176L);
  1563. END (cosh);
  1564. }
  1565. #if 0
  1566. static void
  1567. cpow_test (void)
  1568. {
  1569. errno = 0;
  1570. FUNC(cpow) (BUILD_COMPLEX (1, 0), BUILD_COMPLEX (0, 0));
  1571. if (errno == ENOSYS)
  1572. /* Function not implemented. */
  1573. return;
  1574. START (cpow);
  1575. TEST_cc_c (cpow, 1, 0, 0, 0, 1.0, 0.0);
  1576. TEST_cc_c (cpow, 2, 0, 10, 0, 1024.0, 0.0);
  1577. TEST_cc_c (cpow, M_El, 0, 0, 2 * M_PIl, 1.0, 0.0);
  1578. TEST_cc_c (cpow, 2, 3, 4, 0, -119.0, -120.0);
  1579. TEST_cc_c (cpow, nan_value, nan_value, nan_value, nan_value, nan_value, nan_value);
  1580. TEST_cc_c (cpow, 0.75L, 1.25L, 0.75L, 1.25L, 0.117506293914473555420279832210420483L, 0.346552747708338676483025352060418001L);
  1581. TEST_cc_c (cpow, 0.75L, 1.25L, 1.0L, 1.0L, 0.0846958290317209430433805274189191353L, 0.513285749182902449043287190519090481L);
  1582. TEST_cc_c (cpow, 0.75L, 1.25L, 1.0L, 0.0L, 0.75L, 1.25L);
  1583. TEST_cc_c (cpow, 0.75L, 1.25L, 0.0L, 1.0L, 0.331825439177608832276067945276730566L, 0.131338600281188544930936345230903032L);
  1584. END (cpow, complex);
  1585. }
  1586. static void
  1587. cproj_test (void)
  1588. {
  1589. START (cproj);
  1590. TEST_c_c (cproj, 0.0, 0.0, 0.0, 0.0);
  1591. TEST_c_c (cproj, minus_zero, minus_zero, minus_zero, minus_zero);
  1592. TEST_c_c (cproj, 0.0, minus_zero, 0.0, minus_zero);
  1593. TEST_c_c (cproj, minus_zero, 0.0, minus_zero, 0.0);
  1594. TEST_c_c (cproj, nan_value, nan_value, nan_value, nan_value);
  1595. TEST_c_c (cproj, plus_infty, plus_infty, plus_infty, 0.0);
  1596. TEST_c_c (cproj, plus_infty, minus_infty, plus_infty, minus_zero);
  1597. TEST_c_c (cproj, minus_infty, plus_infty, plus_infty, 0.0);
  1598. TEST_c_c (cproj, minus_infty, minus_infty, plus_infty, minus_zero);
  1599. TEST_c_c (cproj, 1.0, 0.0, 1.0, 0.0);
  1600. TEST_c_c (cproj, 2.0, 3.0, 0.2857142857142857142857142857142857L, 0.42857142857142857142857142857142855L);
  1601. END (cproj, complex);
  1602. }
  1603. static void
  1604. creal_test (void)
  1605. {
  1606. START (creal);
  1607. TEST_c_f (creal, 0.0, 1.0, 0.0);
  1608. TEST_c_f (creal, minus_zero, 1.0, minus_zero);
  1609. TEST_c_f (creal, nan_value, 1.0, nan_value);
  1610. TEST_c_f (creal, nan_value, nan_value, nan_value);
  1611. TEST_c_f (creal, plus_infty, 1.0, plus_infty);
  1612. TEST_c_f (creal, minus_infty, 1.0, minus_infty);
  1613. TEST_c_f (creal, 2.0, 3.0, 2.0);
  1614. END (creal);
  1615. }
  1616. static void
  1617. csin_test (void)
  1618. {
  1619. errno = 0;
  1620. FUNC(csin) (BUILD_COMPLEX (0.7L, 1.2L));
  1621. if (errno == ENOSYS)
  1622. /* Function not implemented. */
  1623. return;
  1624. START (csin);
  1625. TEST_c_c (csin, 0.0, 0.0, 0.0, 0.0);
  1626. TEST_c_c (csin, minus_zero, 0.0, minus_zero, 0.0);
  1627. TEST_c_c (csin, 0.0, minus_zero, 0, minus_zero);
  1628. TEST_c_c (csin, minus_zero, minus_zero, minus_zero, minus_zero);
  1629. TEST_c_c (csin, 0.0, plus_infty, 0.0, plus_infty);
  1630. TEST_c_c (csin, minus_zero, plus_infty, minus_zero, plus_infty);
  1631. TEST_c_c (csin, 0.0, minus_infty, 0.0, minus_infty);
  1632. TEST_c_c (csin, minus_zero, minus_infty, minus_zero, minus_infty);
  1633. TEST_c_c (csin, plus_infty, 0.0, nan_value, 0.0, INVALID_EXCEPTION|IGNORE_ZERO_INF_SIGN);
  1634. TEST_c_c (csin, minus_infty, 0.0, nan_value, 0.0, INVALID_EXCEPTION|IGNORE_ZERO_INF_SIGN);
  1635. TEST_c_c (csin, plus_infty, minus_zero, nan_value, 0.0, INVALID_EXCEPTION|IGNORE_ZERO_INF_SIGN);
  1636. TEST_c_c (csin, minus_infty, minus_zero, nan_value, 0.0, INVALID_EXCEPTION|IGNORE_ZERO_INF_SIGN);
  1637. TEST_c_c (csin, plus_infty, plus_infty, nan_value, plus_infty, INVALID_EXCEPTION|IGNORE_ZERO_INF_SIGN);
  1638. TEST_c_c (csin, minus_infty, plus_infty, nan_value, plus_infty, INVALID_EXCEPTION|IGNORE_ZERO_INF_SIGN);
  1639. TEST_c_c (csin, plus_infty, minus_infty, nan_value, plus_infty, INVALID_EXCEPTION|IGNORE_ZERO_INF_SIGN);
  1640. TEST_c_c (csin, minus_infty, minus_infty, nan_value, plus_infty, INVALID_EXCEPTION|IGNORE_ZERO_INF_SIGN);
  1641. TEST_c_c (csin, plus_infty, 6.75, nan_value, nan_value, INVALID_EXCEPTION);
  1642. TEST_c_c (csin, plus_infty, -6.75, nan_value, nan_value, INVALID_EXCEPTION);
  1643. TEST_c_c (csin, minus_infty, 6.75, nan_value, nan_value, INVALID_EXCEPTION);
  1644. TEST_c_c (csin, minus_infty, -6.75, nan_value, nan_value, INVALID_EXCEPTION);
  1645. TEST_c_c (csin, 4.625, plus_infty, minus_infty, minus_infty);
  1646. TEST_c_c (csin, 4.625, minus_infty, minus_infty, plus_infty);
  1647. TEST_c_c (csin, -4.625, plus_infty, plus_infty, minus_infty);
  1648. TEST_c_c (csin, -4.625, minus_infty, plus_infty, plus_infty);
  1649. TEST_c_c (csin, nan_value, 0.0, nan_value, 0.0, IGNORE_ZERO_INF_SIGN);
  1650. TEST_c_c (csin, nan_value, minus_zero, nan_value, 0.0, IGNORE_ZERO_INF_SIGN);
  1651. TEST_c_c (csin, nan_value, plus_infty, nan_value, plus_infty, IGNORE_ZERO_INF_SIGN);
  1652. TEST_c_c (csin, nan_value, minus_infty, nan_value, plus_infty, IGNORE_ZERO_INF_SIGN);
  1653. TEST_c_c (csin, nan_value, 9.0, nan_value, nan_value, INVALID_EXCEPTION_OK);
  1654. TEST_c_c (csin, nan_value, -9.0, nan_value, nan_value, INVALID_EXCEPTION_OK);
  1655. TEST_c_c (csin, 0.0, nan_value, 0.0, nan_value);
  1656. TEST_c_c (csin, minus_zero, nan_value, minus_zero, nan_value);
  1657. TEST_c_c (csin, 10.0, nan_value, nan_value, nan_value, INVALID_EXCEPTION_OK);
  1658. TEST_c_c (csin, nan_value, -10.0, nan_value, nan_value, INVALID_EXCEPTION_OK);
  1659. TEST_c_c (csin, plus_infty, nan_value, nan_value, nan_value, INVALID_EXCEPTION_OK);
  1660. TEST_c_c (csin, minus_infty, nan_value, nan_value, nan_value, INVALID_EXCEPTION_OK);
  1661. TEST_c_c (csin, nan_value, nan_value, nan_value, nan_value);
  1662. TEST_c_c (csin, 0.75L, 1.25L, 1.28722291002649188575873510790565441L, 1.17210635989270256101081285116138863L);
  1663. TEST_c_c (csin, -2, -3, -9.15449914691142957346729954460983256L, 4.16890695996656435075481305885375484L);
  1664. END (csin, complex);
  1665. }
  1666. static void
  1667. csinh_test (void)
  1668. {
  1669. errno = 0;
  1670. FUNC(csinh) (BUILD_COMPLEX (0.7L, 1.2L));
  1671. if (errno == ENOSYS)
  1672. /* Function not implemented. */
  1673. return;
  1674. START (csinh);
  1675. TEST_c_c (csinh, 0.0, 0.0, 0.0, 0.0);
  1676. TEST_c_c (csinh, minus_zero, 0.0, minus_zero, 0.0);
  1677. TEST_c_c (csinh, 0.0, minus_zero, 0.0, minus_zero);
  1678. TEST_c_c (csinh, minus_zero, minus_zero, minus_zero, minus_zero);
  1679. TEST_c_c (csinh, 0.0, plus_infty, 0.0, nan_value, INVALID_EXCEPTION|IGNORE_ZERO_INF_SIGN);
  1680. TEST_c_c (csinh, minus_zero, plus_infty, 0.0, nan_value, INVALID_EXCEPTION|IGNORE_ZERO_INF_SIGN);
  1681. TEST_c_c (csinh, 0.0, minus_infty, 0.0, nan_value, INVALID_EXCEPTION|IGNORE_ZERO_INF_SIGN);
  1682. TEST_c_c (csinh, minus_zero, minus_infty, 0.0, nan_value, INVALID_EXCEPTION|IGNORE_ZERO_INF_SIGN);
  1683. TEST_c_c (csinh, plus_infty, 0.0, plus_infty, 0.0);
  1684. TEST_c_c (csinh, minus_infty, 0.0, minus_infty, 0.0);
  1685. TEST_c_c (csinh, plus_infty, minus_zero, plus_infty, minus_zero);
  1686. TEST_c_c (csinh, minus_infty, minus_zero, minus_infty, minus_zero);
  1687. TEST_c_c (csinh, plus_infty, plus_infty, plus_infty, nan_value, INVALID_EXCEPTION|IGNORE_ZERO_INF_SIGN);
  1688. TEST_c_c (csinh, minus_infty, plus_infty, plus_infty, nan_value, INVALID_EXCEPTION|IGNORE_ZERO_INF_SIGN);
  1689. TEST_c_c (csinh, plus_infty, minus_infty, plus_infty, nan_value, INVALID_EXCEPTION|IGNORE_ZERO_INF_SIGN);
  1690. TEST_c_c (csinh, minus_infty, minus_infty, plus_infty, nan_value, INVALID_EXCEPTION|IGNORE_ZERO_INF_SIGN);
  1691. TEST_c_c (csinh, plus_infty, 4.625, minus_infty, minus_infty);
  1692. TEST_c_c (csinh, minus_infty, 4.625, plus_infty, minus_infty);
  1693. TEST_c_c (csinh, plus_infty, -4.625, minus_infty, plus_infty);
  1694. TEST_c_c (csinh, minus_infty, -4.625, plus_infty, plus_infty);
  1695. TEST_c_c (csinh, 6.75, plus_infty, nan_value, nan_value, INVALID_EXCEPTION);
  1696. TEST_c_c (csinh, -6.75, plus_infty, nan_value, nan_value, INVALID_EXCEPTION);
  1697. TEST_c_c (csinh, 6.75, minus_infty, nan_value, nan_value, INVALID_EXCEPTION);
  1698. TEST_c_c (csinh, -6.75, minus_infty, nan_value, nan_value, INVALID_EXCEPTION);
  1699. TEST_c_c (csinh, 0.0, nan_value, 0.0, nan_value, IGNORE_ZERO_INF_SIGN);
  1700. TEST_c_c (csinh, minus_zero, nan_value, 0.0, nan_value, IGNORE_ZERO_INF_SIGN);
  1701. TEST_c_c (csinh, plus_infty, nan_value, plus_infty, nan_value, IGNORE_ZERO_INF_SIGN);
  1702. TEST_c_c (csinh, minus_infty, nan_value, plus_infty, nan_value, IGNORE_ZERO_INF_SIGN);
  1703. TEST_c_c (csinh, 9.0, nan_value, nan_value, nan_value, INVALID_EXCEPTION_OK);
  1704. TEST_c_c (csinh, -9.0, nan_value, nan_value, nan_value, INVALID_EXCEPTION_OK);
  1705. TEST_c_c (csinh, nan_value, 0.0, nan_value, 0.0);
  1706. TEST_c_c (csinh, nan_value, minus_zero, nan_value, minus_zero);
  1707. TEST_c_c (csinh, nan_value, 10.0, nan_value, nan_value, INVALID_EXCEPTION_OK);
  1708. TEST_c_c (csinh, nan_value, -10.0, nan_value, nan_value, INVALID_EXCEPTION_OK);
  1709. TEST_c_c (csinh, nan_value, plus_infty, nan_value, nan_value, INVALID_EXCEPTION_OK);
  1710. TEST_c_c (csinh, nan_value, minus_infty, nan_value, nan_value, INVALID_EXCEPTION_OK);
  1711. TEST_c_c (csinh, nan_value, nan_value, nan_value, nan_value);
  1712. TEST_c_c (csinh, 0.75L, 1.25L, 0.259294854551162779153349830618433028L, 1.22863452409509552219214606515777594L);
  1713. TEST_c_c (csinh, -2, -3, 3.59056458998577995201256544779481679L, -0.530921086248519805267040090660676560L);
  1714. END (csinh, complex);
  1715. }
  1716. static void
  1717. csqrt_test (void)
  1718. {
  1719. errno = 0;
  1720. FUNC(csqrt) (BUILD_COMPLEX (-1, 0));
  1721. if (errno == ENOSYS)
  1722. /* Function not implemented. */
  1723. return;
  1724. START (csqrt);
  1725. TEST_c_c (csqrt, 0, 0, 0.0, 0.0);
  1726. TEST_c_c (csqrt, 0, minus_zero, 0, minus_zero);
  1727. TEST_c_c (csqrt, minus_zero, 0, 0.0, 0.0);
  1728. TEST_c_c (csqrt, minus_zero, minus_zero, 0.0, minus_zero);
  1729. TEST_c_c (csqrt, minus_infty, 0, 0.0, plus_infty);
  1730. TEST_c_c (csqrt, minus_infty, 6, 0.0, plus_infty);
  1731. TEST_c_c (csqrt, minus_infty, minus_zero, 0.0, minus_infty);
  1732. TEST_c_c (csqrt, minus_infty, -6, 0.0, minus_infty);
  1733. TEST_c_c (csqrt, plus_infty, 0, plus_infty, 0.0);
  1734. TEST_c_c (csqrt, plus_infty, 6, plus_infty, 0.0);
  1735. TEST_c_c (csqrt, plus_infty, minus_zero, plus_infty, minus_zero);
  1736. TEST_c_c (csqrt, plus_infty, -6, plus_infty, minus_zero);
  1737. TEST_c_c (csqrt, 0, plus_infty, plus_infty, plus_infty);
  1738. TEST_c_c (csqrt, 4, plus_infty, plus_infty, plus_infty);
  1739. TEST_c_c (csqrt, plus_infty, plus_infty, plus_infty, plus_infty);
  1740. TEST_c_c (csqrt, minus_zero, plus_infty, plus_infty, plus_infty);
  1741. TEST_c_c (csqrt, -4, plus_infty, plus_infty, plus_infty);
  1742. TEST_c_c (csqrt, minus_infty, plus_infty, plus_infty, plus_infty);
  1743. TEST_c_c (csqrt, 0, minus_infty, plus_infty, minus_infty);
  1744. TEST_c_c (csqrt, 4, minus_infty, plus_infty, minus_infty);
  1745. TEST_c_c (csqrt, plus_infty, minus_infty, plus_infty, minus_infty);
  1746. TEST_c_c (csqrt, minus_zero, minus_infty, plus_infty, minus_infty);
  1747. TEST_c_c (csqrt, -4, minus_infty, plus_infty, minus_infty);
  1748. TEST_c_c (csqrt, minus_infty, minus_infty, plus_infty, minus_infty);
  1749. TEST_c_c (csqrt, minus_infty, nan_value, nan_value, plus_infty, IGNORE_ZERO_INF_SIGN);
  1750. TEST_c_c (csqrt, plus_infty, nan_value, plus_infty, nan_value);
  1751. TEST_c_c (csqrt, 0, nan_value, nan_value, nan_value, INVALID_EXCEPTION_OK);
  1752. TEST_c_c (csqrt, 1, nan_value, nan_value, nan_value, INVALID_EXCEPTION_OK);
  1753. TEST_c_c (csqrt, minus_zero, nan_value, nan_value, nan_value, INVALID_EXCEPTION_OK);
  1754. TEST_c_c (csqrt, -1, nan_value, nan_value, nan_value, INVALID_EXCEPTION_OK);
  1755. TEST_c_c (csqrt, nan_value, 0, nan_value, nan_value, INVALID_EXCEPTION_OK);
  1756. TEST_c_c (csqrt, nan_value, 8, nan_value, nan_value, INVALID_EXCEPTION_OK);
  1757. TEST_c_c (csqrt, nan_value, minus_zero, nan_value, nan_value, INVALID_EXCEPTION_OK);
  1758. TEST_c_c (csqrt, nan_value, -8, nan_value, nan_value, INVALID_EXCEPTION_OK);
  1759. TEST_c_c (csqrt, nan_value, nan_value, nan_value, nan_value);
  1760. TEST_c_c (csqrt, 16.0, -30.0, 5.0, -3.0);
  1761. TEST_c_c (csqrt, -1, 0, 0.0, 1.0);
  1762. TEST_c_c (csqrt, 0, 2, 1.0, 1.0);
  1763. TEST_c_c (csqrt, 119, 120, 12.0, 5.0);
  1764. TEST_c_c (csqrt, 0.75L, 1.25L, 1.05065169626078392338656675760808326L, 0.594868882070379067881984030639932657L);
  1765. TEST_c_c (csqrt, -2, -3, 0.89597747612983812471573375529004348L, -1.6741492280355400404480393008490519L);
  1766. TEST_c_c (csqrt, -2, 3, 0.89597747612983812471573375529004348L, 1.6741492280355400404480393008490519L);
  1767. /* Principal square root should be returned (i.e., non-negative real
  1768. part). */
  1769. TEST_c_c (csqrt, 0, -1, M_SQRT_2_2, -M_SQRT_2_2);
  1770. END (csqrt, complex);
  1771. }
  1772. static void
  1773. ctan_test (void)
  1774. {
  1775. errno = 0;
  1776. FUNC(ctan) (BUILD_COMPLEX (0.7L, 1.2L));
  1777. if (errno == ENOSYS)
  1778. /* Function not implemented. */
  1779. return;
  1780. START (ctan);
  1781. TEST_c_c (ctan, 0, 0, 0.0, 0.0);
  1782. TEST_c_c (ctan, 0, minus_zero, 0.0, minus_zero);
  1783. TEST_c_c (ctan, minus_zero, 0, minus_zero, 0.0);
  1784. TEST_c_c (ctan, minus_zero, minus_zero, minus_zero, minus_zero);
  1785. TEST_c_c (ctan, 0, plus_infty, 0.0, 1.0);
  1786. TEST_c_c (ctan, 1, plus_infty, 0.0, 1.0);
  1787. TEST_c_c (ctan, minus_zero, plus_infty, minus_zero, 1.0);
  1788. TEST_c_c (ctan, -1, plus_infty, minus_zero, 1.0);
  1789. TEST_c_c (ctan, 0, minus_infty, 0.0, -1.0);
  1790. TEST_c_c (ctan, 1, minus_infty, 0.0, -1.0);
  1791. TEST_c_c (ctan, minus_zero, minus_infty, minus_zero, -1.0);
  1792. TEST_c_c (ctan, -1, minus_infty, minus_zero, -1.0);
  1793. TEST_c_c (ctan, plus_infty, 0, nan_value, nan_value, INVALID_EXCEPTION);
  1794. TEST_c_c (ctan, plus_infty, 2, nan_value, nan_value, INVALID_EXCEPTION);
  1795. TEST_c_c (ctan, minus_infty, 0, nan_value, nan_value, INVALID_EXCEPTION);
  1796. TEST_c_c (ctan, minus_infty, 2, nan_value, nan_value, INVALID_EXCEPTION);
  1797. TEST_c_c (ctan, plus_infty, minus_zero, nan_value, nan_value, INVALID_EXCEPTION);
  1798. TEST_c_c (ctan, plus_infty, -2, nan_value, nan_value, INVALID_EXCEPTION);
  1799. TEST_c_c (ctan, minus_infty, minus_zero, nan_value, nan_value, INVALID_EXCEPTION);
  1800. TEST_c_c (ctan, minus_infty, -2, nan_value, nan_value, INVALID_EXCEPTION);
  1801. TEST_c_c (ctan, nan_value, plus_infty, 0.0, 1.0, IGNORE_ZERO_INF_SIGN);
  1802. TEST_c_c (ctan, nan_value, minus_infty, 0.0, -1.0, IGNORE_ZERO_INF_SIGN);
  1803. TEST_c_c (ctan, 0, nan_value, 0.0, nan_value);
  1804. TEST_c_c (ctan, minus_zero, nan_value, minus_zero, nan_value);
  1805. TEST_c_c (ctan, 0.5, nan_value, nan_value, nan_value, INVALID_EXCEPTION_OK);
  1806. TEST_c_c (ctan, -4.5, nan_value, nan_value, nan_value, INVALID_EXCEPTION_OK);
  1807. TEST_c_c (ctan, nan_value, 0, nan_value, nan_value, INVALID_EXCEPTION_OK);
  1808. TEST_c_c (ctan, nan_value, 5, nan_value, nan_value, INVALID_EXCEPTION_OK);
  1809. TEST_c_c (ctan, nan_value, minus_zero, nan_value, nan_value, INVALID_EXCEPTION_OK);
  1810. TEST_c_c (ctan, nan_value, -0.25, nan_value, nan_value, INVALID_EXCEPTION_OK);
  1811. TEST_c_c (ctan, nan_value, nan_value, nan_value, nan_value);
  1812. TEST_c_c (ctan, 0.75L, 1.25L, 0.160807785916206426725166058173438663L, 0.975363285031235646193581759755216379L);
  1813. TEST_c_c (ctan, -2, -3, 0.376402564150424829275122113032269084e-2L, -1.00323862735360980144635859782192726L);
  1814. END (ctan, complex);
  1815. }
  1816. static void
  1817. ctanh_test (void)
  1818. {
  1819. errno = 0;
  1820. FUNC(ctanh) (BUILD_COMPLEX (0, 0));
  1821. if (errno == ENOSYS)
  1822. /* Function not implemented. */
  1823. return;
  1824. START (ctanh);
  1825. TEST_c_c (ctanh, 0, 0, 0.0, 0.0);
  1826. TEST_c_c (ctanh, 0, minus_zero, 0.0, minus_zero);
  1827. TEST_c_c (ctanh, minus_zero, 0, minus_zero, 0.0);
  1828. TEST_c_c (ctanh, minus_zero, minus_zero, minus_zero, minus_zero);
  1829. TEST_c_c (ctanh, plus_infty, 0, 1.0, 0.0);
  1830. TEST_c_c (ctanh, plus_infty, 1, 1.0, 0.0);
  1831. TEST_c_c (ctanh, plus_infty, minus_zero, 1.0, minus_zero);
  1832. TEST_c_c (ctanh, plus_infty, -1, 1.0, minus_zero);
  1833. TEST_c_c (ctanh, minus_infty, 0, -1.0, 0.0);
  1834. TEST_c_c (ctanh, minus_infty, 1, -1.0, 0.0);
  1835. TEST_c_c (ctanh, minus_infty, minus_zero, -1.0, minus_zero);
  1836. TEST_c_c (ctanh, minus_infty, -1, -1.0, minus_zero);
  1837. TEST_c_c (ctanh, 0, plus_infty, nan_value, nan_value, INVALID_EXCEPTION);
  1838. TEST_c_c (ctanh, 2, plus_infty, nan_value, nan_value, INVALID_EXCEPTION);
  1839. TEST_c_c (ctanh, 0, minus_infty, nan_value, nan_value, INVALID_EXCEPTION);
  1840. TEST_c_c (ctanh, 2, minus_infty, nan_value, nan_value, INVALID_EXCEPTION);
  1841. TEST_c_c (ctanh, minus_zero, plus_infty, nan_value, nan_value, INVALID_EXCEPTION);
  1842. TEST_c_c (ctanh, -2, plus_infty, nan_value, nan_value, INVALID_EXCEPTION);
  1843. TEST_c_c (ctanh, minus_zero, minus_infty, nan_value, nan_value, INVALID_EXCEPTION);
  1844. TEST_c_c (ctanh, -2, minus_infty, nan_value, nan_value, INVALID_EXCEPTION);
  1845. TEST_c_c (ctanh, plus_infty, nan_value, 1.0, 0.0, IGNORE_ZERO_INF_SIGN);
  1846. TEST_c_c (ctanh, minus_infty, nan_value, -1.0, 0.0, IGNORE_ZERO_INF_SIGN);
  1847. TEST_c_c (ctanh, nan_value, 0, nan_value, 0.0);
  1848. TEST_c_c (ctanh, nan_value, minus_zero, nan_value, minus_zero);
  1849. TEST_c_c (ctanh, nan_value, 0.5, nan_value, nan_value, INVALID_EXCEPTION_OK);
  1850. TEST_c_c (ctanh, nan_value, -4.5, nan_value, nan_value, INVALID_EXCEPTION_OK);
  1851. TEST_c_c (ctanh, 0, nan_value, nan_value, nan_value, INVALID_EXCEPTION_OK);
  1852. TEST_c_c (ctanh, 5, nan_value, nan_value, nan_value, INVALID_EXCEPTION_OK);
  1853. TEST_c_c (ctanh, minus_zero, nan_value, nan_value, nan_value, INVALID_EXCEPTION_OK);
  1854. TEST_c_c (ctanh, -0.25, nan_value, nan_value, nan_value, INVALID_EXCEPTION_OK);
  1855. TEST_c_c (ctanh, nan_value, nan_value, nan_value, nan_value);
  1856. TEST_c_c (ctanh, 0, M_PI_4l, 0.0, 1.0);
  1857. TEST_c_c (ctanh, 0.75L, 1.25L, 1.37260757053378320258048606571226857L, 0.385795952609750664177596760720790220L);
  1858. TEST_c_c (ctanh, -2, -3, -0.965385879022133124278480269394560686L, 0.988437503832249372031403430350121098e-2L);
  1859. END (ctanh, complex);
  1860. }
  1861. #endif
  1862. static void
  1863. erf_test (void)
  1864. {
  1865. errno = 0;
  1866. FUNC(erf) (0);
  1867. if (errno == ENOSYS)
  1868. /* Function not implemented. */
  1869. return;
  1870. START (erf);
  1871. TEST_f_f (erf, 0, 0);
  1872. TEST_f_f (erf, minus_zero, minus_zero);
  1873. TEST_f_f (erf, plus_infty, 1);
  1874. TEST_f_f (erf, minus_infty, -1);
  1875. TEST_f_f (erf, nan_value, nan_value);
  1876. TEST_f_f (erf, 0.125L, 0.140316204801333817393029446521623398L);
  1877. TEST_f_f (erf, 0.75L, 0.711155633653515131598937834591410777L);
  1878. TEST_f_f (erf, 1.25L, 0.922900128256458230136523481197281140L);
  1879. TEST_f_f (erf, 2.0L, 0.995322265018952734162069256367252929L);
  1880. TEST_f_f (erf, 4.125L, 0.999999994576599200434933994687765914L);
  1881. TEST_f_f (erf, 27.0L, 1.0L);
  1882. END (erf);
  1883. }
  1884. static void
  1885. erfc_test (void)
  1886. {
  1887. errno = 0;
  1888. FUNC(erfc) (0);
  1889. if (errno == ENOSYS)
  1890. /* Function not implemented. */
  1891. return;
  1892. START (erfc);
  1893. TEST_f_f (erfc, plus_infty, 0.0);
  1894. TEST_f_f (erfc, minus_infty, 2.0);
  1895. TEST_f_f (erfc, 0.0, 1.0);
  1896. TEST_f_f (erfc, minus_zero, 1.0);
  1897. TEST_f_f (erfc, nan_value, nan_value);
  1898. TEST_f_f (erfc, 0.125L, 0.859683795198666182606970553478376602L);
  1899. TEST_f_f (erfc, 0.75L, 0.288844366346484868401062165408589223L);
  1900. TEST_f_f (erfc, 1.25L, 0.0770998717435417698634765188027188596L);
  1901. TEST_f_f (erfc, 2.0L, 0.00467773498104726583793074363274707139L);
  1902. TEST_f_f (erfc, 4.125L, 0.542340079956506600531223408575531062e-8L);
  1903. #ifdef TEST_LDOUBLE
  1904. /* The result can only be represented in long double. */
  1905. # if LDBL_MIN_10_EXP < -319
  1906. TEST_f_f (erfc, 27.0L, 0.523704892378925568501606768284954709e-318L);
  1907. # endif
  1908. #endif
  1909. END (erfc);
  1910. }
  1911. static void
  1912. exp_test (void)
  1913. {
  1914. errno = 0;
  1915. FUNC(exp) (0);
  1916. if (errno == ENOSYS)
  1917. /* Function not implemented. */
  1918. return;
  1919. START (exp);
  1920. TEST_f_f (exp, 0, 1);
  1921. TEST_f_f (exp, minus_zero, 1);
  1922. #ifndef TEST_INLINE
  1923. TEST_f_f (exp, plus_infty, plus_infty);
  1924. TEST_f_f (exp, minus_infty, 0);
  1925. #endif
  1926. TEST_f_f (exp, nan_value, nan_value);
  1927. TEST_f_f (exp, 1, M_El);
  1928. TEST_f_f (exp, 2, M_E2l);
  1929. TEST_f_f (exp, 3, M_E3l);
  1930. TEST_f_f (exp, 0.75L, 2.11700001661267466854536981983709561L);
  1931. TEST_f_f (exp, 50.0L, 5184705528587072464087.45332293348538L);
  1932. #ifdef TEST_LDOUBLE
  1933. /* The result can only be represented in long double. */
  1934. TEST_f_f (exp, 1000.0L, 0.197007111401704699388887935224332313e435L);
  1935. #endif
  1936. END (exp);
  1937. }
  1938. #if 0
  1939. static void
  1940. exp10_test (void)
  1941. {
  1942. errno = 0;
  1943. FUNC(exp10) (0);
  1944. if (errno == ENOSYS)
  1945. /* Function not implemented. */
  1946. return;
  1947. START (exp10);
  1948. TEST_f_f (exp10, 0, 1);
  1949. TEST_f_f (exp10, minus_zero, 1);
  1950. TEST_f_f (exp10, plus_infty, plus_infty);
  1951. TEST_f_f (exp10, minus_infty, 0);
  1952. TEST_f_f (exp10, nan_value, nan_value);
  1953. TEST_f_f (exp10, 3, 1000);
  1954. TEST_f_f (exp10, -1, 0.1L);
  1955. TEST_f_f (exp10, 1e6, plus_infty);
  1956. TEST_f_f (exp10, -1e6, 0);
  1957. TEST_f_f (exp10, 0.75L, 5.62341325190349080394951039776481231L);
  1958. END (exp10);
  1959. }
  1960. static void
  1961. exp2_test (void)
  1962. {
  1963. errno = 0;
  1964. FUNC(exp2) (0);
  1965. if (errno == ENOSYS)
  1966. /* Function not implemented. */
  1967. return;
  1968. START (exp2);
  1969. TEST_f_f (exp2, 0, 1);
  1970. TEST_f_f (exp2, minus_zero, 1);
  1971. TEST_f_f (exp2, plus_infty, plus_infty);
  1972. TEST_f_f (exp2, minus_infty, 0);
  1973. TEST_f_f (exp2, nan_value, nan_value);
  1974. TEST_f_f (exp2, 10, 1024);
  1975. TEST_f_f (exp2, -1, 0.5);
  1976. TEST_f_f (exp2, 1e6, plus_infty);
  1977. TEST_f_f (exp2, -1e6, 0);
  1978. TEST_f_f (exp2, 0.75L, 1.68179283050742908606225095246642979L);
  1979. END (exp2);
  1980. }
  1981. #endif
  1982. static void
  1983. expm1_test (void)
  1984. {
  1985. errno = 0;
  1986. FUNC(expm1) (0);
  1987. if (errno == ENOSYS)
  1988. /* Function not implemented. */
  1989. return;
  1990. START (expm1);
  1991. TEST_f_f (expm1, 0, 0);
  1992. TEST_f_f (expm1, minus_zero, minus_zero);
  1993. #ifndef TEST_INLINE
  1994. TEST_f_f (expm1, plus_infty, plus_infty);
  1995. TEST_f_f (expm1, minus_infty, -1);
  1996. #endif
  1997. TEST_f_f (expm1, nan_value, nan_value);
  1998. TEST_f_f (expm1, 1, M_El - 1.0);
  1999. TEST_f_f (expm1, 0.75L, 1.11700001661267466854536981983709561L);
  2000. END (expm1);
  2001. }
  2002. static void
  2003. fabs_test (void)
  2004. {
  2005. START (fabs);
  2006. TEST_f_f (fabs, 0, 0);
  2007. TEST_f_f (fabs, minus_zero, 0);
  2008. TEST_f_f (fabs, plus_infty, plus_infty);
  2009. TEST_f_f (fabs, minus_infty, plus_infty);
  2010. TEST_f_f (fabs, nan_value, nan_value);
  2011. TEST_f_f (fabs, 38.0, 38.0);
  2012. TEST_f_f (fabs, -M_El, M_El);
  2013. END (fabs);
  2014. }
  2015. #if 0
  2016. static void
  2017. fdim_test (void)
  2018. {
  2019. START (fdim);
  2020. TEST_ff_f (fdim, 0, 0, 0);
  2021. TEST_ff_f (fdim, 9, 0, 9);
  2022. TEST_ff_f (fdim, 0, 9, 0);
  2023. TEST_ff_f (fdim, -9, 0, 0);
  2024. TEST_ff_f (fdim, 0, -9, 9);
  2025. TEST_ff_f (fdim, plus_infty, 9, plus_infty);
  2026. TEST_ff_f (fdim, plus_infty, -9, plus_infty);
  2027. TEST_ff_f (fdim, minus_infty, 9, 0);
  2028. TEST_ff_f (fdim, minus_infty, -9, 0);
  2029. TEST_ff_f (fdim, 9, minus_infty, plus_infty);
  2030. TEST_ff_f (fdim, -9, minus_infty, plus_infty);
  2031. TEST_ff_f (fdim, 9, plus_infty, 0);
  2032. TEST_ff_f (fdim, -9, plus_infty, 0);
  2033. TEST_ff_f (fdim, 0, nan_value, nan_value);
  2034. TEST_ff_f (fdim, 9, nan_value, nan_value);
  2035. TEST_ff_f (fdim, -9, nan_value, nan_value);
  2036. TEST_ff_f (fdim, nan_value, 9, nan_value);
  2037. TEST_ff_f (fdim, nan_value, -9, nan_value);
  2038. TEST_ff_f (fdim, plus_infty, nan_value, nan_value);
  2039. TEST_ff_f (fdim, minus_infty, nan_value, nan_value);
  2040. TEST_ff_f (fdim, nan_value, plus_infty, nan_value);
  2041. TEST_ff_f (fdim, nan_value, minus_infty, nan_value);
  2042. TEST_ff_f (fdim, nan_value, nan_value, nan_value);
  2043. TEST_ff_f (fdim, plus_infty, plus_infty, 0);
  2044. END (fdim);
  2045. }
  2046. #endif
  2047. static void
  2048. floor_test (void)
  2049. {
  2050. START (floor);
  2051. TEST_f_f (floor, 0.0, 0.0);
  2052. TEST_f_f (floor, minus_zero, minus_zero);
  2053. TEST_f_f (floor, plus_infty, plus_infty);
  2054. TEST_f_f (floor, minus_infty, minus_infty);
  2055. TEST_f_f (floor, nan_value, nan_value);
  2056. TEST_f_f (floor, M_PIl, 3.0);
  2057. TEST_f_f (floor, -M_PIl, -4.0);
  2058. TEST_f_f (floor, 0.25, 0.0);
  2059. TEST_f_f (floor, -0.25, -1.0);
  2060. #ifdef TEST_LDOUBLE
  2061. /* The result can only be represented in long double. */
  2062. TEST_f_f (floor, 4503599627370495.5L, 4503599627370495.0L);
  2063. TEST_f_f (floor, 4503599627370496.25L, 4503599627370496.0L);
  2064. TEST_f_f (floor, 4503599627370496.5L, 4503599627370496.0L);
  2065. TEST_f_f (floor, 4503599627370496.75L, 4503599627370496.0L);
  2066. TEST_f_f (floor, 4503599627370497.5L, 4503599627370497.0L);
  2067. TEST_f_f (floor, -4503599627370495.5L, -4503599627370496.0L);
  2068. TEST_f_f (floor, -4503599627370496.25L, -4503599627370497.0L);
  2069. TEST_f_f (floor, -4503599627370496.5L, -4503599627370497.0L);
  2070. TEST_f_f (floor, -4503599627370496.75L, -4503599627370497.0L);
  2071. TEST_f_f (floor, -4503599627370497.5L, -4503599627370498.0L);
  2072. TEST_f_f (floor, 9007199254740991.5L, 9007199254740991.0L);
  2073. TEST_f_f (floor, 9007199254740992.25L, 9007199254740992.0L);
  2074. TEST_f_f (floor, 9007199254740992.5L, 9007199254740992.0L);
  2075. TEST_f_f (floor, 9007199254740992.75L, 9007199254740992.0L);
  2076. TEST_f_f (floor, 9007199254740993.5L, 9007199254740993.0L);
  2077. TEST_f_f (floor, -9007199254740991.5L, -9007199254740992.0L);
  2078. TEST_f_f (floor, -9007199254740992.25L, -9007199254740993.0L);
  2079. TEST_f_f (floor, -9007199254740992.5L, -9007199254740993.0L);
  2080. TEST_f_f (floor, -9007199254740992.75L, -9007199254740993.0L);
  2081. TEST_f_f (floor, -9007199254740993.5L, -9007199254740994.0L);
  2082. TEST_f_f (floor, 72057594037927935.5L, 72057594037927935.0L);
  2083. TEST_f_f (floor, 72057594037927936.25L, 72057594037927936.0L);
  2084. TEST_f_f (floor, 72057594037927936.5L, 72057594037927936.0L);
  2085. TEST_f_f (floor, 72057594037927936.75L, 72057594037927936.0L);
  2086. TEST_f_f (floor, 72057594037927937.5L, 72057594037927937.0L);
  2087. TEST_f_f (floor, -72057594037927935.5L, -72057594037927936.0L);
  2088. TEST_f_f (floor, -72057594037927936.25L, -72057594037927937.0L);
  2089. TEST_f_f (floor, -72057594037927936.5L, -72057594037927937.0L);
  2090. TEST_f_f (floor, -72057594037927936.75L, -72057594037927937.0L);
  2091. TEST_f_f (floor, -72057594037927937.5L, -72057594037927938.0L);
  2092. TEST_f_f (floor, 10141204801825835211973625643007.5L, 10141204801825835211973625643007.0L);
  2093. TEST_f_f (floor, 10141204801825835211973625643008.25L, 10141204801825835211973625643008.0L);
  2094. TEST_f_f (floor, 10141204801825835211973625643008.5L, 10141204801825835211973625643008.0L);
  2095. TEST_f_f (floor, 10141204801825835211973625643008.75L, 10141204801825835211973625643008.0L);
  2096. TEST_f_f (floor, 10141204801825835211973625643009.5L, 10141204801825835211973625643009.0L);
  2097. #endif
  2098. END (floor);
  2099. }
  2100. #if 0
  2101. static void
  2102. fma_test (void)
  2103. {
  2104. START (fma);
  2105. TEST_fff_f (fma, 1.0, 2.0, 3.0, 5.0);
  2106. TEST_fff_f (fma, nan_value, 2.0, 3.0, nan_value);
  2107. TEST_fff_f (fma, 1.0, nan_value, 3.0, nan_value);
  2108. TEST_fff_f (fma, 1.0, 2.0, nan_value, nan_value, INVALID_EXCEPTION_OK);
  2109. TEST_fff_f (fma, plus_infty, 0.0, nan_value, nan_value, INVALID_EXCEPTION_OK);
  2110. TEST_fff_f (fma, minus_infty, 0.0, nan_value, nan_value, INVALID_EXCEPTION_OK);
  2111. TEST_fff_f (fma, 0.0, plus_infty, nan_value, nan_value, INVALID_EXCEPTION_OK);
  2112. TEST_fff_f (fma, 0.0, minus_infty, nan_value, nan_value, INVALID_EXCEPTION_OK);
  2113. TEST_fff_f (fma, plus_infty, 0.0, 1.0, nan_value, INVALID_EXCEPTION);
  2114. TEST_fff_f (fma, minus_infty, 0.0, 1.0, nan_value, INVALID_EXCEPTION);
  2115. TEST_fff_f (fma, 0.0, plus_infty, 1.0, nan_value, INVALID_EXCEPTION);
  2116. TEST_fff_f (fma, 0.0, minus_infty, 1.0, nan_value, INVALID_EXCEPTION);
  2117. TEST_fff_f (fma, plus_infty, plus_infty, minus_infty, nan_value, INVALID_EXCEPTION);
  2118. TEST_fff_f (fma, minus_infty, plus_infty, plus_infty, nan_value, INVALID_EXCEPTION);
  2119. TEST_fff_f (fma, plus_infty, minus_infty, plus_infty, nan_value, INVALID_EXCEPTION);
  2120. TEST_fff_f (fma, minus_infty, minus_infty, minus_infty, nan_value, INVALID_EXCEPTION);
  2121. TEST_fff_f (fma, 1.25L, 0.75L, 0.0625L, 1.0L);
  2122. END (fma);
  2123. }
  2124. static void
  2125. fmax_test (void)
  2126. {
  2127. START (fmax);
  2128. TEST_ff_f (fmax, 0, 0, 0);
  2129. TEST_ff_f (fmax, minus_zero, minus_zero, minus_zero);
  2130. TEST_ff_f (fmax, 9, 0, 9);
  2131. TEST_ff_f (fmax, 0, 9, 9);
  2132. TEST_ff_f (fmax, -9, 0, 0);
  2133. TEST_ff_f (fmax, 0, -9, 0);
  2134. TEST_ff_f (fmax, plus_infty, 9, plus_infty);
  2135. TEST_ff_f (fmax, 0, plus_infty, plus_infty);
  2136. TEST_ff_f (fmax, -9, plus_infty, plus_infty);
  2137. TEST_ff_f (fmax, plus_infty, -9, plus_infty);
  2138. TEST_ff_f (fmax, minus_infty, 9, 9);
  2139. TEST_ff_f (fmax, minus_infty, -9, -9);
  2140. TEST_ff_f (fmax, 9, minus_infty, 9);
  2141. TEST_ff_f (fmax, -9, minus_infty, -9);
  2142. TEST_ff_f (fmax, 0, nan_value, 0);
  2143. TEST_ff_f (fmax, 9, nan_value, 9);
  2144. TEST_ff_f (fmax, -9, nan_value, -9);
  2145. TEST_ff_f (fmax, nan_value, 0, 0);
  2146. TEST_ff_f (fmax, nan_value, 9, 9);
  2147. TEST_ff_f (fmax, nan_value, -9, -9);
  2148. TEST_ff_f (fmax, plus_infty, nan_value, plus_infty);
  2149. TEST_ff_f (fmax, minus_infty, nan_value, minus_infty);
  2150. TEST_ff_f (fmax, nan_value, plus_infty, plus_infty);
  2151. TEST_ff_f (fmax, nan_value, minus_infty, minus_infty);
  2152. TEST_ff_f (fmax, nan_value, nan_value, nan_value);
  2153. END (fmax);
  2154. }
  2155. static void
  2156. fmin_test (void)
  2157. {
  2158. START (fmin);
  2159. TEST_ff_f (fmin, 0, 0, 0);
  2160. TEST_ff_f (fmin, minus_zero, minus_zero, minus_zero);
  2161. TEST_ff_f (fmin, 9, 0, 0);
  2162. TEST_ff_f (fmin, 0, 9, 0);
  2163. TEST_ff_f (fmin, -9, 0, -9);
  2164. TEST_ff_f (fmin, 0, -9, -9);
  2165. TEST_ff_f (fmin, plus_infty, 9, 9);
  2166. TEST_ff_f (fmin, 9, plus_infty, 9);
  2167. TEST_ff_f (fmin, plus_infty, -9, -9);
  2168. TEST_ff_f (fmin, -9, plus_infty, -9);
  2169. TEST_ff_f (fmin, minus_infty, 9, minus_infty);
  2170. TEST_ff_f (fmin, minus_infty, -9, minus_infty);
  2171. TEST_ff_f (fmin, 9, minus_infty, minus_infty);
  2172. TEST_ff_f (fmin, -9, minus_infty, minus_infty);
  2173. TEST_ff_f (fmin, 0, nan_value, 0);
  2174. TEST_ff_f (fmin, 9, nan_value, 9);
  2175. TEST_ff_f (fmin, -9, nan_value, -9);
  2176. TEST_ff_f (fmin, nan_value, 0, 0);
  2177. TEST_ff_f (fmin, nan_value, 9, 9);
  2178. TEST_ff_f (fmin, nan_value, -9, -9);
  2179. TEST_ff_f (fmin, plus_infty, nan_value, plus_infty);
  2180. TEST_ff_f (fmin, minus_infty, nan_value, minus_infty);
  2181. TEST_ff_f (fmin, nan_value, plus_infty, plus_infty);
  2182. TEST_ff_f (fmin, nan_value, minus_infty, minus_infty);
  2183. TEST_ff_f (fmin, nan_value, nan_value, nan_value);
  2184. END (fmin);
  2185. }
  2186. #endif
  2187. static void
  2188. fmod_test (void)
  2189. {
  2190. errno = 0;
  2191. FUNC(fmod) (6.5, 2.3L);
  2192. if (errno == ENOSYS)
  2193. /* Function not implemented. */
  2194. return;
  2195. START (fmod);
  2196. /* fmod (+0, y) == +0 for y != 0. */
  2197. TEST_ff_f (fmod, 0, 3, 0);
  2198. /* fmod (-0, y) == -0 for y != 0. */
  2199. TEST_ff_f (fmod, minus_zero, 3, minus_zero);
  2200. /* fmod (+inf, y) == NaN plus invalid exception. */
  2201. TEST_ff_f (fmod, plus_infty, 3, nan_value, INVALID_EXCEPTION);
  2202. /* fmod (-inf, y) == NaN plus invalid exception. */
  2203. TEST_ff_f (fmod, minus_infty, 3, nan_value, INVALID_EXCEPTION);
  2204. /* fmod (x, +0) == NaN plus invalid exception. */
  2205. TEST_ff_f (fmod, 3, 0, nan_value, INVALID_EXCEPTION);
  2206. /* fmod (x, -0) == NaN plus invalid exception. */
  2207. TEST_ff_f (fmod, 3, minus_zero, nan_value, INVALID_EXCEPTION);
  2208. /* fmod (x, +inf) == x for x not infinite. */
  2209. TEST_ff_f (fmod, 3.0, plus_infty, 3.0);
  2210. /* fmod (x, -inf) == x for x not infinite. */
  2211. TEST_ff_f (fmod, 3.0, minus_infty, 3.0);
  2212. TEST_ff_f (fmod, nan_value, nan_value, nan_value);
  2213. TEST_ff_f (fmod, 6.5, 2.25L, 2.0L);
  2214. TEST_ff_f (fmod, -6.5, 2.25L, -2.0L);
  2215. TEST_ff_f (fmod, 6.5, -2.25L, 2.0L);
  2216. TEST_ff_f (fmod, -6.5, -2.25L, -2.0L);
  2217. END (fmod);
  2218. }
  2219. static void
  2220. fpclassify_test (void)
  2221. {
  2222. START (fpclassify);
  2223. TEST_f_i (fpclassify, nan_value, FP_NAN);
  2224. TEST_f_i (fpclassify, plus_infty, FP_INFINITE);
  2225. TEST_f_i (fpclassify, minus_infty, FP_INFINITE);
  2226. TEST_f_i (fpclassify, plus_zero, FP_ZERO);
  2227. TEST_f_i (fpclassify, minus_zero, FP_ZERO);
  2228. TEST_f_i (fpclassify, 1000, FP_NORMAL);
  2229. END (fpclassify);
  2230. }
  2231. static void
  2232. frexp_test (void)
  2233. {
  2234. int x;
  2235. START (frexp);
  2236. TEST_fI_f1 (frexp, plus_infty, plus_infty, IGNORE);
  2237. TEST_fI_f1 (frexp, minus_infty, minus_infty, IGNORE);
  2238. TEST_fI_f1 (frexp, nan_value, nan_value, IGNORE);
  2239. TEST_fI_f1 (frexp, 0.0, 0.0, 0.0);
  2240. TEST_fI_f1 (frexp, minus_zero, minus_zero, 0.0);
  2241. TEST_fI_f1 (frexp, 12.8L, 0.8L, 4);
  2242. TEST_fI_f1 (frexp, -27.34L, -0.854375L, 5);
  2243. END (frexp);
  2244. }
  2245. static void
  2246. gamma_test (void)
  2247. {
  2248. errno = 0;
  2249. FUNC(gamma) (1);
  2250. if (errno == ENOSYS)
  2251. /* Function not implemented. */
  2252. return;
  2253. feclearexcept (FE_ALL_EXCEPT);
  2254. START (gamma);
  2255. TEST_f_f (gamma, plus_infty, plus_infty);
  2256. TEST_f_f (gamma, 0, plus_infty, DIVIDE_BY_ZERO_EXCEPTION);
  2257. TEST_f_f (gamma, -3, plus_infty, DIVIDE_BY_ZERO_EXCEPTION);
  2258. TEST_f_f (gamma, minus_infty, plus_infty);
  2259. TEST_f_f (gamma, nan_value, nan_value);
  2260. TEST_f_f1 (gamma, 1, 0, 1);
  2261. TEST_f_f1 (gamma, 3, M_LN2l, 1);
  2262. TEST_f_f1 (gamma, 0.5, M_LOG_SQRT_PIl, 1);
  2263. TEST_f_f1 (gamma, -0.5, M_LOG_2_SQRT_PIl, -1);
  2264. END (gamma);
  2265. }
  2266. static void
  2267. hypot_test (void)
  2268. {
  2269. errno = 0;
  2270. FUNC(hypot) (0.7L, 12.4L);
  2271. if (errno == ENOSYS)
  2272. /* Function not implemented. */
  2273. return;
  2274. START (hypot);
  2275. TEST_ff_f (hypot, plus_infty, 1, plus_infty, IGNORE_ZERO_INF_SIGN);
  2276. TEST_ff_f (hypot, minus_infty, 1, plus_infty, IGNORE_ZERO_INF_SIGN);
  2277. #ifndef TEST_INLINE
  2278. TEST_ff_f (hypot, plus_infty, nan_value, plus_infty);
  2279. TEST_ff_f (hypot, minus_infty, nan_value, plus_infty);
  2280. TEST_ff_f (hypot, nan_value, plus_infty, plus_infty);
  2281. TEST_ff_f (hypot, nan_value, minus_infty, plus_infty);
  2282. #endif
  2283. TEST_ff_f (hypot, nan_value, nan_value, nan_value);
  2284. /* hypot (x,y) == hypot (+-x, +-y) */
  2285. TEST_ff_f (hypot, 0.7L, 12.4L, 12.419742348374220601176836866763271L);
  2286. TEST_ff_f (hypot, -0.7L, 12.4L, 12.419742348374220601176836866763271L);
  2287. TEST_ff_f (hypot, 0.7L, -12.4L, 12.419742348374220601176836866763271L);
  2288. TEST_ff_f (hypot, -0.7L, -12.4L, 12.419742348374220601176836866763271L);
  2289. TEST_ff_f (hypot, 12.4L, 0.7L, 12.419742348374220601176836866763271L);
  2290. TEST_ff_f (hypot, -12.4L, 0.7L, 12.419742348374220601176836866763271L);
  2291. TEST_ff_f (hypot, 12.4L, -0.7L, 12.419742348374220601176836866763271L);
  2292. TEST_ff_f (hypot, -12.4L, -0.7L, 12.419742348374220601176836866763271L);
  2293. /* hypot (x,0) == fabs (x) */
  2294. TEST_ff_f (hypot, 0.75L, 0, 0.75L);
  2295. TEST_ff_f (hypot, -0.75L, 0, 0.75L);
  2296. TEST_ff_f (hypot, -5.7e7, 0, 5.7e7L);
  2297. TEST_ff_f (hypot, 0.75L, 1.25L, 1.45773797371132511771853821938639577L);
  2298. END (hypot);
  2299. }
  2300. static void
  2301. ilogb_test (void)
  2302. {
  2303. START (ilogb);
  2304. TEST_f_i (ilogb, 1, 0);
  2305. TEST_f_i (ilogb, M_El, 1);
  2306. TEST_f_i (ilogb, 1024, 10);
  2307. TEST_f_i (ilogb, -2000, 10);
  2308. /* XXX We have a problem here: the standard does not tell us whether
  2309. exceptions are allowed/required. ignore them for now. */
  2310. TEST_f_i (ilogb, 0.0, FP_ILOGB0, EXCEPTIONS_OK);
  2311. TEST_f_i (ilogb, nan_value, FP_ILOGBNAN, EXCEPTIONS_OK);
  2312. TEST_f_i (ilogb, plus_infty, INT_MAX, EXCEPTIONS_OK);
  2313. TEST_f_i (ilogb, minus_infty, INT_MAX, EXCEPTIONS_OK);
  2314. END (ilogb);
  2315. }
  2316. static void
  2317. isfinite_test (void)
  2318. {
  2319. START (isfinite);
  2320. TEST_f_b (isfinite, 0, 1);
  2321. TEST_f_b (isfinite, minus_zero, 1);
  2322. TEST_f_b (isfinite, 10, 1);
  2323. TEST_f_b (isfinite, plus_infty, 0);
  2324. TEST_f_b (isfinite, minus_infty, 0);
  2325. TEST_f_b (isfinite, nan_value, 0);
  2326. END (isfinite);
  2327. }
  2328. static void
  2329. isnormal_test (void)
  2330. {
  2331. START (isnormal);
  2332. TEST_f_b (isnormal, 0, 0);
  2333. TEST_f_b (isnormal, minus_zero, 0);
  2334. TEST_f_b (isnormal, 10, 1);
  2335. TEST_f_b (isnormal, plus_infty, 0);
  2336. TEST_f_b (isnormal, minus_infty, 0);
  2337. TEST_f_b (isnormal, nan_value, 0);
  2338. END (isnormal);
  2339. }
  2340. static void
  2341. j0_test (void)
  2342. {
  2343. errno = 0;
  2344. #if 0
  2345. FLOAT s, c;
  2346. FUNC (sincos) (0, &s, &c);
  2347. if (errno == ENOSYS)
  2348. /* Required function not implemented. */
  2349. return;
  2350. #endif
  2351. FUNC(j0) (0);
  2352. if (errno == ENOSYS)
  2353. /* Function not implemented. */
  2354. return;
  2355. START (j0);
  2356. /* j0 is the Bessel function of the first kind of order 0 */
  2357. TEST_f_f (j0, nan_value, nan_value);
  2358. TEST_f_f (j0, plus_infty, 0);
  2359. TEST_f_f (j0, -1.0, 0.765197686557966551449717526102663221L);
  2360. TEST_f_f (j0, 0.0, 1.0);
  2361. TEST_f_f (j0, 0.125L, 0.996097563041985204620768999453174712L);
  2362. TEST_f_f (j0, 0.75L, 0.864242275166648623555731103820923211L);
  2363. TEST_f_f (j0, 1.0, 0.765197686557966551449717526102663221L);
  2364. TEST_f_f (j0, 1.5, 0.511827671735918128749051744283411720L);
  2365. TEST_f_f (j0, 2.0, 0.223890779141235668051827454649948626L);
  2366. TEST_f_f (j0, 8.0, 0.171650807137553906090869407851972001L);
  2367. TEST_f_f (j0, 10.0, -0.245935764451348335197760862485328754L);
  2368. TEST_f_f (j0, 4.0, -3.9714980986384737228659076845169804197562E-1L);
  2369. TEST_f_f (j0, -4.0, -3.9714980986384737228659076845169804197562E-1L);
  2370. END (j0);
  2371. }
  2372. static void
  2373. j1_test (void)
  2374. {
  2375. errno = 0;
  2376. #if 0
  2377. FLOAT s, c;
  2378. FUNC (sincos) (0, &s, &c);
  2379. if (errno == ENOSYS)
  2380. /* Required function not implemented. */
  2381. return;
  2382. #endif
  2383. FUNC(j1) (0);
  2384. if (errno == ENOSYS)
  2385. /* Function not implemented. */
  2386. return;
  2387. /* j1 is the Bessel function of the first kind of order 1 */
  2388. START (j1);
  2389. TEST_f_f (j1, nan_value, nan_value);
  2390. TEST_f_f (j1, plus_infty, 0);
  2391. TEST_f_f (j1, -1.0, -0.440050585744933515959682203718914913L);
  2392. TEST_f_f (j1, 0.0, 0.0);
  2393. TEST_f_f (j1, 0.125L, 0.0623780091344946810942311355879361177L);
  2394. TEST_f_f (j1, 0.75L, 0.349243602174862192523281016426251335L);
  2395. TEST_f_f (j1, 1.0, 0.440050585744933515959682203718914913L);
  2396. TEST_f_f (j1, 1.5, 0.557936507910099641990121213156089400L);
  2397. TEST_f_f (j1, 2.0, 0.576724807756873387202448242269137087L);
  2398. TEST_f_f (j1, 8.0, 0.234636346853914624381276651590454612L);
  2399. TEST_f_f (j1, 10.0, 0.0434727461688614366697487680258592883L);
  2400. END (j1);
  2401. }
  2402. static void
  2403. jn_test (void)
  2404. {
  2405. errno = 0;
  2406. #if 0
  2407. FLOAT s, c;
  2408. FUNC (sincos) (0, &s, &c);
  2409. if (errno == ENOSYS)
  2410. /* Required function not implemented. */
  2411. return;
  2412. #endif
  2413. FUNC(jn) (1, 1);
  2414. if (errno == ENOSYS)
  2415. /* Function not implemented. */
  2416. return;
  2417. /* jn is the Bessel function of the first kind of order n. */
  2418. START (jn);
  2419. /* jn (0, x) == j0 (x) */
  2420. TEST_ff_f (jn, 0, nan_value, nan_value);
  2421. TEST_ff_f (jn, 0, plus_infty, 0);
  2422. TEST_ff_f (jn, 0, -1.0, 0.765197686557966551449717526102663221L);
  2423. TEST_ff_f (jn, 0, 0.0, 1.0);
  2424. TEST_ff_f (jn, 0, 0.125L, 0.996097563041985204620768999453174712L);
  2425. TEST_ff_f (jn, 0, 0.75L, 0.864242275166648623555731103820923211L);
  2426. TEST_ff_f (jn, 0, 1.0, 0.765197686557966551449717526102663221L);
  2427. TEST_ff_f (jn, 0, 1.5, 0.511827671735918128749051744283411720L);
  2428. TEST_ff_f (jn, 0, 2.0, 0.223890779141235668051827454649948626L);
  2429. TEST_ff_f (jn, 0, 8.0, 0.171650807137553906090869407851972001L);
  2430. TEST_ff_f (jn, 0, 10.0, -0.245935764451348335197760862485328754L);
  2431. TEST_ff_f (jn, 0, 4.0, -3.9714980986384737228659076845169804197562E-1L);
  2432. TEST_ff_f (jn, 0, -4.0, -3.9714980986384737228659076845169804197562E-1L);
  2433. /* jn (1, x) == j1 (x) */
  2434. TEST_ff_f (jn, 1, nan_value, nan_value);
  2435. TEST_ff_f (jn, 1, plus_infty, 0);
  2436. TEST_ff_f (jn, 1, -1.0, -0.440050585744933515959682203718914913L);
  2437. TEST_ff_f (jn, 1, 0.0, 0.0);
  2438. TEST_ff_f (jn, 1, 0.125L, 0.0623780091344946810942311355879361177L);
  2439. TEST_ff_f (jn, 1, 0.75L, 0.349243602174862192523281016426251335L);
  2440. TEST_ff_f (jn, 1, 1.0, 0.440050585744933515959682203718914913L);
  2441. TEST_ff_f (jn, 1, 1.5, 0.557936507910099641990121213156089400L);
  2442. TEST_ff_f (jn, 1, 2.0, 0.576724807756873387202448242269137087L);
  2443. TEST_ff_f (jn, 1, 8.0, 0.234636346853914624381276651590454612L);
  2444. TEST_ff_f (jn, 1, 10.0, 0.0434727461688614366697487680258592883L);
  2445. /* jn (3, x) */
  2446. TEST_ff_f (jn, 3, nan_value, nan_value);
  2447. TEST_ff_f (jn, 3, plus_infty, 0);
  2448. TEST_ff_f (jn, 3, -1.0, -0.0195633539826684059189053216217515083L);
  2449. TEST_ff_f (jn, 3, 0.0, 0.0);
  2450. TEST_ff_f (jn, 3, 0.125L, 0.406503832554912875023029337653442868e-4L);
  2451. TEST_ff_f (jn, 3, 0.75L, 0.848438342327410884392755236884386804e-2L);
  2452. TEST_ff_f (jn, 3, 1.0, 0.0195633539826684059189053216217515083L);
  2453. TEST_ff_f (jn, 3, 2.0, 0.128943249474402051098793332969239835L);
  2454. TEST_ff_f (jn, 3, 10.0, 0.0583793793051868123429354784103409563L);
  2455. /* jn (10, x) */
  2456. TEST_ff_f (jn, 10, nan_value, nan_value);
  2457. TEST_ff_f (jn, 10, plus_infty, 0);
  2458. TEST_ff_f (jn, 10, -1.0, 0.263061512368745320699785368779050294e-9L);
  2459. TEST_ff_f (jn, 10, 0.0, 0.0);
  2460. TEST_ff_f (jn, 10, 0.125L, 0.250543369809369890173993791865771547e-18L);
  2461. TEST_ff_f (jn, 10, 0.75L, 0.149621713117596814698712483621682835e-10L);
  2462. TEST_ff_f (jn, 10, 1.0, 0.263061512368745320699785368779050294e-9L);
  2463. TEST_ff_f (jn, 10, 2.0, 0.251538628271673670963516093751820639e-6L);
  2464. TEST_ff_f (jn, 10, 10.0, 0.207486106633358857697278723518753428L);
  2465. END (jn);
  2466. }
  2467. static void
  2468. ldexp_test (void)
  2469. {
  2470. TEST_ff_f (ldexp, 0, 0, 0);
  2471. TEST_ff_f (ldexp, minus_zero, 0, minus_zero);
  2472. TEST_ff_f (ldexp, plus_infty, 1, plus_infty);
  2473. TEST_ff_f (ldexp, minus_infty, 1, minus_infty);
  2474. TEST_ff_f (ldexp, nan_value, 1, nan_value);
  2475. TEST_ff_f (ldexp, 0.8L, 4, 12.8L);
  2476. TEST_ff_f (ldexp, -0.854375L, 5, -27.34L);
  2477. /* ldexp (x, 0) == x. */
  2478. TEST_ff_f (ldexp, 1.0L, 0L, 1.0L);
  2479. }
  2480. static void
  2481. lgamma_test (void)
  2482. {
  2483. errno = 0;
  2484. FUNC(lgamma) (0);
  2485. if (errno == ENOSYS)
  2486. /* Function not implemented. */
  2487. return;
  2488. feclearexcept (FE_ALL_EXCEPT);
  2489. START (lgamma);
  2490. TEST_f_f (lgamma, plus_infty, plus_infty);
  2491. TEST_f_f (lgamma, 0, plus_infty, DIVIDE_BY_ZERO_EXCEPTION);
  2492. TEST_f_f (lgamma, nan_value, nan_value);
  2493. /* lgamma (x) == +inf plus divide by zero exception for integer x <= 0. */
  2494. TEST_f_f (lgamma, -3, plus_infty, DIVIDE_BY_ZERO_EXCEPTION);
  2495. TEST_f_f (lgamma, minus_infty, plus_infty);
  2496. TEST_f_f1 (lgamma, 1, 0, 1);
  2497. TEST_f_f1 (lgamma, 3, M_LN2l, 1);
  2498. TEST_f_f1 (lgamma, 0.5, M_LOG_SQRT_PIl, 1);
  2499. TEST_f_f1 (lgamma, -0.5, M_LOG_2_SQRT_PIl, -1);
  2500. TEST_f_f1 (lgamma, 0.7L, 0.260867246531666514385732417016759578L, 1);
  2501. TEST_f_f1 (lgamma, 1.2L, -0.853740900033158497197028392998854470e-1L, 1);
  2502. END (lgamma);
  2503. }
  2504. #if 0
  2505. static void
  2506. lrint_test (void)
  2507. {
  2508. /* XXX this test is incomplete. We need to have a way to specifiy
  2509. the rounding method and test the critical cases. So far, only
  2510. unproblematic numbers are tested. */
  2511. START (lrint);
  2512. TEST_f_l (lrint, 0.0, 0);
  2513. TEST_f_l (lrint, minus_zero, 0);
  2514. TEST_f_l (lrint, 0.2L, 0);
  2515. TEST_f_l (lrint, -0.2L, 0);
  2516. TEST_f_l (lrint, 1.4L, 1);
  2517. TEST_f_l (lrint, -1.4L, -1);
  2518. TEST_f_l (lrint, 8388600.3L, 8388600);
  2519. TEST_f_l (lrint, -8388600.3L, -8388600);
  2520. TEST_f_l (lrint, 1071930.0008, 1071930);
  2521. #ifndef TEST_FLOAT
  2522. TEST_f_l (lrint, 1073741824.01, 1073741824);
  2523. # if LONG_MAX > 281474976710656
  2524. TEST_f_l (lrint, 281474976710656.025, 281474976710656);
  2525. # endif
  2526. #endif
  2527. END (lrint);
  2528. }
  2529. static void
  2530. llrint_test (void)
  2531. {
  2532. /* XXX this test is incomplete. We need to have a way to specifiy
  2533. the rounding method and test the critical cases. So far, only
  2534. unproblematic numbers are tested. */
  2535. START (llrint);
  2536. TEST_f_L (llrint, 0.0, 0);
  2537. TEST_f_L (llrint, minus_zero, 0);
  2538. TEST_f_L (llrint, 0.2L, 0);
  2539. TEST_f_L (llrint, -0.2L, 0);
  2540. TEST_f_L (llrint, 1.4L, 1);
  2541. TEST_f_L (llrint, -1.4L, -1);
  2542. TEST_f_L (llrint, 8388600.3L, 8388600);
  2543. TEST_f_L (llrint, -8388600.3L, -8388600);
  2544. TEST_f_l (llrint, 1071930.0008, 1071930);
  2545. /* Test boundary conditions. */
  2546. /* 0x1FFFFF */
  2547. TEST_f_L (llrint, 2097151.0,2097151LL);
  2548. /* 0x800000 */
  2549. TEST_f_L (llrint, 8388608.0, 8388608LL);
  2550. /* 0x1000000 */
  2551. TEST_f_L (llrint, 16777216.0, 16777216LL);
  2552. /* 0x20000000000 */
  2553. TEST_f_L (llrint, 2199023255552.0, 2199023255552LL);
  2554. /* 0x40000000000 */
  2555. TEST_f_L (llrint, 4398046511104.0, 4398046511104LL);
  2556. /* 0x1000000000000 */
  2557. TEST_f_L (llrint, 281474976710656.0, 281474976710656LL);
  2558. /* 0x10000000000000 */
  2559. TEST_f_L (llrint, 4503599627370496.0, 4503599627370496LL);
  2560. /* 0x10000080000000 */
  2561. TEST_f_L (llrint, 4503601774854144.0, 4503601774854144LL);
  2562. /* 0x20000000000000 */
  2563. TEST_f_L (llrint, 9007199254740992.0, 9007199254740992LL);
  2564. /* 0x80000000000000 */
  2565. TEST_f_L (llrint, 36028797018963968.0, 36028797018963968LL);
  2566. /* 0x100000000000000 */
  2567. TEST_f_L (llrint, 72057594037927936.0, 72057594037927936LL);
  2568. #ifdef TEST_LDOUBLE
  2569. /* The input can only be represented in long double. */
  2570. TEST_f_L (llrint, 4503599627370495.5L, 4503599627370496LL);
  2571. TEST_f_L (llrint, 4503599627370496.25L, 4503599627370496LL);
  2572. TEST_f_L (llrint, 4503599627370496.5L, 4503599627370496LL);
  2573. TEST_f_L (llrint, 4503599627370496.75L, 4503599627370497LL);
  2574. TEST_f_L (llrint, 4503599627370497.5L, 4503599627370498LL);
  2575. TEST_f_L (llrint, -4503599627370495.5L, -4503599627370496LL);
  2576. TEST_f_L (llrint, -4503599627370496.25L, -4503599627370496LL);
  2577. TEST_f_L (llrint, -4503599627370496.5L, -4503599627370496LL);
  2578. TEST_f_L (llrint, -4503599627370496.75L, -4503599627370497LL);
  2579. TEST_f_L (llrint, -4503599627370497.5L, -4503599627370498LL);
  2580. TEST_f_L (llrint, 9007199254740991.5L, 9007199254740992LL);
  2581. TEST_f_L (llrint, 9007199254740992.25L, 9007199254740992LL);
  2582. TEST_f_L (llrint, 9007199254740992.5L, 9007199254740992LL);
  2583. TEST_f_L (llrint, 9007199254740992.75L, 9007199254740993LL);
  2584. TEST_f_L (llrint, 9007199254740993.5L, 9007199254740994LL);
  2585. TEST_f_L (llrint, -9007199254740991.5L, -9007199254740992LL);
  2586. TEST_f_L (llrint, -9007199254740992.25L, -9007199254740992LL);
  2587. TEST_f_L (llrint, -9007199254740992.5L, -9007199254740992LL);
  2588. TEST_f_L (llrint, -9007199254740992.75L, -9007199254740993LL);
  2589. TEST_f_L (llrint, -9007199254740993.5L, -9007199254740994LL);
  2590. TEST_f_L (llrint, 72057594037927935.5L, 72057594037927936LL);
  2591. TEST_f_L (llrint, 72057594037927936.25L, 72057594037927936LL);
  2592. TEST_f_L (llrint, 72057594037927936.5L, 72057594037927936LL);
  2593. TEST_f_L (llrint, 72057594037927936.75L, 72057594037927937LL);
  2594. TEST_f_L (llrint, 72057594037927937.5L, 72057594037927938LL);
  2595. TEST_f_L (llrint, -72057594037927935.5L, -72057594037927936LL);
  2596. TEST_f_L (llrint, -72057594037927936.25L, -72057594037927936LL);
  2597. TEST_f_L (llrint, -72057594037927936.5L, -72057594037927936LL);
  2598. TEST_f_L (llrint, -72057594037927936.75L, -72057594037927937LL);
  2599. TEST_f_L (llrint, -72057594037927937.5L, -72057594037927938LL);
  2600. #endif
  2601. END (llrint);
  2602. }
  2603. #endif
  2604. static void
  2605. log_test (void)
  2606. {
  2607. errno = 0;
  2608. FUNC(log) (1);
  2609. if (errno == ENOSYS)
  2610. /* Function not implemented. */
  2611. return;
  2612. START (log);
  2613. TEST_f_f (log, 0, minus_infty, DIVIDE_BY_ZERO_EXCEPTION);
  2614. TEST_f_f (log, minus_zero, minus_infty, DIVIDE_BY_ZERO_EXCEPTION);
  2615. TEST_f_f (log, 1, 0);
  2616. TEST_f_f (log, -1, nan_value, INVALID_EXCEPTION);
  2617. TEST_f_f (log, plus_infty, plus_infty);
  2618. TEST_f_f (log, M_El, 1);
  2619. TEST_f_f (log, 1.0 / M_El, -1);
  2620. TEST_f_f (log, 2, M_LN2l);
  2621. TEST_f_f (log, 10, M_LN10l);
  2622. TEST_f_f (log, 0.75L, -0.287682072451780927439219005993827432L);
  2623. END (log);
  2624. }
  2625. static void
  2626. log10_test (void)
  2627. {
  2628. errno = 0;
  2629. FUNC(log10) (1);
  2630. if (errno == ENOSYS)
  2631. /* Function not implemented. */
  2632. return;
  2633. START (log10);
  2634. TEST_f_f (log10, 0, minus_infty, DIVIDE_BY_ZERO_EXCEPTION);
  2635. TEST_f_f (log10, minus_zero, minus_infty, DIVIDE_BY_ZERO_EXCEPTION);
  2636. TEST_f_f (log10, 1, 0);
  2637. /* log10 (x) == NaN plus invalid exception if x < 0. */
  2638. TEST_f_f (log10, -1, nan_value, INVALID_EXCEPTION);
  2639. TEST_f_f (log10, plus_infty, plus_infty);
  2640. TEST_f_f (log10, nan_value, nan_value);
  2641. TEST_f_f (log10, 0.1L, -1);
  2642. TEST_f_f (log10, 10.0, 1);
  2643. TEST_f_f (log10, 100.0, 2);
  2644. TEST_f_f (log10, 10000.0, 4);
  2645. TEST_f_f (log10, M_El, M_LOG10El);
  2646. TEST_f_f (log10, 0.75L, -0.124938736608299953132449886193870744L);
  2647. END (log10);
  2648. }
  2649. static void
  2650. log1p_test (void)
  2651. {
  2652. errno = 0;
  2653. FUNC(log1p) (0);
  2654. if (errno == ENOSYS)
  2655. /* Function not implemented. */
  2656. return;
  2657. START (log1p);
  2658. TEST_f_f (log1p, 0, 0);
  2659. TEST_f_f (log1p, minus_zero, minus_zero);
  2660. TEST_f_f (log1p, -1, minus_infty, DIVIDE_BY_ZERO_EXCEPTION);
  2661. TEST_f_f (log1p, -2, nan_value, INVALID_EXCEPTION);
  2662. TEST_f_f (log1p, plus_infty, plus_infty);
  2663. TEST_f_f (log1p, nan_value, nan_value);
  2664. TEST_f_f (log1p, M_El - 1.0, 1);
  2665. TEST_f_f (log1p, -0.25L, -0.287682072451780927439219005993827432L);
  2666. TEST_f_f (log1p, -0.875, -2.07944154167983592825169636437452970L);
  2667. END (log1p);
  2668. }
  2669. #if 0
  2670. static void
  2671. log2_test (void)
  2672. {
  2673. errno = 0;
  2674. FUNC(log2) (1);
  2675. if (errno == ENOSYS)
  2676. /* Function not implemented. */
  2677. return;
  2678. START (log2);
  2679. TEST_f_f (log2, 0, minus_infty, DIVIDE_BY_ZERO_EXCEPTION);
  2680. TEST_f_f (log2, minus_zero, minus_infty, DIVIDE_BY_ZERO_EXCEPTION);
  2681. TEST_f_f (log2, 1, 0);
  2682. TEST_f_f (log2, -1, nan_value, INVALID_EXCEPTION);
  2683. TEST_f_f (log2, plus_infty, plus_infty);
  2684. TEST_f_f (log2, nan_value, nan_value);
  2685. TEST_f_f (log2, M_El, M_LOG2El);
  2686. TEST_f_f (log2, 2.0, 1);
  2687. TEST_f_f (log2, 16.0, 4);
  2688. TEST_f_f (log2, 256.0, 8);
  2689. TEST_f_f (log2, 0.75L, -.415037499278843818546261056052183492L);
  2690. END (log2);
  2691. }
  2692. #endif
  2693. static void
  2694. logb_test (void)
  2695. {
  2696. START (logb);
  2697. TEST_f_f (logb, plus_infty, plus_infty);
  2698. TEST_f_f (logb, minus_infty, plus_infty);
  2699. TEST_f_f (logb, 0, minus_infty, DIVIDE_BY_ZERO_EXCEPTION);
  2700. TEST_f_f (logb, minus_zero, minus_infty, DIVIDE_BY_ZERO_EXCEPTION);
  2701. TEST_f_f (logb, nan_value, nan_value);
  2702. TEST_f_f (logb, 1, 0);
  2703. TEST_f_f (logb, M_El, 1);
  2704. TEST_f_f (logb, 1024, 10);
  2705. TEST_f_f (logb, -2000, 10);
  2706. END (logb);
  2707. }
  2708. #if 0
  2709. static void
  2710. lround_test (void)
  2711. {
  2712. START (lround);
  2713. TEST_f_l (lround, 0, 0);
  2714. TEST_f_l (lround, minus_zero, 0);
  2715. TEST_f_l (lround, 0.2L, 0.0);
  2716. TEST_f_l (lround, -0.2L, 0);
  2717. TEST_f_l (lround, 0.5, 1);
  2718. TEST_f_l (lround, -0.5, -1);
  2719. TEST_f_l (lround, 0.8L, 1);
  2720. TEST_f_l (lround, -0.8L, -1);
  2721. TEST_f_l (lround, 1.5, 2);
  2722. TEST_f_l (lround, -1.5, -2);
  2723. TEST_f_l (lround, 22514.5, 22515);
  2724. TEST_f_l (lround, -22514.5, -22515);
  2725. TEST_f_l (lround, 1071930.0008, 1071930);
  2726. #ifndef TEST_FLOAT
  2727. TEST_f_l (lround, 1073741824.01, 1073741824);
  2728. # if LONG_MAX > 281474976710656
  2729. TEST_f_l (lround, 281474976710656.025, 281474976710656);
  2730. # endif
  2731. TEST_f_l (lround, 2097152.5, 2097153);
  2732. TEST_f_l (lround, -2097152.5, -2097153);
  2733. #endif
  2734. END (lround);
  2735. }
  2736. static void
  2737. llround_test (void)
  2738. {
  2739. START (llround);
  2740. TEST_f_L (llround, 0, 0);
  2741. TEST_f_L (llround, minus_zero, 0);
  2742. TEST_f_L (llround, 0.2L, 0.0);
  2743. TEST_f_L (llround, -0.2L, 0);
  2744. TEST_f_L (llround, 0.5, 1);
  2745. TEST_f_L (llround, -0.5, -1);
  2746. TEST_f_L (llround, 0.8L, 1);
  2747. TEST_f_L (llround, -0.8L, -1);
  2748. TEST_f_L (llround, 1.5, 2);
  2749. TEST_f_L (llround, -1.5, -2);
  2750. TEST_f_L (llround, 22514.5, 22515);
  2751. TEST_f_L (llround, -22514.5, -22515);
  2752. TEST_f_l (llround, 1071930.0008, 1071930);
  2753. #ifndef TEST_FLOAT
  2754. TEST_f_L (llround, 2097152.5, 2097153);
  2755. TEST_f_L (llround, -2097152.5, -2097153);
  2756. TEST_f_L (llround, 34359738368.5, 34359738369ll);
  2757. TEST_f_L (llround, -34359738368.5, -34359738369ll);
  2758. #endif
  2759. /* Test boundary conditions. */
  2760. /* 0x1FFFFF */
  2761. TEST_f_L (llround, 2097151.0, 2097151LL);
  2762. /* 0x800000 */
  2763. TEST_f_L (llround, 8388608.0, 8388608LL);
  2764. /* 0x1000000 */
  2765. TEST_f_L (llround, 16777216.0, 16777216LL);
  2766. /* 0x20000000000 */
  2767. TEST_f_L (llround, 2199023255552.0, 2199023255552LL);
  2768. /* 0x40000000000 */
  2769. TEST_f_L (llround, 4398046511104.0, 4398046511104LL);
  2770. /* 0x1000000000000 */
  2771. TEST_f_L (llround, 281474976710656.0, 281474976710656LL);
  2772. /* 0x10000000000000 */
  2773. TEST_f_L (llround, 4503599627370496.0, 4503599627370496LL);
  2774. /* 0x10000080000000 */
  2775. TEST_f_L (llround, 4503601774854144.0, 4503601774854144LL);
  2776. /* 0x20000000000000 */
  2777. TEST_f_L (llround, 9007199254740992.0, 9007199254740992LL);
  2778. /* 0x80000000000000 */
  2779. TEST_f_L (llround, 36028797018963968.0, 36028797018963968LL);
  2780. /* 0x100000000000000 */
  2781. TEST_f_L (llround, 72057594037927936.0, 72057594037927936LL);
  2782. #ifndef TEST_FLOAT
  2783. /* 0x100000000 */
  2784. TEST_f_L (llround, 4294967295.5, 4294967296LL);
  2785. /* 0x200000000 */
  2786. TEST_f_L (llround, 8589934591.5, 8589934592LL);
  2787. #endif
  2788. #ifdef TEST_LDOUBLE
  2789. /* The input can only be represented in long double. */
  2790. TEST_f_L (llround, 4503599627370495.5L, 4503599627370496LL);
  2791. TEST_f_L (llround, 4503599627370496.25L, 4503599627370496LL);
  2792. TEST_f_L (llround, 4503599627370496.5L, 4503599627370497LL);
  2793. TEST_f_L (llround, 4503599627370496.75L, 4503599627370497LL);
  2794. TEST_f_L (llround, 4503599627370497.5L, 4503599627370498LL);
  2795. TEST_f_L (llround, -4503599627370495.5L, -4503599627370496LL);
  2796. TEST_f_L (llround, -4503599627370496.25L, -4503599627370496LL);
  2797. TEST_f_L (llround, -4503599627370496.5L, -4503599627370497LL);
  2798. TEST_f_L (llround, -4503599627370496.75L, -4503599627370497LL);
  2799. TEST_f_L (llround, -4503599627370497.5L, -4503599627370498LL);
  2800. TEST_f_L (llround, 9007199254740991.5L, 9007199254740992LL);
  2801. TEST_f_L (llround, 9007199254740992.25L, 9007199254740992LL);
  2802. TEST_f_L (llround, 9007199254740992.5L, 9007199254740993LL);
  2803. TEST_f_L (llround, 9007199254740992.75L, 9007199254740993LL);
  2804. TEST_f_L (llround, 9007199254740993.5L, 9007199254740994LL);
  2805. TEST_f_L (llround, -9007199254740991.5L, -9007199254740992LL);
  2806. TEST_f_L (llround, -9007199254740992.25L, -9007199254740992LL);
  2807. TEST_f_L (llround, -9007199254740992.5L, -9007199254740993LL);
  2808. TEST_f_L (llround, -9007199254740992.75L, -9007199254740993LL);
  2809. TEST_f_L (llround, -9007199254740993.5L, -9007199254740994LL);
  2810. TEST_f_L (llround, 72057594037927935.5L, 72057594037927936LL);
  2811. TEST_f_L (llround, 72057594037927936.25L, 72057594037927936LL);
  2812. TEST_f_L (llround, 72057594037927936.5L, 72057594037927937LL);
  2813. TEST_f_L (llround, 72057594037927936.75L, 72057594037927937LL);
  2814. TEST_f_L (llround, 72057594037927937.5L, 72057594037927938LL);
  2815. TEST_f_L (llround, -72057594037927935.5L, -72057594037927936LL);
  2816. TEST_f_L (llround, -72057594037927936.25L, -72057594037927936LL);
  2817. TEST_f_L (llround, -72057594037927936.5L, -72057594037927937LL);
  2818. TEST_f_L (llround, -72057594037927936.75L, -72057594037927937LL);
  2819. TEST_f_L (llround, -72057594037927937.5L, -72057594037927938LL);
  2820. TEST_f_L (llround, 9223372036854775806.25L, 9223372036854775806LL);
  2821. TEST_f_L (llround, -9223372036854775806.25L, -9223372036854775806LL);
  2822. TEST_f_L (llround, 9223372036854775806.5L, 9223372036854775807LL);
  2823. TEST_f_L (llround, -9223372036854775806.5L, -9223372036854775807LL);
  2824. TEST_f_L (llround, 9223372036854775807.0L, 9223372036854775807LL);
  2825. TEST_f_L (llround, -9223372036854775807.0L, -9223372036854775807LL);
  2826. #endif
  2827. END (llround);
  2828. }
  2829. #endif
  2830. static void
  2831. modf_test (void)
  2832. {
  2833. FLOAT x;
  2834. START (modf);
  2835. TEST_fF_f1 (modf, plus_infty, 0, plus_infty);
  2836. TEST_fF_f1 (modf, minus_infty, minus_zero, minus_infty);
  2837. TEST_fF_f1 (modf, nan_value, nan_value, nan_value);
  2838. TEST_fF_f1 (modf, 0, 0, 0);
  2839. TEST_fF_f1 (modf, 1.5, 0.5, 1);
  2840. TEST_fF_f1 (modf, 2.5, 0.5, 2);
  2841. TEST_fF_f1 (modf, -2.5, -0.5, -2);
  2842. TEST_fF_f1 (modf, 20, 0, 20);
  2843. TEST_fF_f1 (modf, 21, 0, 21);
  2844. TEST_fF_f1 (modf, 89.5, 0.5, 89);
  2845. END (modf);
  2846. }
  2847. #if 0
  2848. static void
  2849. nearbyint_test (void)
  2850. {
  2851. START (nearbyint);
  2852. TEST_f_f (nearbyint, 0.0, 0.0);
  2853. TEST_f_f (nearbyint, minus_zero, minus_zero);
  2854. TEST_f_f (nearbyint, plus_infty, plus_infty);
  2855. TEST_f_f (nearbyint, minus_infty, minus_infty);
  2856. TEST_f_f (nearbyint, nan_value, nan_value);
  2857. /* Default rounding mode is round to nearest. */
  2858. TEST_f_f (nearbyint, 0.5, 0.0);
  2859. TEST_f_f (nearbyint, 1.5, 2.0);
  2860. TEST_f_f (nearbyint, -0.5, minus_zero);
  2861. TEST_f_f (nearbyint, -1.5, -2.0);
  2862. END (nearbyint);
  2863. }
  2864. static void
  2865. nextafter_test (void)
  2866. {
  2867. START (nextafter);
  2868. TEST_ff_f (nextafter, 0, 0, 0);
  2869. TEST_ff_f (nextafter, minus_zero, 0, 0);
  2870. TEST_ff_f (nextafter, 0, minus_zero, minus_zero);
  2871. TEST_ff_f (nextafter, minus_zero, minus_zero, minus_zero);
  2872. TEST_ff_f (nextafter, 9, 9, 9);
  2873. TEST_ff_f (nextafter, -9, -9, -9);
  2874. TEST_ff_f (nextafter, plus_infty, plus_infty, plus_infty);
  2875. TEST_ff_f (nextafter, minus_infty, minus_infty, minus_infty);
  2876. TEST_ff_f (nextafter, nan_value, 1.1L, nan_value);
  2877. TEST_ff_f (nextafter, 1.1L, nan_value, nan_value);
  2878. TEST_ff_f (nextafter, nan_value, nan_value, nan_value);
  2879. FLOAT fltmax = CHOOSE (LDBL_MAX, DBL_MAX, FLT_MAX,
  2880. LDBL_MAX, DBL_MAX, FLT_MAX);
  2881. TEST_ff_f (nextafter, fltmax, plus_infty, plus_infty);
  2882. TEST_ff_f (nextafter, -fltmax, minus_infty, minus_infty);
  2883. #ifdef TEST_LDOUBLE
  2884. // XXX Enable once gcc is fixed.
  2885. //TEST_ff_f (nextafter, 0x0.00000040000000000000p-16385L, -0.1L, 0x0.0000003ffffffff00000p-16385L);
  2886. #endif
  2887. /* XXX We need the hexadecimal FP number representation here for further
  2888. tests. */
  2889. END (nextafter);
  2890. }
  2891. static void
  2892. nexttoward_test (void)
  2893. {
  2894. START (nexttoward);
  2895. TEST_ff_f (nexttoward, 0, 0, 0);
  2896. TEST_ff_f (nexttoward, minus_zero, 0, 0);
  2897. TEST_ff_f (nexttoward, 0, minus_zero, minus_zero);
  2898. TEST_ff_f (nexttoward, minus_zero, minus_zero, minus_zero);
  2899. TEST_ff_f (nexttoward, 9, 9, 9);
  2900. TEST_ff_f (nexttoward, -9, -9, -9);
  2901. TEST_ff_f (nexttoward, plus_infty, plus_infty, plus_infty);
  2902. TEST_ff_f (nexttoward, minus_infty, minus_infty, minus_infty);
  2903. TEST_ff_f (nexttoward, nan_value, 1.1L, nan_value);
  2904. TEST_ff_f (nexttoward, 1.1L, nan_value, nan_value);
  2905. TEST_ff_f (nexttoward, nan_value, nan_value, nan_value);
  2906. /* XXX We need the hexadecimal FP number representation here for further
  2907. tests. */
  2908. END (nexttoward);
  2909. }
  2910. #endif
  2911. static void
  2912. pow_test (void)
  2913. {
  2914. errno = 0;
  2915. FUNC(pow) (0, 0);
  2916. if (errno == ENOSYS)
  2917. /* Function not implemented. */
  2918. return;
  2919. START (pow);
  2920. TEST_ff_f (pow, 0, 0, 1);
  2921. TEST_ff_f (pow, 0, minus_zero, 1);
  2922. TEST_ff_f (pow, minus_zero, 0, 1);
  2923. TEST_ff_f (pow, minus_zero, minus_zero, 1);
  2924. TEST_ff_f (pow, 10, 0, 1);
  2925. TEST_ff_f (pow, 10, minus_zero, 1);
  2926. TEST_ff_f (pow, -10, 0, 1);
  2927. TEST_ff_f (pow, -10, minus_zero, 1);
  2928. TEST_ff_f (pow, nan_value, 0, 1);
  2929. TEST_ff_f (pow, nan_value, minus_zero, 1);
  2930. #ifndef TEST_INLINE
  2931. TEST_ff_f (pow, 1.1L, plus_infty, plus_infty);
  2932. TEST_ff_f (pow, plus_infty, plus_infty, plus_infty);
  2933. TEST_ff_f (pow, -1.1L, plus_infty, plus_infty);
  2934. TEST_ff_f (pow, minus_infty, plus_infty, plus_infty);
  2935. TEST_ff_f (pow, 0.9L, plus_infty, 0);
  2936. TEST_ff_f (pow, 1e-7L, plus_infty, 0);
  2937. TEST_ff_f (pow, -0.9L, plus_infty, 0);
  2938. TEST_ff_f (pow, -1e-7L, plus_infty, 0);
  2939. TEST_ff_f (pow, 1.1L, minus_infty, 0);
  2940. TEST_ff_f (pow, plus_infty, minus_infty, 0);
  2941. TEST_ff_f (pow, -1.1L, minus_infty, 0);
  2942. TEST_ff_f (pow, minus_infty, minus_infty, 0);
  2943. TEST_ff_f (pow, 0.9L, minus_infty, plus_infty);
  2944. TEST_ff_f (pow, 1e-7L, minus_infty, plus_infty);
  2945. TEST_ff_f (pow, -0.9L, minus_infty, plus_infty);
  2946. TEST_ff_f (pow, -1e-7L, minus_infty, plus_infty);
  2947. TEST_ff_f (pow, plus_infty, 1e-7L, plus_infty);
  2948. TEST_ff_f (pow, plus_infty, 1, plus_infty);
  2949. TEST_ff_f (pow, plus_infty, 1e7L, plus_infty);
  2950. TEST_ff_f (pow, plus_infty, -1e-7L, 0);
  2951. TEST_ff_f (pow, plus_infty, -1, 0);
  2952. TEST_ff_f (pow, plus_infty, -1e7L, 0);
  2953. TEST_ff_f (pow, minus_infty, 1, minus_infty);
  2954. TEST_ff_f (pow, minus_infty, 11, minus_infty);
  2955. TEST_ff_f (pow, minus_infty, 1001, minus_infty);
  2956. TEST_ff_f (pow, minus_infty, 2, plus_infty);
  2957. TEST_ff_f (pow, minus_infty, 12, plus_infty);
  2958. TEST_ff_f (pow, minus_infty, 1002, plus_infty);
  2959. TEST_ff_f (pow, minus_infty, 0.1L, plus_infty);
  2960. TEST_ff_f (pow, minus_infty, 1.1L, plus_infty);
  2961. TEST_ff_f (pow, minus_infty, 11.1L, plus_infty);
  2962. TEST_ff_f (pow, minus_infty, 1001.1L, plus_infty);
  2963. TEST_ff_f (pow, minus_infty, -1, minus_zero);
  2964. TEST_ff_f (pow, minus_infty, -11, minus_zero);
  2965. TEST_ff_f (pow, minus_infty, -1001, minus_zero);
  2966. TEST_ff_f (pow, minus_infty, -2, 0);
  2967. TEST_ff_f (pow, minus_infty, -12, 0);
  2968. TEST_ff_f (pow, minus_infty, -1002, 0);
  2969. TEST_ff_f (pow, minus_infty, -0.1L, 0);
  2970. TEST_ff_f (pow, minus_infty, -1.1L, 0);
  2971. TEST_ff_f (pow, minus_infty, -11.1L, 0);
  2972. TEST_ff_f (pow, minus_infty, -1001.1L, 0);
  2973. #endif
  2974. TEST_ff_f (pow, nan_value, nan_value, nan_value);
  2975. TEST_ff_f (pow, 0, nan_value, nan_value);
  2976. TEST_ff_f (pow, 1, nan_value, 1);
  2977. TEST_ff_f (pow, -1, nan_value, nan_value);
  2978. TEST_ff_f (pow, nan_value, 1, nan_value);
  2979. TEST_ff_f (pow, nan_value, -1, nan_value);
  2980. /* pow (x, NaN) == NaN. */
  2981. TEST_ff_f (pow, 3.0, nan_value, nan_value);
  2982. TEST_ff_f (pow, 1, plus_infty, 1);
  2983. TEST_ff_f (pow, -1, plus_infty, 1);
  2984. TEST_ff_f (pow, 1, minus_infty, 1);
  2985. TEST_ff_f (pow, -1, minus_infty, 1);
  2986. TEST_ff_f (pow, 1, 1, 1);
  2987. TEST_ff_f (pow, 1, -1, 1);
  2988. TEST_ff_f (pow, 1, 1.25, 1);
  2989. TEST_ff_f (pow, 1, -1.25, 1);
  2990. TEST_ff_f (pow, 1, 0x1p62L, 1);
  2991. TEST_ff_f (pow, 1, 0x1p63L, 1);
  2992. TEST_ff_f (pow, 1, 0x1p64L, 1);
  2993. TEST_ff_f (pow, 1, 0x1p72L, 1);
  2994. /* pow (x, +-0) == 1. */
  2995. TEST_ff_f (pow, plus_infty, 0, 1);
  2996. TEST_ff_f (pow, plus_infty, minus_zero, 1);
  2997. TEST_ff_f (pow, minus_infty, 0, 1);
  2998. TEST_ff_f (pow, minus_infty, minus_zero, 1);
  2999. TEST_ff_f (pow, 32.75L, 0, 1);
  3000. TEST_ff_f (pow, 32.75L, minus_zero, 1);
  3001. TEST_ff_f (pow, -32.75L, 0, 1);
  3002. TEST_ff_f (pow, -32.75L, minus_zero, 1);
  3003. TEST_ff_f (pow, 0x1p72L, 0, 1);
  3004. TEST_ff_f (pow, 0x1p72L, minus_zero, 1);
  3005. TEST_ff_f (pow, 0x1p-72L, 0, 1);
  3006. TEST_ff_f (pow, 0x1p-72L, minus_zero, 1);
  3007. TEST_ff_f (pow, -0.1L, 1.1L, nan_value, INVALID_EXCEPTION);
  3008. TEST_ff_f (pow, -0.1L, -1.1L, nan_value, INVALID_EXCEPTION);
  3009. TEST_ff_f (pow, -10.1L, 1.1L, nan_value, INVALID_EXCEPTION);
  3010. TEST_ff_f (pow, -10.1L, -1.1L, nan_value, INVALID_EXCEPTION);
  3011. TEST_ff_f (pow, 0, -1, plus_infty, DIVIDE_BY_ZERO_EXCEPTION);
  3012. TEST_ff_f (pow, 0, -11, plus_infty, DIVIDE_BY_ZERO_EXCEPTION);
  3013. TEST_ff_f (pow, minus_zero, -1, minus_infty, DIVIDE_BY_ZERO_EXCEPTION);
  3014. TEST_ff_f (pow, minus_zero, -11, minus_infty, DIVIDE_BY_ZERO_EXCEPTION);
  3015. TEST_ff_f (pow, 0, -2, plus_infty, DIVIDE_BY_ZERO_EXCEPTION);
  3016. TEST_ff_f (pow, 0, -11.1L, plus_infty, DIVIDE_BY_ZERO_EXCEPTION);
  3017. TEST_ff_f (pow, minus_zero, -2, plus_infty, DIVIDE_BY_ZERO_EXCEPTION);
  3018. TEST_ff_f (pow, minus_zero, -11.1L, plus_infty, DIVIDE_BY_ZERO_EXCEPTION);
  3019. TEST_ff_f (pow, 0x1p72L, 0x1p72L, plus_infty);
  3020. TEST_ff_f (pow, 10, -0x1p72L, 0);
  3021. TEST_ff_f (pow, max_value, max_value, plus_infty);
  3022. TEST_ff_f (pow, 10, -max_value, 0);
  3023. TEST_ff_f (pow, 0, 1, 0);
  3024. TEST_ff_f (pow, 0, 11, 0);
  3025. TEST_ff_f (pow, minus_zero, 1, minus_zero);
  3026. TEST_ff_f (pow, minus_zero, 11, minus_zero);
  3027. TEST_ff_f (pow, 0, 2, 0);
  3028. TEST_ff_f (pow, 0, 11.1L, 0);
  3029. TEST_ff_f (pow, minus_zero, 2, 0);
  3030. TEST_ff_f (pow, minus_zero, 11.1L, 0);
  3031. TEST_ff_f (pow, 0, plus_infty, 0);
  3032. TEST_ff_f (pow, minus_zero, plus_infty, 0);
  3033. #ifndef TEST_INLINE
  3034. /* pow (x, +inf) == +inf for |x| > 1. */
  3035. TEST_ff_f (pow, 1.5, plus_infty, plus_infty);
  3036. /* pow (x, +inf) == +0 for |x| < 1. */
  3037. TEST_ff_f (pow, 0.5, plus_infty, 0.0);
  3038. /* pow (x, -inf) == +0 for |x| > 1. */
  3039. TEST_ff_f (pow, 1.5, minus_infty, 0.0);
  3040. /* pow (x, -inf) == +inf for |x| < 1. */
  3041. TEST_ff_f (pow, 0.5, minus_infty, plus_infty);
  3042. #endif
  3043. /* pow (+inf, y) == +inf for y > 0. */
  3044. TEST_ff_f (pow, plus_infty, 2, plus_infty);
  3045. /* pow (+inf, y) == +0 for y < 0. */
  3046. TEST_ff_f (pow, plus_infty, -1, 0.0);
  3047. /* pow (-inf, y) == -inf for y an odd integer > 0. */
  3048. TEST_ff_f (pow, minus_infty, 27, minus_infty);
  3049. /* pow (-inf, y) == +inf for y > 0 and not an odd integer. */
  3050. TEST_ff_f (pow, minus_infty, 28, plus_infty);
  3051. /* pow (-inf, y) == -0 for y an odd integer < 0. */
  3052. TEST_ff_f (pow, minus_infty, -3, minus_zero);
  3053. /* pow (-inf, y) == +0 for y < 0 and not an odd integer. */
  3054. TEST_ff_f (pow, minus_infty, -2.0, 0.0);
  3055. /* pow (+0, y) == +0 for y an odd integer > 0. */
  3056. TEST_ff_f (pow, 0.0, 27, 0.0);
  3057. /* pow (-0, y) == -0 for y an odd integer > 0. */
  3058. TEST_ff_f (pow, minus_zero, 27, minus_zero);
  3059. /* pow (+0, y) == +0 for y > 0 and not an odd integer. */
  3060. TEST_ff_f (pow, 0.0, 4, 0.0);
  3061. /* pow (-0, y) == +0 for y > 0 and not an odd integer. */
  3062. TEST_ff_f (pow, minus_zero, 4, 0.0);
  3063. TEST_ff_f (pow, 16, 0.25L, 2);
  3064. TEST_ff_f (pow, 0x1p64L, 0.125L, 256);
  3065. TEST_ff_f (pow, 2, 4, 16);
  3066. TEST_ff_f (pow, 256, 8, 0x1p64L);
  3067. TEST_ff_f (pow, 0.75L, 1.25L, 0.697953644326574699205914060237425566L);
  3068. #if defined TEST_DOUBLE || defined TEST_LDOUBLE
  3069. TEST_ff_f (pow, -7.49321e+133, -9.80818e+16, 0);
  3070. #endif
  3071. END (pow);
  3072. }
  3073. static void
  3074. remainder_test (void)
  3075. {
  3076. errno = 0;
  3077. FUNC(remainder) (1.625, 1.0);
  3078. if (errno == ENOSYS)
  3079. /* Function not implemented. */
  3080. return;
  3081. START (remainder);
  3082. TEST_ff_f (remainder, 1, 0, nan_value, INVALID_EXCEPTION);
  3083. TEST_ff_f (remainder, 1, minus_zero, nan_value, INVALID_EXCEPTION);
  3084. TEST_ff_f (remainder, plus_infty, 1, nan_value, INVALID_EXCEPTION);
  3085. TEST_ff_f (remainder, minus_infty, 1, nan_value, INVALID_EXCEPTION);
  3086. TEST_ff_f (remainder, nan_value, nan_value, nan_value);
  3087. TEST_ff_f (remainder, 1.625, 1.0, -0.375);
  3088. TEST_ff_f (remainder, -1.625, 1.0, 0.375);
  3089. TEST_ff_f (remainder, 1.625, -1.0, -0.375);
  3090. TEST_ff_f (remainder, -1.625, -1.0, 0.375);
  3091. TEST_ff_f (remainder, 5.0, 2.0, 1.0);
  3092. TEST_ff_f (remainder, 3.0, 2.0, -1.0);
  3093. END (remainder);
  3094. }
  3095. #if 0
  3096. static void
  3097. remquo_test (void)
  3098. {
  3099. /* x is needed. */
  3100. int x;
  3101. errno = 0;
  3102. FUNC(remquo) (1.625, 1.0, &x);
  3103. if (errno == ENOSYS)
  3104. /* Function not implemented. */
  3105. return;
  3106. START (remquo);
  3107. TEST_ffI_f1 (remquo, 1, 0, nan_value, IGNORE, INVALID_EXCEPTION);
  3108. TEST_ffI_f1 (remquo, 1, minus_zero, nan_value, IGNORE, INVALID_EXCEPTION);
  3109. TEST_ffI_f1 (remquo, plus_infty, 1, nan_value, IGNORE, INVALID_EXCEPTION);
  3110. TEST_ffI_f1 (remquo, minus_infty, 1, nan_value, IGNORE, INVALID_EXCEPTION);
  3111. TEST_ffI_f1 (remquo, nan_value, nan_value, nan_value, IGNORE);
  3112. TEST_ffI_f1 (remquo, 1.625, 1.0, -0.375, 2);
  3113. TEST_ffI_f1 (remquo, -1.625, 1.0, 0.375, -2);
  3114. TEST_ffI_f1 (remquo, 1.625, -1.0, -0.375, -2);
  3115. TEST_ffI_f1 (remquo, -1.625, -1.0, 0.375, 2);
  3116. TEST_ffI_f1 (remquo, 5, 2, 1, 2);
  3117. TEST_ffI_f1 (remquo, 3, 2, -1, 2);
  3118. END (remquo);
  3119. }
  3120. #endif
  3121. static void
  3122. rint_test (void)
  3123. {
  3124. START (rint);
  3125. TEST_f_f (rint, 0.0, 0.0);
  3126. TEST_f_f (rint, minus_zero, minus_zero);
  3127. TEST_f_f (rint, plus_infty, plus_infty);
  3128. TEST_f_f (rint, minus_infty, minus_infty);
  3129. /* Default rounding mode is round to even. */
  3130. TEST_f_f (rint, 0.5, 0.0);
  3131. TEST_f_f (rint, 1.5, 2.0);
  3132. TEST_f_f (rint, 2.5, 2.0);
  3133. TEST_f_f (rint, 3.5, 4.0);
  3134. TEST_f_f (rint, 4.5, 4.0);
  3135. TEST_f_f (rint, -0.5, -0.0);
  3136. TEST_f_f (rint, -1.5, -2.0);
  3137. TEST_f_f (rint, -2.5, -2.0);
  3138. TEST_f_f (rint, -3.5, -4.0);
  3139. TEST_f_f (rint, -4.5, -4.0);
  3140. #ifdef TEST_LDOUBLE
  3141. /* The result can only be represented in long double. */
  3142. TEST_f_f (rint, 4503599627370495.5L, 4503599627370496.0L);
  3143. TEST_f_f (rint, 4503599627370496.25L, 4503599627370496.0L);
  3144. TEST_f_f (rint, 4503599627370496.5L, 4503599627370496.0L);
  3145. TEST_f_f (rint, 4503599627370496.75L, 4503599627370497.0L);
  3146. TEST_f_f (rint, 4503599627370497.5L, 4503599627370498.0L);
  3147. TEST_f_f (rint, -4503599627370495.5L, -4503599627370496.0L);
  3148. TEST_f_f (rint, -4503599627370496.25L, -4503599627370496.0L);
  3149. TEST_f_f (rint, -4503599627370496.5L, -4503599627370496.0L);
  3150. TEST_f_f (rint, -4503599627370496.75L, -4503599627370497.0L);
  3151. TEST_f_f (rint, -4503599627370497.5L, -4503599627370498.0L);
  3152. TEST_f_f (rint, 9007199254740991.5L, 9007199254740992.0L);
  3153. TEST_f_f (rint, 9007199254740992.25L, 9007199254740992.0L);
  3154. TEST_f_f (rint, 9007199254740992.5L, 9007199254740992.0L);
  3155. TEST_f_f (rint, 9007199254740992.75L, 9007199254740993.0L);
  3156. TEST_f_f (rint, 9007199254740993.5L, 9007199254740994.0L);
  3157. TEST_f_f (rint, -9007199254740991.5L, -9007199254740992.0L);
  3158. TEST_f_f (rint, -9007199254740992.25L, -9007199254740992.0L);
  3159. TEST_f_f (rint, -9007199254740992.5L, -9007199254740992.0L);
  3160. TEST_f_f (rint, -9007199254740992.75L, -9007199254740993.0L);
  3161. TEST_f_f (rint, -9007199254740993.5L, -9007199254740994.0L);
  3162. TEST_f_f (rint, 72057594037927935.5L, 72057594037927936.0L);
  3163. TEST_f_f (rint, 72057594037927936.25L, 72057594037927936.0L);
  3164. TEST_f_f (rint, 72057594037927936.5L, 72057594037927936.0L);
  3165. TEST_f_f (rint, 72057594037927936.75L, 72057594037927937.0L);
  3166. TEST_f_f (rint, 72057594037927937.5L, 72057594037927938.0L);
  3167. TEST_f_f (rint, -72057594037927935.5L, -72057594037927936.0L);
  3168. TEST_f_f (rint, -72057594037927936.25L, -72057594037927936.0L);
  3169. TEST_f_f (rint, -72057594037927936.5L, -72057594037927936.0L);
  3170. TEST_f_f (rint, -72057594037927936.75L, -72057594037927937.0L);
  3171. TEST_f_f (rint, -72057594037927937.5L, -72057594037927938.0L);
  3172. TEST_f_f (rint, 10141204801825835211973625643007.5L, 10141204801825835211973625643008.0L);
  3173. TEST_f_f (rint, 10141204801825835211973625643008.25L, 10141204801825835211973625643008.0L);
  3174. TEST_f_f (rint, 10141204801825835211973625643008.5L, 10141204801825835211973625643008.0L);
  3175. TEST_f_f (rint, 10141204801825835211973625643008.75L, 10141204801825835211973625643009.0L);
  3176. TEST_f_f (rint, 10141204801825835211973625643009.5L, 10141204801825835211973625643010.0L);
  3177. #endif
  3178. END (rint);
  3179. }
  3180. #if 0
  3181. static void
  3182. rint_test_tonearest (void)
  3183. {
  3184. int save_round_mode;
  3185. START (rint_tonearest);
  3186. save_round_mode = fegetround();
  3187. if (!fesetround (FE_TONEAREST))
  3188. {
  3189. TEST_f_f (rint, 2.0, 2.0);
  3190. TEST_f_f (rint, 1.5, 2.0);
  3191. TEST_f_f (rint, 1.0, 1.0);
  3192. TEST_f_f (rint, 0.5, 0.0);
  3193. TEST_f_f (rint, 0.0, 0.0);
  3194. TEST_f_f (rint, minus_zero, minus_zero);
  3195. TEST_f_f (rint, -0.5, -0.0);
  3196. TEST_f_f (rint, -1.0, -1.0);
  3197. TEST_f_f (rint, -1.5, -2.0);
  3198. TEST_f_f (rint, -2.0, -2.0);
  3199. }
  3200. fesetround(save_round_mode);
  3201. END (rint_tonearest);
  3202. }
  3203. static void
  3204. rint_test_towardzero (void)
  3205. {
  3206. int save_round_mode;
  3207. START (rint_towardzero);
  3208. save_round_mode = fegetround();
  3209. if (!fesetround (FE_TOWARDZERO))
  3210. {
  3211. TEST_f_f (rint, 2.0, 2.0);
  3212. TEST_f_f (rint, 1.5, 1.0);
  3213. TEST_f_f (rint, 1.0, 1.0);
  3214. TEST_f_f (rint, 0.5, 0.0);
  3215. TEST_f_f (rint, 0.0, 0.0);
  3216. TEST_f_f (rint, minus_zero, minus_zero);
  3217. TEST_f_f (rint, -0.5, -0.0);
  3218. TEST_f_f (rint, -1.0, -1.0);
  3219. TEST_f_f (rint, -1.5, -1.0);
  3220. TEST_f_f (rint, -2.0, -2.0);
  3221. }
  3222. fesetround(save_round_mode);
  3223. END (rint_towardzero);
  3224. }
  3225. static void
  3226. rint_test_downward (void)
  3227. {
  3228. int save_round_mode;
  3229. START (rint_downward);
  3230. save_round_mode = fegetround();
  3231. if (!fesetround (FE_DOWNWARD))
  3232. {
  3233. TEST_f_f (rint, 2.0, 2.0);
  3234. TEST_f_f (rint, 1.5, 1.0);
  3235. TEST_f_f (rint, 1.0, 1.0);
  3236. TEST_f_f (rint, 0.5, 0.0);
  3237. TEST_f_f (rint, 0.0, 0.0);
  3238. TEST_f_f (rint, minus_zero, minus_zero);
  3239. TEST_f_f (rint, -0.5, -1.0);
  3240. TEST_f_f (rint, -1.0, -1.0);
  3241. TEST_f_f (rint, -1.5, -2.0);
  3242. TEST_f_f (rint, -2.0, -2.0);
  3243. }
  3244. fesetround(save_round_mode);
  3245. END (rint_downward);
  3246. }
  3247. static void
  3248. rint_test_upward (void)
  3249. {
  3250. int save_round_mode;
  3251. START (rint_upward);
  3252. save_round_mode = fegetround();
  3253. if (!fesetround (FE_UPWARD))
  3254. {
  3255. TEST_f_f (rint, 2.0, 2.0);
  3256. TEST_f_f (rint, 1.5, 2.0);
  3257. TEST_f_f (rint, 1.0, 1.0);
  3258. TEST_f_f (rint, 0.5, 1.0);
  3259. TEST_f_f (rint, 0.0, 0.0);
  3260. TEST_f_f (rint, minus_zero, minus_zero);
  3261. TEST_f_f (rint, -0.5, -0.0);
  3262. TEST_f_f (rint, -1.0, -1.0);
  3263. TEST_f_f (rint, -1.5, -1.0);
  3264. TEST_f_f (rint, -2.0, -2.0);
  3265. }
  3266. fesetround(save_round_mode);
  3267. END (rint_upward);
  3268. }
  3269. static void
  3270. round_test (void)
  3271. {
  3272. START (round);
  3273. TEST_f_f (round, 0, 0);
  3274. TEST_f_f (round, minus_zero, minus_zero);
  3275. TEST_f_f (round, 0.2L, 0.0);
  3276. TEST_f_f (round, -0.2L, minus_zero);
  3277. TEST_f_f (round, 0.5, 1.0);
  3278. TEST_f_f (round, -0.5, -1.0);
  3279. TEST_f_f (round, 0.8L, 1.0);
  3280. TEST_f_f (round, -0.8L, -1.0);
  3281. TEST_f_f (round, 1.5, 2.0);
  3282. TEST_f_f (round, -1.5, -2.0);
  3283. TEST_f_f (round, 2097152.5, 2097153);
  3284. TEST_f_f (round, -2097152.5, -2097153);
  3285. #ifdef TEST_LDOUBLE
  3286. /* The result can only be represented in long double. */
  3287. TEST_f_f (round, 4503599627370495.5L, 4503599627370496.0L);
  3288. TEST_f_f (round, 4503599627370496.25L, 4503599627370496.0L);
  3289. TEST_f_f (round, 4503599627370496.5L, 4503599627370497.0L);
  3290. TEST_f_f (round, 4503599627370496.75L, 4503599627370497.0L);
  3291. TEST_f_f (round, 4503599627370497.5L, 4503599627370498.0L);
  3292. TEST_f_f (round, -4503599627370495.5L, -4503599627370496.0L);
  3293. TEST_f_f (round, -4503599627370496.25L, -4503599627370496.0L);
  3294. TEST_f_f (round, -4503599627370496.5L, -4503599627370497.0L);
  3295. TEST_f_f (round, -4503599627370496.75L, -4503599627370497.0L);
  3296. TEST_f_f (round, -4503599627370497.5L, -4503599627370498.0L);
  3297. TEST_f_f (round, 9007199254740991.5L, 9007199254740992.0L);
  3298. TEST_f_f (round, 9007199254740992.25L, 9007199254740992.0L);
  3299. TEST_f_f (round, 9007199254740992.5L, 9007199254740993.0L);
  3300. TEST_f_f (round, 9007199254740992.75L, 9007199254740993.0L);
  3301. TEST_f_f (round, 9007199254740993.5L, 9007199254740994.0L);
  3302. TEST_f_f (round, -9007199254740991.5L, -9007199254740992.0L);
  3303. TEST_f_f (round, -9007199254740992.25L, -9007199254740992.0L);
  3304. TEST_f_f (round, -9007199254740992.5L, -9007199254740993.0L);
  3305. TEST_f_f (round, -9007199254740992.75L, -9007199254740993.0L);
  3306. TEST_f_f (round, -9007199254740993.5L, -9007199254740994.0L);
  3307. TEST_f_f (round, 72057594037927935.5L, 72057594037927936.0L);
  3308. TEST_f_f (round, 72057594037927936.25L, 72057594037927936.0L);
  3309. TEST_f_f (round, 72057594037927936.5L, 72057594037927937.0L);
  3310. TEST_f_f (round, 72057594037927936.75L, 72057594037927937.0L);
  3311. TEST_f_f (round, 72057594037927937.5L, 72057594037927938.0L);
  3312. TEST_f_f (round, -72057594037927935.5L, -72057594037927936.0L);
  3313. TEST_f_f (round, -72057594037927936.25L, -72057594037927936.0L);
  3314. TEST_f_f (round, -72057594037927936.5L, -72057594037927937.0L);
  3315. TEST_f_f (round, -72057594037927936.75L, -72057594037927937.0L);
  3316. TEST_f_f (round, -72057594037927937.5L, -72057594037927938.0L);
  3317. TEST_f_f (round, 10141204801825835211973625643007.5L, 10141204801825835211973625643008.0L);
  3318. TEST_f_f (round, 10141204801825835211973625643008.25L, 10141204801825835211973625643008.0L);
  3319. TEST_f_f (round, 10141204801825835211973625643008.5L, 10141204801825835211973625643009.0L);
  3320. TEST_f_f (round, 10141204801825835211973625643008.75L, 10141204801825835211973625643009.0L);
  3321. TEST_f_f (round, 10141204801825835211973625643009.5L, 10141204801825835211973625643010.0L);
  3322. #endif
  3323. END (round);
  3324. }
  3325. #endif
  3326. static void
  3327. scalb_test (void)
  3328. {
  3329. START (scalb);
  3330. #ifndef TEST_LDOUBLE /* uclibc doesn't have scalbl */
  3331. #ifdef __UCLIBC_SUSV3_LEGACY__ /* scalbf is susv3 legacy */
  3332. TEST_ff_f (scalb, 2.0, 0.5, nan_value, INVALID_EXCEPTION);
  3333. TEST_ff_f (scalb, 3.0, -2.5, nan_value, INVALID_EXCEPTION);
  3334. TEST_ff_f (scalb, 0, nan_value, nan_value);
  3335. TEST_ff_f (scalb, 1, nan_value, nan_value);
  3336. TEST_ff_f (scalb, 1, 0, 1);
  3337. TEST_ff_f (scalb, -1, 0, -1);
  3338. TEST_ff_f (scalb, 0, plus_infty, nan_value, INVALID_EXCEPTION);
  3339. TEST_ff_f (scalb, minus_zero, plus_infty, nan_value, INVALID_EXCEPTION);
  3340. TEST_ff_f (scalb, 0, 2, 0);
  3341. TEST_ff_f (scalb, minus_zero, -4, minus_zero);
  3342. TEST_ff_f (scalb, 0, 0, 0);
  3343. TEST_ff_f (scalb, minus_zero, 0, minus_zero);
  3344. TEST_ff_f (scalb, 0, -1, 0);
  3345. TEST_ff_f (scalb, minus_zero, -10, minus_zero);
  3346. TEST_ff_f (scalb, 0, minus_infty, 0);
  3347. TEST_ff_f (scalb, minus_zero, minus_infty, minus_zero);
  3348. TEST_ff_f (scalb, plus_infty, -1, plus_infty);
  3349. TEST_ff_f (scalb, minus_infty, -10, minus_infty);
  3350. TEST_ff_f (scalb, plus_infty, 0, plus_infty);
  3351. TEST_ff_f (scalb, minus_infty, 0, minus_infty);
  3352. TEST_ff_f (scalb, plus_infty, 2, plus_infty);
  3353. TEST_ff_f (scalb, minus_infty, 100, minus_infty);
  3354. TEST_ff_f (scalb, 0.1L, minus_infty, 0.0);
  3355. TEST_ff_f (scalb, -0.1L, minus_infty, minus_zero);
  3356. TEST_ff_f (scalb, 1, plus_infty, plus_infty);
  3357. TEST_ff_f (scalb, -1, plus_infty, minus_infty);
  3358. TEST_ff_f (scalb, plus_infty, plus_infty, plus_infty);
  3359. TEST_ff_f (scalb, minus_infty, plus_infty, minus_infty);
  3360. TEST_ff_f (scalb, plus_infty, minus_infty, nan_value, INVALID_EXCEPTION);
  3361. TEST_ff_f (scalb, minus_infty, minus_infty, nan_value, INVALID_EXCEPTION);
  3362. TEST_ff_f (scalb, nan_value, 1, nan_value);
  3363. TEST_ff_f (scalb, 1, nan_value, nan_value);
  3364. TEST_ff_f (scalb, nan_value, 0, nan_value);
  3365. TEST_ff_f (scalb, 0, nan_value, nan_value);
  3366. TEST_ff_f (scalb, nan_value, plus_infty, nan_value);
  3367. TEST_ff_f (scalb, plus_infty, nan_value, nan_value);
  3368. TEST_ff_f (scalb, nan_value, nan_value, nan_value);
  3369. TEST_ff_f (scalb, 0.8L, 4, 12.8L);
  3370. TEST_ff_f (scalb, -0.854375L, 5, -27.34L);
  3371. #endif /* __UCLIBC_SUSV3_LEGACY__ */
  3372. #endif /* TEST_LDOUBLE */
  3373. END (scalb);
  3374. }
  3375. static void
  3376. scalbn_test (void)
  3377. {
  3378. START (scalbn);
  3379. TEST_fi_f (scalbn, 0, 0, 0);
  3380. TEST_fi_f (scalbn, minus_zero, 0, minus_zero);
  3381. TEST_fi_f (scalbn, plus_infty, 1, plus_infty);
  3382. TEST_fi_f (scalbn, minus_infty, 1, minus_infty);
  3383. TEST_fi_f (scalbn, nan_value, 1, nan_value);
  3384. TEST_fi_f (scalbn, 0.8L, 4, 12.8L);
  3385. TEST_fi_f (scalbn, -0.854375L, 5, -27.34L);
  3386. TEST_fi_f (scalbn, 1, 0L, 1);
  3387. END (scalbn);
  3388. }
  3389. #if 0
  3390. static void
  3391. scalbln_test (void)
  3392. {
  3393. START (scalbln);
  3394. TEST_fl_f (scalbln, 0, 0, 0);
  3395. TEST_fl_f (scalbln, minus_zero, 0, minus_zero);
  3396. TEST_fl_f (scalbln, plus_infty, 1, plus_infty);
  3397. TEST_fl_f (scalbln, minus_infty, 1, minus_infty);
  3398. TEST_fl_f (scalbln, nan_value, 1, nan_value);
  3399. TEST_fl_f (scalbln, 0.8L, 4, 12.8L);
  3400. TEST_fl_f (scalbln, -0.854375L, 5, -27.34L);
  3401. TEST_fl_f (scalbln, 1, 0L, 1);
  3402. END (scalbn);
  3403. }
  3404. #endif
  3405. static void
  3406. signbit_test (void)
  3407. {
  3408. START (signbit);
  3409. TEST_f_b (signbit, 0, 0);
  3410. TEST_f_b (signbit, minus_zero, 1);
  3411. TEST_f_b (signbit, plus_infty, 0);
  3412. TEST_f_b (signbit, minus_infty, 1);
  3413. /* signbit (x) != 0 for x < 0. */
  3414. TEST_f_b (signbit, -1, 1);
  3415. /* signbit (x) == 0 for x >= 0. */
  3416. TEST_f_b (signbit, 1, 0);
  3417. END (signbit);
  3418. }
  3419. static void
  3420. sin_test (void)
  3421. {
  3422. errno = 0;
  3423. FUNC(sin) (0);
  3424. if (errno == ENOSYS)
  3425. /* Function not implemented. */
  3426. return;
  3427. START (sin);
  3428. TEST_f_f (sin, 0, 0);
  3429. TEST_f_f (sin, minus_zero, minus_zero);
  3430. TEST_f_f (sin, plus_infty, nan_value, INVALID_EXCEPTION);
  3431. TEST_f_f (sin, minus_infty, nan_value, INVALID_EXCEPTION);
  3432. TEST_f_f (sin, nan_value, nan_value);
  3433. TEST_f_f (sin, M_PI_6l, 0.5);
  3434. TEST_f_f (sin, -M_PI_6l, -0.5);
  3435. TEST_f_f (sin, M_PI_2l, 1);
  3436. TEST_f_f (sin, -M_PI_2l, -1);
  3437. TEST_f_f (sin, 0.75L, 0.681638760023334166733241952779893935L);
  3438. #ifdef TEST_DOUBLE
  3439. TEST_f_f (sin, 0.80190127184058835, 0.71867942238767868);
  3440. #endif
  3441. END (sin);
  3442. }
  3443. #if 0
  3444. static void
  3445. sincos_test (void)
  3446. {
  3447. FLOAT sin_res, cos_res;
  3448. errno = 0;
  3449. FUNC(sincos) (0, &sin_res, &cos_res);
  3450. if (errno == ENOSYS)
  3451. /* Function not implemented. */
  3452. return;
  3453. START (sincos);
  3454. /* sincos is treated differently because it returns void. */
  3455. TEST_extra (sincos, 0, 0, 1);
  3456. TEST_extra (sincos, minus_zero, minus_zero, 1);
  3457. TEST_extra (sincos, plus_infty, nan_value, nan_value, INVALID_EXCEPTION);
  3458. TEST_extra (sincos, minus_infty, nan_value, nan_value, INVALID_EXCEPTION);
  3459. TEST_extra (sincos, nan_value, nan_value, nan_value);
  3460. TEST_extra (sincos, M_PI_2l, 1, 0);
  3461. TEST_extra (sincos, M_PI_6l, 0.5, 0.86602540378443864676372317075293616L);
  3462. TEST_extra (sincos, M_PI_6l*2.0, 0.86602540378443864676372317075293616L, 0.5);
  3463. TEST_extra (sincos, 0.75L, 0.681638760023334166733241952779893935L, 0.731688868873820886311838753000084544L);
  3464. #ifdef TEST_DOUBLE
  3465. TEST_extra (sincos, 0.80190127184058835, 0.71867942238767868, 0.69534156199418473);
  3466. #endif
  3467. END (sincos);
  3468. }
  3469. #endif
  3470. static void
  3471. sinh_test (void)
  3472. {
  3473. errno = 0;
  3474. FUNC(sinh) (0.7L);
  3475. if (errno == ENOSYS)
  3476. /* Function not implemented. */
  3477. return;
  3478. START (sinh);
  3479. TEST_f_f (sinh, 0, 0);
  3480. TEST_f_f (sinh, minus_zero, minus_zero);
  3481. #ifndef TEST_INLINE
  3482. TEST_f_f (sinh, plus_infty, plus_infty);
  3483. TEST_f_f (sinh, minus_infty, minus_infty);
  3484. #endif
  3485. TEST_f_f (sinh, nan_value, nan_value);
  3486. TEST_f_f (sinh, 0.75L, 0.822316731935829980703661634446913849L);
  3487. TEST_f_f (sinh, 0x8p-32L, 1.86264514923095703232705808926175479e-9L);
  3488. END (sinh);
  3489. }
  3490. static void
  3491. sqrt_test (void)
  3492. {
  3493. errno = 0;
  3494. FUNC(sqrt) (1);
  3495. if (errno == ENOSYS)
  3496. /* Function not implemented. */
  3497. return;
  3498. START (sqrt);
  3499. TEST_f_f (sqrt, 0, 0);
  3500. TEST_f_f (sqrt, nan_value, nan_value);
  3501. TEST_f_f (sqrt, plus_infty, plus_infty);
  3502. TEST_f_f (sqrt, minus_zero, minus_zero);
  3503. /* sqrt (x) == NaN plus invalid exception for x < 0. */
  3504. TEST_f_f (sqrt, -1, nan_value, INVALID_EXCEPTION);
  3505. TEST_f_f (sqrt, minus_infty, nan_value, INVALID_EXCEPTION);
  3506. TEST_f_f (sqrt, nan_value, nan_value);
  3507. TEST_f_f (sqrt, 2209, 47);
  3508. TEST_f_f (sqrt, 4, 2);
  3509. TEST_f_f (sqrt, 2, M_SQRT2l);
  3510. TEST_f_f (sqrt, 0.25, 0.5);
  3511. TEST_f_f (sqrt, 6642.25, 81.5);
  3512. TEST_f_f (sqrt, 15190.5625L, 123.25L);
  3513. TEST_f_f (sqrt, 0.75L, 0.866025403784438646763723170752936183L);
  3514. END (sqrt);
  3515. }
  3516. static void
  3517. tan_test (void)
  3518. {
  3519. errno = 0;
  3520. FUNC(tan) (0);
  3521. if (errno == ENOSYS)
  3522. /* Function not implemented. */
  3523. return;
  3524. START (tan);
  3525. TEST_f_f (tan, 0, 0);
  3526. TEST_f_f (tan, minus_zero, minus_zero);
  3527. TEST_f_f (tan, plus_infty, nan_value, INVALID_EXCEPTION);
  3528. TEST_f_f (tan, minus_infty, nan_value, INVALID_EXCEPTION);
  3529. TEST_f_f (tan, nan_value, nan_value);
  3530. TEST_f_f (tan, M_PI_4l, 1);
  3531. TEST_f_f (tan, 0.75L, 0.931596459944072461165202756573936428L);
  3532. END (tan);
  3533. }
  3534. static void
  3535. tanh_test (void)
  3536. {
  3537. errno = 0;
  3538. FUNC(tanh) (0.7L);
  3539. if (errno == ENOSYS)
  3540. /* Function not implemented. */
  3541. return;
  3542. START (tanh);
  3543. TEST_f_f (tanh, 0, 0);
  3544. /* vda: uclibc: added IGNORE_ZERO_INF_SIGN to treat -0 as ok */
  3545. TEST_f_f (tanh, minus_zero, minus_zero, IGNORE_ZERO_INF_SIGN);
  3546. #ifndef TEST_INLINE
  3547. TEST_f_f (tanh, plus_infty, 1);
  3548. TEST_f_f (tanh, minus_infty, -1);
  3549. #endif
  3550. TEST_f_f (tanh, nan_value, nan_value);
  3551. TEST_f_f (tanh, 0.75L, 0.635148952387287319214434357312496495L);
  3552. TEST_f_f (tanh, -0.75L, -0.635148952387287319214434357312496495L);
  3553. TEST_f_f (tanh, 1.0L, 0.7615941559557648881194582826047935904L);
  3554. TEST_f_f (tanh, -1.0L, -0.7615941559557648881194582826047935904L);
  3555. /* 2^-57 */
  3556. TEST_f_f (tanh, 0x1p-57L, 6.938893903907228377647697925567626953125e-18L);
  3557. END (tanh);
  3558. }
  3559. #if 0
  3560. static void
  3561. tgamma_test (void)
  3562. {
  3563. errno = 0;
  3564. FUNC(tgamma) (1);
  3565. if (errno == ENOSYS)
  3566. /* Function not implemented. */
  3567. return;
  3568. feclearexcept (FE_ALL_EXCEPT);
  3569. START (tgamma);
  3570. TEST_f_f (tgamma, plus_infty, plus_infty);
  3571. TEST_f_f (tgamma, 0, plus_infty, DIVIDE_BY_ZERO_EXCEPTION);
  3572. TEST_f_f (tgamma, minus_zero, minus_infty, DIVIDE_BY_ZERO_EXCEPTION);
  3573. /* tgamma (x) == NaN plus invalid exception for integer x <= 0. */
  3574. TEST_f_f (tgamma, -2, nan_value, INVALID_EXCEPTION);
  3575. TEST_f_f (tgamma, minus_infty, nan_value, INVALID_EXCEPTION);
  3576. TEST_f_f (tgamma, nan_value, nan_value);
  3577. TEST_f_f (tgamma, 0.5, M_SQRT_PIl);
  3578. TEST_f_f (tgamma, -0.5, -M_2_SQRT_PIl);
  3579. TEST_f_f (tgamma, 1, 1);
  3580. TEST_f_f (tgamma, 4, 6);
  3581. TEST_f_f (tgamma, 0.7L, 1.29805533264755778568117117915281162L);
  3582. TEST_f_f (tgamma, 1.2L, 0.918168742399760610640951655185830401L);
  3583. END (tgamma);
  3584. }
  3585. #endif
  3586. #if 0
  3587. static void
  3588. trunc_test (void)
  3589. {
  3590. START (trunc);
  3591. TEST_f_f (trunc, plus_infty, plus_infty);
  3592. TEST_f_f (trunc, minus_infty, minus_infty);
  3593. TEST_f_f (trunc, nan_value, nan_value);
  3594. TEST_f_f (trunc, 0, 0);
  3595. TEST_f_f (trunc, minus_zero, minus_zero);
  3596. TEST_f_f (trunc, 0.625, 0);
  3597. TEST_f_f (trunc, -0.625, minus_zero);
  3598. TEST_f_f (trunc, 1, 1);
  3599. TEST_f_f (trunc, -1, -1);
  3600. TEST_f_f (trunc, 1.625, 1);
  3601. TEST_f_f (trunc, -1.625, -1);
  3602. TEST_f_f (trunc, 1048580.625L, 1048580L);
  3603. TEST_f_f (trunc, -1048580.625L, -1048580L);
  3604. TEST_f_f (trunc, 8388610.125L, 8388610.0L);
  3605. TEST_f_f (trunc, -8388610.125L, -8388610.0L);
  3606. TEST_f_f (trunc, 4294967296.625L, 4294967296.0L);
  3607. TEST_f_f (trunc, -4294967296.625L, -4294967296.0L);
  3608. #ifdef TEST_LDOUBLE
  3609. /* The result can only be represented in long double. */
  3610. TEST_f_f (trunc, 4503599627370495.5L, 4503599627370495.0L);
  3611. TEST_f_f (trunc, 4503599627370496.25L, 4503599627370496.0L);
  3612. TEST_f_f (trunc, 4503599627370496.5L, 4503599627370496.0L);
  3613. TEST_f_f (trunc, 4503599627370496.75L, 4503599627370496.0L);
  3614. TEST_f_f (trunc, 4503599627370497.5L, 4503599627370497.0L);
  3615. TEST_f_f (trunc, -4503599627370495.5L, -4503599627370495.0L);
  3616. TEST_f_f (trunc, -4503599627370496.25L, -4503599627370496.0L);
  3617. TEST_f_f (trunc, -4503599627370496.5L, -4503599627370496.0L);
  3618. TEST_f_f (trunc, -4503599627370496.75L, -4503599627370496.0L);
  3619. TEST_f_f (trunc, -4503599627370497.5L, -4503599627370497.0L);
  3620. TEST_f_f (trunc, 9007199254740991.5L, 9007199254740991.0L);
  3621. TEST_f_f (trunc, 9007199254740992.25L, 9007199254740992.0L);
  3622. TEST_f_f (trunc, 9007199254740992.5L, 9007199254740992.0L);
  3623. TEST_f_f (trunc, 9007199254740992.75L, 9007199254740992.0L);
  3624. TEST_f_f (trunc, 9007199254740993.5L, 9007199254740993.0L);
  3625. TEST_f_f (trunc, -9007199254740991.5L, -9007199254740991.0L);
  3626. TEST_f_f (trunc, -9007199254740992.25L, -9007199254740992.0L);
  3627. TEST_f_f (trunc, -9007199254740992.5L, -9007199254740992.0L);
  3628. TEST_f_f (trunc, -9007199254740992.75L, -9007199254740992.0L);
  3629. TEST_f_f (trunc, -9007199254740993.5L, -9007199254740993.0L);
  3630. TEST_f_f (trunc, 72057594037927935.5L, 72057594037927935.0L);
  3631. TEST_f_f (trunc, 72057594037927936.25L, 72057594037927936.0L);
  3632. TEST_f_f (trunc, 72057594037927936.5L, 72057594037927936.0L);
  3633. TEST_f_f (trunc, 72057594037927936.75L, 72057594037927936.0L);
  3634. TEST_f_f (trunc, 72057594037927937.5L, 72057594037927937.0L);
  3635. TEST_f_f (trunc, -72057594037927935.5L, -72057594037927935.0L);
  3636. TEST_f_f (trunc, -72057594037927936.25L, -72057594037927936.0L);
  3637. TEST_f_f (trunc, -72057594037927936.5L, -72057594037927936.0L);
  3638. TEST_f_f (trunc, -72057594037927936.75L, -72057594037927936.0L);
  3639. TEST_f_f (trunc, -72057594037927937.5L, -72057594037927937.0L);
  3640. TEST_f_f (trunc, 10141204801825835211973625643007.5L, 10141204801825835211973625643007.0L);
  3641. TEST_f_f (trunc, 10141204801825835211973625643008.25L, 10141204801825835211973625643008.0L);
  3642. TEST_f_f (trunc, 10141204801825835211973625643008.5L, 10141204801825835211973625643008.0L);
  3643. TEST_f_f (trunc, 10141204801825835211973625643008.75L, 10141204801825835211973625643008.0L);
  3644. TEST_f_f (trunc, 10141204801825835211973625643009.5L, 10141204801825835211973625643009.0L);
  3645. #endif
  3646. END (trunc);
  3647. }
  3648. #endif
  3649. static void
  3650. y0_test (void)
  3651. {
  3652. errno = 0;
  3653. #if 0
  3654. FLOAT s, c;
  3655. FUNC (sincos) (0, &s, &c);
  3656. if (errno == ENOSYS)
  3657. /* Required function not implemented. */
  3658. return;
  3659. #endif
  3660. FUNC(y0) (1);
  3661. if (errno == ENOSYS)
  3662. /* Function not implemented. */
  3663. return;
  3664. /* y0 is the Bessel function of the second kind of order 0 */
  3665. START (y0);
  3666. TEST_f_f (y0, -1.0, minus_infty, INVALID_EXCEPTION);
  3667. TEST_f_f (y0, 0.0, minus_infty);
  3668. TEST_f_f (y0, nan_value, nan_value);
  3669. TEST_f_f (y0, plus_infty, 0);
  3670. TEST_f_f (y0, 0.125L, -1.38968062514384052915582277745018693L);
  3671. TEST_f_f (y0, 0.75L, -0.137172769385772397522814379396581855L);
  3672. TEST_f_f (y0, 1.0, 0.0882569642156769579829267660235151628L);
  3673. TEST_f_f (y0, 1.5, 0.382448923797758843955068554978089862L);
  3674. TEST_f_f (y0, 2.0, 0.510375672649745119596606592727157873L);
  3675. TEST_f_f (y0, 8.0, 0.223521489387566220527323400498620359L);
  3676. TEST_f_f (y0, 10.0, 0.0556711672835993914244598774101900481L);
  3677. END (y0);
  3678. }
  3679. static void
  3680. y1_test (void)
  3681. {
  3682. errno = 0;
  3683. #if 0
  3684. FLOAT s, c;
  3685. FUNC (sincos) (0, &s, &c);
  3686. if (errno == ENOSYS)
  3687. /* Required function not implemented. */
  3688. return;
  3689. #endif
  3690. FUNC(y1) (1);
  3691. if (errno == ENOSYS)
  3692. /* Function not implemented. */
  3693. return;
  3694. /* y1 is the Bessel function of the second kind of order 1 */
  3695. START (y1);
  3696. TEST_f_f (y1, -1.0, minus_infty, INVALID_EXCEPTION);
  3697. TEST_f_f (y1, 0.0, minus_infty);
  3698. TEST_f_f (y1, plus_infty, 0);
  3699. TEST_f_f (y1, nan_value, nan_value);
  3700. TEST_f_f (y1, 0.125L, -5.19993611253477499595928744876579921L);
  3701. TEST_f_f (y1, 0.75L, -1.03759455076928541973767132140642198L);
  3702. TEST_f_f (y1, 1.0, -0.781212821300288716547150000047964821L);
  3703. TEST_f_f (y1, 1.5, -0.412308626973911295952829820633445323L);
  3704. TEST_f_f (y1, 2.0, -0.107032431540937546888370772277476637L);
  3705. TEST_f_f (y1, 8.0, -0.158060461731247494255555266187483550L);
  3706. TEST_f_f (y1, 10.0, 0.249015424206953883923283474663222803L);
  3707. END (y1);
  3708. }
  3709. static void
  3710. yn_test (void)
  3711. {
  3712. errno = 0;
  3713. #if 0
  3714. FLOAT s, c;
  3715. FUNC (sincos) (0, &s, &c);
  3716. if (errno == ENOSYS)
  3717. /* Required function not implemented. */
  3718. return;
  3719. #endif
  3720. FUNC(yn) (1, 1);
  3721. if (errno == ENOSYS)
  3722. /* Function not implemented. */
  3723. return;
  3724. /* yn is the Bessel function of the second kind of order n */
  3725. START (yn);
  3726. /* yn (0, x) == y0 (x) */
  3727. TEST_ff_f (yn, 0, -1.0, minus_infty, INVALID_EXCEPTION);
  3728. TEST_ff_f (yn, 0, 0.0, minus_infty);
  3729. TEST_ff_f (yn, 0, nan_value, nan_value);
  3730. TEST_ff_f (yn, 0, plus_infty, 0);
  3731. TEST_ff_f (yn, 0, 0.125L, -1.38968062514384052915582277745018693L);
  3732. TEST_ff_f (yn, 0, 0.75L, -0.137172769385772397522814379396581855L);
  3733. TEST_ff_f (yn, 0, 1.0, 0.0882569642156769579829267660235151628L);
  3734. TEST_ff_f (yn, 0, 1.5, 0.382448923797758843955068554978089862L);
  3735. TEST_ff_f (yn, 0, 2.0, 0.510375672649745119596606592727157873L);
  3736. TEST_ff_f (yn, 0, 8.0, 0.223521489387566220527323400498620359L);
  3737. TEST_ff_f (yn, 0, 10.0, 0.0556711672835993914244598774101900481L);
  3738. /* yn (1, x) == y1 (x) */
  3739. TEST_ff_f (yn, 1, -1.0, minus_infty, INVALID_EXCEPTION);
  3740. TEST_ff_f (yn, 1, 0.0, minus_infty);
  3741. TEST_ff_f (yn, 1, plus_infty, 0);
  3742. TEST_ff_f (yn, 1, nan_value, nan_value);
  3743. TEST_ff_f (yn, 1, 0.125L, -5.19993611253477499595928744876579921L);
  3744. TEST_ff_f (yn, 1, 0.75L, -1.03759455076928541973767132140642198L);
  3745. TEST_ff_f (yn, 1, 1.0, -0.781212821300288716547150000047964821L);
  3746. TEST_ff_f (yn, 1, 1.5, -0.412308626973911295952829820633445323L);
  3747. TEST_ff_f (yn, 1, 2.0, -0.107032431540937546888370772277476637L);
  3748. TEST_ff_f (yn, 1, 8.0, -0.158060461731247494255555266187483550L);
  3749. TEST_ff_f (yn, 1, 10.0, 0.249015424206953883923283474663222803L);
  3750. /* yn (3, x) */
  3751. TEST_ff_f (yn, 3, plus_infty, 0);
  3752. TEST_ff_f (yn, 3, nan_value, nan_value);
  3753. TEST_ff_f (yn, 3, 0.125L, -2612.69757350066712600220955744091741L);
  3754. TEST_ff_f (yn, 3, 0.75L, -12.9877176234475433186319774484809207L);
  3755. TEST_ff_f (yn, 3, 1.0, -5.82151760596472884776175706442981440L);
  3756. TEST_ff_f (yn, 3, 2.0, -1.12778377684042778608158395773179238L);
  3757. TEST_ff_f (yn, 3, 10.0, -0.251362657183837329779204747654240998L);
  3758. /* yn (10, x) */
  3759. TEST_ff_f (yn, 10, plus_infty, 0);
  3760. TEST_ff_f (yn, 10, nan_value, nan_value);
  3761. TEST_ff_f (yn, 10, 0.125L, -127057845771019398.252538486899753195L);
  3762. TEST_ff_f (yn, 10, 0.75L, -2133501638.90573424452445412893839236L);
  3763. TEST_ff_f (yn, 10, 1.0, -121618014.278689189288130426667971145L);
  3764. TEST_ff_f (yn, 10, 2.0, -129184.542208039282635913145923304214L);
  3765. TEST_ff_f (yn, 10, 10.0, -0.359814152183402722051986577343560609L);
  3766. END (yn);
  3767. }
  3768. static void
  3769. significand_test (void)
  3770. {
  3771. /* significand returns the mantissa of the exponential representation. */
  3772. START (significand);
  3773. TEST_f_f (significand, 4.0, 1.0);
  3774. TEST_f_f (significand, 6.0, 1.5);
  3775. TEST_f_f (significand, 8.0, 1.0);
  3776. END (significand);
  3777. }
  3778. static void
  3779. initialize (void)
  3780. {
  3781. fpstack_test ("start *init*");
  3782. plus_zero = 0.0;
  3783. nan_value = plus_zero / plus_zero; /* Suppress GCC warning */
  3784. minus_zero = FUNC(copysign) (0.0, -1.0);
  3785. plus_infty = CHOOSE (HUGE_VALL, HUGE_VAL, HUGE_VALF,
  3786. HUGE_VALL, HUGE_VAL, HUGE_VALF);
  3787. minus_infty = CHOOSE (-HUGE_VALL, -HUGE_VAL, -HUGE_VALF,
  3788. -HUGE_VALL, -HUGE_VAL, -HUGE_VALF);
  3789. max_value = CHOOSE (LDBL_MAX, DBL_MAX, FLT_MAX,
  3790. LDBL_MAX, DBL_MAX, FLT_MAX);
  3791. min_value = CHOOSE (LDBL_MIN, DBL_MIN, FLT_MIN,
  3792. LDBL_MIN, DBL_MIN, FLT_MIN);
  3793. (void) &plus_zero;
  3794. (void) &nan_value;
  3795. (void) &minus_zero;
  3796. (void) &plus_infty;
  3797. (void) &minus_infty;
  3798. (void) &max_value;
  3799. (void) &min_value;
  3800. /* Clear all exceptions. From now on we must not get random exceptions. */
  3801. feclearexcept (FE_ALL_EXCEPT);
  3802. /* Test to make sure we start correctly. */
  3803. fpstack_test ("end *init*");
  3804. }
  3805. #if 0
  3806. /* function to check our ulp calculation. */
  3807. void
  3808. check_ulp (void)
  3809. {
  3810. int i;
  3811. FLOAT u, diff, ulp;
  3812. /* This gives one ulp. */
  3813. u = FUNC(nextafter) (10, 20);
  3814. check_equal (10.0, u, 1, &diff, &ulp);
  3815. printf ("One ulp: % .4" PRINTF_NEXPR "\n", ulp);
  3816. /* This gives one more ulp. */
  3817. u = FUNC(nextafter) (u, 20);
  3818. check_equal (10.0, u, 2, &diff, &ulp);
  3819. printf ("two ulp: % .4" PRINTF_NEXPR "\n", ulp);
  3820. /* And now calculate 100 ulp. */
  3821. for (i = 2; i < 100; i++)
  3822. u = FUNC(nextafter) (u, 20);
  3823. check_equal (10.0, u, 100, &diff, &ulp);
  3824. printf ("100 ulp: % .4" PRINTF_NEXPR "\n", ulp);
  3825. }
  3826. #endif
  3827. int
  3828. main (int argc, char **argv)
  3829. {
  3830. int key;
  3831. verbose = 1;
  3832. output_ulps = 0;
  3833. output_max_error = 1;
  3834. output_points = 1;
  3835. /* XXX set to 0 for releases. */
  3836. ignore_max_ulp = 0;
  3837. /* Parse and process arguments. */
  3838. while ((key = getopt(argc, argv, "fi:puv")) > 0) {
  3839. switch (key)
  3840. {
  3841. case 'f':
  3842. output_max_error = 0;
  3843. break;
  3844. case 'i':
  3845. if (strcmp (optarg, "yes") == 0)
  3846. ignore_max_ulp = 1;
  3847. else if (strcmp (optarg, "no") == 0)
  3848. ignore_max_ulp = 0;
  3849. break;
  3850. case 'p':
  3851. output_points = 0;
  3852. break;
  3853. case 'u':
  3854. output_ulps = 1;
  3855. break;
  3856. case 'v':
  3857. verbose = 3;
  3858. break;
  3859. default:
  3860. fprintf (stderr, "Unknown argument: %c", key);
  3861. exit (EXIT_FAILURE);
  3862. }
  3863. }
  3864. if (optind != argc)
  3865. {
  3866. fprintf (stderr, "wrong number of arguments");
  3867. exit (EXIT_FAILURE);
  3868. }
  3869. if (output_ulps)
  3870. {
  3871. ulps_file = fopen ("ULPs", "a");
  3872. if (ulps_file == NULL)
  3873. {
  3874. perror ("can't open file `ULPs' for writing: ");
  3875. exit (1);
  3876. }
  3877. }
  3878. initialize ();
  3879. printf (TEST_MSG);
  3880. #if 0
  3881. check_ulp ();
  3882. #endif
  3883. /* Keep the tests a wee bit ordered (according to ISO C99). */
  3884. /* Classification macros: */
  3885. fpclassify_test ();
  3886. isfinite_test ();
  3887. isnormal_test ();
  3888. signbit_test ();
  3889. /* Trigonometric functions: */
  3890. acos_test ();
  3891. asin_test ();
  3892. atan_test ();
  3893. atan2_test ();
  3894. cos_test ();
  3895. sin_test ();
  3896. #if 0
  3897. sincos_test ();
  3898. #endif
  3899. tan_test ();
  3900. /* Hyperbolic functions: */
  3901. acosh_test ();
  3902. asinh_test ();
  3903. atanh_test ();
  3904. cosh_test ();
  3905. sinh_test ();
  3906. tanh_test ();
  3907. /* Exponential and logarithmic functions: */
  3908. exp_test ();
  3909. #if 0
  3910. exp10_test ();
  3911. exp2_test ();
  3912. #endif
  3913. expm1_test ();
  3914. frexp_test ();
  3915. ldexp_test ();
  3916. log_test ();
  3917. log10_test ();
  3918. log1p_test ();
  3919. #if 0
  3920. log2_test ();
  3921. #endif
  3922. logb_test ();
  3923. modf_test ();
  3924. ilogb_test ();
  3925. scalb_test ();
  3926. scalbn_test ();
  3927. #if 0
  3928. scalbln_test ();
  3929. #endif
  3930. significand_test ();
  3931. /* Power and absolute value functions: */
  3932. cbrt_test ();
  3933. fabs_test ();
  3934. hypot_test ();
  3935. pow_test ();
  3936. sqrt_test ();
  3937. /* Error and gamma functions: */
  3938. erf_test ();
  3939. erfc_test ();
  3940. gamma_test ();
  3941. lgamma_test ();
  3942. #if 0
  3943. tgamma_test ();
  3944. #endif
  3945. /* Nearest integer functions: */
  3946. ceil_test ();
  3947. floor_test ();
  3948. #if 0
  3949. nearbyint_test ();
  3950. #endif
  3951. rint_test ();
  3952. #if 0
  3953. rint_test_tonearest ();
  3954. rint_test_towardzero ();
  3955. rint_test_downward ();
  3956. rint_test_upward ();
  3957. lrint_test ();
  3958. llrint_test ();
  3959. round_test ();
  3960. lround_test ();
  3961. llround_test ();
  3962. trunc_test ();
  3963. #endif
  3964. /* Remainder functions: */
  3965. fmod_test ();
  3966. remainder_test ();
  3967. #if 0
  3968. remquo_test ();
  3969. #endif
  3970. /* Manipulation functions: */
  3971. copysign_test ();
  3972. #if 0
  3973. nextafter_test ();
  3974. nexttoward_test ();
  3975. /* maximum, minimum and positive difference functions */
  3976. fdim_test ();
  3977. fmax_test ();
  3978. fmin_test ();
  3979. /* Multiply and add: */
  3980. fma_test ();
  3981. /* Complex functions: */
  3982. cabs_test ();
  3983. cacos_test ();
  3984. cacosh_test ();
  3985. carg_test ();
  3986. casin_test ();
  3987. casinh_test ();
  3988. catan_test ();
  3989. catanh_test ();
  3990. ccos_test ();
  3991. ccosh_test ();
  3992. cexp_test ();
  3993. cimag_test ();
  3994. clog10_test ();
  3995. clog_test ();
  3996. #if 0
  3997. conj_test ();
  3998. #endif
  3999. cpow_test ();
  4000. cproj_test ();
  4001. creal_test ();
  4002. csin_test ();
  4003. csinh_test ();
  4004. csqrt_test ();
  4005. ctan_test ();
  4006. ctanh_test ();
  4007. #endif
  4008. /* Bessel functions: */
  4009. #if 0
  4010. j0_test ();
  4011. j1_test ();
  4012. jn_test ();
  4013. y0_test ();
  4014. y1_test ();
  4015. yn_test ();
  4016. #endif
  4017. if (output_ulps)
  4018. fclose (ulps_file);
  4019. printf ("\nTest suite completed:\n");
  4020. printf (" %d test cases plus %d tests for exception flags executed.\n",
  4021. noTests, noExcTests);
  4022. if (noXFails)
  4023. printf (" %d expected failures occurred.\n", noXFails);
  4024. if (noXPasses)
  4025. printf (" %d unexpected passes occurred.\n", noXPasses);
  4026. if (noErrors)
  4027. {
  4028. printf (" %d errors occurred.\n", noErrors);
  4029. return 1;
  4030. }
  4031. printf (" All tests passed successfully.\n");
  4032. return 0;
  4033. }
  4034. /*
  4035. * Local Variables:
  4036. * mode:c
  4037. * End:
  4038. */