"{ Encoding: utf8 }"

"{ Package: 'stx:goodies/regression' }"

"{ NameSpace: RegressionTests }"

TestCase subclass:#QDoubleTests
	instanceVariableNames:''
	classVariableNames:''
	poolDictionaries:''
	category:'tests-Regression-Numbers'
!

!QDoubleTests class methodsFor:'documentation'!

documentation
"
    documentation to be added.

    [author:]
        cg

    [instance variables:]

    [class variables:]

    [see also:]

"
! !

!QDoubleTests methodsFor:'tests'!

test99_misc
     self skip:'unfinished work in progress'.
     
     self assert:( 1 asQDouble / (3 factorial) - (QDouble invFact at:1) ) = 0.0. 
     self assert:( 1 asQDouble / (4 factorial) - (QDouble invFact at:2) ) = 0.0. 
     self assert:( 1 asQDouble / (5 factorial) - (QDouble invFact at:3) ) = 0.0. 
     self assert:( 1 asQDouble / (6 factorial) - (QDouble invFact at:4) ) = 0.0. 
     self assert:( 1 asQDouble / (7 factorial) - (QDouble invFact at:5) ) = 0.0. 
     self assert:( 1 asQDouble / (8 factorial) - (QDouble invFact at:6) ) = 0.0. 
     self assert:( 1 asQDouble / (9 factorial) - (QDouble invFact at:7) ) = 0.0. 
     self assert:( 1 asQDouble / (10 factorial) - (QDouble invFact at:8) ) = 0.0. 
     self assert:( 1 asQDouble / (11 factorial) - (QDouble invFact at:9) ) = 0.0. 
     self assert:( 1 asQDouble / (12 factorial) - (QDouble invFact at:10) ) = 0.0. 
     self assert:( 1 asQDouble / (13 factorial) - (QDouble invFact at:11) ) = 0.0. 
     self assert:( 1 asQDouble / (14 factorial) - (QDouble invFact at:12) ) = 0.0. 
     self assert:( 1 asQDouble / (15 factorial) - (QDouble invFact at:13) ) = 0.0. 
     self assert:( 1 asQDouble / (16 factorial) - (QDouble invFact at:14) ) = 0.0. 
     self assert:( 1 asQDouble / (17 factorial) - (QDouble invFact at:15) ) = 0.0.

    "Created: / 20-06-2017 / 13:06:04 / cg"
    "Modified: / 28-03-2019 / 12:20:51 / Claus Gittinger"
!

test_01_instance_creation
    |q|

    q := QDouble fromFloat:1.0.
    self assert:(q d0 = 1.0).
    self assert:(q d1 = 0.0).
    self assert:(q d2 = 0.0).
    self assert:(q d3 = 0.0).

    q := 1.0 asQDouble.
    self assert:(q d0 = 1.0).
    self assert:(q d1 = 0.0).
    self assert:(q d2 = 0.0).
    self assert:(q d3 = 0.0).

    q := 1e40 asQDouble + 1.0 asQDouble.
    self assert:(q d0 = 1e40).
    self assert:(q d1 = 1.0).
    self assert:(q d2 = 0.0).
    self assert:(q d3 = 0.0).

    q := 1e40 asQDouble + 1.0.
    self assert:(q d0 = 1e40).
    self assert:(q d1 = 1.0).
    self assert:(q d2 = 0.0).
    self assert:(q d3 = 0.0).

    q := 1e120 asQDouble + 1e40 + 1.0.
    self assert:(q d0 = 1e120).
    self assert:(q d1 = 1e40).
    self assert:(q d2 = 1.0).
    self assert:(q d3 = 0.0).

    q := 1e200 asQDouble + 1e120 + 1e40 + 1.0.
    self assert:(q d0 = 1e200).
    self assert:(q d1 = 1e120).
    self assert:(q d2 = 1e40).
    self assert:(q d3 = 1.0).

    "
     self run:#test_01_instance_creation
     self new test_01_instance_creation
    "

    "Created: / 12-06-2017 / 17:03:46 / cg"
    "Modified: / 20-06-2017 / 07:48:51 / cg"
    "Modified: / 05-06-2019 / 20:15:30 / Claus Gittinger"
