| 1183 | } |
| 1184 | |
| 1185 | static void |
| 1186 | nc_atan( npy_intp n, npy_complex64 *x, npy_intp sb1, npy_complex64 *r) |
| 1187 | { |
| 1188 | npy_complex64 a; |
| 1189 | if( sb1 == sizeof(npy_complex64) ) { // Aligned |
| 1190 | for( npy_intp I = 0; I < n; I++ ) { |
| 1191 | _inline_sub( Z_J, x[I], a ); |
| 1192 | _inline_add( Z_J, x[I], r[I] ); |
| 1193 | _inline_div( r[I], a, r[I] ); |
| 1194 | _inline_log( r[I], r[I] ); |
| 1195 | _inline_mul( Z_J2, r[I], r[I] ); |
| 1196 | } |
| 1197 | } |
| 1198 | else { |
| 1199 | sb1 /= sizeof(npy_complex64); |
| 1200 | for( npy_intp I = 0; I < n; I++ ) { |
| 1201 | _inline_sub( Z_J, x[I*sb1], a ); |
| 1202 | _inline_add( Z_J, x[I*sb1], r[I] ); |
| 1203 | _inline_div( r[I], a, r[I] ); |
| 1204 | _inline_log( r[I], r[I] ); |
| 1205 | _inline_mul( Z_J2, r[I], r[I] ); |
| 1206 | } |
| 1207 | } |
| 1208 | } |
| 1209 | |
| 1210 | static void |
| 1211 | nc_atan( npy_intp n, npy_complex128 *x, npy_intp sb1, npy_complex128 *r) |
nothing calls this directly
no test coverage detected
searching dependent graphs…