| 1067 | } |
| 1068 | |
| 1069 | static void |
| 1070 | nc_asin( npy_intp n, npy_complex64 *x, npy_intp sb1, npy_complex64 *r) |
| 1071 | { |
| 1072 | npy_complex64 a; |
| 1073 | if( sb1 == sizeof(npy_complex64) ) { // Aligned |
| 1074 | for( npy_intp I = 0; I < n; I++ ) { |
| 1075 | _inline_muli( x[I], a ); |
| 1076 | _inline_mul( x[I], x[I], r[I] ); |
| 1077 | _inline_sub( Z_1, r[I], r[I] ); |
| 1078 | _inline_sqrt( r[I], r[I] ); |
| 1079 | _inline_add( a, r[I], r[I] ); |
| 1080 | _inline_log( r[I], r[I] ); |
| 1081 | _inline_muli( r[I], r[I] ); |
| 1082 | _inline_neg( r[I], r[I] ); |
| 1083 | } |
| 1084 | } |
| 1085 | else { |
| 1086 | sb1 /= sizeof(npy_complex64); |
| 1087 | for( npy_intp I = 0; I < n; I++ ) { |
| 1088 | _inline_muli( x[I*sb1], a ); |
| 1089 | _inline_mul( x[I*sb1], x[I*sb1], r[I] ); |
| 1090 | _inline_sub( Z_1, r[I], r[I] ); |
| 1091 | _inline_sqrt( r[I], r[I] ); |
| 1092 | _inline_add( a, r[I], r[I] ); |
| 1093 | _inline_log( r[I], r[I] ); |
| 1094 | _inline_muli( r[I], r[I] ); |
| 1095 | _inline_neg( r[I], r[I] ); |
| 1096 | } |
| 1097 | } |
| 1098 | } |
| 1099 | |
| 1100 | static void |
| 1101 | nc_asin( npy_intp n, npy_complex128 *x, npy_intp sb1, npy_complex128 *r) |
nothing calls this directly
no test coverage detected
searching dependent graphs…