| 1633 | } |
| 1634 | |
| 1635 | static void |
| 1636 | nc_angle(npy_intp n, npy_complex64 *a, npy_intp sb1, npy_float32 *r) { |
| 1637 | |
| 1638 | if( sb1 == sizeof(npy_complex64) ) { // Aligned |
| 1639 | for( npy_intp I = 0; I < n; I++ ) { |
| 1640 | r[I] = atan2f(a[I].imag, a[I].real); |
| 1641 | } |
| 1642 | } |
| 1643 | else { |
| 1644 | sb1 /= sizeof(npy_complex64); |
| 1645 | for( npy_intp I = 0; I < n; I++ ) { |
| 1646 | // Take product of A and complex conjugate of B |
| 1647 | r[I*sb1] = atan2f(a[I*sb1].imag, a[I*sb1].real); |
| 1648 | } |
| 1649 | } |
| 1650 | } |
| 1651 | |
| 1652 | static void |
| 1653 | nc_angle(npy_intp n, npy_complex128 *a, npy_intp sb1, npy_float64 *r) { |
nothing calls this directly
no outgoing calls
no test coverage detected
searching dependent graphs…