!

test_02_addition
    |d1 d2 l1 l2 sum_l q1 q2 sum_d sum_q|

    "/ self skip:'unfinished work in progress'.

    d1 := 1.0.
    d2 := 1.0e-16.
    sum_d := d1 + d2.
    "/ precision is lost with doubles!!
    self assert:( sum_d - 1.0 = 0.0 ).
    
    l1 := 1.0 asLongFloat.
    l2 := 1.0e-16 asLongFloat.
    sum_l := l1 + l2.
    "/ no precision is lost with long doubles!!
    self assert:( (sum_l - (1.0 asLongFloat)) > 0.0 ).

    l1 := 1.0 asLongFloat.
    l2 := 1.0e-30 asLongFloat.
    sum_l := l1 + l2.
    "/ precision is lost!!
    self assert:( (sum_l - (1.0 asLongFloat)) = 0.0 ).

    q1 := QDouble fromFloat:1.0.
    q2 := QDouble fromFloat:1.0e-16.
    sum_q := q1 + q2.
    
    "/ no precision is lost
    self assert:( (sum_q - (QDouble fromFloat:1.0)) > 0.0 ).

    q1 := QDouble fromFloat:1.0.
    q2 := QDouble fromFloat:2.0.
    sum_q := q1 + q2.
    self assert:( sum_q = 3.0 ).
    
    q1 := QDouble fromFloat:1.0.
    q2 := QDouble fromFloat:2.0.
    sum_q := q1 + 2.0.
    self assert:( sum_q = 3.0 ).

    q1 := QDouble fromFloat:1.0.
    q2 := QDouble fromFloat:2.0.
    sum_q := 1.0 + q2.
    self assert:( sum_q = 3.0 ).

    q1 := 1e200 asQDouble + 1e120 + 1e40 + 1.0.
    q2 := 1e120.
    sum_q := q1 + q2.
    self assert:(sum_q d0 = 1e200).
    self assert:(sum_q d1 = 2e120).
    self assert:(sum_q d2 = 1e40).
    self assert:(sum_q d3 = 1.0).

    "
     self run:#test_02_addition
     self new test_02_addition
    "

    "Created: / 12-06-2017 / 17:05:07 / cg"
    "Modified: / 20-06-2017 / 07:49:32 / cg"
    "Modified: / 05-06-2019 / 20:15:38 / Claus Gittinger"
!

test_03_subtract
    |d1 d2 l1 l2 diff_l q1 q2 diff_d diff_q|

    "/ self skip:'unfinished work in progress'.

    d1 := 1.0.
    d2 := 0.99999.
    diff_d := d1 - d2.
    "/ precision is lost with doubles!!
    
    l1 := 1.0 asLongFloat.
    l2 := 0.99999 asLongFloat.
    diff_l := l1 - l2.
    "/ no precision is lost with long doubles!!

    q1 := QDouble fromFloat:1.0.
    q2 := QDouble fromFloat:0.99999.

    diff_q := q1 - q2.
    
    "
     self run:#test_03_subtract
     self new test_03_subtract
    "

    "Created: / 12-06-2017 / 23:12:00 / cg"
    "Modified: / 05-06-2019 / 20:16:12 / Claus Gittinger"
!

