Complex conjugate pairs
Sort complex numbers into complex conjugate pairs ordered by increasing real part.
cplxpair(z, tol = 100 * .Machine$double.eps, MARGIN = 2)
z |
Vector, matrix, or array of complex numbers. |
tol |
Weighting factor |
MARGIN |
Vector giving the subscripts which the function will be applied over. E.g., for a matrix 1 indicates rows, 2 indicates columns, c(1, 2) indicates rows and columns. Where X has named dimnames, it can be a character vector selecting dimension names. Default: 2 (columns). |
The negative imaginary complex numbers are placed first within each pair. All
real numbers (those with abs(Im (z) / z) < tol) are placed after the
complex pairs.
An error is signaled if some complex numbers could not be paired and if all
complex numbers are not exact conjugates (to within tol).
Vector, matrix or array containing ordered complex conjugate pairs by increasing real parts.
There is no defined order for pairs with identical real parts but differing imaginary parts.
Geert van Boxtel, G.J.M.vanBoxtel@gmail.com.
r <- rbind(t(cplxpair(exp(2i * pi * 0:4 / 5))),
t(exp(2i * pi *c(3, 2, 4, 1, 0) / 5)))Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.