| 1005 | } |
| 1006 | |
| 1007 | static void |
| 1008 | nc_acosh( npy_intp n, npy_complex64 *x, npy_intp sb1, npy_complex64 *r) |
| 1009 | { |
| 1010 | npy_complex64 t, a; |
| 1011 | if( sb1 == sizeof(npy_complex64) ) { // Aligned |
| 1012 | for( npy_intp I = 0; I < n; I++ ) { |
| 1013 | a = x[I]; |
| 1014 | _inline_add( Z_1, x[I], t ); |
| 1015 | _inline_sqrt( t, t ); |
| 1016 | _inline_sub( x[I], Z_1, r[I] ); |
| 1017 | _inline_sqrt( r[I] , r[I] ); |
| 1018 | _inline_mul( t, r[I], r[I] ); |
| 1019 | _inline_add( a, r[I], r[I] ); |
| 1020 | _inline_log( r[I], r[I] ); |
| 1021 | } |
| 1022 | } |
| 1023 | else { |
| 1024 | sb1 /= sizeof(npy_complex64); |
| 1025 | for( npy_intp I = 0; I < n; I++ ) { |
| 1026 | a = x[I*sb1]; |
| 1027 | _inline_add( Z_1, x[I*sb1], t ); |
| 1028 | _inline_sqrt( t, t ); |
| 1029 | _inline_sub( x[I*sb1], Z_1, r[I] ); |
| 1030 | _inline_sqrt( r[I] , r[I] ); |
| 1031 | _inline_mul( t, r[I], r[I] ); |
| 1032 | _inline_add( a, r[I], r[I] ); |
| 1033 | _inline_log( r[I], r[I] ); |
| 1034 | } |
| 1035 | } |
| 1036 | } |
| 1037 | |
| 1038 | static void |
| 1039 | nc_acosh( npy_intp n, npy_complex128 *x, npy_intp sb1, npy_complex128 *r) |
nothing calls this directly
no test coverage detected
searching dependent graphs…