DQFUN quick check of operations and functions

Test data:

t1 = pi:
3.141592653589793238462643383279502884197169399375105820974944592307791e0       
t2 = -log(2):
-6.931471805599453094172321214581765680755001343602552541206800094933957e-1     
z1 = (0.5*pi, exp(0.5)):
1.570796326794896619231321691639751442098584699687552910487472296153897e0       
1.648721270700128146848650787814163571653776100710148011575079311639945e0       
z2 = (-Gamma(0.5), Cos(1)):
-1.772453850905516027298167483341145182797549456122387128213807789852874e0      
5.403023058681397174009366074429766037323104206179222276700972553812593e-1      
e1 = 3141/8192:
    3.83422851562500000000000000000000000D-01
e2 = 6931 / 8192:
    8.46069335937500000000000000000000000D-01
ec1 = (3141/8192, 6931/8192)
    3.83422851562500000000000000000000000D-01
    8.46069335937500000000000000000000000D-01
ec2 = (6931/8192, 3141/8192):
   -8.46069335937500000000000000000000000D-01
    3.83422851562500000000000000000000000D-01

Real data operations:

addition: t1+t2 =
2.448445473029847929045411261821326316121669265014850566854264582814371e0       
addition: t1+e2 =
3.987661989527293238462643383279502884197169399375105820974944592307794e0       
addition: e1+t2 =
-3.097243289974453094172321214581765680755001343602552541206800094933956e-1     
subtraction: t1-t2 =
3.834739834149738547879875504737679452272669533735361075095624601801167e0       
subtraction: t1-e2 =
2.295523317652293238462643383279502884197169399375105820974944592307791e0       
subtraction: e1-t2 =
1.076570032122445309417232121458176568075500134360255254120680009493394e0       
multiplication: t1*t2 =
-2.177586090303602130500688898237613947338583700369286294325795253194295e0      
multiplication: t1*e2 =
2.658005210208844840794016270692167296187815076546491509421062130039701e0       
multiplication: e1*t2 =
-2.657684685228012960058015250854654053131281643097609562003242077415467e-1     
division: t1/t2 =
-4.532360141827193809627682945716666810171861467723795584186016547940537e0      
division: t1/e2 =
3.713162172588022826357809060139328758814487335114827136838370523760693e0       
division: e1/t2 =
-5.531622465127238845846586209749686527901727224114415062401668197229706e-1     
exponentiation: t1**i1 =
3.060196847852814532627413100434356064803007066280749905534924436862213e2       
exponentiation: t1**t2 =
4.522743596675148667436794168649302701633236076842417709449509543432488e-1      
equal test: t1 == t2, e1 == t2, t1 == e2
   F   F   F
not-equal test: t1 /= t2, e1 /= t2, t1 =/ e2
   T   T   T
less-than-or-equal test: t1 <= t2, e1 <= t2, t1 <= e2
   F   F   F
greater-than-or-equal test: t1 >= t2, e1 >= t2, t1 >= e2
   T   T   T
less-than test: t1 < t2, e1 < t2, t1 < e2
   F   F   F
greater-than test: t1 > t2, e1 > t2, t1 > e2
   T   T   T
