message: Adding checks which ensure the integration loop can succeed
diff: @@ -659,6 +659,9 @@ void MixedScalarVectorIntegrator::AssembleElementMatrix2(
  MFEM_ASSERT(this->VerifyFiniteElementTypes(trial_fe, test_fe),
  this->FiniteElementTypeFailureMessage());
  + MFEM_VERIFY(VQ, "MixedScalarVectorIntegrator: "
  + "VectorCoefficient must be set");
  +
  const FiniteElement * vec_fe = transpose?&trial_fe:&test_fe;
  const FiniteElement * sca_fe = transpose?&test_fe:&trial_fe;
  @@ -669,13 +672,17 @@ void MixedScalarVectorIntegrator::AssembleElementMatrix2(
  int vdim = GetVDim(*vec_fe);
  double vtmp;
  + MFEM_VERIFY(VQ->GetVDim() == vdim, "MixedScalarVectorIntegrator: "
  + "Dimensions of VectorCoefficient and Vector-valued basis "
  + "functions must match");
  +
  #ifdef MFEM_THREAD_SAFE
  - Vector V(VQ ? VQ->GetVDim() : 0);
  + Vector V(vdim);
  DenseMatrix vshape(vec_nd, vdim);
  Vector shape(sca_nd);
  Vector vshape_tmp(vec_nd);
  #else
  - V.SetSize(VQ ? VQ->GetVDim() : 0);
  + V.SetSize(vdim);
  vshape.SetSize(vec_nd, vdim);
  shape.SetSize(sca_nd);
  vshape_tmp.SetSize(vec_nd);
  