Qelement_face_coordinate_translation_schemes.h
Go to the documentation of this file.
1 // LIC// ====================================================================
2 // LIC// This file forms part of oomph-lib, the object-oriented,
3 // LIC// multi-physics finite-element library, available
4 // LIC// at http://www.oomph-lib.org.
5 // LIC//
6 // LIC// Copyright (C) 2006-2022 Matthias Heil and Andrew Hazel
7 // LIC//
8 // LIC// This library is free software; you can redistribute it and/or
9 // LIC// modify it under the terms of the GNU Lesser General Public
10 // LIC// License as published by the Free Software Foundation; either
11 // LIC// version 2.1 of the License, or (at your option) any later version.
12 // LIC//
13 // LIC// This library is distributed in the hope that it will be useful,
14 // LIC// but WITHOUT ANY WARRANTY; without even the implied warranty of
15 // LIC// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
16 // LIC// Lesser General Public License for more details.
17 // LIC//
18 // LIC// You should have received a copy of the GNU Lesser General Public
19 // LIC// License along with this library; if not, write to the Free Software
20 // LIC// Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA
21 // LIC// 02110-1301 USA.
22 // LIC//
23 // LIC// The authors may be contacted at oomph-lib@maths.man.ac.uk.
24 // LIC//
25 // LIC//====================================================================
26 // Header containing functions that map face coordinates on
27 // Qelements to bulk coordinates in those elements
28 
29 // Include guards
30 #ifndef OOMPH_QELEMENT_FACE_COORDINATE_TRANSLATION_HEADER
31 #define OOMPH_QELEMENT_FACE_COORDINATE_TRANSLATION_HEADER
32 
33 // Config header generated by autoconfig
34 #ifdef HAVE_CONFIG_H
35 #include <oomph-lib-config.h>
36 #endif
37 
38 #ifdef OOMPH_HAS_MPI
39 #include "mpi.h"
40 #endif
41 
42 // oomph-lib headers
43 #include "Vector.h"
44 #include "matrices.h"
45 
46 namespace oomph
47 {
48  //=============================================================
51  //==============================================================
52  namespace QElement1FaceToBulkCoordinates
53  {
55  void face0(const Vector<double>& s, Vector<double>& s_bulk);
56 
58  void face1(const Vector<double>& s, Vector<double>& s_bulk);
59  } // namespace QElement1FaceToBulkCoordinates
60 
61 
62  //=============================================================
66  //=============================================================
67  namespace QElement1BulkCoordinateDerivatives
68  {
70  void faces0(const Vector<double>& s,
71  DenseMatrix<double>& dsbulk_dsface,
72  unsigned& interior_direction);
73  } // namespace QElement1BulkCoordinateDerivatives
74 
75 
76  //===================================================================
79  //==================================================================
80  namespace QElement2FaceToBulkCoordinates
81  {
83  void face0(const Vector<double>& s, Vector<double>& s_bulk);
84 
86  void face1(const Vector<double>& s, Vector<double>& s_bulk);
87 
89  void face2(const Vector<double>& s, Vector<double>& s_bulk);
90 
92  void face3(const Vector<double>& s, Vector<double>& s_bulk);
93  } // namespace QElement2FaceToBulkCoordinates
94 
95 
96  //=============================================================
100  //=============================================================
101  namespace QElement2BulkCoordinateDerivatives
102  {
104  void faces0(const Vector<double>& s,
105  DenseMatrix<double>& dsbulk_dsface,
106  unsigned& interior_direction);
107 
109  void faces1(const Vector<double>& s,
110  DenseMatrix<double>& dsbulk_dsface,
111  unsigned& interior_direction);
112  } // namespace QElement2BulkCoordinateDerivatives
113 
114 
115  //===================================================================
118  //==================================================================
119  namespace QElement3FaceToBulkCoordinates
120  {
122  void face0(const Vector<double>& s, Vector<double>& s_bulk);
123 
125  void face1(const Vector<double>& s, Vector<double>& s_bulk);
126 
128  void face2(const Vector<double>& s, Vector<double>& s_bulk);
129 
131  void face3(const Vector<double>& s, Vector<double>& s_bulk);
132 
134  void face4(const Vector<double>& s, Vector<double>& s_bulk);
135 
137  void face5(const Vector<double>& s, Vector<double>& s_bulk);
138  } // namespace QElement3FaceToBulkCoordinates
139 
140 
141  //=============================================================
145  //=============================================================
146  namespace QElement3BulkCoordinateDerivatives
147  {
149  void faces0(const Vector<double>& s,
150  DenseMatrix<double>& dsbulk_dsface,
151  unsigned& interior_direction);
152 
154  void faces1(const Vector<double>& s,
155  DenseMatrix<double>& dsbulk_dsface,
156  unsigned& interior_direction);
157 
159  void faces2(const Vector<double>& s,
160  DenseMatrix<double>& dsbulk_dsface,
161  unsigned& interior_direction);
162  } // namespace QElement3BulkCoordinateDerivatives
163 
164 } // namespace oomph
165 
166 #endif
Matrix< Scalar, Dynamic, Dynamic > DenseMatrix
Definition: BenchSparseUtil.h:23
RealScalar s
Definition: level1_cplx_impl.h:130
void faces0(const Vector< double > &s, DenseMatrix< double > &dsbulk_dsface, unsigned &interior_direction)
Function for both faces – the bulk coordinate is fixed on both.
Definition: Qelement_face_coordinate_translation_schemes.cc:59
void face1(const Vector< double > &s, Vector< double > &s_bulk)
The translation scheme for the face s0 = 1.0.
Definition: Qelement_face_coordinate_translation_schemes.cc:44
void face0(const Vector< double > &s, Vector< double > &s_bulk)
The translation scheme for the face s0 = -1.0.
Definition: Qelement_face_coordinate_translation_schemes.cc:38
void faces1(const Vector< double > &s, DenseMatrix< double > &dsbulk_dsface, unsigned &interior_direction)
Function for the north and south faces, along which s1 is fixed.
Definition: Qelement_face_coordinate_translation_schemes.cc:131
void faces0(const Vector< double > &s, DenseMatrix< double > &dsbulk_dsface, unsigned &interior_direction)
Function for the east and west faces, along which s0 is fixed.
Definition: Qelement_face_coordinate_translation_schemes.cc:116
void face2(const Vector< double > &s, Vector< double > &s_bulk)
The translation scheme for the east face (s0 = 1.0)
Definition: Qelement_face_coordinate_translation_schemes.cc:93
void face1(const Vector< double > &s, Vector< double > &s_bulk)
The translation scheme for the south face (s1 = -1.0)
Definition: Qelement_face_coordinate_translation_schemes.cc:86
void face0(const Vector< double > &s, Vector< double > &s_bulk)
The translation scheme for the west face (s0 = -1.0)
Definition: Qelement_face_coordinate_translation_schemes.cc:79
void face3(const Vector< double > &s, Vector< double > &s_bulk)
The translation scheme for the north face (s1 = 1.0)
Definition: Qelement_face_coordinate_translation_schemes.cc:100
void faces2(const Vector< double > &s, DenseMatrix< double > &dsbulk_dsface, unsigned &interior_direction)
Function for the left and right faces, along which s2 is fixed.
Definition: Qelement_face_coordinate_translation_schemes.cc:249
void faces0(const Vector< double > &s, DenseMatrix< double > &dsbulk_dsface, unsigned &interior_direction)
Function for the back and front faces, along which s0 is fixed.
Definition: Qelement_face_coordinate_translation_schemes.cc:210
void faces1(const Vector< double > &s, DenseMatrix< double > &dsbulk_dsface, unsigned &interior_direction)
Function for the up and down faces, along which s1 is fixed.
Definition: Qelement_face_coordinate_translation_schemes.cc:230
void face4(const Vector< double > &s, Vector< double > &s_bulk)
The translation scheme for the up face (s1 = 1.0)
Definition: Qelement_face_coordinate_translation_schemes.cc:185
void face1(const Vector< double > &s, Vector< double > &s_bulk)
The translation scheme for the down face (s1 = -1.0)
Definition: Qelement_face_coordinate_translation_schemes.cc:161
void face2(const Vector< double > &s, Vector< double > &s_bulk)
The translation scheme for the back face (s2 = -1.0)
Definition: Qelement_face_coordinate_translation_schemes.cc:169
void face0(const Vector< double > &s, Vector< double > &s_bulk)
The translation scheme for the left face (s0 = -1.0)
Definition: Qelement_face_coordinate_translation_schemes.cc:153
void face5(const Vector< double > &s, Vector< double > &s_bulk)
The translation scheme for the front face (s2 = 1.0)
Definition: Qelement_face_coordinate_translation_schemes.cc:193
void face3(const Vector< double > &s, Vector< double > &s_bulk)
The translation scheme for the right face (s0 = 1.0)
Definition: Qelement_face_coordinate_translation_schemes.cc:177
DRAIG: Change all instances of (SPATIAL_DIM) to (DIM-1).
Definition: AnisotropicHookean.h:10