Introduce mad(), as low-accuracy version of fma()

1 files changed, 8 insertions, 0 deletions

diff --git a/vec_base.h b/vec_base.h
index 3db2e24..702a405 100644
--- a/vec_base.h
+++ b/vec_base.h
@@ -519,6 +519,14 @@ namespace vecmathlib {
   }
   
   template<typename real_t, int size>
+  inline realvec<real_t, size> mad(realvec<real_t, size> x,
+                                   realvec<real_t, size> y,
+                                   realvec<real_t, size> z)
+  {
+    return x.mad(y, z);
+  }
+  
+  template<typename real_t, int size>
   inline realvec<real_t, size> nextafter(realvec<real_t, size> x,
                                          realvec<real_t, size> y)
   {