test_04_relops
    "/ self skip:'unfinished work in progress'.

    self assert:(1.0 = 1.0 asQDouble).
    self assert:(1.0 asQDouble = 1.0).
    self assert:(1.0 asQDouble = 1.0 asQDouble).

    self assert:(1 = 1 asQDouble).
    self assert:(1 asQDouble = 1).
    self assert:(1 asQDouble = 1 asQDouble).

    self assert:(1 = 1.0 asQDouble).
    self assert:(1.0 asQDouble = 1).
    self assert:(1.0 asQDouble = 1 asQDouble).
    
    self assert:(1.0 = 1 asQDouble).
    self assert:(1 asQDouble = 1.0).
    self assert:(1 asQDouble = 1.0 asQDouble).

    
    self deny:(1.0 = 1.1 asQDouble).
    self deny:(1.0 asQDouble = 1.1).
    self deny:(1.0 asQDouble = 1.1 asQDouble).

    self deny:(1 = 2 asQDouble).
    self deny:(1 asQDouble = 2).
    self deny:(1 asQDouble = 2 asQDouble).

    self deny:(1 = 2.0 asQDouble).
    self deny:(1.0 asQDouble = 2).
    self deny:(1.0 asQDouble = 2 asQDouble).

    self deny:(2.0 = 1 asQDouble).
    self deny:(2 asQDouble = 1.0).
    self deny:(2 asQDouble = 1.0 asQDouble).


    
    self assert:(1.0 < 1.1 asQDouble).
    self assert:(1.0 asQDouble < 1.1).
    self assert:(1.0 asQDouble < 1.1 asQDouble).

    self assert:(1 < 2 asQDouble).
    self assert:(1 asQDouble < 2).
    self assert:(1 asQDouble < 2 asQDouble).

    self assert:(1 < 2.0 asQDouble).
    self assert:(1.0 asQDouble < 2).
    self assert:(1.0 asQDouble < 2 asQDouble).

    self assert:(1.0 < 2 asQDouble).
    self assert:(1 asQDouble < 2.0).
    self assert:(1 asQDouble < 2.0 asQDouble).

    
    self deny:(1.0 < 1.0 asQDouble).
    self deny:(1.0 asQDouble < 1.0).
    self deny:(1.0 asQDouble < 1.0 asQDouble).

    self deny:(1 < 1 asQDouble).
    self deny:(1 asQDouble < 1).
    self deny:(1 asQDouble < 1 asQDouble).

    self deny:(1 < 1.0 asQDouble).
    self deny:(1.0 asQDouble < 1).
    self deny:(1.0 asQDouble < 1 asQDouble).

    self deny:(1.0 < 1 asQDouble).
    self deny:(1 asQDouble < 1.0).
    self deny:(1 asQDouble < 1.0 asQDouble).



    self assert:(2.0 > 1.1 asQDouble).
    self assert:(2.0 asQDouble > 1.1).
    self assert:(2.0 asQDouble > 1.1 asQDouble).

    self assert:(3 > 2 asQDouble).
    self assert:(3 asQDouble > 2).
    self assert:(3 asQDouble > 2 asQDouble).

    self assert:(3 > 2.0 asQDouble).
    self assert:(3.0 asQDouble > 2).
    self assert:(3.0 asQDouble > 2 asQDouble).

    self assert:(3.0 > 2 asQDouble).
    self assert:(3 asQDouble > 2.0).
    self assert:(3 asQDouble > 2.0 asQDouble).


    self deny:(1.0 > 3.0 asQDouble).
    self deny:(1.0 asQDouble > 3.0).
    self deny:(1.0 asQDouble > 3.0 asQDouble).

    self deny:(1 > 3 asQDouble).
    self deny:(1 asQDouble > 3).
    self deny:(1 asQDouble > 3 asQDouble).

    self deny:(1 > 3.0 asQDouble).
    self deny:(1.0 asQDouble > 3).
    self deny:(1.0 asQDouble > 3 asQDouble).

    self deny:(1.0 > 3 asQDouble).
    self deny:(1 asQDouble > 3.0).
    self deny:(1 asQDouble > 3.0 asQDouble).


    self assert:(1.0 <= 3.0 asQDouble).
    self assert:(1.0 asQDouble <= 3.0).
    self assert:(1.0 asQDouble <= 3.0 asQDouble).

    self assert:(1.0 <= 1.0 asQDouble).
    self assert:(1.0 asQDouble <= 1.0).
    self assert:(1.0 asQDouble <= 1.0 asQDouble).
    
    self assert:(2.0 >= 1.0 asQDouble).
    self assert:(2.0 asQDouble >= 1.0).
    self assert:(2.0 asQDouble >= 1.0 asQDouble).

    self assert:(1.0 >= 1.0 asQDouble).
    self assert:(1.0 asQDouble >= 1.0).
    self assert:(1.0 asQDouble >= 1.0 asQDouble).


    
    self assert: ((1.0 asQDouble + 1e-40) > 1.0 asQDouble).
    self assert: (1.0 asQDouble < (1.0 asQDouble + 1e-40)).
    self assert: (1.0 asQDouble < (1.0 asQDouble + 1e-80)).
    self assert: (1.0 asQDouble < (1.0 asQDouble + 1e-120)).
    self assert: (1.0 asQDouble < (1.0 asQDouble + 1e-300)).
    self assert: (1.0 asQDouble < (1.0 asQDouble + 1e-300)).
    
    self assert: ((1.0 asQDouble + 1e-40) > 1.0).
    self assert: (1.0 < (1.0 asQDouble + 1e-40)).
    self assert: (1.0 < (1.0 asQDouble + 1e-80)).
    self assert: (1.0 < (1.0 asQDouble + 1e-120)).
    self assert: (1.0 < (1.0 asQDouble + 1e-300)).
    self assert: (1.0 < (1.0 asQDouble + 1e-300)).

    self assert: (1e10 asQDouble + 1.0) = 10000000001.0.
    "/ that's the problem with Floats...
    "/ (1e20 asQDouble + 1.0) ~= 100000000000000000001.0. 
    "/ (1e20 + 1.0) - 1e20 ~= 1.0. 
    self assert: (1e20 asQDouble + 1.0) - 1e20 = 1.0. 

    "/ here the precision of QDOubles is no longer enough:
    "/ self assert:  1.0 asQDouble < (1.0 asQDouble + 1e-310).

    "
     self run:#test_03_subtract
     self new test_03_subtract
    "

    "Created: / 13-06-2017 / 17:11:37 / cg"
    "Modified (comment): / 20-06-2017 / 09:18:50 / cg"
    "Modified: / 05-06-2019 / 20:16:21 / Claus Gittinger"
