message: Fixing with_coef logic in flux calculation.
diff: @@ -1012,9 +1012,11 @@ void DiffusionIntegrator::ComputeElementFlux
  CalcInverse(Trans.Jacobian(), invdfdx);
  invdfdx.MultTranspose(vec, vecdxt);
  + if (with_coef)
  + {
  if (!MQ && !VQ)
  {
  - if (Q && with_coef)
  + if (Q)
  {
  vecdxt *= Q->Eval(Trans,ip);
  }
  @@ -1037,7 +1039,6 @@ void DiffusionIntegrator::ComputeElementFlux
  {
  pointflux[j] = D[j] * vecdxt[j];
  }
  -
  }
  for (j = 0; j < spaceDim; j++)
  {
  @@ -1045,6 +1046,14 @@ void DiffusionIntegrator::ComputeElementFlux
  }
  }
  }
  + else
  + {
  + for (j = 0; j < spaceDim; j++)
  + {
  + flux(fnd*j+i) = vecdxt(j);
  + }
  + }
  + }
  }
  double DiffusionIntegrator::ComputeFluxEnergy
  