File: math_domain_special.py

package info (click to toggle)
giac 1.9.0.93%2Bdfsg2-3
  • links: PTS, VCS
  • area: main
  • in suites: forky, sid, trixie
  • size: 117,732 kB
  • sloc: cpp: 404,272; ansic: 205,462; python: 30,548; javascript: 28,788; makefile: 17,997; yacc: 2,690; lex: 2,464; sh: 705; perl: 314; lisp: 216; asm: 62; java: 41; xml: 36; sed: 16; csh: 7; pascal: 6
file content (36 lines) | stat: -rw-r--r-- 938 bytes parent folder | download | duplicates (3)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
# Tests domain errors in special math functions

try:
    import math
    math.erf
except (ImportError, AttributeError):
    print("SKIP")
    raise SystemExit

inf = float('inf')
nan = float('nan')

# single argument functions
for name, f, args in (
    ('expm1', math.exp, ()),
    ('log2', math.log2, (-1, 0)),
    ('log10', math.log10, (-1, 0)),
    ('sinh', math.sinh, ()),
    ('cosh', math.cosh, ()),
    ('tanh', math.tanh, ()),
    ('asinh', math.asinh, ()),
    ('acosh', math.acosh, (-1, 0.9, 1)),
    ('atanh', math.atanh, (-1, 1)),
    ('erf', math.erf, ()),
    ('erfc', math.erfc, ()),
    ('gamma', math.gamma, (-2, -1, 0, 1)),
    ('lgamma', math.lgamma, (-2, -1, 0, 1)),
    ):
    for x in args + (inf, -inf, nan):
        try:
            ans = f(x)
            print('%.4f' % ans)
        except ValueError:
            print(name, 'ValueError')
        except OverflowError:
            print(name, 'OverflowError')