abs(t2) =
6.931471805599453094172321214581765680755001343602552541206800094933957e-1      
acos(t2) =
2.336642521613976834675716408789519705045173801836318809655685251968249e0       
acosh(t1) =
1.811526272460853107021852049305420510220702081057922474861595622974743e0       
agm(t1,abs(t2)) =
1.689291197815604730290885986081005314786426904408367845571534212726961e0       
aint(t1) =
3.000000000000000000000000000000000000000000000000000000000000000000000e0       
anint(t1) =
3.000000000000000000000000000000000000000000000000000000000000000000000e0       
asin(t2) =
-7.658461948190802154443947171497682629465891021487658991682129558148067e-1     
asinh(t1) =
1.862295743310848219888361325182620574902674184961554765612879514424698e0       
atan(t1) =
1.262627255678911683444322083605698343508947670424383596973809952252131e0       
atan2(t1,t2) =
1.787952808146344472915858377805388991973691256720279064632274601883306e0       
atanh(t2) =
-8.539880479975240067710412491943246423317581329392131056994132774145449e-1     
bessel_i(t1,-t2) =
5.130642193812675854289557060873826757735015116248008989545595298990872e-3      
bessel_in(3,t1) =
1.157485470370689730634534574426664942067969051104256539648672369942761e0       
bessel_in(3,50*t1) =
5.123106847141514880796133624509529692484608379983241653530211773542147e66      
bessel_j(t1,-t2) =
4.841508727944457944998908303542887700696401281506000138937950029198399e-3      
bessel_jn(3,t1) =
3.334583362029895353902185817398867483022440044539132183547815705509442e-1      
bessel_jn(3,50*t1) =
4.374859995867834281590427180474704719443236913408894823507678618262503e-2      
bessel_k(t1,-t2) =
3.022519056795630839310461514145783446622819281661286085524319185156057e1       
bessel_kn(3,t1) =
9.897437621437765384270398796055839340999519968489535071859127597785593e-2      
bessel_kn(3,50*t1) =
6.212118292148577739193828918353542133962529838069671627703602406114359e-70     
bessel_y(t1,-t2) =
-2.153400705916899924604164743623541376239146851380688212271783491830878e1      
bessel_yn(3,t1) =
-4.860704562984321414013025441338746300148774425680023496916586377635847e-1     
bessel_yn(3,50*t1) =
4.625609271829441376706267446468629294193265850060158080327671133141398e-2      
cos(t2) =
7.692389013639721265783299936612707014408959949119638531698715074289018e-1      
cosh(t1) =
1.159195327552152062775175205256013769577091717620542253821288304846468e1       
digamma_be(nbe,be,t1) =
9.772133079420067332920694864061823436408346099943256380095232865324421e-1      
dqcssh(t1) =
1.159195327552152062775175205256013769577091717620542253821288304846468e1       
1.154873935725774837797733431538840968449518906639478945523216336106371e1       
dqcssn(t2) =
7.692389013639721265783299936612707014408959949119638531698715074289018e-1      
-6.389612763136348011500329114647017842572305378305797294955869566465273e-1     
dqeform(t1,n1,n2,chr1) =
3.141592653589793238462643383279502884197169399375105820974944592307791e0       
dqfform(t1,n1,n2,chr1) =
        3.1415926535897932384626433832795028841971693993751058209749445923077911
dqnrt(t1,i1) =
1.257274115669185059384522114110448293906166310039658173534181622627157e0       
dqreal (chr120) =
3.141592653589793238462643383279502884197169399375105820974944592307791e0       
erf(t1) =
9.999911238536323583947316207812029447123820815128765990475863916464051e-1      
erfc(t1) =
8.876146367641605268379218797055287617918487123400952413608353592671352e-6      
exp(t1) =
2.314069263277926900572908636794854738026610624260021199344504640952834e1       
expint(t1) =
1.092837438933141034863844590690753517156633883505578841248057416030997e1       
gamma(t1) =
2.288037795340032417959588909060233922889688153356222441199380745467604e0       
hurwitz_zetan(3,1/t1) =
3.158409420922361280724292999398432297405053423434495301095473674960161e1       
hurwitz_zetan_be(nbe,be,5,t1) =
4.607288802411928060528967862644377638954673528307916219793433944913618e-3      
hypergeom_pfq(np,nq,aa,bb,t1) =
2.490587773704512509448250304178393529239999768813747978841470938888019e0       
incgamma(t1,t2+2) =
2.005126486991562070180754490480130943426440840883982772740395204028099e0       
incgamma(-3,t1) =
2.096769170449034722999282570976305875330320834572594492656318729484928e-4      
log(t1) =
1.144729885849400174143427351353058711647294812915311571513623071472723e0       
log10(t1) =
4.971498726941338543512682882908988736516783243804424461340534999252671e-1      
max(t1,t2) =
3.141592653589793238462643383279502884197169399375105820974944592307791e0       
min(t1,t2) =
-6.931471805599453094172321214581765680755001343602552541206800094933957e-1     
mod(t1,t2) =
3.690039313500120007937148974467966118951688619340848044922245543342143e-1      
polygamma(3,1/t1) =
5.867329744355181201226190645218983678756215473427307405990288143669559e2       
polygamma_be(nbe,be,5,t1) =
1.591352138903679163757538458364354595126832313926580413273307986658081e-1      
polylog_neg (-10, rr, -t1) =
-8.025408854473784508947588772003746749247905347018584022236183194583041e0      
polylog_pos (10, 1/t1) =
3.184093889782650585151509188778627801898787499438544547844038745175928e-1      
sign(t1,t2) =
-3.141592653589793238462643383279502884197169399375105820974944592307791e0      
sin(t2) =
-6.389612763136348011500329114647017842572305378305797294955869566465273e-1     
sinh(t1) =
1.154873935725774837797733431538840968449518906639478945523216336106371e1       
sqrt(t1) =
1.772453850905516027298167483341145182797549456122387128213807789852874e0       
struve_hn(10,t1) =
1.181983397769752114529000299689581613912341555707830352880505317267653e-5      
tan(t2) =
-8.306408778607839470304590233006621122265435930965995358779049026972790e-1     
tanh(t1) =
9.962720762207499442646905800125367118968991908045876143626124159785444e-1      
zeta(t1) =
1.176241738382582758872150451938052091169738990021655834960508346219137e0       
zeta_be(nbe,be,t1) =
1.176241738382582758872150451938052091169738990021655834960508346230422e0       
zeta_int(10) =
1.000994575127818085337145958900319017006019531564477517257788994636262e0       

