adds C++ function for FMS

dune/codegen/common/muladd_workarounds.hh

@@ -13,4 +13,10 @@ inline T mul_add(T op1, T op2, T op3)
   return op1 * op2 + op3;
 }
 
+template<typename T>
+inline T mul_sub(T op1, T op2, T op3)
+{
+  return op1 * op2 - op3;
+}
+
 #endif