!

test_05_multiply
    |q1 q2 p|

    "/ self skip:'unfinished work in progress'.

    q1 := QDouble fromFloat:3.0. "/ DoubleArray(2.0 0.0 0.0 0.0)
    p := 2.0 * q1.

    self assert:(p d0 = 6.0).
    self assert:(p d1 = 0.0).
    self assert:(p d2 = 0.0).
    self assert:(p d3 = 0.0).
    self assert:(p = 6.0).
    self assert:(6.0 = p).

    q1 := QDouble fromFloat:2.0. "/ DoubleArray(2.0 0.0 0.0 0.0)
    p := q1 * 3.0.

    self assert:(p d0 = 6.0).
    self assert:(p d1 = 0.0).
    self assert:(p d2 = 0.0).
    self assert:(p d3 = 0.0).
    self assert:(p = 6.0).
    self assert:(6.0 = p).

    q1 := QDouble fromFloat:2.0. "/ DoubleArray(2.0 0.0 0.0 0.0)
    q2 := QDouble fromFloat:3.0. "/ DoubleArray(3.0 0.0 0.0 0.0)
    p := q1 * q2.
    
    self assert:(p d0 = 6.0).
    self assert:(p d1 = 0.0).
    self assert:(p d2 = 0.0).
    self assert:(p d3 = 0.0).
    self assert:(p = 6.0).
    self assert:(6.0 = p).

    "/ the 1 is not forgotten!!
    q1 := 1 asQDouble + 1e40.
    q2 := 2.0.
    p := q1 * q2.

    self assert:(p d0 = 2e40).
    self assert:(p d1 = 2.0).
    self assert:(p d2 = 0.0).
    self assert:(p d3 = 0.0).
    self assert:(p > 2e40).
    self assert:(p - 2e40) = 2.0.

    "/ nothing is forgotten!!
    q1 := 1 asQDouble + 1e40 + 1e120.
    q2 := 2.0.
    p := q1 * q2.

    self assert:(p d0 = 2e120).
    self assert:(p d1 = 2e40).
    self assert:(p d2 = 2.0).
    self assert:(p d3 = 0.0).
    self assert:(p > 2e120).
    self assert:(p > (2e120 asQDouble + 2)).
    self assert:(p > (2e120 asQDouble + 1e40)).
    self assert:(p - 2e120) = (2 asQDouble + 2e40).
    self assert:(p - 2e120 - 2e40) = 2.0.

    "
     self run:#test_05_multiply
     self new test_05_multiply
    "

    "Created: / 19-06-2017 / 16:58:29 / cg"
    "Modified (comment): / 20-06-2017 / 07:42:52 / cg"
    "Modified: / 05-06-2019 / 20:16:30 / Claus Gittinger"
