99 lines
2.1 KiB
Go
99 lines
2.1 KiB
Go
// Copyright (c) 2019 Alexander Medvednikov. All rights reserved.
|
|
// Use of this source code is governed by an MIT license
|
|
// that can be found in the LICENSE file.
|
|
|
|
module math
|
|
|
|
struct Complex {
|
|
re f64
|
|
im f64
|
|
}
|
|
|
|
pub fn complex(re f64,im f64) Complex {
|
|
return Complex{re,im}
|
|
}
|
|
|
|
// To String method
|
|
pub fn (c Complex) str() string {
|
|
mut out := '$c.re'
|
|
out += if c.im >= 0 {
|
|
'+$c.im'
|
|
}
|
|
else {
|
|
'$c.im'
|
|
}
|
|
out += 'i'
|
|
return out
|
|
}
|
|
|
|
// Complex Angle
|
|
pub fn (c Complex) angle() f64 {
|
|
return atan2(c.im, c.re)
|
|
}
|
|
|
|
// Complex Addition c1 + c2
|
|
pub fn (c1 Complex) + (c2 Complex) Complex {
|
|
return Complex{c1.re+c2.re,c1.im+c2.im}
|
|
}
|
|
|
|
// Complex Substraction c1 - c2
|
|
pub fn (c1 Complex) - (c2 Complex) Complex {
|
|
return Complex{c1.re-c2.re,c1.im-c2.im}
|
|
}
|
|
|
|
// Complex Multiplication c1 * c2
|
|
// Currently Not Supported
|
|
// pub fn (c1 Complex) * (c2 Complex) Complex {
|
|
// return Complex{
|
|
// (c1.re * c2.re) + ((c1.im * c2.im) * -1),
|
|
// (c1.re * c2.im) + (c1.im * c2.re)
|
|
// }
|
|
// }
|
|
|
|
// Complex Division c1 / c2
|
|
// Currently Not Supported
|
|
// pub fn (c1 Complex) / (c2 Complex) Complex {
|
|
// denom := (c2.re * c2.re) + (c2.im * c2.im)
|
|
// return Complex {
|
|
// ((c1.re * c2.re) + ((c1.im * -c2.im) * -1))/denom,
|
|
// ((c1.re * -c2.im) + (c1.im * c2.re))/denom
|
|
// }
|
|
// }
|
|
|
|
// Complex Addition c1.add(c2)
|
|
pub fn (c1 Complex) add(c2 Complex) Complex {
|
|
return c1 + c2
|
|
}
|
|
|
|
// Complex Subtraction c1.subtract(c2)
|
|
pub fn (c1 Complex) subtract(c2 Complex) Complex {
|
|
return c1 - c2
|
|
}
|
|
|
|
// Complex Multiplication c1.multiply(c2)
|
|
pub fn (c1 Complex) multiply(c2 Complex) Complex {
|
|
return Complex{
|
|
(c1.re * c2.re) + ((c1.im * c2.im) * -1),
|
|
(c1.re * c2.im) + (c1.im * c2.re)
|
|
}
|
|
}
|
|
|
|
// Complex Division c1.divide(c2)
|
|
pub fn (c1 Complex) divide(c2 Complex) Complex {
|
|
denom := (c2.re * c2.re) + (c2.im * c2.im)
|
|
return Complex {
|
|
((c1.re * c2.re) + ((c1.im * -c2.im) * -1))/denom,
|
|
((c1.re * -c2.im) + (c1.im * c2.re))/denom
|
|
}
|
|
}
|
|
|
|
// Complex Conjugate
|
|
pub fn (c1 Complex) conjugate() Complex{
|
|
return Complex{c1.re,-c1.im}
|
|
}
|
|
|
|
// Complex Equals
|
|
pub fn (c1 Complex) equals(c2 Complex) bool {
|
|
return (c1.re == c2.re) && (c1.im == c2.im)
|
|
}
|