19 Sep 2026 [krishna@varuthini ddfun-v04-ftn95]$ WINEDEBUG=-all wine testddfun.exe ddknd = 2 ddpicon(1) = 0.31415926535897931E+01 ddpicon(2) = 0.12246467991473532E-15 DDFUN quick check of operations and functions Test data: -- t1 = pi: 3.1415926535897932384626433832795e0 FTN95 3.1415926535897932384626433832795e0 gfortran 3.1415926535897932384626433832795e0 ref -- t2 = -log(2): -6.9314718055994530941723212145817e-1 FTN95 -6.9314718055994530941723212145817e-1 gfortran -6.9314718055994530941723212145817e-1 ref -- z1 = (0.5*pi, exp(0.5)): 1.5707963267948966192313216916397e0 FTN95 1.6487212707001281468486507878140e0 FTN95 1.5707963267948966192313216916397e0 gfortran 1.6487212707001281468486507878140e0 gfortran 1.5707963267948966192313216916397e0 ref 1.6487212707001281468486507878140e0 ref -- z2 = (-Gamma(0.5), Cos(1)): -1.7724538509055160272981674833411e0 FTN95 5.4030230586813971740093660744290e-1 FTN95 -1.7724538509055160272981674833411e0 gfortran 5.4030230586813971740093660744290e-1 gfortran -1.7724538509055160272981674833411e0 ref 5.4030230586813971740093660744290e-1 ref -- e1 = 3141/8192: 3.83422851562500000000000000000000000D-01 FTN95 3.83422851562500000000000000000000000D-01 gfortran 3.83422851562500000000000000000000000D-01 ref -- e2 = 6931 / 8192: 8.46069335937500000000000000000000000D-01 FTN95 8.46069335937500000000000000000000000D-01 gfortran 8.46069335937500000000000000000000000D-01 ref -- ec1 = (3141/8192, 6931/8192) 3.83422851562500000000000000000000000D-01 FTN95 8.46069335937500000000000000000000000D-01 FTN95 3.83422851562500000000000000000000000D-01 gfortran 8.46069335937500000000000000000000000D-01 gfortran 3.83422851562500000000000000000000000D-01 ref 8.46069335937500000000000000000000000D-01 ref -- ec2 = (6931/8192, 3141/8192): -8.46069335937500000000000000000000000D-01 FTN95 3.83422851562500000000000000000000000D-01 FTN95 -8.46069335937500000000000000000000000D-01 gfortran 3.83422851562500000000000000000000000D-01 gfortran -8.46069335937500000000000000000000000D-01 ref 3.83422851562500000000000000000000000D-01 ref -- Real data operations: addition: t1+t2 = 2.4484454730298479290454112618213e0 FTN95 2.4484454730298479290454112618213e0 gfortran addition: t1+e2 = 3.9876619895272932384626433832795e0 FTN95 3.9876619895272932384626433832795e0 gfortran addition: e1+t2 = -3.0972432899744530941723212145817e-1 FTN95 -3.0972432899744530941723212145817e-1 gfortran subtraction: t1-t2 = 3.8347398341497385478798755047376e0 FTN95 3.8347398341497385478798755047376e0 gfortran subtraction: t1-e2 = 2.2955233176522932384626433832795e0 FTN95 2.2955233176522932384626433832795e0 gfortran subtraction: e1-t2 = 1.0765700321224453094172321214581e0 FTN95 1.0765700321224453094172321214581e0 gfortran multiplication: t1*t2 = -2.1775860903036021305006888982375e0 FTN95 -2.1775860903036021305006888982375e0 gfortran multiplication: t1*e2 = 2.6580052102088448407940162706921e0 FTN95 2.6580052102088448407940162706921e0 gfortran multiplication: e1*t2 = -2.6576846852280129600580152508546e-1 FTN95 -2.6576846852280129600580152508546e-1 gfortran division: t1/t2 = -4.5323601418271938096276829457166e0 FTN95 -4.5323601418271938096276829457166e0 gfortran division: t1/e2 = 3.7131621725880228263578090601393e0 FTN95 3.7131621725880228263578090601393e0 gfortran division: e1/t2 = -5.5316224651272388458465862097496e-1 FTN95 -5.5316224651272388458465862097496e-1 gfortran exponentiation: t1**i1 = 3.0601968478528145326274131004343e2 FT95 3.0601968478528145326274131004343e2 gfortran exponentiation: t1**t2 = 4.5227435966751486674367941686477e-1 FTN95 4.5227435966751486674367941686475e-1 gfortran 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.9314718055994530941723212145817e-1 FTN95 6.9314718055994530941723212145817e-1 gfortran acos(t2) = 2.3366425216139768346757164087890e0 FTN95 2.3366425216139768346757164087890e0 gfortran acosh(t1) = 1.8115262724608531070218520493050e0 FTN95 ERROR: 7.484598D-31 1.8115262724608531070218520493064e0 gfortran agm(t1,abs(t2)) = 1.6892911978156047302908859860809e0 FTN95 1.6892911978156047302908859860809e0 gfortran aint(t1) = 3.0000000000000000000000000000000e0 FTN95 3.0000000000000000000000000000000e0 gfortran anint(t1) = 3.0000000000000000000000000000000e0 FTN95 3.0000000000000000000000000000000e0 gfortran asin(t2) = -7.6584619481908021544439471714795e-1 FTN95 ERROR: 5.492266D-31 -7.6584619481908021544439471714795e-1 gfortran ERROR: 5.492266D-31 asinh(t1) = 1.8622957433108482198883613251825e0 FTN95 1.8622957433108482198883613251826e0 gfortran atan(t1) = 1.2626272556789116834443220836060e0 FTN95 1.2626272556789116834443220836060e0 gfortran atan2(t1,t2) = 1.7879528081463444729158583778055e0 FTN95 1.7879528081463444729158583778055e0 gfortran atanh(t2) = -8.5398804799752400677104124919429e-1 FTN95 ERROR: 1.876342D-31 -8.5398804799752400677104124919447e-1 gfortran bessel_i(t1,-t2) = 5.1306421938126758542895570608681e-3 FTN95 ERROR: 1.501512D-30 5.1306421938126758542895570608606e-3 gfortran bessel_in(3,t1) = 1.1574854703706897306345345744266e0 FTN95 1.1574854703706897306345345744266e0 gfortran bessel_in(3,50*t1) = 5.1231068471415148807961336245119e66 FTN95 5.1231068471415148807961336245119e66 gfortran bessel_j(t1,-t2) = 4.8415087279444579449989083035375e-3 FTN95 ERROR: 1.541457D-30 4.8415087279444579449989083035302e-3 gfortran bessel_jn(3,t1) = 3.3345833620298953539021858173988e-1 FTN95 3.3345833620298953539021858173988e-1 gfortran bessel_jn(3,50*t1) = 4.3748599958678342815904271804742e-2 FTN95 ERROR: 1.646448D-30 4.3748599958678342815904271804671e-2 bessel_k(t1,-t2) = 3.0225190567956308393104615141416e1 FTN95 ERROR: 4.567404D-31 3.0225190567956308393104615141430e1 gfortran bessel_kn(3,t1) = 9.8974376214377653842703987960065e-2 FTN95 ERROR: 1.879727D-29 9.8974376214377653842703987961926e-2 gfortran bessel_kn(3,50*t1) = 6.2121182921485777391938289183505e-70 FTN95 6.2121182921485777391938289183505e-70 bessel_y(t1,-t2) = -2.1534007059168999246041647436205e1 FTN95 ERROR: 4.602053D-31 -2.1534007059168999246041647436215e1 gfortran bessel_yn(3,t1) = -4.8607045629843214140130254413396e-1 FTN95 ERROR: 7.353926D-31 -4.8607045629843214140130254413362e-1 gfortran bessel_yn(3,50*t1) = 4.6256092718294413767062674464689e-2 FTN95 ERROR: 1.482250D-30 4.6256092718294413767062674464758e-2 gfortran cos(t2) = 7.6923890136397212657832999366221e-1 FTN95 ERROR: 1.906805D-30 7.6923890136397212657832999366076e-1 gfortran cosh(t1) = 1.1591953275521520627751752052560e1 FTN95 ERROR: 3.636553D-31 1.1591953275521520627751752052565e1 gfortran digamma_be(nbe,be,t1) = 9.7721330794200673329206948640576e-1 FTN95 9.7721330794200673329206948640578e-1 gfortran ddcssh(t1) = 1.1591953275521520627751752052560e1 FTN95 ERROR: 3.636553D-31 1.1548739357257748377977334315388e1 FTN95 ERROR: 3.500739D-31 1.1591953275521520627751752052565e1 gfortran 1.1548739357257748377977334315393e1 gfortran ddcssn(t2) = 7.6923890136397212657832999366221e-1 FTN95 ERROR: 1.906805D-30 -6.3896127631363480115003291146355e-1 FTN95 ERROR: 2.719976D-30 7.6923890136397212657832999366076e-1 gfortran -6.3896127631363480115003291146531e-1 gfortran ddeform(t1,n1,n2,chr1) = 3.1415926535897932384626433832795e0 FTN95 3.1415926535897932384626433832795e0 gfortran ddfform(t1,n1,n2,chr1) = 3.14159265358979323846264338327950 FTN95 3.14159265358979323846264338327950 gfortran ddnrt(t1,i1) = 1.2572741156691850593845221141104e0 FTN95 1.2572741156691850593845221141104e0 gfortran ddreal (chr120) = 3.1415926535897932384626433832794e0 FTN95 3.1415926535897932384626433832794e0 gfortran erf(t1) = 9.9999112385363235839473162078226e-1 FTN95 ERROR: 9.121285D-31 9.9999112385363235839473162078137e-1 gfortran erfc(t1) = 8.8761463676416052683792177343467e-6 FTN95 ERROR: 9.998354D-26 8.8761463676416052683792186218152e-6 gfortran exp(t1) = 2.3140692632779269005729086367949e1 FTN95 ERROR: 3.920324D-31 2.3140692632779269005729086367959e1 gfortran expint(t1) = 1.0928374389331410348638445906920e1 FTN95 ERROR: 1.624155D-31 1.0928374389331410348638445906919e1 gfortran gamma(t1) = 2.2880377953400324179595889090664e0 FTN95 2.2880377953400324179595889090663e0 gfortran hurwitz_zetan(3,1/t1) = 3.1584094209223612807242929993908e1 FTN95 ERROR: 1.373334D-28 3.1584094209223612807242929989571e1 gfortran hurwitz_zetan_be(nbe,be,5,t1) = 4.6072888024119280605289678626442e-3 FTN95 4.6072888024119280605289678626444e-3 gfortran hypergeom_pfq(np,nq,aa,bb,t1) = 2.4905877737045125094482503041784e0 FTN95 2.4905877737045125094482503041784e0 gfortran incgamma(t1,t2+2) = 2.0051264869915620701807544904875e0 FTN95 ERROR: 1.967110D-31 2.0051264869915620701807544904879e0 gfortran incgamma(-3,t1) = 2.0967691704490347229992825696212e-4 FTN95 ERROR: 4.337149D-29 2.0967691704490347229992825695303e-4 gfortran log(t1) = 1.1447298858494001741434273513530e0 FTN95 1.1447298858494001741434273513530e0 gfortran log10(t1) = 4.9714987269413385435126828829091e-1 FTN95 ERROR: 6.126640D-31 4.9714987269413385435126828829123e-1 gfortran max(t1,t2) = 3.1415926535897932384626433832795e0 FTN95 3.1415926535897932384626433832795e0 gfortran min(t1,t2) = -6.9314718055994530941723212145817e-1 FTN95 -6.9314718055994530941723212145817e-1 gfortran mod(t1,t2) = 3.6900393135001200079371489744680e-1 FTN95 3.6900393135001200079371489744680e-1 gfortran polygamma(3,1/t1) = 5.8673297443551812012261906460730e2 FTN95 ERROR: 2.670064D-28 5.8673297443551812012261906445065e2 gfortran polygamma_be(nbe,be,5,t1) = 1.5913521389036791637575384583643e-1 FTN95 1.5913521389036791637575384583643e-1 gfortran polylog_neg (-10, rr, -t1) = -8.0254088544737845089475887720048e0 FTN95 -8.0254088544737845089475887720048e0 gfortran polylog_pos (10, 1/t1) = 3.1840938897826505851515091887786e-1 FTN95 3.1840938897826505851515091887786e-1 gfortran sign(t1,t2) = -3.1415926535897932384626433832795e0 FTN95 -3.1415926535897932384626433832795e0 gfortran sin(t2) = -6.3896127631363480115003291146355e-1 FTN95 ERROR: 2.719976D-30 -6.3896127631363480115003291146531e-1 gfortran sinh(t1) = 1.1548739357257748377977334315388e1 FTN95 ERROR: 3.500739D-31 1.1548739357257748377977334315393e1 gfortran sqrt(t1) = 1.7724538509055160272981674833411e0 FTN95 1.7724538509055160272981674833411e0 gfortran struve_hn(10,t1) = 1.1819833977697521145290002996895e-5 FTN95 1.1819833977697521145290002996895e-5 gfortran tan(t2) = -8.3064087786078394703045902329814e-1 FTN95 ERROR: 4.629795D-30 -8.3064087786078394703045902330201e-1 gfortran tanh(t1) = 9.9627207622074994426469058001252e-1 FTN95 9.9627207622074994426469058001253e-1 gfortran zeta(t1) = 1.1762417383825827588721504519380e0 FTN95 1.1762417383825827588721504519380e0 gfortran zeta_be(nbe,be,t1) = 1.1762417383825827588721504519381e0 FTN95 ERROR: 1.100305D-31 1.1762417383825827588721504519380e0 gfortran zeta_int(10) = 1.0009945751278180853371459589003e0 FTN95 1.0009945751278180853371459589003e0 gfortran Complex data operations: addition: z1+z2 = -2.0165752411061940806684579170139e-1 FTN95 2.1890235765682678642495873952570e0 FTN95 addition: z1+e2 = 2.4168656627323966192313216916397e0 FTN95 1.6487212707001281468486507878140e0 FTN95 addition: e1+z2 = -1.3890309993430160272981674833411e0 FTN95 5.4030230586813971740093660744290e-1 FTN95 addition: z1+ec2 = 7.2472699085739661923132169163975e-1 FTN95 2.0321441222626281468486507878140e0 FTN95 addition: ec1+z2 = -1.3890309993430160272981674833411e0 FTN95 1.3863716418056397174009366074429e0 FTN95 addition: z1+t2 = 8.7764914623495130981408957018157e-1 FTN95 1.6487212707001281468486507878140e0 FTN95 addition: t1+z2 = 1.3691388026842772111644758999383e0 FTN95 5.4030230586813971740093660744290e-1 FTN95 subtraction: z1-z2 = 3.3432501777004126465294891749808e0 FTN95 1.1084189648319884294477141803711e0 FTN95 subtraction: z1-e2 = 7.2472699085739661923132169163975e-1 FTN95 1.6487212707001281468486507878140e0 FTN95 subtraction: e1-z2 = 2.1558767024680160272981674833411e0 FTN95 -5.4030230586813971740093660744290e-1 FTN95 subtraction: z1-ec2 = 2.4168656627323966192313216916397e0 FTN95 1.2652984191376281468486507878140e0 FTN95 subtraction: ec1-z2 = 2.1558767024680160272981674833411e0 FTN95 3.0576703006936028259906339255709e-1 FTN95 subtraction: z1-t2 = 2.2639435073548419286485538130979e0 FTN95 1.6487212707001281468486507878140e0 FTN95 subtraction: t1-z2 = 4.9140465044953092657608108666206e0 FTN95 -5.4030230586813971740093660744290e-1 FTN95 multiplication: z1*z2 = -3.6749719027089825421986358520178e0 FTN95 -2.0735774879057912796960818196054e0 FTN95 multiplication: z1*e2 = 1.3290026051044224203970081353460e0 FTN95 1.3949325106472885969003904553637e0 FTN95 multiplication: e1*z2 = -6.7959930977712717794720996889337e-1 FTN95 2.0716425082175620756303001513405e-1 FTN95 multiplication: z1*ec2 = -1.9611600161480140352958865685155e0 FTN95 -7.9265330380368870908232631553945e-1 FTN95 multiplication: ec1*z2 = -1.1367325229084841580870893177931e0 FTN95 -1.2924546017937383707088936698070e0 FTN95 multiplication: z1*t2 = -1.0887930451518010652503444491187e0 FTN95 -1.1428065003150041926917235553723e0 FTN95 multiplication: t1*z2 = -5.5683279968317078452848179821188e0 FTN95 1.6974097548329731696910825130042e0 FTN95 multiplication: z1*e2 = 1.3290026051044224203970081353460e0 FTN95 1.3949325106472885969003904553637e0 FTN95 multiplication: e1*z2 = -6.7959930977712717794720996889337e-1 FTN95 2.0716425082175620756303001513405e-1 FTN95 division: z1/z2 = -5.5143308202505235309319492678347e-1 FTN95 -1.0982863308151707791148976578178e0 FTN95 division: z1/e2 = 1.8565810862940114131789045300696e0 FTN95 1.9486834006024310747344030087682e0 FTN95 division: e1/z2 = -1.9793082933305223971873848218311e-1 FTN95 -6.0335835224373422321537161536479e-2 FTN95 division: z1/ec2 = -8.0761209709718786022689577180586e-1 FTN95 -2.3146781483996128766349483296804e0 FTN95 division: ec1/z2 = -6.4792442214258164548227795583516e-2 FTN95 -4.9709437648746959344238288583166e-1 FTN95 division: z1/t2 = -2.2661800709135969048138414728583e0 FTN95 -2.3786020010472250834113107924484e0 FTN95 division: t1/z2 = -1.6217552939728537940734027673102e0 FTN95 -4.9436442276888664087117814236064e-1 FTN95 exponentiation: z1**i1 = -3.7758122042342049433937797240400e1 FTN95 -4.8211113863064078910484120179802e1 FTN95 exponentiation: z1**z2 = 8.2348324476658808357864024273515e-2 FTN95 -1.2556301651409895918663118275737e-1 FTN95 exponentiation: t1**z2 = 1.0711458726885084804417651222829e-1 FTN95 7.6227798774539096243706577479256e-2 FTN95 ERROR: 1.498392D-30 exponentiation: z1**t2 = 4.7858855946366079834403314748766e-1 FTN95 -3.0083157713034823132511536931490e-1 FTN95 ERROR: 1.573438D-30 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.8529757783943701742247845600847e0 FTN95 aimag(z1) = 1.6487212707001281468486507878140e0 FTN95 conjg(z1) = 1.5707963267948966192313216916397e0 FTN95 -1.6487212707001281468486507878140e0 FTN95 cos(z2) = -2.3024722091021678526847749091911e-1 FTN95 5.5548753833677987141958048746338e-1 FTN95 ERROR: 1.470263D-30 exp(z1) = -3.7447692223345517631657155895988e-1 FTN95 4.7958794467223451852222400975742e0 FTN95 ERROR: 1.533165D-31 log(z1) = 8.2295077269804973154300843933861e-1 FTN95 8.0959735783473218711648994416608e-1 FTN95 ERROR: 1.899561D-31 ddcmplx(t1,t2) = 3.1415926535897932384626433832795e0 FTN95 -6.9314718055994530941723212145817e-1 FTN95 ddreal(z1) = 1.5707963267948966192313216916397e0 FTN95 sin(z2) = -1.1262543685294698363059631197625e0 FTN95 -1.1356178988170950050303278823902e-1 FTN95 ERROR: 4.842711D-31 sqrt(z1) = 1.3870843138895871220472192078626e0 FTN95 5.9431184326382898915869848949756e-1 FTN95 sqrt(z2) = 2.0065134872316974651058426007592e-1 FTN95 1.3463709795780444928274633384031e0 FTN95 Max relative error = 9.998354D-26 ONE OR MORE TESTS FAILED