diff options
Diffstat (limited to 'mathfuncs_log.h')
| -rw-r--r-- | mathfuncs_log.h | 46 |
1 files changed, 23 insertions, 23 deletions
diff --git a/mathfuncs_log.h b/mathfuncs_log.h index 0a40eec..53b5f94 100644 --- a/mathfuncs_log.h +++ b/mathfuncs_log.h @@ -34,41 +34,41 @@ namespace vecmathlib { #ifdef VML_HAVE_FP_CONTRACT // float, error=5.98355642684398209498469870525e-9 r = RV(0.410981538282433293325329456838); - r = fma(r, y2, RV(0.402155483172044562892705980539)); - r = fma(r, y2, RV(0.57755014627178237959721643293)); - r = fma(r, y2, RV(0.96178780600659929206930296869)); - r = fma(r, y2, RV(2.88539012786343587248965772685)); + r = mad(r, y2, RV(0.402155483172044562892705980539)); + r = mad(r, y2, RV(0.57755014627178237959721643293)); + r = mad(r, y2, RV(0.96178780600659929206930296869)); + r = mad(r, y2, RV(2.88539012786343587248965772685)); #else //flaot, error=2.25468184051947656525068987795e-7 r = RV(0.498866687070343238590910977481); - r = fma(r, y2, RV(0.57002741193682764193895550312)); - r = fma(r, y2, RV(0.96200215034262628756932169194)); - r = fma(r, y2, RV(2.88538850388042106595516956395)); + r = mad(r, y2, RV(0.57002741193682764193895550312)); + r = mad(r, y2, RV(0.96200215034262628756932169194)); + r = mad(r, y2, RV(2.88538850388042106595516956395)); #endif break; case 8: #ifdef VML_HAVE_FP_CONTRACT // double, error=9.45037202901655672811489051683e-17 r = RV(0.259935726478127940817401224248); - r = fma(r, y2, RV(0.140676370079882918464564658472)); - r = fma(r, y2, RV(0.196513478841924000569879320851)); - r = fma(r, y2, RV(0.221596471338300882039273355617)); - r = fma(r, y2, RV(0.262327298560598641020007602127)); - r = fma(r, y2, RV(0.320598261015170101859472461613)); - r = fma(r, y2, RV(0.412198595799726905825871956187)); - r = fma(r, y2, RV(0.57707801621733949207376840932)); - r = fma(r, y2, RV(0.96179669392666302667713134701)); - r = fma(r, y2, RV(2.88539008177792581277410991327)); + r = mad(r, y2, RV(0.140676370079882918464564658472)); + r = mad(r, y2, RV(0.196513478841924000569879320851)); + r = mad(r, y2, RV(0.221596471338300882039273355617)); + r = mad(r, y2, RV(0.262327298560598641020007602127)); + r = mad(r, y2, RV(0.320598261015170101859472461613)); + r = mad(r, y2, RV(0.412198595799726905825871956187)); + r = mad(r, y2, RV(0.57707801621733949207376840932)); + r = mad(r, y2, RV(0.96179669392666302667713134701)); + r = mad(r, y2, RV(2.88539008177792581277410991327)); #else // double, error=1.21820548287702216975532695788e-13 r = RV(0.293251364683280430617251942017); - r = fma(r, y2, RV(0.201364223624519571276587631354)); - r = fma(r, y2, RV(0.264443947645547871780098560836)); - r = fma(r, y2, RV(0.320475051320227723946459855458)); - r = fma(r, y2, RV(0.412202612052105347480086431555)); - r = fma(r, y2, RV(0.57707794741938820005328259256)); - r = fma(r, y2, RV(0.96179669445173881282808321929)); - r = fma(r, y2, RV(2.88539008177676567117601117274)); + r = mad(r, y2, RV(0.201364223624519571276587631354)); + r = mad(r, y2, RV(0.264443947645547871780098560836)); + r = mad(r, y2, RV(0.320475051320227723946459855458)); + r = mad(r, y2, RV(0.412202612052105347480086431555)); + r = mad(r, y2, RV(0.57707794741938820005328259256)); + r = mad(r, y2, RV(0.96179669445173881282808321929)); + r = mad(r, y2, RV(2.88539008177676567117601117274)); #endif break; default: |
