message: enable the use of custom allocators for dg trace integrator
diff: @@ -136,6 +136,9 @@ static void PADGTraceSetup(const int dim,
  void DGTraceIntegrator::SetupPA(const FiniteElementSpace &fes, FaceType type)
  {
  + const MemoryType mt = (pa_mt == MemoryType::DEFAULT) ?
  + Device::GetDeviceMemoryType() : pa_mt;
  +
  nf = fes.GetNFbyType(type);
  if (nf==0) { return; }
  // Assumes tensor-product elements
  @@ -153,7 +156,7 @@ void DGTraceIntegrator::SetupPA(const FiniteElementSpace &fes, FaceType type)
  geom = mesh->GetFaceGeometricFactors(
  *ir,
  FaceGeometricFactors::DETERMINANTS |
  - FaceGeometricFactors::NORMALS, type);
  + FaceGeometricFactors::NORMALS, type, mt);
  maps = &el.GetDofToQuad(*ir, DofToQuad::TENSOR);
  dofs1D = maps->ndof;
  quad1D = maps->nqpt;
  