2019-07-10 21:40:29 +02:00
|
|
|
import math
|
|
|
|
import cmath
|
2019-07-06 20:19:09 +02:00
|
|
|
|
|
|
|
// Tests are based on and verified from practice examples of Khan Academy
|
|
|
|
// https://www.khanacademy.org/math/precalculus/imaginary-and-complex-numbers
|
|
|
|
|
|
|
|
fn test_complex_addition() {
|
2019-07-10 21:40:29 +02:00
|
|
|
mut c1 := cmath.complex(0,-10)
|
|
|
|
mut c2 := cmath.complex(-40,8)
|
2019-07-06 20:19:09 +02:00
|
|
|
mut result := c1 + c2
|
2019-07-10 21:40:29 +02:00
|
|
|
assert result.equals(cmath.complex(-40,-2))
|
|
|
|
c1 = cmath.complex(-71,2)
|
|
|
|
c2 = cmath.complex(88,-12)
|
2019-07-06 20:19:09 +02:00
|
|
|
result = c1 + c2
|
2019-07-10 21:40:29 +02:00
|
|
|
assert result.equals(cmath.complex(17,-10))
|
|
|
|
c1 = cmath.complex(0,-30)
|
|
|
|
c2 = cmath.complex(52,-30)
|
2019-07-06 20:19:09 +02:00
|
|
|
result = c1 + c2
|
2019-07-10 21:40:29 +02:00
|
|
|
assert result.equals(cmath.complex(52,-60))
|
|
|
|
c1 = cmath.complex(12,-9)
|
|
|
|
c2 = cmath.complex(32,-6)
|
2019-07-06 20:19:09 +02:00
|
|
|
result = c1 + c2
|
2019-07-10 21:40:29 +02:00
|
|
|
assert result.equals(cmath.complex(44,-15))
|
2019-07-06 20:19:09 +02:00
|
|
|
}
|
|
|
|
|
|
|
|
fn test_complex_subtraction() {
|
2019-07-10 21:40:29 +02:00
|
|
|
mut c1 := cmath.complex(-8,0)
|
|
|
|
mut c2 := cmath.complex(6,30)
|
2019-07-06 20:19:09 +02:00
|
|
|
mut result := c1 - c2
|
2019-07-10 21:40:29 +02:00
|
|
|
assert result.equals(cmath.complex(-14,-30))
|
|
|
|
c1 = cmath.complex(-19,7)
|
|
|
|
c2 = cmath.complex(29,32)
|
2019-07-06 20:19:09 +02:00
|
|
|
result = c1 - c2
|
2019-07-10 21:40:29 +02:00
|
|
|
assert result.equals(cmath.complex(-48,-25))
|
|
|
|
c1 = cmath.complex(12,0)
|
|
|
|
c2 = cmath.complex(23,13)
|
2019-07-06 20:19:09 +02:00
|
|
|
result = c1 - c2
|
2019-07-10 21:40:29 +02:00
|
|
|
assert result.equals(cmath.complex(-11,-13))
|
|
|
|
c1 = cmath.complex(-14,3)
|
|
|
|
c2 = cmath.complex(0,14)
|
2019-07-06 20:19:09 +02:00
|
|
|
result = c1 - c2
|
2019-07-10 21:40:29 +02:00
|
|
|
assert result.equals(cmath.complex(-14,-11))
|
2019-07-06 20:19:09 +02:00
|
|
|
}
|
|
|
|
|
|
|
|
fn test_complex_multiplication() {
|
2019-07-10 21:40:29 +02:00
|
|
|
mut c1 := cmath.complex(1,2)
|
|
|
|
mut c2 := cmath.complex(1,-4)
|
2019-07-06 20:19:09 +02:00
|
|
|
mut result := c1.multiply(c2)
|
2019-07-10 21:40:29 +02:00
|
|
|
assert result.equals(cmath.complex(9,-2))
|
|
|
|
c1 = cmath.complex(-4,-4)
|
|
|
|
c2 = cmath.complex(-5,-3)
|
2019-07-06 20:19:09 +02:00
|
|
|
result = c1.multiply(c2)
|
2019-07-10 21:40:29 +02:00
|
|
|
assert result.equals(cmath.complex(8,32))
|
|
|
|
c1 = cmath.complex(4,4)
|
|
|
|
c2 = cmath.complex(-2,-5)
|
2019-07-06 20:19:09 +02:00
|
|
|
result = c1.multiply(c2)
|
2019-07-10 21:40:29 +02:00
|
|
|
assert result.equals(cmath.complex(12,-28))
|
|
|
|
c1 = cmath.complex(2,-2)
|
|
|
|
c2 = cmath.complex(4,-4)
|
2019-07-06 20:19:09 +02:00
|
|
|
result = c1.multiply(c2)
|
2019-07-10 21:40:29 +02:00
|
|
|
assert result.equals(cmath.complex(0,-16))
|
2019-07-06 20:19:09 +02:00
|
|
|
}
|
|
|
|
|
|
|
|
fn test_complex_division() {
|
2019-07-10 21:40:29 +02:00
|
|
|
mut c1 := cmath.complex(-9,-6)
|
|
|
|
mut c2 := cmath.complex(-3,-2)
|
2019-07-06 20:19:09 +02:00
|
|
|
mut result := c1.divide(c2)
|
2019-07-10 21:40:29 +02:00
|
|
|
assert result.equals(cmath.complex(3,0))
|
|
|
|
c1 = cmath.complex(-23,11)
|
|
|
|
c2 = cmath.complex(5,1)
|
2019-07-06 20:19:09 +02:00
|
|
|
result = c1.divide(c2)
|
2019-07-10 21:40:29 +02:00
|
|
|
assert result.equals(cmath.complex(-4,3))
|
|
|
|
c1 = cmath.complex(8,-2)
|
|
|
|
c2 = cmath.complex(-4,1)
|
2019-07-06 20:19:09 +02:00
|
|
|
result = c1.divide(c2)
|
2019-07-10 21:40:29 +02:00
|
|
|
assert result.equals(cmath.complex(-2,0))
|
|
|
|
c1 = cmath.complex(11,24)
|
|
|
|
c2 = cmath.complex(-4,-1)
|
2019-07-06 20:19:09 +02:00
|
|
|
result = c1.divide(c2)
|
2019-07-10 21:40:29 +02:00
|
|
|
assert result.equals(cmath.complex(-4,-5))
|
2019-07-06 20:19:09 +02:00
|
|
|
}
|
|
|
|
|
|
|
|
fn test_complex_conjugate() {
|
2019-07-10 21:40:29 +02:00
|
|
|
mut c1 := cmath.complex(0,8)
|
2019-07-06 20:19:09 +02:00
|
|
|
mut result := c1.conjugate()
|
2019-07-10 21:40:29 +02:00
|
|
|
assert result.equals(cmath.complex(0,-8))
|
|
|
|
c1 = cmath.complex(7,3)
|
2019-07-06 20:19:09 +02:00
|
|
|
result = c1.conjugate()
|
2019-07-10 21:40:29 +02:00
|
|
|
assert result.equals(cmath.complex(7,-3))
|
|
|
|
c1 = cmath.complex(2,2)
|
2019-07-06 20:19:09 +02:00
|
|
|
result = c1.conjugate()
|
2019-07-10 21:40:29 +02:00
|
|
|
assert result.equals(cmath.complex(2,-2))
|
|
|
|
c1 = cmath.complex(7,0)
|
2019-07-06 20:19:09 +02:00
|
|
|
result = c1.conjugate()
|
2019-07-10 21:40:29 +02:00
|
|
|
assert result.equals(cmath.complex(7,0))
|
2019-07-06 20:19:09 +02:00
|
|
|
}
|
|
|
|
|
|
|
|
fn test_complex_equals() {
|
2019-07-10 21:40:29 +02:00
|
|
|
mut c1 := cmath.complex(0,8)
|
|
|
|
mut c2 := cmath.complex(0,8)
|
2019-07-06 20:19:09 +02:00
|
|
|
assert c1.equals(c2)
|
2019-07-10 21:40:29 +02:00
|
|
|
c1 = cmath.complex(-3,19)
|
|
|
|
c2 = cmath.complex(-3,19)
|
2019-07-06 20:19:09 +02:00
|
|
|
assert c1.equals(c2)
|
2019-07-07 16:13:17 +02:00
|
|
|
}
|
|
|
|
|
2019-07-09 21:12:52 +02:00
|
|
|
fn test_complex_abs() {
|
2019-07-10 21:40:29 +02:00
|
|
|
mut c1 := cmath.complex(3,4)
|
2019-07-09 21:12:52 +02:00
|
|
|
assert c1.abs() == 5
|
2019-07-10 21:40:29 +02:00
|
|
|
c1 = cmath.complex(1,2)
|
2019-07-09 21:12:52 +02:00
|
|
|
assert c1.abs() == math.sqrt(5)
|
|
|
|
assert c1.abs() == c1.conjugate().abs()
|
2019-07-10 21:40:29 +02:00
|
|
|
c1 = cmath.complex(7,0)
|
2019-07-09 21:12:52 +02:00
|
|
|
assert c1.abs() == 7
|
|
|
|
}
|
|
|
|
|
2019-07-07 16:13:17 +02:00
|
|
|
fn test_complex_angle(){
|
2019-07-10 21:40:29 +02:00
|
|
|
mut c := cmath.complex(1, 0)
|
2019-07-07 16:13:17 +02:00
|
|
|
assert c.angle() * 180 / math.Pi == 0
|
2019-07-10 21:40:29 +02:00
|
|
|
c = cmath.complex(1, 1)
|
2019-07-07 16:13:17 +02:00
|
|
|
assert c.angle() * 180 / math.Pi == 45
|
2019-07-10 21:40:29 +02:00
|
|
|
c = cmath.complex(0, 1)
|
2019-07-07 16:13:17 +02:00
|
|
|
assert c.angle() * 180 / math.Pi == 90
|
2019-07-10 21:40:29 +02:00
|
|
|
c = cmath.complex(-1, 1)
|
2019-07-07 16:13:17 +02:00
|
|
|
assert c.angle() * 180 / math.Pi == 135
|
2019-07-10 21:40:29 +02:00
|
|
|
c = cmath.complex(-1, -1)
|
2019-07-07 16:13:17 +02:00
|
|
|
assert c.angle() * 180 / math.Pi == -135
|
|
|
|
mut cc := c.conjugate()
|
|
|
|
assert cc.angle() + c.angle() == 0
|
|
|
|
}
|
2019-07-08 19:41:37 +02:00
|
|
|
|
|
|
|
|
|
|
|
fn test_complex_addinv() {
|
|
|
|
// Tests were also verified on Wolfram Alpha
|
2019-07-10 21:40:29 +02:00
|
|
|
mut c1 := cmath.complex(5,7)
|
|
|
|
mut c2 := cmath.complex(-5,-7)
|
2019-07-08 19:41:37 +02:00
|
|
|
mut result := c1.addinv()
|
|
|
|
assert result.equals(c2)
|
2019-07-10 21:40:29 +02:00
|
|
|
c1 = cmath.complex(-3,4)
|
|
|
|
c2 = cmath.complex(3,-4)
|
2019-07-08 19:41:37 +02:00
|
|
|
result = c1.addinv()
|
|
|
|
assert result.equals(c2)
|
2019-07-10 21:40:29 +02:00
|
|
|
c1 = cmath.complex(-1,-2)
|
|
|
|
c2 = cmath.complex(1,2)
|
2019-07-08 19:41:37 +02:00
|
|
|
result = c1.addinv()
|
|
|
|
assert result.equals(c2)
|
|
|
|
}
|
|
|
|
|
|
|
|
fn test_complex_mulinv() {
|
|
|
|
// Tests were also verified on Wolfram Alpha
|
2019-07-10 21:40:29 +02:00
|
|
|
mut c1 := cmath.complex(5,7)
|
|
|
|
mut c2 := cmath.complex(0.067568,-0.094595)
|
2019-07-08 19:41:37 +02:00
|
|
|
mut result := c1.mulinv()
|
|
|
|
// Some issue with precision comparison in f64 using == operator hence serializing to string
|
|
|
|
assert result.str().eq(c2.str())
|
2019-07-10 21:40:29 +02:00
|
|
|
c1 = cmath.complex(-3,4)
|
|
|
|
c2 = cmath.complex(-0.12,-0.16)
|
2019-07-08 19:41:37 +02:00
|
|
|
result = c1.mulinv()
|
|
|
|
assert result.str().eq(c2.str())
|
2019-07-10 21:40:29 +02:00
|
|
|
c1 = cmath.complex(-1,-2)
|
|
|
|
c2 = cmath.complex(-0.2,0.4)
|
2019-07-08 19:41:37 +02:00
|
|
|
result = c1.mulinv()
|
|
|
|
assert result.equals(c2)
|
|
|
|
}
|
|
|
|
|
|
|
|
fn test_complex_mod() {
|
|
|
|
// Tests were also verified on Wolfram Alpha
|
2019-07-10 21:40:29 +02:00
|
|
|
mut c1 := cmath.complex(5,7)
|
2019-07-08 19:41:37 +02:00
|
|
|
mut result := c1.mod()
|
|
|
|
// Some issue with precision comparison in f64 using == operator hence serializing to string
|
|
|
|
assert result.str().eq('8.602325')
|
2019-07-10 21:40:29 +02:00
|
|
|
c1 = cmath.complex(-3,4)
|
2019-07-08 19:41:37 +02:00
|
|
|
result = c1.mod()
|
|
|
|
assert result == 5
|
2019-07-10 21:40:29 +02:00
|
|
|
c1 = cmath.complex(-1,-2)
|
2019-07-08 19:41:37 +02:00
|
|
|
result = c1.mod()
|
|
|
|
// Some issue with precision comparison in f64 using == operator hence serializing to string
|
|
|
|
assert result.str().eq('2.236068')
|
|
|
|
}
|
|
|
|
|
|
|
|
fn test_complex_pow() {
|
|
|
|
// Tests were also verified on Wolfram Alpha
|
2019-07-10 21:40:29 +02:00
|
|
|
mut c1 := cmath.complex(5,7)
|
|
|
|
mut c2 := cmath.complex(-24.0,70.0)
|
2019-07-08 19:41:37 +02:00
|
|
|
mut result := c1.pow(2)
|
|
|
|
// Some issue with precision comparison in f64 using == operator hence serializing to string
|
|
|
|
assert result.str().eq(c2.str())
|
2019-07-10 21:40:29 +02:00
|
|
|
c1 = cmath.complex(-3,4)
|
|
|
|
c2 = cmath.complex(117,44)
|
2019-07-08 19:41:37 +02:00
|
|
|
result = c1.pow(3)
|
|
|
|
// Some issue with precision comparison in f64 using == operator hence serializing to string
|
|
|
|
assert result.str().eq(c2.str())
|
2019-07-10 21:40:29 +02:00
|
|
|
c1 = cmath.complex(-1,-2)
|
|
|
|
c2 = cmath.complex(-7,-24)
|
2019-07-08 19:41:37 +02:00
|
|
|
result = c1.pow(4)
|
|
|
|
// Some issue with precision comparison in f64 using == operator hence serializing to string
|
|
|
|
assert result.str().eq(c2.str())
|
|
|
|
}
|
|
|
|
|
|
|
|
fn test_complex_root() {
|
|
|
|
// Tests were also verified on Wolfram Alpha
|
2019-07-10 21:40:29 +02:00
|
|
|
mut c1 := cmath.complex(5,7)
|
|
|
|
mut c2 := cmath.complex(2.607904,1.342074)
|
2019-07-08 19:41:37 +02:00
|
|
|
mut result := c1.root(2)
|
|
|
|
// Some issue with precision comparison in f64 using == operator hence serializing to string
|
|
|
|
assert result.str().eq(c2.str())
|
2019-07-10 21:40:29 +02:00
|
|
|
c1 = cmath.complex(-3,4)
|
|
|
|
c2 = cmath.complex(1.264953,1.150614)
|
2019-07-08 19:41:37 +02:00
|
|
|
result = c1.root(3)
|
|
|
|
// Some issue with precision comparison in f64 using == operator hence serializing to string
|
|
|
|
assert result.str().eq(c2.str())
|
2019-07-10 21:40:29 +02:00
|
|
|
c1 = cmath.complex(-1,-2)
|
|
|
|
c2 = cmath.complex(1.068059,-0.595482)
|
2019-07-08 19:41:37 +02:00
|
|
|
result = c1.root(4)
|
|
|
|
// Some issue with precision comparison in f64 using == operator hence serializing to string
|
|
|
|
assert result.str().eq(c2.str())
|
|
|
|
}
|
|
|
|
|
|
|
|
fn test_complex_exp() {
|
|
|
|
// Tests were also verified on Wolfram Alpha
|
2019-07-10 21:40:29 +02:00
|
|
|
mut c1 := cmath.complex(5,7)
|
|
|
|
mut c2 := cmath.complex(111.889015,97.505457)
|
2019-07-08 19:41:37 +02:00
|
|
|
mut result := c1.exp()
|
|
|
|
// Some issue with precision comparison in f64 using == operator hence serializing to string
|
|
|
|
assert result.str().eq(c2.str())
|
2019-07-10 21:40:29 +02:00
|
|
|
c1 = cmath.complex(-3,4)
|
|
|
|
c2 = cmath.complex(-0.032543,-0.037679)
|
2019-07-08 19:41:37 +02:00
|
|
|
result = c1.exp()
|
|
|
|
// Some issue with precision comparison in f64 using == operator hence serializing to string
|
|
|
|
assert result.str().eq(c2.str())
|
2019-07-10 21:40:29 +02:00
|
|
|
c1 = cmath.complex(-1,-2)
|
|
|
|
c2 = cmath.complex(-0.153092,-0.334512)
|
2019-07-08 19:41:37 +02:00
|
|
|
result = c1.exp()
|
|
|
|
// Some issue with precision comparison in f64 using == operator hence serializing to string
|
|
|
|
assert result.str().eq(c2.str())
|
|
|
|
}
|
|
|
|
|
|
|
|
fn test_complex_ln() {
|
|
|
|
// Tests were also verified on Wolfram Alpha
|
2019-07-10 21:40:29 +02:00
|
|
|
mut c1 := cmath.complex(5,7)
|
|
|
|
mut c2 := cmath.complex(2.152033,0.950547)
|
2019-07-08 19:41:37 +02:00
|
|
|
mut result := c1.ln()
|
|
|
|
// Some issue with precision comparison in f64 using == operator hence serializing to string
|
|
|
|
assert result.str().eq(c2.str())
|
2019-07-10 21:40:29 +02:00
|
|
|
c1 = cmath.complex(-3,4)
|
|
|
|
c2 = cmath.complex(1.609438,2.214297)
|
2019-07-08 19:41:37 +02:00
|
|
|
result = c1.ln()
|
|
|
|
// Some issue with precision comparison in f64 using == operator hence serializing to string
|
|
|
|
assert result.str().eq(c2.str())
|
2019-07-10 21:40:29 +02:00
|
|
|
c1 = cmath.complex(-1,-2)
|
|
|
|
c2 = cmath.complex(0.804719,-2.034444)
|
2019-07-08 19:41:37 +02:00
|
|
|
result = c1.ln()
|
|
|
|
// Some issue with precision comparison in f64 using == operator hence serializing to string
|
|
|
|
assert result.str().eq(c2.str())
|
|
|
|
}
|
2019-07-10 19:28:19 +02:00
|
|
|
|
|
|
|
fn test_complex_sin() {
|
|
|
|
// Tests were also verified on Wolfram Alpha
|
2019-07-10 21:40:29 +02:00
|
|
|
mut c1 := cmath.complex(5,7)
|
|
|
|
mut c2 := cmath.complex(-525.794515,155.536550)
|
2019-07-10 19:28:19 +02:00
|
|
|
mut result := c1.sin()
|
|
|
|
// Some issue with precision comparison in f64 using == operator hence serializing to string
|
|
|
|
assert result.str().eq(c2.str())
|
2019-07-10 21:40:29 +02:00
|
|
|
c1 = cmath.complex(-3,4)
|
|
|
|
c2 = cmath.complex(-3.853738,-27.016813)
|
2019-07-10 19:28:19 +02:00
|
|
|
result = c1.sin()
|
|
|
|
// Some issue with precision comparison in f64 using == operator hence serializing to string
|
|
|
|
assert result.str().eq(c2.str())
|
2019-07-10 21:40:29 +02:00
|
|
|
c1 = cmath.complex(-1,-2)
|
|
|
|
c2 = cmath.complex(-3.165779,-1.959601)
|
2019-07-10 19:28:19 +02:00
|
|
|
result = c1.sin()
|
|
|
|
// Some issue with precision comparison in f64 using == operator hence serializing to string
|
|
|
|
assert result.str().eq(c2.str())
|
|
|
|
}
|
|
|
|
|
|
|
|
fn test_complex_cos() {
|
|
|
|
// Tests were also verified on Wolfram Alpha
|
2019-07-10 21:40:29 +02:00
|
|
|
mut c1 := cmath.complex(5,7)
|
|
|
|
mut c2 := cmath.complex(155.536809,525.793641)
|
2019-07-10 19:28:19 +02:00
|
|
|
mut result := c1.cos()
|
|
|
|
// Some issue with precision comparison in f64 using == operator hence serializing to string
|
|
|
|
assert result.str().eq(c2.str())
|
2019-07-10 21:40:29 +02:00
|
|
|
c1 = cmath.complex(-3,4)
|
|
|
|
c2 = cmath.complex(-27.034946,3.851153)
|
2019-07-10 19:28:19 +02:00
|
|
|
result = c1.cos()
|
|
|
|
// Some issue with precision comparison in f64 using == operator hence serializing to string
|
|
|
|
assert result.str().eq(c2.str())
|
2019-07-10 21:40:29 +02:00
|
|
|
c1 = cmath.complex(-1,-2)
|
|
|
|
c2 = cmath.complex(2.032723,-3.051898)
|
2019-07-10 19:28:19 +02:00
|
|
|
result = c1.cos()
|
|
|
|
// Some issue with precision comparison in f64 using == operator hence serializing to string
|
|
|
|
assert result.str().eq(c2.str())
|
|
|
|
}
|
|
|
|
|
|
|
|
fn test_complex_tan() {
|
|
|
|
// Tests were also verified on Wolfram Alpha
|
2019-07-10 21:40:29 +02:00
|
|
|
mut c1 := cmath.complex(5,7)
|
|
|
|
mut c2 := cmath.complex(-0.000001,1.000001)
|
2019-07-10 19:28:19 +02:00
|
|
|
mut result := c1.tan()
|
|
|
|
// Some issue with precision comparison in f64 using == operator hence serializing to string
|
|
|
|
assert result.str().eq(c2.str())
|
2019-07-10 21:40:29 +02:00
|
|
|
c1 = cmath.complex(-3,4)
|
|
|
|
c2 = cmath.complex(0.000187,0.999356)
|
2019-07-10 19:28:19 +02:00
|
|
|
result = c1.tan()
|
|
|
|
// Some issue with precision comparison in f64 using == operator hence serializing to string
|
|
|
|
assert result.str().eq(c2.str())
|
2019-07-10 21:40:29 +02:00
|
|
|
c1 = cmath.complex(-1,-2)
|
|
|
|
c2 = cmath.complex(-0.033813,-1.014794)
|
2019-07-10 19:28:19 +02:00
|
|
|
result = c1.tan()
|
|
|
|
// Some issue with precision comparison in f64 using == operator hence serializing to string
|
|
|
|
assert result.str().eq(c2.str())
|
|
|
|
}
|
|
|
|
|
2019-07-12 12:56:21 +02:00
|
|
|
fn test_complex_cot() {
|
|
|
|
// Tests were also verified on Wolfram Alpha
|
|
|
|
mut c1 := cmath.complex(5,7)
|
|
|
|
mut c2 := cmath.complex(-0.000001,-0.999999)
|
|
|
|
mut result := c1.cot()
|
|
|
|
// Some issue with precision comparison in f64 using == operator hence serializing to string
|
|
|
|
assert result.str().eq(c2.str())
|
|
|
|
c1 = cmath.complex(-3,4)
|
|
|
|
c2 = cmath.complex(0.000188,-1.000644)
|
|
|
|
result = c1.cot()
|
|
|
|
// Some issue with precision comparison in f64 using == operator hence serializing to string
|
|
|
|
assert result.str().eq(c2.str())
|
|
|
|
c1 = cmath.complex(-1,-2)
|
|
|
|
c2 = cmath.complex(-0.032798,0.984329)
|
|
|
|
result = c1.cot()
|
|
|
|
// Some issue with precision comparison in f64 using == operator hence serializing to string
|
|
|
|
assert result.str().eq(c2.str())
|
|
|
|
}
|
|
|
|
|
|
|
|
fn test_complex_sec() {
|
|
|
|
// Tests were also verified on Wolfram Alpha
|
|
|
|
mut c1 := cmath.complex(5,7)
|
|
|
|
mut c2 := cmath.complex(0.000517,-0.001749)
|
|
|
|
mut result := c1.sec()
|
|
|
|
// Some issue with precision comparison in f64 using == operator hence serializing to string
|
|
|
|
assert result.str().eq(c2.str())
|
|
|
|
c1 = cmath.complex(-3,4)
|
|
|
|
c2 = cmath.complex(-0.036253,-0.005164)
|
|
|
|
result = c1.sec()
|
|
|
|
// Some issue with precision comparison in f64 using == operator hence serializing to string
|
|
|
|
assert result.str().eq(c2.str())
|
|
|
|
c1 = cmath.complex(-1,-2)
|
|
|
|
c2 = cmath.complex(0.151176,0.226974)
|
|
|
|
result = c1.sec()
|
|
|
|
// Some issue with precision comparison in f64 using == operator hence serializing to string
|
|
|
|
assert result.str().eq(c2.str())
|
|
|
|
}
|
|
|
|
|
|
|
|
fn test_complex_csc() {
|
|
|
|
// Tests were also verified on Wolfram Alpha
|
|
|
|
mut c1 := cmath.complex(5,7)
|
|
|
|
mut c2 := cmath.complex(-0.001749,-0.000517)
|
|
|
|
mut result := c1.csc()
|
|
|
|
// Some issue with precision comparison in f64 using == operator hence serializing to string
|
|
|
|
assert result.str().eq(c2.str())
|
|
|
|
c1 = cmath.complex(-3,4)
|
|
|
|
c2 = cmath.complex(-0.005174,0.036276)
|
|
|
|
result = c1.csc()
|
|
|
|
// Some issue with precision comparison in f64 using == operator hence serializing to string
|
|
|
|
assert result.str().eq(c2.str())
|
|
|
|
c1 = cmath.complex(-1,-2)
|
|
|
|
c2 = cmath.complex(-0.228375,0.141363)
|
|
|
|
result = c1.csc()
|
|
|
|
// Some issue with precision comparison in f64 using == operator hence serializing to string
|
|
|
|
assert result.str().eq(c2.str())
|
|
|
|
}
|
|
|
|
|
2019-07-11 15:35:06 +02:00
|
|
|
fn test_complex_asin() {
|
|
|
|
// Tests were also verified on Wolfram Alpha
|
|
|
|
mut c1 := cmath.complex(5,7)
|
|
|
|
mut c2 := cmath.complex(0.617064,2.846289)
|
|
|
|
mut result := c1.asin()
|
|
|
|
// Some issue with precision comparison in f64 using == operator hence serializing to string
|
|
|
|
assert result.str().eq(c2.str())
|
|
|
|
c1 = cmath.complex(-3,4)
|
|
|
|
c2 = cmath.complex(-0.633984,2.305509)
|
|
|
|
result = c1.asin()
|
|
|
|
// Some issue with precision comparison in f64 using == operator hence serializing to string
|
|
|
|
assert result.str().eq(c2.str())
|
|
|
|
c1 = cmath.complex(-1,-2)
|
|
|
|
c2 = cmath.complex(-0.427079,-1.528571)
|
|
|
|
result = c1.asin()
|
|
|
|
// Some issue with precision comparison in f64 using == operator hence serializing to string
|
|
|
|
assert result.str().eq(c2.str())
|
|
|
|
}
|
|
|
|
|
|
|
|
fn test_complex_acos() {
|
|
|
|
// Tests were also verified on Wolfram Alpha
|
|
|
|
mut c1 := cmath.complex(5,7)
|
|
|
|
mut c2 := cmath.complex(0.953732,-2.846289)
|
|
|
|
mut result := c1.acos()
|
|
|
|
// Some issue with precision comparison in f64 using == operator hence serializing to string
|
|
|
|
assert result.str().eq(c2.str())
|
|
|
|
c1 = cmath.complex(-3,4)
|
|
|
|
c2 = cmath.complex(2.204780,-2.305509)
|
|
|
|
result = c1.acos()
|
|
|
|
// Some issue with precision comparison in f64 using == operator hence serializing to string
|
|
|
|
assert result.str().eq(c2.str())
|
|
|
|
c1 = cmath.complex(-1,-2)
|
|
|
|
c2 = cmath.complex(1.997875,1.528571)
|
|
|
|
result = c1.acos()
|
|
|
|
// Some issue with precision comparison in f64 using == operator hence serializing to string
|
|
|
|
assert result.str().eq(c2.str())
|
|
|
|
}
|
|
|
|
|
|
|
|
fn test_complex_atan() {
|
|
|
|
// Tests were also verified on Wolfram Alpha
|
|
|
|
mut c1 := cmath.complex(5,7)
|
|
|
|
mut c2 := cmath.complex(1.502727,0.094441)
|
|
|
|
mut result := c1.atan()
|
|
|
|
// Some issue with precision comparison in f64 using == operator hence serializing to string
|
|
|
|
assert result.str().eq(c2.str())
|
|
|
|
c1 = cmath.complex(-3,4)
|
|
|
|
c2 = cmath.complex(-1.448307,0.158997)
|
|
|
|
result = c1.atan()
|
|
|
|
// Some issue with precision comparison in f64 using == operator hence serializing to string
|
|
|
|
assert result.str().eq(c2.str())
|
|
|
|
c1 = cmath.complex(-1,-2)
|
|
|
|
c2 = cmath.complex(-1.338973,-0.402359)
|
|
|
|
result = c1.atan()
|
|
|
|
// Some issue with precision comparison in f64 using == operator hence serializing to string
|
|
|
|
assert result.str().eq(c2.str())
|
|
|
|
}
|
|
|
|
|
2019-07-12 12:56:21 +02:00
|
|
|
fn test_complex_acot() {
|
|
|
|
// Tests were also verified on Wolfram Alpha
|
|
|
|
mut c1 := cmath.complex(5,7)
|
|
|
|
mut c2 := cmath.complex(0.068069,-0.094441)
|
|
|
|
mut result := c1.acot()
|
|
|
|
// Some issue with precision comparison in f64 using == operator hence serializing to string
|
|
|
|
assert result.str().eq(c2.str())
|
|
|
|
c1 = cmath.complex(-3,4)
|
|
|
|
c2 = cmath.complex(-0.122489,-0.158997)
|
|
|
|
result = c1.acot()
|
|
|
|
// Some issue with precision comparison in f64 using == operator hence serializing to string
|
|
|
|
assert result.str().eq(c2.str())
|
|
|
|
c1 = cmath.complex(-1,-2)
|
|
|
|
c2 = cmath.complex(-0.231824,0.402359)
|
|
|
|
result = c1.acot()
|
|
|
|
// Some issue with precision comparison in f64 using == operator hence serializing to string
|
|
|
|
assert result.str().eq(c2.str())
|
|
|
|
}
|
|
|
|
|
|
|
|
fn test_complex_asec() {
|
|
|
|
// Tests were also verified on Wolfram Alpha
|
|
|
|
mut c1 := cmath.complex(5,7)
|
|
|
|
mut c2 := cmath.complex(1.503480,0.094668)
|
|
|
|
mut result := c1.asec()
|
|
|
|
// Some issue with precision comparison in f64 using == operator hence serializing to string
|
|
|
|
assert result.str().eq(c2.str())
|
|
|
|
c1 = cmath.complex(-3,4)
|
|
|
|
c2 = cmath.complex(1.689547,0.160446)
|
|
|
|
result = c1.asec()
|
|
|
|
// Some issue with precision comparison in f64 using == operator hence serializing to string
|
|
|
|
assert result.str().eq(c2.str())
|
|
|
|
c1 = cmath.complex(-1,-2)
|
|
|
|
c2 = cmath.complex(1.757114,-0.396568)
|
|
|
|
result = c1.asec()
|
|
|
|
// Some issue with precision comparison in f64 using == operator hence serializing to string
|
|
|
|
assert result.str().eq(c2.str())
|
|
|
|
}
|
|
|
|
|
|
|
|
fn test_complex_acsc() {
|
|
|
|
// Tests were also verified on Wolfram Alpha
|
|
|
|
mut c1 := cmath.complex(5,7)
|
|
|
|
mut c2 := cmath.complex(0.067317,-0.094668)
|
|
|
|
mut result := c1.acsc()
|
|
|
|
// Some issue with precision comparison in f64 using == operator hence serializing to string
|
|
|
|
assert result.str().eq(c2.str())
|
|
|
|
c1 = cmath.complex(-3,4)
|
|
|
|
c2 = cmath.complex(-0.118751,-0.160446)
|
|
|
|
result = c1.acsc()
|
|
|
|
// Some issue with precision comparison in f64 using == operator hence serializing to string
|
|
|
|
assert result.str().eq(c2.str())
|
|
|
|
c1 = cmath.complex(-1,-2)
|
|
|
|
c2 = cmath.complex(-0.186318,0.396568)
|
|
|
|
result = c1.acsc()
|
|
|
|
// Some issue with precision comparison in f64 using == operator hence serializing to string
|
|
|
|
assert result.str().eq(c2.str())
|
|
|
|
}
|
|
|
|
|
2019-07-10 19:28:19 +02:00
|
|
|
fn test_complex_sinh() {
|
|
|
|
// Tests were also verified on Wolfram Alpha
|
2019-07-10 21:40:29 +02:00
|
|
|
mut c1 := cmath.complex(5,7)
|
|
|
|
mut c2 := cmath.complex(55.941968,48.754942)
|
2019-07-10 19:28:19 +02:00
|
|
|
mut result := c1.sinh()
|
|
|
|
// Some issue with precision comparison in f64 using == operator hence serializing to string
|
|
|
|
assert result.str().eq(c2.str())
|
2019-07-10 21:40:29 +02:00
|
|
|
c1 = cmath.complex(-3,4)
|
|
|
|
c2 = cmath.complex(6.548120,-7.619232)
|
2019-07-10 19:28:19 +02:00
|
|
|
result = c1.sinh()
|
|
|
|
// Some issue with precision comparison in f64 using == operator hence serializing to string
|
|
|
|
assert result.str().eq(c2.str())
|
2019-07-10 21:40:29 +02:00
|
|
|
c1 = cmath.complex(-1,-2)
|
|
|
|
c2 = cmath.complex(0.489056,-1.403119)
|
2019-07-10 19:28:19 +02:00
|
|
|
result = c1.sinh()
|
|
|
|
// Some issue with precision comparison in f64 using == operator hence serializing to string
|
|
|
|
assert result.str().eq(c2.str())
|
|
|
|
}
|
|
|
|
|
|
|
|
fn test_complex_cosh() {
|
|
|
|
// Tests were also verified on Wolfram Alpha
|
2019-07-10 21:40:29 +02:00
|
|
|
mut c1 := cmath.complex(5,7)
|
|
|
|
mut c2 := cmath.complex(55.947047,48.750515)
|
2019-07-10 19:28:19 +02:00
|
|
|
mut result := c1.cosh()
|
|
|
|
// Some issue with precision comparison in f64 using == operator hence serializing to string
|
|
|
|
assert result.str().eq(c2.str())
|
2019-07-10 21:40:29 +02:00
|
|
|
c1 = cmath.complex(-3,4)
|
|
|
|
c2 = cmath.complex(-6.580663,7.581553)
|
2019-07-10 19:28:19 +02:00
|
|
|
result = c1.cosh()
|
|
|
|
// Some issue with precision comparison in f64 using == operator hence serializing to string
|
|
|
|
assert result.str().eq(c2.str())
|
2019-07-10 21:40:29 +02:00
|
|
|
c1 = cmath.complex(-1,-2)
|
|
|
|
c2 = cmath.complex(-0.642148,1.068607)
|
2019-07-10 19:28:19 +02:00
|
|
|
result = c1.cosh()
|
|
|
|
// Some issue with precision comparison in f64 using == operator hence serializing to string
|
|
|
|
assert result.str().eq(c2.str())
|
|
|
|
}
|
|
|
|
|
|
|
|
fn test_complex_tanh() {
|
|
|
|
// Tests were also verified on Wolfram Alpha
|
2019-07-10 21:40:29 +02:00
|
|
|
mut c1 := cmath.complex(5,7)
|
|
|
|
mut c2 := cmath.complex(0.999988,0.000090)
|
2019-07-10 19:28:19 +02:00
|
|
|
mut result := c1.tanh()
|
|
|
|
// Some issue with precision comparison in f64 using == operator hence serializing to string
|
|
|
|
assert result.str().eq(c2.str())
|
2019-07-10 21:40:29 +02:00
|
|
|
c1 = cmath.complex(-3,4)
|
|
|
|
c2 = cmath.complex(-1.000710,0.004908)
|
2019-07-10 19:28:19 +02:00
|
|
|
result = c1.tanh()
|
|
|
|
// Some issue with precision comparison in f64 using == operator hence serializing to string
|
|
|
|
assert result.str().eq(c2.str())
|
2019-07-10 21:40:29 +02:00
|
|
|
c1 = cmath.complex(-1,-2)
|
|
|
|
c2 = cmath.complex(-1.166736,0.243458)
|
2019-07-10 19:28:19 +02:00
|
|
|
result = c1.tanh()
|
|
|
|
// Some issue with precision comparison in f64 using == operator hence serializing to string
|
|
|
|
assert result.str().eq(c2.str())
|
2019-07-10 21:40:29 +02:00
|
|
|
}
|
2019-07-12 12:56:21 +02:00
|
|
|
|
|
|
|
fn test_complex_coth() {
|
|
|
|
// Tests were also verified on Wolfram Alpha
|
|
|
|
mut c1 := cmath.complex(5,7)
|
|
|
|
mut c2 := cmath.complex(1.000012,-0.000090)
|
|
|
|
mut result := c1.coth()
|
|
|
|
// Some issue with precision comparison in f64 using == operator hence serializing to string
|
|
|
|
assert result.str().eq(c2.str())
|
|
|
|
c1 = cmath.complex(-3,4)
|
|
|
|
c2 = cmath.complex(-0.999267,-0.004901)
|
|
|
|
result = c1.coth()
|
|
|
|
// Some issue with precision comparison in f64 using == operator hence serializing to string
|
|
|
|
assert result.str().eq(c2.str())
|
|
|
|
c1 = cmath.complex(-1,-2)
|
|
|
|
c2 = cmath.complex(-0.821330,-0.171384)
|
|
|
|
result = c1.coth()
|
|
|
|
// Some issue with precision comparison in f64 using == operator hence serializing to string
|
|
|
|
assert result.str().eq(c2.str())
|
|
|
|
}
|
2019-07-10 21:40:29 +02:00
|
|
|
|
2019-07-12 12:56:21 +02:00
|
|
|
fn test_complex_sech() {
|
|
|
|
// Tests were also verified on Wolfram Alpha
|
|
|
|
mut c1 := cmath.complex(5,7)
|
|
|
|
mut c2 := cmath.complex(0.010160,-0.008853)
|
|
|
|
mut result := c1.sech()
|
|
|
|
// Some issue with precision comparison in f64 using == operator hence serializing to string
|
|
|
|
assert result.str().eq(c2.str())
|
|
|
|
c1 = cmath.complex(-3,4)
|
|
|
|
c2 = cmath.complex(-0.065294,-0.075225)
|
|
|
|
result = c1.sech()
|
|
|
|
// Some issue with precision comparison in f64 using == operator hence serializing to string
|
|
|
|
assert result.str().eq(c2.str())
|
|
|
|
c1 = cmath.complex(-1,-2)
|
|
|
|
c2 = cmath.complex(-0.413149,-0.687527)
|
|
|
|
result = c1.sech()
|
|
|
|
// Some issue with precision comparison in f64 using == operator hence serializing to string
|
|
|
|
assert result.str().eq(c2.str())
|
|
|
|
}
|
|
|
|
|
|
|
|
fn test_complex_csch() {
|
|
|
|
// Tests were also verified on Wolfram Alpha
|
|
|
|
mut c1 := cmath.complex(5,7)
|
|
|
|
mut c2 := cmath.complex(0.010159,-0.008854)
|
|
|
|
mut result := c1.csch()
|
|
|
|
// Some issue with precision comparison in f64 using == operator hence serializing to string
|
|
|
|
assert result.str().eq(c2.str())
|
|
|
|
c1 = cmath.complex(-3,4)
|
|
|
|
c2 = cmath.complex(0.064877,0.075490)
|
|
|
|
result = c1.csch()
|
|
|
|
// Some issue with precision comparison in f64 using == operator hence serializing to string
|
|
|
|
assert result.str().eq(c2.str())
|
|
|
|
c1 = cmath.complex(-1,-2)
|
|
|
|
c2 = cmath.complex(0.221501,0.635494)
|
|
|
|
result = c1.csch()
|
|
|
|
// Some issue with precision comparison in f64 using == operator hence serializing to string
|
|
|
|
assert result.str().eq(c2.str())
|
|
|
|
}
|
|
|
|
|
2019-07-11 15:35:06 +02:00
|
|
|
fn test_complex_asinh() {
|
|
|
|
// Tests were also verified on Wolfram Alpha
|
|
|
|
mut c1 := cmath.complex(5,7)
|
|
|
|
mut c2 := cmath.complex(2.844098,0.947341)
|
|
|
|
mut result := c1.asinh()
|
|
|
|
// Some issue with precision comparison in f64 using == operator hence serializing to string
|
|
|
|
assert result.str().eq(c2.str())
|
|
|
|
c1 = cmath.complex(-3,4)
|
|
|
|
c2 = cmath.complex(-2.299914,0.917617)
|
|
|
|
result = c1.asinh()
|
|
|
|
// Some issue with precision comparison in f64 using == operator hence serializing to string
|
|
|
|
assert result.str().eq(c2.str())
|
|
|
|
c1 = cmath.complex(-1,-2)
|
|
|
|
c2 = cmath.complex(-1.469352,-1.063440)
|
|
|
|
result = c1.asinh()
|
|
|
|
// Some issue with precision comparison in f64 using == operator hence serializing to string
|
|
|
|
assert result.str().eq(c2.str())
|
|
|
|
}
|
|
|
|
|
|
|
|
fn test_complex_acosh() {
|
|
|
|
// Tests were also verified on Wolfram Alpha
|
|
|
|
mut c1 := cmath.complex(5,7)
|
|
|
|
mut c2 := cmath.complex(2.846289,0.953732)
|
|
|
|
mut result := c1.acosh()
|
|
|
|
// Some issue with precision comparison in f64 using == operator hence serializing to string
|
|
|
|
assert result.str().eq(c2.str())
|
|
|
|
c1 = cmath.complex(-3,4)
|
|
|
|
c2 = cmath.complex(2.305509,2.204780)
|
|
|
|
result = c1.acosh()
|
|
|
|
// Some issue with precision comparison in f64 using == operator hence serializing to string
|
|
|
|
assert result.str().eq(c2.str())
|
|
|
|
c1 = cmath.complex(-1,-2)
|
|
|
|
c2 = cmath.complex(1.528571,-1.997875)
|
|
|
|
result = c1.acosh()
|
|
|
|
// Some issue with precision comparison in f64 using == operator hence serializing to string
|
|
|
|
assert result.str().eq(c2.str())
|
|
|
|
}
|
2019-07-10 21:40:29 +02:00
|
|
|
|
2019-07-11 15:35:06 +02:00
|
|
|
fn test_complex_atanh() {
|
|
|
|
// Tests were also verified on Wolfram Alpha
|
|
|
|
mut c1 := cmath.complex(5,7)
|
|
|
|
mut c2 := cmath.complex(0.067066,1.476056)
|
|
|
|
mut result := c1.atanh()
|
|
|
|
// Some issue with precision comparison in f64 using == operator hence serializing to string
|
|
|
|
assert result.str().eq(c2.str())
|
|
|
|
c1 = cmath.complex(-3,4)
|
|
|
|
c2 = cmath.complex(-0.117501,1.409921)
|
|
|
|
result = c1.atanh()
|
|
|
|
// Some issue with precision comparison in f64 using == operator hence serializing to string
|
|
|
|
assert result.str().eq(c2.str())
|
|
|
|
c1 = cmath.complex(-1,-2)
|
|
|
|
c2 = cmath.complex(-0.173287,-1.178097)
|
|
|
|
result = c1.atanh()
|
|
|
|
// Some issue with precision comparison in f64 using == operator hence serializing to string
|
|
|
|
assert result.str().eq(c2.str())
|
2019-07-12 12:56:21 +02:00
|
|
|
}
|
|
|
|
|
|
|
|
fn test_complex_acoth() {
|
|
|
|
// Tests were also verified on Wolfram Alpha
|
|
|
|
mut c1 := cmath.complex(5,7)
|
|
|
|
mut c2 := cmath.complex(0.067066,-0.094740)
|
|
|
|
mut result := c1.acoth()
|
|
|
|
// Some issue with precision comparison in f64 using == operator hence serializing to string
|
|
|
|
assert result.str().eq(c2.str())
|
|
|
|
c1 = cmath.complex(-3,4)
|
|
|
|
c2 = cmath.complex(-0.117501,-0.160875)
|
|
|
|
result = c1.acoth()
|
|
|
|
// Some issue with precision comparison in f64 using == operator hence serializing to string
|
|
|
|
assert result.str().eq(c2.str())
|
|
|
|
c1 = cmath.complex(-1,-2)
|
|
|
|
c2 = cmath.complex(-0.173287,0.392699)
|
|
|
|
result = c1.acoth()
|
|
|
|
// Some issue with precision comparison in f64 using == operator hence serializing to string
|
|
|
|
assert result.str().eq(c2.str())
|
|
|
|
}
|
|
|
|
|
|
|
|
// fn test_complex_asech() {
|
|
|
|
// // Tests were also verified on Wolfram Alpha
|
|
|
|
// mut c1 := cmath.complex(5,7)
|
|
|
|
// mut c2 := cmath.complex(0.094668,-1.503480)
|
|
|
|
// mut result := c1.asech()
|
|
|
|
// // Some issue with precision comparison in f64 using == operator hence serializing to string
|
|
|
|
// assert result.str().eq(c2.str())
|
|
|
|
// c1 = cmath.complex(-3,4)
|
|
|
|
// c2 = cmath.complex(0.160446,-1.689547)
|
|
|
|
// result = c1.asech()
|
|
|
|
// // Some issue with precision comparison in f64 using == operator hence serializing to string
|
|
|
|
// assert result.str().eq(c2.str())
|
|
|
|
// c1 = cmath.complex(-1,-2)
|
|
|
|
// c2 = cmath.complex(0.396568,1.757114)
|
|
|
|
// result = c1.asech()
|
|
|
|
// // Some issue with precision comparison in f64 using == operator hence serializing to string
|
|
|
|
// assert result.str().eq(c2.str())
|
|
|
|
// }
|
|
|
|
|
|
|
|
fn test_complex_acsch() {
|
|
|
|
// Tests were also verified on Wolfram Alpha
|
|
|
|
mut c1 := cmath.complex(5,7)
|
|
|
|
mut c2 := cmath.complex(0.067819,-0.094518)
|
|
|
|
mut result := c1.acsch()
|
|
|
|
// Some issue with precision comparison in f64 using == operator hence serializing to string
|
|
|
|
assert result.str().eq(c2.str())
|
|
|
|
c1 = cmath.complex(-3,4)
|
|
|
|
c2 = cmath.complex(-0.121246,-0.159507)
|
|
|
|
result = c1.acsch()
|
|
|
|
// Some issue with precision comparison in f64 using == operator hence serializing to string
|
|
|
|
assert result.str().eq(c2.str())
|
|
|
|
c1 = cmath.complex(-1,-2)
|
|
|
|
c2 = cmath.complex(-0.215612,0.401586)
|
|
|
|
result = c1.acsch()
|
|
|
|
// Some issue with precision comparison in f64 using == operator hence serializing to string
|
|
|
|
assert result.str().eq(c2.str())
|
2019-07-11 15:35:06 +02:00
|
|
|
}
|