FEDEM Solver  R8.0
Source code of the dynamics solver
FFa3DLocation.H
Go to the documentation of this file.
1 // SPDX-FileCopyrightText: 2023 SAP SE
2 //
3 // SPDX-License-Identifier: Apache-2.0
4 //
5 // This file is part of FEDEM - https://openfedem.org
7 
8 #ifndef FFA_3D_LOCATION_H
9 #define FFA_3D_LOCATION_H
10 
13 #include <array>
14 
15 class FaMat33;
16 class FaMat34;
17 
18 
20 {
21 public:
22 
23  enum PosType {
31  };
32 
33  enum RotType {
37  DIR_EX_EXY
38  };
39 
40 private:
41 
43  FFaEnumEntry(CART_X_Y_Z , "CART_X_Y_Z");
44  FFaEnumEntry(CYL_R_YR_X , "CYL_R_YR_X");
45  FFaEnumEntry(CYL_R_ZR_Y , "CYL_R_ZR_Y");
46  FFaEnumEntry(CYL_R_XR_Z , "CYL_R_XR_Z");
47  FFaEnumEntry(SPH_R_XR_YZ, "SPH_R_XR_YZ");
48  FFaEnumEntry(SPH_R_YR_ZX, "SPH_R_YR_ZX");
49  FFaEnumEntry(SPH_R_ZR_XY, "SPH_R_ZR_XY");
51  };
52 
54  FFaEnumEntry(EUL_Z_Y_X , "EUL_Z_Y_X" );
55  FFaEnumEntry(PNT_PZ_PXZ, "PNT_PZ_PXZ");
56  FFaEnumEntry(PNT_PX_PXY, "PNT_PX_PXY");
57  FFaEnumEntry(DIR_EX_EXY, "DIR_EX_EXY");
59  };
60 
61  PosTypeEnum myPosTyp = CART_X_Y_Z;
62  RotTypeEnum myRotTyp = EUL_Z_Y_X;
63 
64  std::array<FaVec3,3> myL;
65 
66  bool saveNumericalData = false;
67 
68 public:
72  FFa3DLocation(PosType t, const FaVec3& v0,
73  RotType r, const FaVec3& v1, const FaVec3& v2);
75  FFa3DLocation(PosType t, const FaVec3& v0, RotType r, const FaVec3& v1);
77  FFa3DLocation(const FaMat34& m);
79  FFa3DLocation(const FFa3DLocation& m) { *this = m; }
80 
83 
87  bool operator==(const FFa3DLocation& m) const;
88 
90  const FaVec3& operator[](int i) const;
92  FaVec3& operator[](int i);
93 
94  PosType getPosType() const { return myPosTyp; }
95  RotType getRotType() const { return myRotTyp; }
96  int getNumFields() const { return myRotTyp == EUL_Z_Y_X ? 6 : 9; }
97 
98  // Converting and setting
99 
100  bool changePosType(PosType newType);
101  bool changeRotType(RotType newType);
102 
103  FFa3DLocation& changePosRefCS(const FaMat34& newRef, const FaMat34& oldRef);
104  FFa3DLocation& changeRotRefCS(const FaMat34& newRef, const FaMat34& oldRef);
105 
106  FFa3DLocation& set(PosType p, const FaMat34& posRelMx,
107  RotType r, const FaMat34& rotRelMx, const FaMat34& glbPos);
108  FFa3DLocation& set(PosType p, RotType r, const FaMat34& glbPos);
109 
110  FFa3DLocation& setPos(PosType p, const FaVec3& cartPos);
111  FFa3DLocation& setRot(RotType r, const FaMat33& rotMat);
112 
113  FaVec3 translation() const;
114  FaMat33 direction() const;
115  FaMat34 getMatrix() const;
116  FaMat34 getMatrix(const FaMat34& posRelMx, const FaMat34& rotRelMx) const;
117  bool isValid() const;
118 
119  // Global operators
120 
121  friend std::ostream& operator<<(std::ostream& s, const FFa3DLocation& m);
122  friend std::istream& operator>>(std::istream& s, FFa3DLocation& m);
123 };
124 
125 
126 inline const FaVec3& FFa3DLocation::operator[] (int i) const
127 {
128 #ifdef FFA_INDEXCHECK
129  if (i < 0 || i > 2)
130  std::cerr <<"FFa3DLocation::operator[]: index i="<< i
131  <<" is out of range [0,2]"<< std::endl;
132 #endif
133  return myL[i];
134 }
135 
137 {
138 #ifdef FFA_INDEXCHECK
139  if (i < 0 || i > 2)
140  std::cerr <<"FFa3DLocation::operator[]: index i="<< i
141  << " is out of range [0,2]"<< std::endl;
142 #endif
143  return myL[i];
144 }
145 
146 #endif
Enum variables with text representations.
#define FFaEnumEntryEnd
Defines the end of an enum mapping definition.
Definition: FFaEnum.H:237
#define FFaEnumEntry(Value, Text)
Adds one entry to the current enum mapping definition.
Definition: FFaEnum.H:234
Point vectors in 3D space.
Definition: FFa3DLocation.H:20
FFa3DLocation & setPos(PosType p, const FaVec3 &cartPos)
Definition: FFa3DLocation.C:201
void saveMyNumericalData()
Let this instance save its numerical data.
Definition: FFa3DLocation.H:82
PosType getPosType() const
Definition: FFa3DLocation.H:94
std::array< FaVec3, 3 > myL
Definition: FFa3DLocation.H:64
FFa3DLocation & set(PosType p, const FaMat34 &posRelMx, RotType r, const FaMat34 &rotRelMx, const FaMat34 &glbPos)
Definition: FFa3DLocation.C:159
RotType
Definition: FFa3DLocation.H:33
@ EUL_Z_Y_X
Euler ZYX-angles. Supposed to be entered [ X, Y, Z].
Definition: FFa3DLocation.H:34
@ PNT_PX_PXY
Point on X-axis and point in XY-plane.
Definition: FFa3DLocation.H:35
@ PNT_PZ_PXZ
Point on Z-axis and point in XZ-plane.
Definition: FFa3DLocation.H:36
@ DIR_EX_EXY
Direction vector along X-axis, direction vector in XY-plane.
Definition: FFa3DLocation.H:37
FaVec3 translation() const
Definition: FFa3DLocation.C:289
bool isValid() const
Definition: FFa3DLocation.C:391
const FaVec3 & operator[](int i) const
Indexing operator.
Definition: FFa3DLocation.H:126
FFa3DLocation & changeRotRefCS(const FaMat34 &newRef, const FaMat34 &oldRef)
Definition: FFa3DLocation.C:132
FFa3DLocation & operator=(const FFa3DLocation &m)
Assignment operator.
Definition: FFa3DLocation.C:52
int getNumFields() const
Definition: FFa3DLocation.H:96
FFa3DLocation()
Default constructor.
Definition: FFa3DLocation.H:70
FaMat33 direction() const
Definition: FFa3DLocation.C:330
friend std::ostream & operator<<(std::ostream &s, const FFa3DLocation &m)
Definition: FFa3DLocation.C:430
bool operator==(const FFa3DLocation &m) const
Comparison operator.
Definition: FFa3DLocation.C:68
PosType
Definition: FFa3DLocation.H:23
@ CART_X_Y_Z
Cartesian coordinates.
Definition: FFa3DLocation.H:24
@ CYL_R_YR_X
Cylindric coordinates, angle from Y-axis about X-axis.
Definition: FFa3DLocation.H:25
@ SPH_R_ZR_XY
Spherical coordinates, asimuth angle from Z-axis.
Definition: FFa3DLocation.H:30
@ CYL_R_ZR_Y
Cylindric coordinates, angle from Z-axis about Y-axis.
Definition: FFa3DLocation.H:26
@ SPH_R_XR_YZ
Spherical coordinates, asimuth angle from X-axis.
Definition: FFa3DLocation.H:28
@ CYL_R_XR_Z
Cylindric coordinates, angle from X-axis about Z-axis.
Definition: FFa3DLocation.H:27
@ SPH_R_YR_ZX
Spherical coordinates, asimuth angle from Y-axis.
Definition: FFa3DLocation.H:29
friend std::istream & operator>>(std::istream &s, FFa3DLocation &m)
Definition: FFa3DLocation.C:453
FFaEnumMapping(PosType)
Definition: FFa3DLocation.H:42
FFa3DLocation(const FFa3DLocation &m)
Copy constructor.
Definition: FFa3DLocation.H:79
PosTypeEnum myPosTyp
Definition: FFa3DLocation.H:61
FaMat34 getMatrix() const
Definition: FFa3DLocation.C:366
RotTypeEnum myRotTyp
Definition: FFa3DLocation.H:62
FFa3DLocation & setRot(RotType r, const FaMat33 &rotMat)
Definition: FFa3DLocation.C:245
bool changeRotType(RotType newType)
Definition: FFa3DLocation.C:92
RotType getRotType() const
Definition: FFa3DLocation.H:95
bool saveNumericalData
If false, only the type enums are saved.
Definition: FFa3DLocation.H:66
FFaEnumMapping(RotType)
Definition: FFa3DLocation.H:53
FFa3DLocation & changePosRefCS(const FaMat34 &newRef, const FaMat34 &oldRef)
Definition: FFa3DLocation.C:107
bool changePosType(PosType newType)
Definition: FFa3DLocation.C:78
Definition: FFaMat33.H:15
Definition: FFaMat34.H:16
Class for point vectors in 3D space.
Definition: FFaVec3.H:40
real(dp), dimension(:,:), allocatable m
Definition: inverseModule.f90:38
real(dp), dimension(:), allocatable v0
Definition: inverseModule.f90:26
real(dp), dimension(:), allocatable v1
Definition: inverseModule.f90:27
real(dp), dimension(:), allocatable v2
Definition: inverseModule.f90:28