void matrixProduct(const double *A, int n1, int p1, const double *B, int n2, int p2, double *C)
{
void matrixProduct(const double *A, int n1, int p1, const double *B, int n2, int p2, double *C)
{
+ if(p1!=n2)
+ {
+ std::ostringstream oss; oss << "matrixProduct : the size of input matrix are not coherent the nb of cols of input matrix #0 is " << p1 << " whereas the number of rows of input matrix #1 is " << n2 << " !";