libm-test.inc 158 KB

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