Complex data operations:

addition: z1+z2 =
-2.016575241106194080668457917013937406989647564348342177263354936989816e-1     
2.189023576568267864249587395257140175386086521328070239245176567021200e0       
addition: z1+e2 =
2.416865662732396619231321691639751442098584699687552910487472296153897e0       
1.648721270700128146848650787814163571653776100710148011575079311639945e0       
addition: e1+z2 =
-1.389030999343016027298167483341145182797549456122387128213807789852878e0      
5.403023058681397174009366074429766037323104206179222276700972553812593e-1      
addition: z1+ec2 =
7.247269908573966192313216916397514420985846996875529104874722961538985e-1      
2.032144122262628146848650787814163571653776100710148011575079311639944e0       
addition: ec1+z2 =
-1.389030999343016027298167483341145182797549456122387128213807789852878e0      
1.386371641805639717400936607442976603732310420617922227670097255381259e0       
addition: z1+t2 =
8.776491462349513098140895701815748740230845653272976563667922866604953e-1      
1.648721270700128146848650787814163571653776100710148011575079311639945e0       
addition: t1+z2 =
1.369138802684277211164475899938357701399619943252718692761136802454914e0       
5.403023058681397174009366074429766037323104206179222276700972553812593e-1      
subtraction: z1-z2 =
3.343250177700412646529489174980896624896134155809940038701280086006768e0       
1.108418964831988429447714180371186967921465680092225783904982056258682e0       
subtraction: z1-e2 =
7.247269908573966192313216916397514420985846996875529104874722961538985e-1      
1.648721270700128146848650787814163571653776100710148011575079311639945e0       
subtraction: e1-z2 =
2.155876702468016027298167483341145182797549456122387128213807789852878e0       
-5.403023058681397174009366074429766037323104206179222276700972553812593e-1     
subtraction: z1-ec2 =
2.416865662732396619231321691639751442098584699687552910487472296153897e0       
1.265298419137628146848650787814163571653776100710148011575079311639945e0       
subtraction: ec1-z2 =
2.155876702468016027298167483341145182797549456122387128213807789852878e0       
3.057670300693602825990633925570233962676895793820777723299027446187406e-1      
subtraction: z1-t2 =
2.263943507354841928648553813097928010174084834047808164608152305647283e0       
1.648721270700128146848650787814163571653776100710148011575079311639945e0       
subtraction: t1-z2 =
4.914046504495309265760810866620648066994718855497492949188752382160662e0       
-5.403023058681397174009366074429766037323104206179222276700972553812593e-1     
multiplication: z1*z2 =
-3.674971902708982542198635852017981319841852850774359696725603445591010e0      
-2.073577487905791279696081819605453953259704351198766508148485917823680e0      
multiplication: z1*e2 =
1.329002605104422420397008135346083648093907538273245754710531065019850e0       
1.394932510647288596900390455363765590226113544192143050320663416623106e0       
multiplication: e1*z2 =
-6.795993097771271779472099688933761009725467336035666466942834799716652e-1     
2.071642508217562075630300151340807510160140418897575338271210301699865e-1      
multiplication: z1*ec2 =
-1.961160016148014035295886568515552371063208164778442886590087231750941e0      
-7.926533038036887090823263155394907758179525893925088105939604768990666e-1     
multiplication: ec1*z2 =
-1.136732522908484158087089317793189435990691695188383414270143352658122e0      
-1.292454601793738370708893669807017547564285674954024837468032997109100e0      
multiplication: z1*t2 =
-1.088793045151801065250344449118806973669291850184643147162897626597147e0      
-1.142806500315004192691723555372366060570162479396662094822552632300805e0      
multiplication: t1*z2 =
-5.568327996831707845284817982118835702013624390283243910753675818829579e0      
1.697409754832973169691082513004456724276406402895600701738508187293934e0       
multiplication: z1*e2 =
1.329002605104422420397008135346083648093907538273245754710531065019850e0       
1.394932510647288596900390455363765590226113544192143050320663416623106e0       
multiplication: e1*z2 =
-6.795993097771271779472099688933761009725467336035666466942834799716652e-1     
2.071642508217562075630300151340807510160140418897575338271210301699865e-1      
division: z1/z2 =
-5.514330820250523530931949267834197475150248203763202584252541652219403e-1     
-1.098286330815170779114897657817863265978541105684486762844034968131209e0      
division: z1/e2 =
1.856581086294011413178904530069664379407243667557413568419185261880350e0       
1.948683400602431074734403008768378008799269054540114342926424718071622e0       
division: e1/z2 =
-1.979308293330522397187384821831121789177743189239603821511400268951152e-1     
-6.033583522437342232153716153648717126869711198333166271431973341291785e-2     
division: z1/ec2 =
-8.076120970971878602268957718058416627557671161372861017766276741107253e-1     
-2.314678148399612876634948329680533348968921992325024982903395937864141e0      
division: ec1/z2 =
-6.479244221425816454822779558349658386418002310846475837751695389768576e-2     
-4.970943764874695934423828858316640296192522869155091248886436832410103e-1     
division: z1/t2 =
-2.266180070913596904813841472858333405085930733861897792093008273970268e0      
-2.378602001047225083411310792448537021398121087403685302094245105621295e0      
division: t1/z2 =
-1.621755293972853794073402767310224041921312741308016934712892864246883e0      
-4.943644227688866408711781423607018856899370201434641400233363405503943e-1     
exponentiation: z1**i1 =
-3.775812204234204943393779724040022886179695100824356159019756726468250e1      
-4.821111386306407891048412017981088073717203203291274845888293013796447e1      
exponentiation: z1**z2 =
8.234832447665880835786402427344714596914545008614885681605697787613370e-2      
-1.255630165140989591866311827573673007364534853816634275905725883238808e-1     
exponentiation: t1**z2 =
1.071145872688508480441765122284336425817904314402935676387561158603185e-1      
7.622779877453909624370657747910955969219249825585782434207429682069564e-2      
exponentiation: z1**t2 =
4.785885594636607983440331474875715798354128945336830837561568670005892e-1      
-3.008315771303482313251153693152478986309470691567529325616472853231124e-1     
equal test: z1 == z2, e1 == z2, z1 == e2, ec1 == z2, z1 == ec2
   F   F   F   F   F
