Martin Rubey wrote:
>
> Dear all,
>
> whoever has some cpu cycles to spare
>
> PLEASE
>
> build the new revision of FriCAS 448, and report the output of elemnum.input,
> together with lisp variant (CLISP anybody?) and hardware architecture.
>
64-bit Fedora 9 on Intel Core-2
Using clisp-2.47, given
log(0.0::Complex DFLOAT)
FriCAS apparently goes into infinite loop. Up to this point the
results are essentially the same as sbcl.
Closure-CL 1.2 has essentialy the same results (modulo some extra
messages and lowerest digits) as sbcl, except for two atan tests:
testLibraryError("atan(%i::Complex DFLOAT)")
testLibraryError("atan(-%i::Complex DFLOAT)")
which using Closure CL give:
+Daly Bug
+ >> System error:
+ DIVISION-BY-ZERO detected
+performing LOG on (0.0)
+
sbcl 1.22 results in attachement.
BTW: most failures seem to be since we do not simplify special
values of inverse functions. I am not sure if we really want
such simplifications (at least as default): they may be expensive.
--
Waldek Hebisch
[EMAIL PROTECTED]
--~--~---------~--~----~------------~-------~--~----~
You received this message because you are subscribed to the Google Groups
"FriCAS - computer algebra system" group.
To post to this group, send email to [email protected]
To unsubscribe from this group, send email to [EMAIL PROTECTED]
For more options, visit this group at
http://groups.google.com/group/fricas-devel?hl=en
-~----------~----~----~----~------~----~------~--~---
=============================================================================
Testsuite: inverses
failed (total): 0 (6)
=============================================================================
Testsuite: branch cuts
failed (total): 4 (5)
failed testcases were:
---------------------------------------------------------------------------
Testcase: log
failed (total): 1 (198)
failed tests were:
100: EQLC: (valOld, valNew)
Output1: 0.6931471805599453 - 3.141592653589793 %i
Output2: 0.6931471805599453 + 3.0787608005179976 %i
Output3: abs err: 38.69279709403072
Output4: rel err: 3.885122820389911
---------------------------------------------------------------------------
Testcase: acos
failed (total): 4 (396)
failed tests were:
1: EQLC: (valOld, valNew)
Output1: - 1.3169578969 248167086 %i
Output2: 3.1053226275 946066269 - 1.3171476554 947782944 %i
Output3: abs err: 9.6430286574 593868287
Output4: rel err: 0.8475222221 2187734941
100: EQLC: (valOld, valNew)
Output1: 3.141592653589793 + 1.3169578969248166 %i
Output2: 3.105322627594607 - 1.317147655494778 %i
Output3: abs err: 6.939827576073428
Output4: rel err: 0.6099388788880831
199: EQLC: (valOld, valNew)
Output1: 3.1415926535 897932385 + 1.3169578969 248167086 %i
Output2: 0.0362700259 951866115 + 1.3171476554 947782944 %i
Output3: abs err: 9.6430286574 593868289
Output4: rel err: 5.5541210545 802845491
298: EQLC: (valOld, valNew)
Output1: - 1.3169578969248166 %i
Output2: 0.03627002599518661 + 1.317147655494778 %i
Output3: abs err: 6.939827576073428
Output4: rel err: 3.9971510844375646
---------------------------------------------------------------------------
Testcase: asin
failed (total): 4 (396)
failed tests were:
1: EQLC: (valOld, valNew)
Output1: 1.5707963267 948966192 + 1.3169578969 248167086 %i
Output2: - 1.5345263007 997100077 + 1.3171476554 947782944 %i
Output3: abs err: 9.6430286574 593868287
Output4: rel err: 2.3579111213 988348043
100: EQLC: (valOld, valNew)
Output1: - 1.5707963267948966 - 1.3169578969248166 %i
Output2: - 1.5345263007997103 + 1.317147655494778 %i
Output3: abs err: 6.939827576073428
Output4: rel err: 1.6969250225711858
199: EQLC: (valOld, valNew)
Output1: - 1.5707963267 948966192 - 1.3169578969 248167086 %i
Output2: 1.5345263007 997100077 - 1.3171476554 947782944 %i
Output3: abs err: 9.6430286574 593868289
Output4: rel err: 2.3579111213 988348043
298: EQLC: (valOld, valNew)
Output1: 1.5707963267948966 + 1.3169578969248166 %i
Output2: 1.53452630079971 - 1.317147655494778 %i
Output3: abs err: 6.939827576073428
Output4: rel err: 1.6969250225711854
---------------------------------------------------------------------------
Testcase: atan
failed (total): 4 (396)
failed tests were:
1: EQLC: (valOld, valNew)
Output1: - 1.5707963267 948966192 - 0.5493061443 340548457 %i
Output2: 1.5498680512 236059035 - 0.5487583583 8648399359 %i
Output3: abs err: 9.7385464603 03051567
Output4: rel err: 3.6025637127 298097619
100: EQLC: (valOld, valNew)
Output1: - 1.5707963267948966 - 0.5493061443340549 %i
Output2: 1.5498680512236058 - 0.5487583583864839 %i
Output3: abs err: 9.738546460303052
Output4: rel err: 3.60256371272981
199: EQLC: (valOld, valNew)
Output1: 1.5707963267 948966192 + 0.5493061443 340548457 %i
Output2: - 1.5498680512 236059035 + 0.5487583583 8648399359 %i
Output3: abs err: 9.7385464603 030515671
Output4: rel err: 3.6025637127 298097619
298: EQLC: (valOld, valNew)
Output1: 1.5707963267948966 + 0.5493061443340549 %i
Output2: - 1.5498680512236058 + 0.548758358386484 %i
Output3: abs err: 9.738546460303052
Output4: rel err: 3.6025637127298094
=============================================================================
Testsuite: singularities
failed (total): 1 (2)
failed testcases were:
---------------------------------------------------------------------------
Testcase: atan
failed (total): 4 (6)
failed tests were:
3: ERROR: (atan(%i::Complex DFLOAT))
Output1: 0.7853981633974483 + 177.09910463306602 %i
4: ERROR: (atan(-%i::Complex DFLOAT))
Output1: - 0.7853981633974483 - 177.09910463306602 %i
5: ERROR: (atan(%i))
Output1: atan(%i)
6: ERROR: (atan(-%i))
Output1: - atan(%i)
=============================================================================
Testsuite: specific values
failed (total): 4 (4)
failed testcases were:
---------------------------------------------------------------------------
Testcase: log
failed (total): 12 (16)
failed tests were:
2: EQUL: (log(-%e^0), 0+ %i*%pi)
Output1: log(- 1)
Output2: %i %pi
3: EQUL: (log(%i*%e^0), 0+ %i*%pi/2)
Output1: log(%i)
%i %pi
Output2: ------
2
4: EQUL: (log(-%i*%e^0), 0- %i*%pi/2)
Output1: log(- %i)
%i %pi
Output2: - ------
2
6: EQUL: (log(-%e^1), 1+ %i*%pi)
Output1: log(- %e)
Output2: %i %pi + 1
7: EQUL: (log(%i*%e^1), 1+ %i*%pi/2)
Output1: log(%i %e)
%i %pi + 2
Output2: ----------
2
8: EQUL: (log(-%i*%e^1), 1- %i*%pi/2)
Output1: log(- %i %e)
- %i %pi + 2
Output2: ------------
2
10: EQUL: (log(-%e^2), 2+ %i*%pi)
2
Output1: log(- %e )
Output2: %i %pi + 2
11: EQUL: (log(%i*%e^2), 2+ %i*%pi/2)
2
Output1: log(%i %e )
%i %pi + 4
Output2: ----------
2
12: EQUL: (log(-%i*%e^2), 2- %i*%pi/2)
2
Output1: log(- %i %e )
- %i %pi + 4
Output2: ------------
2
14: EQUL: (log(-%e^3), 3+ %i*%pi)
3
Output1: log(- %e )
Output2: %i %pi + 3
15: EQUL: (log(%i*%e^3), 3+ %i*%pi/2)
3
Output1: log(%i %e )
%i %pi + 6
Output2: ----------
2
16: EQUL: (log(-%i*%e^3), 3- %i*%pi/2)
3
Output1: log(- %i %e )
- %i %pi + 6
Output2: ------------
2
---------------------------------------------------------------------------
Testcase: acos
failed (total): 11 (24)
failed tests were:
3: EQUL: (acos(1/sqrt(1+(2-sqrt 3)^2)), %pi/12)
+-----------+
+-+ | +-+
(\|3 + 2)\|- 4\|3 + 8
Output1: acos(------------------------)
4
%pi
Output2: ---
12
5: EQUL: (acos(1/sqrt(1+(sqrt(1-2/sqrt 5))^2)), %pi/10)
+------------+
+-+ | +-+
(\|5 + 1)\|- 2\|5 + 10
Output1: acos(-------------------------)
8
%pi
Output2: ---
10
7: EQUL: (acos(1/sqrt(1+(sqrt 2 - 1)^2)), %pi/8)
+-----------+
+-+ | +-+
(\|2 + 2)\|- 2\|2 + 4
Output1: acos(------------------------)
4
%pi
Output2: ---
8
9: EQUL: (acos(1/sqrt(1+(1/sqrt 3)^2)), %pi/6)
+-+
\|3
Output1: acos(----)
2
%pi
Output2: ---
6
11: EQUL: (acos(1/sqrt(1+(sqrt(5-2*sqrt 5))^2)), %pi/5)
+-----------+
+-+ | +-+
(\|5 + 3)\|- 2\|5 + 6
Output1: acos(------------------------)
8
%pi
Output2: ---
5
13: EQUL: (acos(1/sqrt(1+(1)^2)), %pi/4)
+-+
\|2
Output1: acos(----)
2
%pi
Output2: ---
4
15: EQUL: (acos(1/sqrt(1+(sqrt(1+2/sqrt 5))^2)), 3*%pi/10)
+----------+
+-+ | +-+
(\|5 - 1)\|2\|5 + 10
Output1: acos(-----------------------)
8
3%pi
Output2: ----
10
17: EQUL: (acos(1/sqrt(1+(sqrt 3)^2)), %pi/3)
1
Output1: acos(-)
2
%pi
Output2: ---
3
19: EQUL: (acos(1/sqrt(1+(1+sqrt 2)^2)), 3*%pi/8)
+---------+
+-+ | +-+
(- \|2 + 2)\|2\|2 + 4
Output1: acos(------------------------)
4
3%pi
Output2: ----
8
21: EQUL: (acos(1/sqrt(1+(sqrt(5+2*sqrt 5))^2)), 2*%pi/5)
+---------+
+-+ | +-+
(- \|5 + 3)\|2\|5 + 6
Output1: acos(------------------------)
8
2%pi
Output2: ----
5
23: EQUL: (acos(1/sqrt(1+(2+sqrt 3)^2)), 5*%pi/12)
+---------+
+-+ | +-+
(- \|3 + 2)\|4\|3 + 8
Output1: acos(------------------------)
4
5%pi
Output2: ----
12
---------------------------------------------------------------------------
Testcase: asin
failed (total): 11 (24)
failed tests were:
3: EQUL: (asin((2-sqrt 3)/sqrt(1+(2-sqrt 3)^2)), %pi/12)
+-----------+
| +-+
\|- 4\|3 + 8
Output1: asin(--------------)
4
%pi
Output2: ---
12
5: EQUL: (asin((sqrt(1-2/sqrt 5))/sqrt(1+(sqrt(1-2/sqrt 5))^2)), %pi/10
)
+-----------+ +------------+
+-+ | +-+ | +-+
(\|5 + 5)\|- 2\|5 + 5 \|- 2\|5 + 10
Output1: asin(---------------------------------------)
40
%pi
Output2: ---
10
7: EQUL: (asin((sqrt 2 - 1)/sqrt(1+(sqrt 2 - 1)^2)), %pi/8)
+-----------+
+-+ | +-+
\|2 \|- 2\|2 + 4
Output1: asin(------------------)
4
%pi
Output2: ---
8
9: EQUL: (asin((1/sqrt 3)/sqrt(1+(1/sqrt 3)^2)), %pi/6)
1
Output1: asin(-)
2
%pi
Output2: ---
6
11: EQUL: (asin((sqrt(5-2*sqrt 5))/sqrt(1+(sqrt(5-2*sqrt 5))^2)), %pi/5
)
+-----------+ +-----------+
+-+ | +-+ | +-+
(\|5 + 3)\|- 2\|5 + 5 \|- 2\|5 + 6
Output1: asin(--------------------------------------)
8
%pi
Output2: ---
5
13: EQUL: (asin((1)/sqrt(1+(1)^2)), %pi/4)
+-+
\|2
Output1: asin(----)
2
%pi
Output2: ---
4
15: EQUL: (asin((sqrt(1+2/sqrt 5))/sqrt(1+(sqrt(1+2/sqrt 5))^2)), 3*%pi
/10)
+---------+ +----------+
+-+ | +-+ | +-+
(\|5 - 5)\|2\|5 + 5 \|2\|5 + 10
Output1: - asin(-----------------------------------)
40
3%pi
Output2: ----
10
17: EQUL: (asin((sqrt 3)/sqrt(1+(sqrt 3)^2)), %pi/3)
+-+
\|3
Output1: asin(----)
2
%pi
Output2: ---
3
19: EQUL: (asin((1+sqrt 2)/sqrt(1+(1+sqrt 2)^2)), 3*%pi/8)
+---------+
+-+ | +-+
\|2 \|2\|2 + 4
Output1: asin(----------------)
4
3%pi
Output2: ----
8
21: EQUL: (asin((sqrt(5+2*sqrt 5))/sqrt(1+(sqrt(5+2*sqrt 5))^2)), 2*%pi
/5)
+---------+ +---------+
+-+ | +-+ | +-+
(\|5 - 3)\|2\|5 + 5 \|2\|5 + 6
Output1: - asin(----------------------------------)
8
2%pi
Output2: ----
5
23: EQUL: (asin((2+sqrt 3)/sqrt(1+(2+sqrt 3)^2)), 5*%pi/12)
+---------+
| +-+
\|4\|3 + 8
Output1: asin(------------)
4
5%pi
Output2: ----
12
---------------------------------------------------------------------------
Testcase: atan
failed (total): 16 (48)
failed tests were:
5: EQUL: (atan(2-sqrt 3), %pi/12)
+-+
Output1: - atan(\|3 - 2)
%pi
Output2: ---
12
6: EQUL: (atan(-(2-sqrt 3)), -(%pi/12))
+-+
Output1: atan(\|3 - 2)
%pi
Output2: - ---
12
9: EQUL: (atan(sqrt(1-2/sqrt 5)), %pi/10)
+-----------+
| +-+
\|- 2\|5 + 5
Output1: atan(--------------)
+-+
\|5
%pi
Output2: ---
10
10: EQUL: (atan(-(sqrt(1-2/sqrt 5))), -(%pi/10))
+-----------+
| +-+
\|- 2\|5 + 5
Output1: - atan(--------------)
+-+
\|5
%pi
Output2: - ---
10
13: EQUL: (atan(sqrt 2 - 1), %pi/8)
+-+
Output1: atan(\|2 - 1)
%pi
Output2: ---
8
14: EQUL: (atan(-(sqrt 2 - 1)), -(%pi/8))
+-+
Output1: - atan(\|2 - 1)
%pi
Output2: - ---
8
21: EQUL: (atan(sqrt(5-2*sqrt 5)), %pi/5)
+-----------+
| +-+
Output1: atan(\|- 2\|5 + 5 )
%pi
Output2: ---
5
22: EQUL: (atan(-(sqrt(5-2*sqrt 5))), -(%pi/5))
+-----------+
| +-+
Output1: - atan(\|- 2\|5 + 5 )
%pi
Output2: - ---
5
29: EQUL: (atan(sqrt(1+2/sqrt 5)), 3*%pi/10)
+---------+
| +-+
\|2\|5 + 5
Output1: atan(------------)
+-+
\|5
3%pi
Output2: ----
10
30: EQUL: (atan(-(sqrt(1+2/sqrt 5))), -(3*%pi/10))
+---------+
| +-+
\|2\|5 + 5
Output1: - atan(------------)
+-+
\|5
3%pi
Output2: - ----
10
37: EQUL: (atan(1+sqrt 2), 3*%pi/8)
+-+
Output1: atan(\|2 + 1)
3%pi
Output2: ----
8
38: EQUL: (atan(-(1+sqrt 2)), -(3*%pi/8))
+-+
Output1: - atan(\|2 + 1)
3%pi
Output2: - ----
8
41: EQUL: (atan(sqrt(5+2*sqrt 5)), 2*%pi/5)
+---------+
| +-+
Output1: atan(\|2\|5 + 5 )
2%pi
Output2: ----
5
42: EQUL: (atan(-(sqrt(5+2*sqrt 5))), -(2*%pi/5))
+---------+
| +-+
Output1: - atan(\|2\|5 + 5 )
2%pi
Output2: - ----
5
45: EQUL: (atan(2+sqrt 3), 5*%pi/12)
+-+
Output1: atan(\|3 + 2)
5%pi
Output2: ----
12
46: EQUL: (atan(-(2+sqrt 3)), -(5*%pi/12))
+-+
Output1: - atan(\|3 + 2)
5%pi
Output2: - ----
12
=============================================================================
Testsuite: simplification
failed (total): 1 (2)
failed testcases were:
---------------------------------------------------------------------------
Testcase: acos
failed (total): 10 (20)
failed tests were:
2: EQLC: (val, 0.0)
Output1: - 3.1415926535 897932385 - 0.1 E -19 %i
Output2: 0.0
Output3: abs err: 9.8696044010 893586189
3: EQLC: (val, 0.0)
Output1: - 3.1415926535 897932385 + 0.7 E -20 %i
Output2: 0.0
Output3: abs err: 9.8696044010 893586188
4: EQLC: (val, 0.0)
Output1: - 3.1415926535 897932385 - 0.7 E -20 %i
Output2: 0.0
Output3: abs err: 9.8696044010 893586189
5: EQLC: (val, 0.0)
Output1: - 3.1415926535 897932385
Output2: 0.0
Output3: abs err: 9.8696044010 893586189
6: EQLC: (val, 0.0)
Output1: - 3.1415926535 897932385
Output2: 0.0
Output3: abs err: 9.8696044010 893586188
11: EQLC: (val, 0.0)
Output1: - 3.1415926535 897932385 + 0.1 E -19 %i
Output2: 0.0
Output3: abs err: 9.8696044010 893586188
17: EQLC: (val, 0.0)
Output1: - 3.1415926535 897932385 + 0.2 E -19 %i
Output2: 0.0
Output3: abs err: 9.8696044010 893586188
18: EQLC: (val, 0.0)
Output1: - 3.1415926535 897932385 - 0.7 E -20 %i
Output2: 0.0
Output3: abs err: 9.8696044010 893586188
19: EQLC: (val, 0.0)
Output1: - 3.1415926535 897932385 - 0.7 E -20 %i
Output2: 0.0
Output3: abs err: 9.8696044010 893586188
20: EQLC: (val, 0.0)
Output1: - 3.1415926535 897932385 + 0.1 E -19 %i
Output2: 0.0
Output3: abs err: 9.8696044010 893586189
=============================================================================
testsuite | testcases: failed (total) | tests: failed (total)
inverses 0 (6) 0 (312)
branch cuts 4 (5) 13 (1584)
singularities 1 (2) 4 (10)
specific values 4 (4) 50 (112)
simplification 1 (2) 10 (40)
Type: Void
