MFEM  v4.5.1
Finite element discretization library
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Pages
eval_by_nodes.cpp
Go to the documentation of this file.
1 // Copyright (c) 2010-2022, 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 "../quadinterpolator.hpp"
13 #include "dispatch.hpp"
14 #include "eval.hpp"
15 
16 namespace mfem
17 {
18 
19 namespace internal
20 {
21 
22 namespace quadrature_interpolator
23 {
24 
25 // Tensor-product evaluation of quadrature point values: dispatch function.
26 // Instantiation for the case QVectorLayout::byNODES.
27 template<>
28 void TensorValues<QVectorLayout::byNODES>(const int NE,
29  const int vdim,
30  const DofToQuad &maps,
31  const Vector &e_vec,
32  Vector &q_val)
33 {
34  if (NE == 0) { return; }
35  const int dim = maps.FE->GetDim();
36  const int D1D = maps.ndof;
37  const int Q1D = maps.nqpt;
38  const double *B = maps.B.Read();
39  const double *X = e_vec.Read();
40  double *Y = q_val.Write();
41 
42  constexpr QVectorLayout L = QVectorLayout::byNODES;
43 
44  const int id = (vdim<<8) | (D1D<<4) | Q1D;
45 
46  if (dim == 1)
47  {
48  MFEM_VERIFY(D1D <= MAX_D1D, "Orders higher than " << MAX_D1D-1
49  << " are not supported!");
50  MFEM_VERIFY(Q1D <= MAX_Q1D, "Quadrature rules with more than "
51  << MAX_Q1D << " 1D points are not supported!");
52  Values1D<L>(NE, B, X, Y, vdim, D1D, Q1D);
53  return;
54  }
55  if (dim == 2)
56  {
57  switch (id)
58  {
59  case 0x133: return Values2D<L,1,3,3>(NE,B,X,Y);
60  case 0x124: return Values2D<L,1,2,4>(NE,B,X,Y);
61  case 0x132: return Values2D<L,1,3,2>(NE,B,X,Y);
62  case 0x134: return Values2D<L,1,3,4>(NE,B,X,Y);
63  case 0x143: return Values2D<L,1,4,3>(NE,B,X,Y);
64  case 0x144: return Values2D<L,1,4,4>(NE,B,X,Y);
65 
66  case 0x222: return Values2D<L,2,2,2>(NE,B,X,Y);
67  case 0x223: return Values2D<L,2,2,3>(NE,B,X,Y);
68  case 0x224: return Values2D<L,2,2,4>(NE,B,X,Y);
69  case 0x225: return Values2D<L,2,2,5>(NE,B,X,Y);
70  case 0x226: return Values2D<L,2,2,6>(NE,B,X,Y);
71 
72  case 0x233: return Values2D<L,2,3,3>(NE,B,X,Y);
73  case 0x234: return Values2D<L,2,3,4>(NE,B,X,Y);
74  case 0x236: return Values2D<L,2,3,6>(NE,B,X,Y);
75 
76  case 0x243: return Values2D<L,2,4,3>(NE,B,X,Y);
77  case 0x244: return Values2D<L,2,4,4>(NE,B,X,Y);
78  case 0x245: return Values2D<L,2,4,5>(NE,B,X,Y);
79  case 0x246: return Values2D<L,2,4,6>(NE,B,X,Y);
80  case 0x247: return Values2D<L,2,4,7>(NE,B,X,Y);
81 
82  case 0x256: return Values2D<L,2,5,6>(NE,B,X,Y);
83 
84  default:
85  {
86  constexpr int MD = MAX_D1D;
87  constexpr int MQ = MAX_Q1D;
88  MFEM_VERIFY(D1D <= MD, "Orders higher than " << MD-1
89  << " are not supported!");
90  MFEM_VERIFY(Q1D <= MQ, "Quadrature rules with more than "
91  << MQ << " 1D points are not supported!");
92  Values2D<L,0,0,0,0,MD,MQ>(NE,B,X,Y,vdim,D1D,Q1D);
93  return;
94  }
95  }
96  }
97  if (dim == 3)
98  {
99  switch (id)
100  {
101  case 0x124: return Values3D<L,1,2,4>(NE,B,X,Y);
102  case 0x133: return Values3D<L,1,3,3>(NE,B,X,Y);
103  case 0x134: return Values3D<L,1,3,4>(NE,B,X,Y);
104  case 0x136: return Values3D<L,1,3,6>(NE,B,X,Y);
105  case 0x143: return Values3D<L,1,4,3>(NE,B,X,Y);
106  case 0x144: return Values3D<L,1,4,4>(NE,B,X,Y);
107  case 0x148: return Values3D<L,1,4,8>(NE,B,X,Y);
108 
109  case 0x222: return Values3D<L,2,2,2>(NE,B,X,Y);
110  case 0x223: return Values3D<L,2,2,3>(NE,B,X,Y);
111  case 0x234: return Values3D<L,2,3,4>(NE,B,X,Y);
112 
113  case 0x323: return Values3D<L,3,2,3>(NE,B,X,Y);
114  case 0x324: return Values3D<L,3,2,4>(NE,B,X,Y);
115  case 0x325: return Values3D<L,3,2,5>(NE,B,X,Y);
116  case 0x326: return Values3D<L,3,2,6>(NE,B,X,Y);
117 
118  case 0x333: return Values3D<L,3,3,3>(NE,B,X,Y);
119  case 0x334: return Values3D<L,3,3,4>(NE,B,X,Y);
120  case 0x335: return Values3D<L,3,3,5>(NE,B,X,Y);
121  case 0x336: return Values3D<L,3,3,6>(NE,B,X,Y);
122 
123  case 0x343: return Values3D<L,3,4,3>(NE,B,X,Y);
124  case 0x344: return Values3D<L,3,4,4>(NE,B,X,Y);
125  case 0x346: return Values3D<L,3,4,6>(NE,B,X,Y);
126  case 0x347: return Values3D<L,3,4,7>(NE,B,X,Y);
127  case 0x348: return Values3D<L,3,4,8>(NE,B,X,Y);
128 
129  default:
130  {
131  constexpr int MD = 8;
132  constexpr int MQ = 8;
133  MFEM_VERIFY(D1D <= MD, "Orders higher than " << MD-1
134  << " are not supported!");
135  MFEM_VERIFY(Q1D <= MQ, "Quadrature rules with more than "
136  << MQ << " 1D points are not supported!");
137  Values3D<L,0,0,0,MD,MQ>(NE,B,X,Y,vdim,D1D,Q1D);
138  return;
139  }
140  }
141  }
142  mfem::out << "Unknown kernel 0x" << std::hex << id << std::endl;
143  MFEM_ABORT("Kernel not supported yet");
144 }
145 
146 } // namespace quadrature_interpolator
147 
148 } // namespace internal
149 
150 } // namespace mfem
mfem::MAX_Q1D
const int MAX_Q1D
Definition: forall.hpp:29
mfem::QVectorLayout
QVectorLayout
Type describing possible layouts for Q-vectors.
Definition: fespace.hpp:52
dim
int dim
Definition: ex24.cpp:53
mfem::MAX_D1D
const int MAX_D1D
Definition: forall.hpp:28
mfem::QVectorLayout::byNODES
NQPT x VDIM x NE (values) / NQPT x VDIM x DIM x NE (grads)
mfem::out
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
Generated on Wed Dec 7 2022 10:41:46 for MFEM by   doxygen 1.8.5