GitHub-Proxy/eigen
1
0
Fork 0
mirror of https://gitlab.com/libeigen/eigen.git synced 2025-07-23 05:14:26 +08:00
Code Issues Packages Projects Releases Wiki Activity

Fix unit test for SIMD engine not supporting sqrt

Browse Source
This commit is contained in:
Gael Guennebaud 2018-04-26 10:47:39 +02:00
parent 8810baaed4
commit 097dd4616d
1 changed files with 31 additions and 25 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

56
test/packetmath.cpp
View File

@ -451,35 +451,41 @@ template<typename Scalar> void packetmath_real()
{ {
data1[0] = std::numeric_limits<Scalar>::quiet_NaN(); data1[0] = std::numeric_limits<Scalar>::quiet_NaN();
data1[1] = std::numeric_limits<Scalar>::epsilon(); data1[1] = std::numeric_limits<Scalar>::epsilon();
packet_helper<PacketTraits::HasLog,Packet> h; {
h.store(data2, internal::plog(h.load(data1))); packet_helper<PacketTraits::HasLog,Packet> h;
VERIFY((numext::isnan)(data2[0])); h.store(data2, internal::plog(h.load(data1)));
VERIFY_IS_EQUAL(std::log(std::numeric_limits<Scalar>::epsilon()), data2[1]); VERIFY((numext::isnan)(data2[0]));
VERIFY_IS_EQUAL(std::log(std::numeric_limits<Scalar>::epsilon()), data2[1]);
data1[0] = -std::numeric_limits<Scalar>::epsilon(); data1[0] = -std::numeric_limits<Scalar>::epsilon();
data1[1] = 0; data1[1] = 0;
h.store(data2, internal::plog(h.load(data1))); h.store(data2, internal::plog(h.load(data1)));
VERIFY((numext::isnan)(data2[0])); VERIFY((numext::isnan)(data2[0]));
VERIFY_IS_EQUAL(std::log(Scalar(0)), data2[1]); VERIFY_IS_EQUAL(std::log(Scalar(0)), data2[1]);
data1[0] = (std::numeric_limits<Scalar>::min)(); data1[0] = (std::numeric_limits<Scalar>::min)();
data1[1] = -(std::numeric_limits<Scalar>::min)(); data1[1] = -(std::numeric_limits<Scalar>::min)();
h.store(data2, internal::plog(h.load(data1))); h.store(data2, internal::plog(h.load(data1)));
VERIFY_IS_EQUAL(std::log((std::numeric_limits<Scalar>::min)()), data2[0]); VERIFY_IS_EQUAL(std::log((std::numeric_limits<Scalar>::min)()), data2[0]);
VERIFY((numext::isnan)(data2[1])); VERIFY((numext::isnan)(data2[1]));
data1[0] = std::numeric_limits<Scalar>::denorm_min(); data1[0] = std::numeric_limits<Scalar>::denorm_min();
data1[1] = -std::numeric_limits<Scalar>::denorm_min(); data1[1] = -std::numeric_limits<Scalar>::denorm_min();
h.store(data2, internal::plog(h.load(data1))); h.store(data2, internal::plog(h.load(data1)));
// VERIFY_IS_EQUAL(std::log(std::numeric_limits<Scalar>::denorm_min()), data2[0]); // VERIFY_IS_EQUAL(std::log(std::numeric_limits<Scalar>::denorm_min()), data2[0]);
VERIFY((numext::isnan)(data2[1])); VERIFY((numext::isnan)(data2[1]));
data1[0] = Scalar(-1.0f); data1[0] = Scalar(-1.0f);
h.store(data2, internal::plog(h.load(data1))); h.store(data2, internal::plog(h.load(data1)));
VERIFY((numext::isnan)(data2[0])); VERIFY((numext::isnan)(data2[0]));
h.store(data2, internal::psqrt(h.load(data1))); }
VERIFY((numext::isnan)(data2[0])); {
VERIFY((numext::isnan)(data2[1])); packet_helper<PacketTraits::HasSqrt,Packet> h;
data1[0] = Scalar(-1.0f);
h.store(data2, internal::psqrt(h.load(data1)));
VERIFY((numext::isnan)(data2[0]));
VERIFY((numext::isnan)(data2[1]));
}
} }
} }