我想简单地将三个矩阵相乘：
boost::numeric::ublas::matrix<double> bTP = trans(b1);
boost::numeric::ublas::matrix<double> tmp3 = prod((M1*elem[e_].alpha_por*Kf), trans(M1));
boost::numeric::ublas::matrix<double> tmp1 = prod(bTP, tmp3);
boost::numeric::ublas::matrix<double> tmp2 = prod(tmp1, b1);
我的tmp3矩阵是对称的，但是每乘以((b1)^T * tmp3 * b1)都会得到一个非对称矩阵，而不是对称矩阵。
boost库的prod函数在处理大型双精度矩阵乘法时是否存在问题？