message: Fixing some unsupported coefficients in ComputeFluxEnergy.
diff: @@ -984,6 +984,8 @@ void DiffusionIntegrator::ComputeElementFlux
  "Unexpected height for MatrixCoefficient");
  }
  + MFEM_VERIFY(!SMQ, "SymmetricMatrixCoefficient not supported here");
  +
  #ifdef MFEM_THREAD_SAFE
  DenseMatrix dshape(nd,dim), invdfdx(dim, spaceDim);
  DenseMatrix M(MQ ? spaceDim : 0);
  @@ -1066,8 +1068,13 @@ double DiffusionIntegrator::ComputeFluxEnergy
  #ifdef MFEM_THREAD_SAFE
  DenseMatrix M;
  + Vector D(VQ ? VQ->GetVDim() : 0);
  +#else
  + D.SetSize(VQ ? VQ->GetVDim() : 0);
  #endif
  + MFEM_VERIFY(!SMQ, "SymmetricMatrixCoefficient not supported here");
  +
  shape.SetSize(nd);
  pointflux.SetSize(spaceDim);
  if (d_energy) { vec.SetSize(spaceDim); }
  @@ -1096,16 +1103,22 @@ double DiffusionIntegrator::ComputeFluxEnergy
  Trans.SetIntPoint(&ip);
  double w = Trans.Weight() * ip.weight;
  - if (!MQ)
  + if (MQ)
  {
  - double e = (pointflux * pointflux);
  - if (Q) { e *= Q->Eval(Trans, ip); }
  - energy += w * e;
  + MQ->Eval(M, Trans, ip);
  + energy += w * M.InnerProduct(pointflux, pointflux);
  + }
  + else if (VQ)
  + {
  + VQ->Eval(D, Trans, ip);
  + D *= pointflux;
  + energy += w * (D * pointflux);
  }
  else
  {
  - MQ->Eval(M, Trans, ip);
  - energy += w * M.InnerProduct(pointflux, pointflux);
  + double e = (pointflux * pointflux);
  + if (Q) { e *= Q->Eval(Trans, ip); }
  + energy += w * e;
  }
  if (d_energy)
  