mirror/eigen
2
0
Fork 0
mirror of https://gitlab.com/libeigen/eigen.git synced 2024-12-27 07:29:52 +08:00
Code Issues Packages Projects Releases Wiki Activity

Fix implicit conversion warnings and use pnegate to negate packets

Browse Source
This commit is contained in:
Christoph Hertzberg 2019-09-23 16:07:43 +02:00
parent ba0736fa8e
commit e4c1b3c1d2
1 changed files with 10 additions and 10 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

20
unsupported/Eigen/src/SpecialFunctions/BesselFunctionsImpl.h
View File

@ -311,7 +311,7 @@ struct generic_i1e<T, float> {
// [-8, 8] and evaluate a branch based off of that. It's possible
// in practice most elements are in this region.
y = pselect(pcmp_le(y, pset1<T>(8.0f)), y_le_eight, y_gt_eight);
return pselect(pcmp_lt(x, pset1<T>(0.0f)), -y, y);
return pselect(pcmp_lt(x, pset1<T>(0.0f)), pnegate(y), y);
}
};
@ -392,7 +392,7 @@ struct generic_i1e<T, double> {
// [-8, 8] and evaluate a branch based off of that. It's possible
// in practice most elements are in this region.
y = pselect(pcmp_le(y, pset1<T>(8.0)), y_le_eight, y_gt_eight);
return pselect(pcmp_lt(x, pset1<T>(0.0f)), -y, y);
return pselect(pcmp_lt(x, pset1<T>(0.0)), pnegate(y), y);
}
};
@ -491,8 +491,8 @@ struct generic_k0e<T, float> {
T x_le_two = internal::pchebevl<T, 7>::run(
pmadd(x, x, pset1<T>(-2.0)), A);
x_le_two = pmadd(
generic_i0<T, float>::run(x), pmul(
pset1<T>(-1.0), plog(pmul(pset1<T>(0.5), x))), x_le_two);
generic_i0<T, float>::run(x), pnegate(
plog(pmul(pset1<T>(0.5), x))), x_le_two);
x_le_two = pmul(pexp(x), x_le_two);
T x_gt_two = pmul(
internal::pchebevl<T, 10>::run(
@ -663,8 +663,8 @@ struct generic_k0<T, float> {
T x_le_two = internal::pchebevl<T, 7>::run(
pmadd(x, x, pset1<T>(-2.0)), A);
x_le_two = pmadd(
generic_i0<T, float>::run(x), pmul(
pset1<T>(-1.0), plog(pmul(pset1<T>(0.5), x))), x_le_two);
generic_i0<T, float>::run(x), pnegate(
plog(pmul(pset1<T>(0.5), x))), x_le_two);
x_le_two = pselect(pcmp_le(x, pset1<T>(0.0)), MAXNUM, x_le_two);
T x_gt_two = pmul(
pmul(
@ -741,8 +741,8 @@ struct generic_k0<T, double> {
T x_le_two = internal::pchebevl<T, 10>::run(
pmadd(x, x, pset1<T>(-2.0)), A);
x_le_two = pmadd(
generic_i0<T, double>::run(x), pmul(
pset1<T>(-1.0), plog(pmul(pset1<T>(0.5), x))), x_le_two);
generic_i0<T, double>::run(x), pnegate(
plog(pmul(pset1<T>(0.5), x))), x_le_two);
x_le_two = pselect(pcmp_le(x, pset1<T>(0.0)), MAXNUM, x_le_two);
T x_gt_two = pmul(
pmul(
@ -1572,7 +1572,7 @@ struct generic_j1<T, float> {
// j1 is an odd function. This implementation differs from cephes to
// take this fact in to account. Cephes returns -j1(x) for y > 2 range.
y_gt_two = pselect(
pcmp_lt(x, pset1<T>(0.0f)), pmul(pset1<T>(-1.0f), y_gt_two), y_gt_two);
pcmp_lt(x, pset1<T>(0.0f)), pnegate(y_gt_two), y_gt_two);
return pselect(pcmp_le(y, pset1<T>(2.0f)), y_le_two, y_gt_two);
}
};
@ -1660,7 +1660,7 @@ struct generic_j1<T, double> {
// j1 is an odd function. This implementation differs from cephes to
// take this fact in to account. Cephes returns -j1(x) for y > 5 range.
y_gt_five = pselect(
pcmp_lt(x, pset1<T>(0.0f)), pmul(pset1<T>(-1.0), y_gt_five), y_gt_five);
pcmp_lt(x, pset1<T>(0.0)), pnegate(y_gt_five), y_gt_five);
return pselect(pcmp_le(y, pset1<T>(5.0)), y_le_five, y_gt_five);
}
};