!

test_06_exp
    |f1 q1|

    "/ self skip:'unfinished work in progress'.

    f1 := 2.0 exp.
    q1 := 2.0 asQDouble exp.

    "/ the first 113 digits are:
    "/                                            7.3890560989306502272304274605750078131803155705518473240871278225225737960790577633843124850791217947737531612654...
    "/ 2.0 exp             
    "/ '%30.28f' printf:{ 2.0 exp }           -> '7.38905609893065040694182243896'
    "/ '%30.28f' printf:{ 2.0 exp asQDouble } -> '7.38905609893065040694182243897'
    "/ '%30.28f' printf:{ 2.0 asQDouble exp } -> '7.38905609893065022723042746058'
    "/ '%40.38f' printf:{ 2.0 asQDouble exp } -> '7.389056098930650227230427460575007822092' (10 iters in exp)
    "/ '%40.38f' printf:{ 2.0 asQDouble exp }    
    "/ 2.0 asQDouble exp   -> 7.3890560989306502272304274605750078* 220914983921672'
    "/ 1.0 asQDouble exp       
    self assert:(q1 d0 = f1).

    "
     self run:#test_06_exp
     self new test_06_exp
    "

    "Created: / 19-06-2017 / 17:58:08 / cg"
    "Modified (comment): / 21-06-2017 / 13:45:56 / cg"
    "Modified: / 05-06-2019 / 20:16:35 / Claus Gittinger"
!

test_06b_log10
    |q1 q2 d3a q3a d3b q3b q4 i|

    "/ self skip:'unfinished work in progress'.

    q1 := QDouble d0: 5.55111512312578E-18 
                  d1: -3.08148791101958E-34 
                  d2: 1.71056941445901E-50 
                  d3: 0.0.

    q2 := 10 / q1.
    self assert:(q2 > 1e18).

    d3a := q2 asFloat ln.
    self assert:(d3a isNaN not).
    
    q3a := q2 ln.
    self assert:(q3a isNaN not).
    
    d3b := q2 asFloat log10.
    self assert:(d3b isNaN not).

    q3b := q2 log10.
    self assert:(q3b isNaN not).

    q4 := q3b floor.
    i := q4 asInteger.


    "
     self run:#test_06b_log10
     self new test_06b_log10
    "

    "Created: / 03-07-2017 / 10:31:53 / cg"
    "Modified: / 05-06-2019 / 20:18:26 / Claus Gittinger"
!

test_07_reading
    |q1 s|

    "/ self skip:'unfinished work in progress'.

    q1 := QDouble readFrom: '3.1415926535897'.
    s := '%55.54f' printf:{ q1 }.
    self assert:(s = '                                        3.1415926535897').

    q1 := QDouble readFrom: '3.141592653589793238462643383'.
    s := '%55.54f' printf:{ q1 }.
    self assert:(s = '                          3.141592653589793238462643383').

    q1 := QDouble readFrom: '3.141592653589793238462643383279502884197169399375105820'.
    s := '%55.54f' printf:{ q1 }.
    self assert:(s = '3.14159265358979323846264338327950288419716939937510582').

    "
     self run:#test_07_reading
     self new test_07_reading
    "

    "Created: / 19-06-2017 / 23:29:34 / cg"
    "Modified: / 05-06-2019 / 20:18:51 / Claus Gittinger"
!

test_08_conversion
    |q i|
    
    "/ self skip:'unfinished work in progress'.

    q := 1e10 asQDouble + 1.0.
    i := q asInteger.

    self assert:(i = 10000000001).

    "
     self run:#test_08_conversion
     self new test_08_conversion
    "

    "Created: / 20-06-2017 / 07:51:50 / cg"
    "Modified: / 05-06-2019 / 20:18:48 / Claus Gittinger"
! !

!QDoubleTests class methodsFor:'documentation'!

version
    ^ '$Header$'
!

version_CVS
    ^ '$Header$'
! !