#include <Eigen/src/Core/util/BlasUtil.h>
combine_scalar_factors extracts and multiplies factors from GEMM and GEMV products. There is a specialization for booleans
combine_scalar_factors