diff --git a/src/mat/impls/dense/seq/dense.c b/src/mat/impls/dense/seq/dense.c index b35238e..6028e9a 100644 --- a/src/mat/impls/dense/seq/dense.c +++ b/src/mat/impls/dense/seq/dense.c @@ -483,7 +483,7 @@ PetscErrorCode MatCholeskyFactor_SeqDense(Mat A,IS perm,const MatFactorInfo *fac A->ops->solvetransposeadd = MatSolveTransposeAdd_SeqDense; A->factortype = MAT_FACTOR_CHOLESKY; - ierr = PetscLogFlops((A->cmap->n*A->cmap->n*A->cmap->n)/3.0);CHKERRQ(ierr); + ierr = PetscLogFlops((1.0*A->cmap->n*A->cmap->n*A->cmap->n)/3.0);CHKERRQ(ierr); #endif PetscFunctionReturn(0); } @@ -1298,7 +1298,7 @@ PetscErrorCode MatDiagonalScale_SeqDense(Mat A,Vec ll,Vec rr) for (j=0; j *nrm) *nrm = sum; } - ierr = PetscLogFlops(A->cmap->n*A->rmap->n);CHKERRQ(ierr); + ierr = PetscLogFlops(1.0*A->cmap->n*A->rmap->n);CHKERRQ(ierr); } else if (type == NORM_INFINITY) { *nrm = 0.0; for (j=0; jrmap->n; j++) { @@ -1362,7 +1362,7 @@ PetscErrorCode MatNorm_SeqDense(Mat A,NormType type,PetscReal *nrm) } if (sum > *nrm) *nrm = sum; } - ierr = PetscLogFlops(A->cmap->n*A->rmap->n);CHKERRQ(ierr); + ierr = PetscLogFlops(1.0*A->cmap->n*A->rmap->n);CHKERRQ(ierr); } else SETERRQ(PETSC_COMM_SELF,PETSC_ERR_SUP,"No two norm"); PetscFunctionReturn(0); }