diff options
| author | Erik Schnetter <schnetter@gmail.com> | 2013-09-02 18:24:44 -0400 |
|---|---|---|
| committer | Erik Schnetter <schnetter@gmail.com> | 2013-09-02 22:15:56 -0400 |
| commit | 18978efbbb1236ec67c22c7a48ef699010d0e435 (patch) | |
| tree | baca3bb7852ba97a02883d00512fcdf34da327c8 | |
| parent | fad83c572e8ef9674e3fe5a217195b5934972730 (diff) | |
| download | vecmathlib-18978efbbb1236ec67c22c7a48ef699010d0e435.zip vecmathlib-18978efbbb1236ec67c22c7a48ef699010d0e435.tar.gz | |
rsqrt: Do not require strictly positive arguments
| -rw-r--r-- | mathfuncs_sqrt.h | 2 |
1 files changed, 1 insertions, 1 deletions
diff --git a/mathfuncs_sqrt.h b/mathfuncs_sqrt.h index 97307d0..a9ab172 100644 --- a/mathfuncs_sqrt.h +++ b/mathfuncs_sqrt.h @@ -73,7 +73,7 @@ namespace vecmathlib { realvec_t mathfuncs<realvec_t>::vml_rsqrt(realvec_t x) { // Initial guess - VML_ASSERT(all(x > RV(0.0))); + // VML_ASSERT(all(x > RV(0.0))); intvec_t ilogb_x = ilogb(x); realvec_t s = ifthen(convert_bool(ilogb_x & IV(I(1))), RV(R(0.583)), RV(R(0.824))); |
