message: Fix for surface meshes
diff: @@ -810,7 +810,6 @@ void MassIntegrator::AssembleElementMatrix2(
  {
  int tr_nd = trial_fe.GetDof();
  int te_nd = test_fe.GetDof();
  - // int dim = trial_fe.GetDim();
  double w;
  #ifdef MFEM_THREAD_SAFE
  @@ -1846,18 +1845,19 @@ void VectorFEMassIntegrator::AssembleElementMatrix2(
  {
  // assume test_fe is scalar FE and trial_fe is vector FE
  int dim = test_fe.GetDim();
  + int spaceDim = Trans.GetSpaceDim();
  int trial_dof = trial_fe.GetDof();
  int test_dof = test_fe.GetDof();
  double w;
  double Kv;
  #ifdef MFEM_THREAD_SAFE
  - DenseMatrix trial_vshape(trial_dof, dim);
  + DenseMatrix trial_vshape(trial_dof, spaceDim);
  Vector shape(test_dof);
  Vector D(VQ ? VQ->GetVDim() : 0);
  DenseMatrix K(MQ ? MQ->GetVDim() : 0, MQ ? MQ->GetVDim() : 0);
  #else
  - trial_vshape.SetSize(trial_dof, dim);
  + trial_vshape.SetSize(trial_dof, spaceDim);
  shape.SetSize(test_dof);
  D.SetSize(VQ ? VQ->GetVDim() : 0);
  K.SetSize(MQ ? MQ->GetVDim() : 0, MQ ? MQ->GetVDim() : 0);
  @@ -1943,19 +1943,19 @@ void VectorFEMassIntegrator::AssembleElementMatrix2(
  && trial_fe.GetRangeType() == FiniteElement::VECTOR)
  {
  // assume both test_fe and trial_fe are vector FE
  - int dim = test_fe.GetDim();
  + int spaceDim = Trans.GetSpaceDim();
  int trial_dof = trial_fe.GetDof();
  int test_dof = test_fe.GetDof();
  double w;
  #ifdef MFEM_THREAD_SAFE
  - DenseMatrix trial_vshape(trial_dof, dim);
  - DenseMatrix test_vshape(test_dof,dim);
  + DenseMatrix trial_vshape(trial_dof,spaceDim);
  + DenseMatrix test_vshape(test_dof,spaceDim);
  Vector D(VQ ? VQ->GetVDim() : 0);
  DenseMatrix K(MQ ? MQ->GetVDim() : 0, MQ ? MQ->GetVDim() : 0);
  #else
  - trial_vshape.SetSize(trial_dof, dim);
  - test_vshape.SetSize(test_dof,dim);
  + trial_vshape.SetSize(trial_dof,spaceDim);
  + test_vshape.SetSize(test_dof,spaceDim);
  D.SetSize(VQ ? VQ->GetVDim() : 0);
  K.SetSize(MQ ? MQ->GetVDim() : 0, MQ ? MQ->GetVDim() : 0);
  #endif
  