MFEM  v4.6.0
Finite element discretization library
grad_phys_by_vdim.cpp
Go to the documentation of this file.
1 // Copyright (c) 2010-2023, Lawrence Livermore National Security, LLC. Produced
2 // at the Lawrence Livermore National Laboratory. All Rights reserved. See files
3 // LICENSE and NOTICE for details. LLNL-CODE-806117.
4 //
5 // This file is part of the MFEM library. For more information and source code
6 // availability visit https://mfem.org.
7 //
8 // MFEM is free software; you can redistribute it and/or modify it under the
9 // terms of the BSD-3 license. We welcome feedback and contributions, see file
10 // CONTRIBUTING.md for details.
11 
12 #include "dispatch.hpp"
13 #include "grad.hpp"
14 
15 namespace mfem
16 {
17 
18 namespace internal
19 {
20 
21 namespace quadrature_interpolator
22 {
23 
24 // Tensor-product evaluation of quadrature point physical derivatives: dispatch
25 // function.
26 // Instantiation for the case QVectorLayout::byVDIM.
27 template<>
28 void TensorPhysDerivatives<QVectorLayout::byVDIM>(const int NE,
29  const int vdim,
30  const DofToQuad &maps,
31  const GeometricFactors &geom,
32  const Vector &e_vec,
33  Vector &q_der)
34 {
35  if (NE == 0) { return; }
36  const int dim = maps.FE->GetDim();
37  const int D1D = maps.ndof;
38  const int Q1D = maps.nqpt;
39 
40  const int sdim = geom.mesh->SpaceDimension();
41 
42  const double *B = maps.B.Read();
43  const double *G = maps.G.Read();
44  const double *J = geom.J.Read();
45  const double *X = e_vec.Read();
46  double *Y = q_der.Write();
47 
49  constexpr bool P = true; // GRAD_PHYS
50 
51  const int id = (vdim<<8) | (D1D<<4) | Q1D;
52 
53  if (dim == 1)
54  {
55  return Derivatives1D<L,P>(NE,G,J,X,Y,sdim,vdim,D1D,Q1D);
56  }
57  if (dim == 2)
58  {
59  switch (id)
60  {
61  case 0x134: return Derivatives2D<L,P,1,3,4,8>(NE,B,G,J,X,Y,sdim);
62  case 0x146: return Derivatives2D<L,P,1,4,6,4>(NE,B,G,J,X,Y,sdim);
63  case 0x158: return Derivatives2D<L,P,1,5,8,2>(NE,B,G,J,X,Y,sdim);
64 
65  case 0x233: return Derivatives2D<L,P,2,3,3,8>(NE,B,G,J,X,Y,sdim);
66  case 0x234: return Derivatives2D<L,P,2,3,4,8>(NE,B,G,J,X,Y,sdim);
67  case 0x246: return Derivatives2D<L,P,2,4,6,4>(NE,B,G,J,X,Y,sdim);
68  case 0x258: return Derivatives2D<L,P,2,5,8,2>(NE,B,G,J,X,Y,sdim);
69  default:
70  {
71  const int MD = DeviceDofQuadLimits::Get().MAX_D1D;
72  const int MQ = DeviceDofQuadLimits::Get().MAX_Q1D;
73  MFEM_VERIFY(D1D <= MD, "Orders higher than " << MD-1
74  << " are not supported!");
75  MFEM_VERIFY(Q1D <= MQ, "Quadrature rules with more than "
76  << MQ << " 1D points are not supported!");
77  Derivatives2D<L,P>(NE,B,G,J,X,Y,sdim,vdim,D1D,Q1D);
78  return;
79  }
80  }
81  }
82  if (dim == 3)
83  {
84  switch (id)
85  {
86  case 0x134: return Derivatives3D<L,P,1,3,4>(NE,B,G,J,X,Y);
87  case 0x146: return Derivatives3D<L,P,1,4,6>(NE,B,G,J,X,Y);
88  case 0x158: return Derivatives3D<L,P,1,5,8>(NE,B,G,J,X,Y);
89 
90  case 0x334: return Derivatives3D<L,P,3,3,4>(NE,B,G,J,X,Y);
91  case 0x346: return Derivatives3D<L,P,3,4,6>(NE,B,G,J,X,Y);
92  case 0x358: return Derivatives3D<L,P,3,5,8>(NE,B,G,J,X,Y);
93  default:
94  {
95  const int MD = DeviceDofQuadLimits::Get().MAX_INTERP_1D;
96  const int MQ = DeviceDofQuadLimits::Get().MAX_INTERP_1D;
97  MFEM_VERIFY(D1D <= MD, "Orders higher than " << MD-1
98  << " are not supported!");
99  MFEM_VERIFY(Q1D <= MQ, "Quadrature rules with more than "
100  << MQ << " 1D points are not supported!");
101  Derivatives3D<L,P>(NE,B,G,J,X,Y,vdim,D1D,Q1D);
102  return;
103  }
104  }
105  }
106  mfem::out << "Unknown kernel 0x" << std::hex << id << std::endl;
107  MFEM_ABORT("Unknown kernel");
108 }
109 
110 } // namespace quadrature_interpolator
111 
112 } // namespace internal
113 
114 } // namespace mfem
int MAX_Q1D
Maximum number of 1D quadrature points.
Definition: forall.hpp:113
int MAX_D1D
Maximum number of 1D nodal points.
Definition: forall.hpp:112
OutStream out(std::cout)
Global stream used by the library for standard output. Initially it uses the same std::streambuf as s...
Definition: globals.hpp:66
static const DeviceDofQuadLimits & Get()
Return a const reference to the DeviceDofQuadLimits singleton.
Definition: forall.hpp:122
QVectorLayout
Type describing possible layouts for Q-vectors.
Definition: fespace.hpp:52
int dim
Definition: ex24.cpp:53
int MAX_INTERP_1D
Maximum number of points for use in QuadratureInterpolator.
Definition: forall.hpp:118
VDIM x NQPT x NE (values) / VDIM x DIM x NQPT x NE (grads)