not-equal test: z1 /= z2, e1 /= z2, z1 /= e2, ec1 /= z2, z1 /= ec2
   T   T   T   T   T
abs(z2) =
1.852975778394370174224784560084798599325600563420540889070602430996008e0       
aimag(z1) =
1.648721270700128146848650787814163571653776100710148011575079311639945e0       
conjg(z1) =
1.570796326794896619231321691639751442098584699687552910487472296153897e0       
-1.648721270700128146848650787814163571653776100710148011575079311639945e0      
cos(z2) =
-2.302472209102167852684774909189701960708700876369202705063348330117822e-1     
5.554875383367798714195804874637118062469751300894303666552491331747987e-1      
exp(z1) =
-3.744769222334551763165715589581777316716078231959876588248279456022400e-1     
4.795879446722345185222240097567601959373325781485449528593383999312460e0       
log(z1) =
8.229507726980497315430084393387342882323535416361858417229737278108838e-1      
8.095973578347321871164899441671339362041949905939811919425943859405850e-1      
dqcmplx(t1,t2) =
3.141592653589793238462643383279502884197169399375105820974944592307791e0       
-6.931471805599453094172321214581765680755001343602552541206800094933957e-1     
dqreal(z1) =
1.570796326794896619231321691639751442098584699687552910487472296153897e0       
sin(z2) =
-1.126254368529469836305963119762696606310488214252259580019015219515298e0      
-1.135617898817095005030327882390028597272728730558821608520216111174161e-1     
sqrt(z1) =
1.387084313889587122047219207862690134580459173107124175034011813693765e0       
5.943118432638289891586984894975830748864270140914484635893821370046578e-1      
sqrt(z2) =
2.006513487231697465105842600759502813354721536377801153996939721287827e-1      
1.346370979578044492827463338403112065365298334792736019513456788069612e0       

Max relative error =   4.248079D-62
ALL TESTS PASSED
