diff --git a/src/language/utils/EmbeddedDiscreteFunctionMathFunctions.cpp b/src/language/utils/EmbeddedDiscreteFunctionMathFunctions.cpp
index eb7cebd57f9b0e4d387a50deef77f2462fd8d810..b53e19ac62e5aab59c181311953b89864c435ab1 100644
--- a/src/language/utils/EmbeddedDiscreteFunctionMathFunctions.cpp
+++ b/src/language/utils/EmbeddedDiscreteFunctionMathFunctions.cpp
@@ -509,7 +509,7 @@ sum_of_Vh_components(const std::shared_ptr<const DiscreteFunctionVariant>& f)
template <size_t Dimension>
void
vectorize_to(const std::vector<std::shared_ptr<const DiscreteFunctionVariant>>& discrete_function_list,
- const Mesh<Connectivity<Dimension>>& mesh,
+ const Mesh<Dimension>& mesh,
DiscreteFunctionP0Vector<Dimension, double>& discrete_vector_function)
{
if (hasSameMesh(discrete_function_list)) {
diff --git a/src/mesh/CartesianMeshBuilder.cpp b/src/mesh/CartesianMeshBuilder.cpp
index 7b98ea861a5361a2d9a555dac3e0673181bc1090..4ae7410ce7b896e7cf1f8147a0b64365cebc9352 100644
--- a/src/mesh/CartesianMeshBuilder.cpp
+++ b/src/mesh/CartesianMeshBuilder.cpp
@@ -117,7 +117,7 @@ CartesianMeshBuilder::_buildCartesianMesh(const TinyVector<Dimension>& a,
NodeValue<TinyVector<Dimension>> xr = _getNodeCoordinates(a, b, cell_size, connectivity);
- m_mesh = std::make_shared<MeshVariant>(std::make_shared<const Mesh<ConnectivityType>>(p_connectivity, xr));
+ m_mesh = std::make_shared<MeshVariant>(std::make_shared<const Mesh<Dimension>>(p_connectivity, xr));
}
template <size_t Dimension>
diff --git a/src/mesh/Dual1DMeshBuilder.cpp b/src/mesh/Dual1DMeshBuilder.cpp
index c2d0f63422784bb0e8aebc57e6dda1b57d2cacc8..495bbbe862c20078e496241cdcbbd412e4626837 100644
--- a/src/mesh/Dual1DMeshBuilder.cpp
+++ b/src/mesh/Dual1DMeshBuilder.cpp
@@ -1,6 +1,5 @@
#include <mesh/Dual1DMeshBuilder.hpp>
-#include <mesh/Connectivity.hpp>
#include <mesh/DualConnectivityManager.hpp>
#include <mesh/ItemValueUtils.hpp>
#include <mesh/Mesh.hpp>
@@ -39,11 +38,11 @@ Dual1DMeshBuilder::_buildDual1DMeshFrom(const MeshType& primal_mesh)
NodeValue<TinyVector<1>> dual_xr{dual_connectivity};
primal_to_dual_1d_connectivity_data_mapper->toDualNode(primal_xr, primal_xj, dual_xr);
- m_mesh = std::make_shared<MeshVariant>(std::make_shared<const Mesh<Connectivity<1>>>(p_dual_connectivity, dual_xr));
+ m_mesh = std::make_shared<MeshVariant>(std::make_shared<const Mesh<1>>(p_dual_connectivity, dual_xr));
}
Dual1DMeshBuilder::Dual1DMeshBuilder(const std::shared_ptr<const MeshVariant>& mesh_v)
{
std::cout << "building Dual1DMesh\n";
- this->_buildDual1DMeshFrom(*(mesh_v->get<Mesh<Connectivity<1>>>()));
+ this->_buildDual1DMeshFrom(*(mesh_v->get<Mesh<1>>()));
}
diff --git a/src/mesh/DualMeshManager.hpp b/src/mesh/DualMeshManager.hpp
index 040dfd8d6ff3edaa51d93c97c7dc3a69f6ad98ef..46ad9c9f8d62cbff40a7d33efcc72c02e79f1ac7 100644
--- a/src/mesh/DualMeshManager.hpp
+++ b/src/mesh/DualMeshManager.hpp
@@ -10,12 +10,6 @@
class MeshVariant;
-template <size_t Dimension>
-class Connectivity;
-
-template <typename ConnectivityType>
-class Mesh;
-
class DualMeshManager
{
private:
diff --git a/src/mesh/GmshReader.cpp b/src/mesh/GmshReader.cpp
index b8ddd2910b9daed61a33aa285da3538a5321db17..50ecb1af18a15b53e193214d3606ae901e2341a2 100644
--- a/src/mesh/GmshReader.cpp
+++ b/src/mesh/GmshReader.cpp
@@ -1502,7 +1502,7 @@ GmshReader::__proceedData()
std::dynamic_pointer_cast<const Connectivity3D>(connectivity_builder.connectivity());
const Connectivity3D& connectivity = *p_connectivity;
- using MeshType = Mesh<Connectivity3D>;
+ using MeshType = Mesh<3>;
using Rd = TinyVector<3, double>;
NodeValue<Rd> xr(connectivity);
@@ -1523,7 +1523,7 @@ GmshReader::__proceedData()
std::dynamic_pointer_cast<const Connectivity2D>(connectivity_builder.connectivity());
const Connectivity2D& connectivity = *p_connectivity;
- using MeshType = Mesh<Connectivity2D>;
+ using MeshType = Mesh<2>;
using Rd = TinyVector<2, double>;
NodeValue<Rd> xr(connectivity);
@@ -1543,7 +1543,7 @@ GmshReader::__proceedData()
std::dynamic_pointer_cast<const Connectivity1D>(connectivity_builder.connectivity());
const Connectivity1D& connectivity = *p_connectivity;
- using MeshType = Mesh<Connectivity1D>;
+ using MeshType = Mesh<1>;
using Rd = TinyVector<1, double>;
NodeValue<Rd> xr(connectivity);
diff --git a/src/mesh/Mesh.cpp b/src/mesh/Mesh.cpp
index b622bfb975d6662bbc9742d9c50be6b2075d311f..baa0d00d71e7371514aad465870b7fafb0d3f460 100644
--- a/src/mesh/Mesh.cpp
+++ b/src/mesh/Mesh.cpp
@@ -4,13 +4,13 @@
#include <mesh/DualMeshManager.hpp>
#include <mesh/MeshDataManager.hpp>
-template <typename ConnectivityType>
-Mesh<ConnectivityType>::~Mesh()
+template <size_t Dimension>
+Mesh<Dimension>::~Mesh()
{
MeshDataManager::instance().deleteMeshData(this->id());
DualMeshManager::instance().deleteMesh(this->id());
}
-template Mesh<Connectivity<1>>::~Mesh();
-template Mesh<Connectivity<2>>::~Mesh();
-template Mesh<Connectivity<3>>::~Mesh();
+template Mesh<1>::~Mesh();
+template Mesh<2>::~Mesh();
+template Mesh<3>::~Mesh();
diff --git a/src/mesh/Mesh.hpp b/src/mesh/Mesh.hpp
index 451d53a31cda6421c34a6d8c7afb68391b81671f..f260dd1fc8b53eacdd7c0362ba587517d9763b67 100644
--- a/src/mesh/Mesh.hpp
+++ b/src/mesh/Mesh.hpp
@@ -2,20 +2,21 @@
#define MESH_HPP
#include <algebra/TinyVector.hpp>
+#include <mesh/Connectivity.hpp>
#include <mesh/ItemValue.hpp>
#include <mesh/MeshVariant.hpp>
#include <utils/GlobalVariableManager.hpp>
#include <memory>
-template <typename ConnectivityType>
-class Mesh : public std::enable_shared_from_this<Mesh<ConnectivityType>>
+template <size_t MeshDimension>
+class Mesh : public std::enable_shared_from_this<Mesh<MeshDimension>>
{
public:
- using Connectivity = ConnectivityType;
+ static constexpr size_t Dimension = MeshDimension;
+ using Connectivity = Connectivity<Dimension>;
- static constexpr size_t Dimension = ConnectivityType::Dimension;
- using Rd = TinyVector<Dimension>;
+ using Rd = TinyVector<Dimension>;
private:
const size_t m_id;
@@ -44,7 +45,7 @@ class Mesh : public std::enable_shared_from_this<Mesh<ConnectivityType>>
}
PUGS_INLINE
- std::shared_ptr<const Mesh<ConnectivityType>>
+ std::shared_ptr<const Mesh<Dimension>>
shared_ptr() const
{
return this->shared_from_this();
diff --git a/src/mesh/MeshBuilderBase.cpp b/src/mesh/MeshBuilderBase.cpp
index 6787e4bdc5aca6d6f728fd9a1c7356cd27be0d1b..4537f83b2177210f2bb0bf7e8e55b3219b57e6f8 100644
--- a/src/mesh/MeshBuilderBase.cpp
+++ b/src/mesh/MeshBuilderBase.cpp
@@ -23,7 +23,7 @@ MeshBuilderBase::_dispatch()
using ConnectivityType = Connectivity<Dimension>;
using Rd = TinyVector<Dimension>;
- using MeshType = Mesh<ConnectivityType>;
+ using MeshType = Mesh<Dimension>;
if (not m_mesh) {
ConnectivityDescriptor descriptor;
@@ -49,8 +49,8 @@ template <size_t Dimension>
void
MeshBuilderBase::_checkMesh() const
{
- using ConnectivityType = Connectivity<Dimension>;
- using MeshType = Mesh<ConnectivityType>;
+ using MeshType = Mesh<Dimension>;
+ using ConnectivityType = typename MeshType::Connectivity;
if (not m_mesh) {
throw UnexpectedError("mesh is not built yet");
diff --git a/src/mesh/MeshCellZone.cpp b/src/mesh/MeshCellZone.cpp
index ac63998b6edc658f39b9b9a1a11cf72f3bcf1c0a..08bdf64ba40c72449cabf2670a9c8c1112d0f2e8 100644
--- a/src/mesh/MeshCellZone.cpp
+++ b/src/mesh/MeshCellZone.cpp
@@ -5,13 +5,13 @@
#include <utils/Messenger.hpp>
template <size_t Dimension>
-MeshCellZone<Dimension>::MeshCellZone(const Mesh<Connectivity<Dimension>>&, const RefCellList& ref_cell_list)
+MeshCellZone<Dimension>::MeshCellZone(const Mesh<Dimension>&, const RefCellList& ref_cell_list)
: m_cell_list(ref_cell_list.list()), m_zone_name(ref_cell_list.refId().tagName())
{}
template <size_t Dimension>
MeshCellZone<Dimension>
-getMeshCellZone(const Mesh<Connectivity<Dimension>>& mesh, const IZoneDescriptor& zone_descriptor)
+getMeshCellZone(const Mesh<Dimension>& mesh, const IZoneDescriptor& zone_descriptor)
{
for (size_t i_ref_cell_list = 0; i_ref_cell_list < mesh.connectivity().template numberOfRefItemList<ItemType::cell>();
++i_ref_cell_list) {
@@ -28,6 +28,6 @@ getMeshCellZone(const Mesh<Connectivity<Dimension>>& mesh, const IZoneDescriptor
throw NormalError(ost.str());
}
-template MeshCellZone<1> getMeshCellZone(const Mesh<Connectivity<1>>&, const IZoneDescriptor&);
-template MeshCellZone<2> getMeshCellZone(const Mesh<Connectivity<2>>&, const IZoneDescriptor&);
-template MeshCellZone<3> getMeshCellZone(const Mesh<Connectivity<3>>&, const IZoneDescriptor&);
+template MeshCellZone<1> getMeshCellZone(const Mesh<1>&, const IZoneDescriptor&);
+template MeshCellZone<2> getMeshCellZone(const Mesh<2>&, const IZoneDescriptor&);
+template MeshCellZone<3> getMeshCellZone(const Mesh<3>&, const IZoneDescriptor&);
diff --git a/src/mesh/MeshCellZone.hpp b/src/mesh/MeshCellZone.hpp
index de4903689b328bddbac3f0bacf363497374f8fe4..1f8827b1b2998678f9f73ddbbdd2c18a18ed753d 100644
--- a/src/mesh/MeshCellZone.hpp
+++ b/src/mesh/MeshCellZone.hpp
@@ -6,9 +6,6 @@
#include <utils/Array.hpp>
template <size_t Dimension>
-class Connectivity;
-
-template <typename ConnectivityType>
class Mesh;
template <size_t Dimension>
@@ -20,30 +17,30 @@ class [[nodiscard]] MeshCellZone // clazy:exclude=copyable-polymorphic
public:
template <size_t MeshDimension>
- friend MeshCellZone<MeshDimension> getMeshCellZone(const Mesh<Connectivity<MeshDimension>>& mesh,
+ friend MeshCellZone<MeshDimension> getMeshCellZone(const Mesh<MeshDimension>& mesh,
const IZoneDescriptor& zone_descriptor);
MeshCellZone& operator=(const MeshCellZone&) = default;
- MeshCellZone& operator=(MeshCellZone&&) = default;
+ MeshCellZone& operator=(MeshCellZone&&) = default;
- const Array<const CellId>& cellList() const
+ const Array<const CellId>&
+ cellList() const
{
return m_cell_list;
}
protected:
- MeshCellZone(const Mesh<Connectivity<Dimension>>& mesh, const RefCellList& ref_cell_list);
+ MeshCellZone(const Mesh<Dimension>& mesh, const RefCellList& ref_cell_list);
public:
MeshCellZone(const MeshCellZone&) = default;
- MeshCellZone(MeshCellZone &&) = default;
+ MeshCellZone(MeshCellZone&&) = default;
MeshCellZone() = default;
virtual ~MeshCellZone() = default;
};
template <size_t Dimension>
-MeshCellZone<Dimension> getMeshCellZone(const Mesh<Connectivity<Dimension>>& mesh,
- const IZoneDescriptor& zone_descriptor);
+MeshCellZone<Dimension> getMeshCellZone(const Mesh<Dimension>& mesh, const IZoneDescriptor& zone_descriptor);
#endif // MESH_CELL_ZONE_HPP
diff --git a/src/mesh/MeshData.hpp b/src/mesh/MeshData.hpp
index 78b9b81e69daf5e3fc0010bd8c737b7303759353..a198aab4526215b1e88f25311a4d77de15513d1d 100644
--- a/src/mesh/MeshData.hpp
+++ b/src/mesh/MeshData.hpp
@@ -9,9 +9,6 @@
#include <utils/PugsUtils.hpp>
template <size_t Dimension>
-class Connectivity;
-
-template <typename ConnectivityType>
class Mesh;
template <size_t Dimension>
@@ -21,7 +18,7 @@ class MeshData : public IMeshData
static_assert(Dimension > 0, "dimension must be strictly positive");
static_assert((Dimension <= 3), "only 1d, 2d and 3d are implemented");
- using MeshType = Mesh<Connectivity<Dimension>>;
+ using MeshType = Mesh<Dimension>;
using Rd = TinyVector<Dimension>;
diff --git a/src/mesh/MeshDataManager.cpp b/src/mesh/MeshDataManager.cpp
index def1d1b617fd21953ec577e38e65572e8c127c07..3607fcad6b476f7b48c922891722b11a639d84a7 100644
--- a/src/mesh/MeshDataManager.cpp
+++ b/src/mesh/MeshDataManager.cpp
@@ -42,7 +42,7 @@ MeshDataManager::deleteMeshData(const size_t mesh_id)
template <size_t Dimension>
MeshData<Dimension>&
-MeshDataManager::getMeshData(const Mesh<Connectivity<Dimension>>& mesh)
+MeshDataManager::getMeshData(const Mesh<Dimension>& mesh)
{
if (auto i_mesh_data = m_mesh_id_mesh_data_map.find(mesh.id()); i_mesh_data != m_mesh_id_mesh_data_map.end()) {
return dynamic_cast<MeshData<Dimension>&>(*i_mesh_data->second);
@@ -55,6 +55,6 @@ MeshDataManager::getMeshData(const Mesh<Connectivity<Dimension>>& mesh)
}
}
-template MeshData<1>& MeshDataManager::getMeshData(const Mesh<Connectivity<1>>&);
-template MeshData<2>& MeshDataManager::getMeshData(const Mesh<Connectivity<2>>&);
-template MeshData<3>& MeshDataManager::getMeshData(const Mesh<Connectivity<3>>&);
+template MeshData<1>& MeshDataManager::getMeshData(const Mesh<1>&);
+template MeshData<2>& MeshDataManager::getMeshData(const Mesh<2>&);
+template MeshData<3>& MeshDataManager::getMeshData(const Mesh<3>&);
diff --git a/src/mesh/MeshDataManager.hpp b/src/mesh/MeshDataManager.hpp
index 8b1672b522671b8ff0a4e43b42004d96bf369137..0d2bf964b02daf341bbcd6e1b30c9cfe8bfee64e 100644
--- a/src/mesh/MeshDataManager.hpp
+++ b/src/mesh/MeshDataManager.hpp
@@ -8,10 +8,7 @@
#include <memory>
#include <unordered_map>
-template <size_t>
-class Connectivity;
-
-template <typename ConnectivityType>
+template <size_t Dimension>
class Mesh;
template <size_t Dimension>
@@ -45,7 +42,7 @@ class MeshDataManager
void deleteMeshData(const size_t mesh_id);
template <size_t Dimension>
- MeshData<Dimension>& getMeshData(const Mesh<Connectivity<Dimension>>&);
+ MeshData<Dimension>& getMeshData(const Mesh<Dimension>&);
};
#endif // MESH_DATA_MANAGER_HPP
diff --git a/src/mesh/MeshEdgeBoundary.cpp b/src/mesh/MeshEdgeBoundary.cpp
index 0a4aed04bb9e54e0f35dda08e3b9ce5e8f584559..0ad15dbad0c229461ba1de4d3d7eaf3ca797610f 100644
--- a/src/mesh/MeshEdgeBoundary.cpp
+++ b/src/mesh/MeshEdgeBoundary.cpp
@@ -9,9 +9,9 @@ template <typename MeshType>
MeshEdgeBoundary::MeshEdgeBoundary(const MeshType&, const RefEdgeList& ref_edge_list) : m_ref_edge_list(ref_edge_list)
{}
-template MeshEdgeBoundary::MeshEdgeBoundary(const Mesh<Connectivity<1>>&, const RefEdgeList&);
-template MeshEdgeBoundary::MeshEdgeBoundary(const Mesh<Connectivity<2>>&, const RefEdgeList&);
-template MeshEdgeBoundary::MeshEdgeBoundary(const Mesh<Connectivity<3>>&, const RefEdgeList&);
+template MeshEdgeBoundary::MeshEdgeBoundary(const Mesh<1>&, const RefEdgeList&);
+template MeshEdgeBoundary::MeshEdgeBoundary(const Mesh<2>&, const RefEdgeList&);
+template MeshEdgeBoundary::MeshEdgeBoundary(const Mesh<3>&, const RefEdgeList&);
template <typename MeshType>
MeshEdgeBoundary::MeshEdgeBoundary(const MeshType& mesh, const RefFaceList& ref_face_list)
@@ -55,8 +55,8 @@ MeshEdgeBoundary::MeshEdgeBoundary(const MeshType& mesh, const RefFaceList& ref_
const_cast<Connectivity<Dimension>&>(mesh.connectivity()).addRefItemList(m_ref_edge_list);
}
-template MeshEdgeBoundary::MeshEdgeBoundary(const Mesh<Connectivity<2>>&, const RefFaceList&);
-template MeshEdgeBoundary::MeshEdgeBoundary(const Mesh<Connectivity<3>>&, const RefFaceList&);
+template MeshEdgeBoundary::MeshEdgeBoundary(const Mesh<2>&, const RefFaceList&);
+template MeshEdgeBoundary::MeshEdgeBoundary(const Mesh<3>&, const RefFaceList&);
template <typename MeshType>
MeshEdgeBoundary
@@ -107,6 +107,6 @@ getMeshEdgeBoundary(const MeshType& mesh, const IBoundaryDescriptor& boundary_de
throw NormalError(ost.str());
}
-template MeshEdgeBoundary getMeshEdgeBoundary(const Mesh<Connectivity<1>>&, const IBoundaryDescriptor&);
-template MeshEdgeBoundary getMeshEdgeBoundary(const Mesh<Connectivity<2>>&, const IBoundaryDescriptor&);
-template MeshEdgeBoundary getMeshEdgeBoundary(const Mesh<Connectivity<3>>&, const IBoundaryDescriptor&);
+template MeshEdgeBoundary getMeshEdgeBoundary(const Mesh<1>&, const IBoundaryDescriptor&);
+template MeshEdgeBoundary getMeshEdgeBoundary(const Mesh<2>&, const IBoundaryDescriptor&);
+template MeshEdgeBoundary getMeshEdgeBoundary(const Mesh<3>&, const IBoundaryDescriptor&);
diff --git a/src/mesh/MeshEdgeInterface.cpp b/src/mesh/MeshEdgeInterface.cpp
index da4d0bc9f1ec42fd68150fab5944d3d446e42dcb..55235d8c7b6ff2ff8e4a5a1d54e3beb46d21c4f4 100644
--- a/src/mesh/MeshEdgeInterface.cpp
+++ b/src/mesh/MeshEdgeInterface.cpp
@@ -5,19 +5,19 @@
#include <mesh/Mesh.hpp>
#include <utils/Messenger.hpp>
-template <size_t Dimension>
-MeshEdgeInterface<Dimension>::MeshEdgeInterface(const Mesh<Connectivity<Dimension>>&, const RefEdgeList& ref_edge_list)
- : m_ref_edge_list(ref_edge_list)
+template <typename MeshType>
+MeshEdgeInterface::MeshEdgeInterface(const MeshType&, const RefEdgeList& ref_edge_list) : m_ref_edge_list(ref_edge_list)
{}
-template MeshEdgeInterface<1>::MeshEdgeInterface(const Mesh<Connectivity<1>>&, const RefEdgeList&);
-template MeshEdgeInterface<2>::MeshEdgeInterface(const Mesh<Connectivity<2>>&, const RefEdgeList&);
-template MeshEdgeInterface<3>::MeshEdgeInterface(const Mesh<Connectivity<3>>&, const RefEdgeList&);
+template MeshEdgeInterface::MeshEdgeInterface(const Mesh<1>&, const RefEdgeList&);
+template MeshEdgeInterface::MeshEdgeInterface(const Mesh<2>&, const RefEdgeList&);
+template MeshEdgeInterface::MeshEdgeInterface(const Mesh<3>&, const RefEdgeList&);
-template <size_t Dimension>
-MeshEdgeInterface<Dimension>::MeshEdgeInterface(const Mesh<Connectivity<Dimension>>& mesh,
- const RefFaceList& ref_face_list)
+template <typename MeshType>
+MeshEdgeInterface::MeshEdgeInterface(const MeshType& mesh, const RefFaceList& ref_face_list)
{
+ constexpr size_t Dimension = MeshType::Dimension;
+
const Array<const FaceId>& face_list = ref_face_list.list();
static_assert(Dimension > 1, "conversion from to edge from face is valid in dimension > 1");
@@ -55,13 +55,15 @@ MeshEdgeInterface<Dimension>::MeshEdgeInterface(const Mesh<Connectivity<Dimensio
const_cast<Connectivity<Dimension>&>(mesh.connectivity()).addRefItemList(m_ref_edge_list);
}
-template MeshEdgeInterface<2>::MeshEdgeInterface(const Mesh<Connectivity<2>>&, const RefFaceList&);
-template MeshEdgeInterface<3>::MeshEdgeInterface(const Mesh<Connectivity<3>>&, const RefFaceList&);
+template MeshEdgeInterface::MeshEdgeInterface(const Mesh<2>&, const RefFaceList&);
+template MeshEdgeInterface::MeshEdgeInterface(const Mesh<3>&, const RefFaceList&);
-template <size_t Dimension>
-MeshEdgeInterface<Dimension>
-getMeshEdgeInterface(const Mesh<Connectivity<Dimension>>& mesh, const IInterfaceDescriptor& interface_descriptor)
+template <typename MeshType>
+MeshEdgeInterface
+getMeshEdgeInterface(const MeshType& mesh, const IInterfaceDescriptor& interface_descriptor)
{
+ constexpr size_t Dimension = MeshType::Dimension;
+
for (size_t i_ref_edge_list = 0; i_ref_edge_list < mesh.connectivity().template numberOfRefItemList<ItemType::edge>();
++i_ref_edge_list) {
const auto& ref_edge_list = mesh.connectivity().template refItemList<ItemType::edge>(i_ref_edge_list);
@@ -76,7 +78,7 @@ getMeshEdgeInterface(const Mesh<Connectivity<Dimension>>& mesh, const IInterface
throw NormalError(ost.str());
}
- return MeshEdgeInterface<Dimension>{mesh, ref_edge_list};
+ return MeshEdgeInterface{mesh, ref_edge_list};
}
}
if constexpr (Dimension > 1) {
@@ -95,7 +97,7 @@ getMeshEdgeInterface(const Mesh<Connectivity<Dimension>>& mesh, const IInterface
throw NormalError(ost.str());
}
- return MeshEdgeInterface<Dimension>{mesh, ref_face_list};
+ return MeshEdgeInterface{mesh, ref_face_list};
}
}
}
@@ -106,6 +108,6 @@ getMeshEdgeInterface(const Mesh<Connectivity<Dimension>>& mesh, const IInterface
throw NormalError(ost.str());
}
-template MeshEdgeInterface<1> getMeshEdgeInterface(const Mesh<Connectivity<1>>&, const IInterfaceDescriptor&);
-template MeshEdgeInterface<2> getMeshEdgeInterface(const Mesh<Connectivity<2>>&, const IInterfaceDescriptor&);
-template MeshEdgeInterface<3> getMeshEdgeInterface(const Mesh<Connectivity<3>>&, const IInterfaceDescriptor&);
+template MeshEdgeInterface getMeshEdgeInterface(const Mesh<1>&, const IInterfaceDescriptor&);
+template MeshEdgeInterface getMeshEdgeInterface(const Mesh<2>&, const IInterfaceDescriptor&);
+template MeshEdgeInterface getMeshEdgeInterface(const Mesh<3>&, const IInterfaceDescriptor&);
diff --git a/src/mesh/MeshEdgeInterface.hpp b/src/mesh/MeshEdgeInterface.hpp
index 25226264c2b906066d2940e6a92a8a65b53042cc..45bcdbe84099303a48741697767887e5a9ca524f 100644
--- a/src/mesh/MeshEdgeInterface.hpp
+++ b/src/mesh/MeshEdgeInterface.hpp
@@ -6,52 +6,47 @@
#include <mesh/RefItemList.hpp>
#include <utils/Array.hpp>
-template <size_t Dimension>
-class Connectivity;
-
-template <typename ConnectivityType>
-class Mesh;
-
-template <size_t Dimension>
class [[nodiscard]] MeshEdgeInterface // clazy:exclude=copyable-polymorphic
{
protected:
RefEdgeList m_ref_edge_list;
public:
- template <size_t MeshDimension>
- friend MeshEdgeInterface<MeshDimension> getMeshEdgeInterface(const Mesh<Connectivity<MeshDimension>>& mesh,
- const IInterfaceDescriptor& interface_descriptor);
+ template <typename MeshType>
+ friend MeshEdgeInterface getMeshEdgeInterface(const MeshType& mesh, const IInterfaceDescriptor& interface_descriptor);
MeshEdgeInterface& operator=(const MeshEdgeInterface&) = default;
- MeshEdgeInterface& operator=(MeshEdgeInterface&&) = default;
+ MeshEdgeInterface& operator=(MeshEdgeInterface&&) = default;
PUGS_INLINE
- const RefEdgeList& refEdgeList() const
+ const RefEdgeList&
+ refEdgeList() const
{
return m_ref_edge_list;
}
PUGS_INLINE
- const Array<const EdgeId>& edgeList() const
+ const Array<const EdgeId>&
+ edgeList() const
{
return m_ref_edge_list.list();
}
protected:
- MeshEdgeInterface(const Mesh<Connectivity<Dimension>>& mesh, const RefEdgeList& ref_edge_list);
- MeshEdgeInterface(const Mesh<Connectivity<Dimension>>& mesh, const RefFaceList& ref_face_list);
+ template <typename MeshType>
+ MeshEdgeInterface(const MeshType& mesh, const RefEdgeList& ref_edge_list);
+ template <typename MeshType>
+ MeshEdgeInterface(const MeshType& mesh, const RefFaceList& ref_face_list);
public:
MeshEdgeInterface(const MeshEdgeInterface&) = default; // LCOV_EXCL_LINE
- MeshEdgeInterface(MeshEdgeInterface &&) = default; // LCOV_EXCL_LINE
+ MeshEdgeInterface(MeshEdgeInterface&&) = default; // LCOV_EXCL_LINE
MeshEdgeInterface() = default;
virtual ~MeshEdgeInterface() = default;
};
-template <size_t Dimension>
-MeshEdgeInterface<Dimension> getMeshEdgeInterface(const Mesh<Connectivity<Dimension>>& mesh,
- const IInterfaceDescriptor& interface_descriptor);
+template <typename MeshType>
+MeshEdgeInterface getMeshEdgeInterface(const MeshType& mesh, const IInterfaceDescriptor& interface_descriptor);
#endif // MESH_EDGE_INTERFACE_HPP
diff --git a/src/mesh/MeshFaceBoundary.cpp b/src/mesh/MeshFaceBoundary.cpp
index 3a578bb7893040ca3838581b0b7c5ad5d297d9a3..d883ba63af186bf6cf63bcca36a8e89db9fc06cd 100644
--- a/src/mesh/MeshFaceBoundary.cpp
+++ b/src/mesh/MeshFaceBoundary.cpp
@@ -8,9 +8,9 @@ template <typename MeshType>
MeshFaceBoundary::MeshFaceBoundary(const MeshType&, const RefFaceList& ref_face_list) : m_ref_face_list(ref_face_list)
{}
-template MeshFaceBoundary::MeshFaceBoundary(const Mesh<Connectivity<1>>&, const RefFaceList&);
-template MeshFaceBoundary::MeshFaceBoundary(const Mesh<Connectivity<2>>&, const RefFaceList&);
-template MeshFaceBoundary::MeshFaceBoundary(const Mesh<Connectivity<3>>&, const RefFaceList&);
+template MeshFaceBoundary::MeshFaceBoundary(const Mesh<1>&, const RefFaceList&);
+template MeshFaceBoundary::MeshFaceBoundary(const Mesh<2>&, const RefFaceList&);
+template MeshFaceBoundary::MeshFaceBoundary(const Mesh<3>&, const RefFaceList&);
template <typename MeshType>
MeshFaceBoundary
@@ -40,6 +40,6 @@ getMeshFaceBoundary(const MeshType& mesh, const IBoundaryDescriptor& boundary_de
throw NormalError(ost.str());
}
-template MeshFaceBoundary getMeshFaceBoundary(const Mesh<Connectivity<1>>&, const IBoundaryDescriptor&);
-template MeshFaceBoundary getMeshFaceBoundary(const Mesh<Connectivity<2>>&, const IBoundaryDescriptor&);
-template MeshFaceBoundary getMeshFaceBoundary(const Mesh<Connectivity<3>>&, const IBoundaryDescriptor&);
+template MeshFaceBoundary getMeshFaceBoundary(const Mesh<1>&, const IBoundaryDescriptor&);
+template MeshFaceBoundary getMeshFaceBoundary(const Mesh<2>&, const IBoundaryDescriptor&);
+template MeshFaceBoundary getMeshFaceBoundary(const Mesh<3>&, const IBoundaryDescriptor&);
diff --git a/src/mesh/MeshFaceInterface.cpp b/src/mesh/MeshFaceInterface.cpp
index afb73bf1872ece258788d1a4c6818f9d588250b6..663ac8d44027767f5318a424e83c9ca08f26cda8 100644
--- a/src/mesh/MeshFaceInterface.cpp
+++ b/src/mesh/MeshFaceInterface.cpp
@@ -4,18 +4,17 @@
#include <mesh/Mesh.hpp>
#include <utils/Messenger.hpp>
-template <size_t Dimension>
-MeshFaceInterface<Dimension>::MeshFaceInterface(const Mesh<Connectivity<Dimension>>&, const RefFaceList& ref_face_list)
- : m_ref_face_list(ref_face_list)
+template <typename MeshType>
+MeshFaceInterface::MeshFaceInterface(const MeshType&, const RefFaceList& ref_face_list) : m_ref_face_list(ref_face_list)
{}
-template MeshFaceInterface<1>::MeshFaceInterface(const Mesh<Connectivity<1>>&, const RefFaceList&);
-template MeshFaceInterface<2>::MeshFaceInterface(const Mesh<Connectivity<2>>&, const RefFaceList&);
-template MeshFaceInterface<3>::MeshFaceInterface(const Mesh<Connectivity<3>>&, const RefFaceList&);
+template MeshFaceInterface::MeshFaceInterface(const Mesh<1>&, const RefFaceList&);
+template MeshFaceInterface::MeshFaceInterface(const Mesh<2>&, const RefFaceList&);
+template MeshFaceInterface::MeshFaceInterface(const Mesh<3>&, const RefFaceList&);
-template <size_t Dimension>
-MeshFaceInterface<Dimension>
-getMeshFaceInterface(const Mesh<Connectivity<Dimension>>& mesh, const IInterfaceDescriptor& interface_descriptor)
+template <typename MeshType>
+MeshFaceInterface
+getMeshFaceInterface(const MeshType& mesh, const IInterfaceDescriptor& interface_descriptor)
{
for (size_t i_ref_face_list = 0; i_ref_face_list < mesh.connectivity().template numberOfRefItemList<ItemType::face>();
++i_ref_face_list) {
@@ -30,7 +29,7 @@ getMeshFaceInterface(const Mesh<Connectivity<Dimension>>& mesh, const IInterface
throw NormalError(ost.str());
}
- return MeshFaceInterface<Dimension>{mesh, ref_face_list};
+ return MeshFaceInterface{mesh, ref_face_list};
}
}
@@ -40,6 +39,6 @@ getMeshFaceInterface(const Mesh<Connectivity<Dimension>>& mesh, const IInterface
throw NormalError(ost.str());
}
-template MeshFaceInterface<1> getMeshFaceInterface(const Mesh<Connectivity<1>>&, const IInterfaceDescriptor&);
-template MeshFaceInterface<2> getMeshFaceInterface(const Mesh<Connectivity<2>>&, const IInterfaceDescriptor&);
-template MeshFaceInterface<3> getMeshFaceInterface(const Mesh<Connectivity<3>>&, const IInterfaceDescriptor&);
+template MeshFaceInterface getMeshFaceInterface(const Mesh<1>&, const IInterfaceDescriptor&);
+template MeshFaceInterface getMeshFaceInterface(const Mesh<2>&, const IInterfaceDescriptor&);
+template MeshFaceInterface getMeshFaceInterface(const Mesh<3>&, const IInterfaceDescriptor&);
diff --git a/src/mesh/MeshFaceInterface.hpp b/src/mesh/MeshFaceInterface.hpp
index 87d836a3705199aab546396d24052025b59b5729..0d8277e0bc9ee7a0504e01235a2ecfcfbcc749c3 100644
--- a/src/mesh/MeshFaceInterface.hpp
+++ b/src/mesh/MeshFaceInterface.hpp
@@ -6,51 +6,45 @@
#include <mesh/RefItemList.hpp>
#include <utils/Array.hpp>
-template <size_t Dimension>
-class Connectivity;
-
-template <typename ConnectivityType>
-class Mesh;
-
-template <size_t Dimension>
class [[nodiscard]] MeshFaceInterface // clazy:exclude=copyable-polymorphic
{
protected:
RefFaceList m_ref_face_list;
public:
- template <size_t MeshDimension>
- friend MeshFaceInterface<MeshDimension> getMeshFaceInterface(const Mesh<Connectivity<MeshDimension>>& mesh,
- const IInterfaceDescriptor& interface_descriptor);
+ template <typename MeshType>
+ friend MeshFaceInterface getMeshFaceInterface(const MeshType& mesh, const IInterfaceDescriptor& interface_descriptor);
MeshFaceInterface& operator=(const MeshFaceInterface&) = default;
- MeshFaceInterface& operator=(MeshFaceInterface&&) = default;
+ MeshFaceInterface& operator=(MeshFaceInterface&&) = default;
PUGS_INLINE
- const RefFaceList& refFaceList() const
+ const RefFaceList&
+ refFaceList() const
{
return m_ref_face_list;
}
PUGS_INLINE
- const Array<const FaceId>& faceList() const
+ const Array<const FaceId>&
+ faceList() const
{
return m_ref_face_list.list();
}
protected:
- MeshFaceInterface(const Mesh<Connectivity<Dimension>>& mesh, const RefFaceList& ref_face_list);
+ template <typename MeshType>
+ MeshFaceInterface(const MeshType& mesh, const RefFaceList& ref_face_list);
public:
MeshFaceInterface(const MeshFaceInterface&) = default; // LCOV_EXCL_LINE
- MeshFaceInterface(MeshFaceInterface &&) = default; // LCOV_EXCL_LINE
+ MeshFaceInterface(MeshFaceInterface&&) = default; // LCOV_EXCL_LINE
MeshFaceInterface() = default;
virtual ~MeshFaceInterface() = default;
};
-template <size_t Dimension>
-MeshFaceInterface<Dimension> getMeshFaceInterface(const Mesh<Connectivity<Dimension>>& mesh,
- const IInterfaceDescriptor& interface_descriptor);
+template <typename MeshType>
+MeshFaceInterface getMeshFaceInterface(const MeshType& mesh, const IInterfaceDescriptor& interface_descriptor);
#endif // MESH_FACE_INTERFACE_HPP
diff --git a/src/mesh/MeshFlatEdgeBoundary.cpp b/src/mesh/MeshFlatEdgeBoundary.cpp
index 3783049ff2422e899f90210c1ffba722f71a4767..8b5669a340902d5a68cfc0c5ee05bb60bee20400 100644
--- a/src/mesh/MeshFlatEdgeBoundary.cpp
+++ b/src/mesh/MeshFlatEdgeBoundary.cpp
@@ -15,9 +15,6 @@ getMeshFlatEdgeBoundary(const MeshType& mesh, const IBoundaryDescriptor& boundar
mesh_flat_node_boundary.outgoingNormal()};
}
-template MeshFlatEdgeBoundary<const Mesh<Connectivity<1>>> getMeshFlatEdgeBoundary(const Mesh<Connectivity<1>>&,
- const IBoundaryDescriptor&);
-template MeshFlatEdgeBoundary<const Mesh<Connectivity<2>>> getMeshFlatEdgeBoundary(const Mesh<Connectivity<2>>&,
- const IBoundaryDescriptor&);
-template MeshFlatEdgeBoundary<const Mesh<Connectivity<3>>> getMeshFlatEdgeBoundary(const Mesh<Connectivity<3>>&,
- const IBoundaryDescriptor&);
+template MeshFlatEdgeBoundary<const Mesh<1>> getMeshFlatEdgeBoundary(const Mesh<1>&, const IBoundaryDescriptor&);
+template MeshFlatEdgeBoundary<const Mesh<2>> getMeshFlatEdgeBoundary(const Mesh<2>&, const IBoundaryDescriptor&);
+template MeshFlatEdgeBoundary<const Mesh<3>> getMeshFlatEdgeBoundary(const Mesh<3>&, const IBoundaryDescriptor&);
diff --git a/src/mesh/MeshFlatFaceBoundary.cpp b/src/mesh/MeshFlatFaceBoundary.cpp
index e05373897dce3421af35fa94a06bb89dd760ec35..12d7f439f19d91dc03614e68a916809fcb9c0bff 100644
--- a/src/mesh/MeshFlatFaceBoundary.cpp
+++ b/src/mesh/MeshFlatFaceBoundary.cpp
@@ -15,9 +15,6 @@ getMeshFlatFaceBoundary(const MeshType& mesh, const IBoundaryDescriptor& boundar
mesh_flat_node_boundary.outgoingNormal()};
}
-template MeshFlatFaceBoundary<const Mesh<Connectivity<1>>> getMeshFlatFaceBoundary(const Mesh<Connectivity<1>>&,
- const IBoundaryDescriptor&);
-template MeshFlatFaceBoundary<const Mesh<Connectivity<2>>> getMeshFlatFaceBoundary(const Mesh<Connectivity<2>>&,
- const IBoundaryDescriptor&);
-template MeshFlatFaceBoundary<const Mesh<Connectivity<3>>> getMeshFlatFaceBoundary(const Mesh<Connectivity<3>>&,
- const IBoundaryDescriptor&);
+template MeshFlatFaceBoundary<const Mesh<1>> getMeshFlatFaceBoundary(const Mesh<1>&, const IBoundaryDescriptor&);
+template MeshFlatFaceBoundary<const Mesh<2>> getMeshFlatFaceBoundary(const Mesh<2>&, const IBoundaryDescriptor&);
+template MeshFlatFaceBoundary<const Mesh<3>> getMeshFlatFaceBoundary(const Mesh<3>&, const IBoundaryDescriptor&);
diff --git a/src/mesh/MeshFlatNodeBoundary.cpp b/src/mesh/MeshFlatNodeBoundary.cpp
index 3ee564332b57c1bdf5a1f219af6417d9b9c0844c..df763198866961715a193dfde9f193fdfa67d649 100644
--- a/src/mesh/MeshFlatNodeBoundary.cpp
+++ b/src/mesh/MeshFlatNodeBoundary.cpp
@@ -34,7 +34,7 @@ MeshFlatNodeBoundary<MeshType>::_checkBoundaryIsFlat(const TinyVector<MeshType::
template <>
TinyVector<1, double>
-MeshFlatNodeBoundary<const Mesh<Connectivity<1>>>::_getNormal(const Mesh<Connectivity<1>>& mesh)
+MeshFlatNodeBoundary<const Mesh<1>>::_getNormal(const Mesh<1>& mesh)
{
using R = TinyVector<1, double>;
@@ -60,7 +60,7 @@ MeshFlatNodeBoundary<const Mesh<Connectivity<1>>>::_getNormal(const Mesh<Connect
template <>
TinyVector<2, double>
-MeshFlatNodeBoundary<const Mesh<Connectivity<2>>>::_getNormal(const Mesh<Connectivity<2>>& mesh)
+MeshFlatNodeBoundary<const Mesh<2>>::_getNormal(const Mesh<2>& mesh)
{
using R2 = TinyVector<2, double>;
@@ -88,9 +88,7 @@ MeshFlatNodeBoundary<const Mesh<Connectivity<2>>>::_getNormal(const Mesh<Connect
template <>
TinyVector<3, double>
-MeshFlatNodeBoundary<const Mesh<Connectivity<3>>>::_getFarestNode(const Mesh<Connectivity<3>>& mesh,
- const Rd& x0,
- const Rd& x1)
+MeshFlatNodeBoundary<const Mesh<3>>::_getFarestNode(const Mesh<3>& mesh, const Rd& x0, const Rd& x1)
{
const NodeValue<const Rd>& xr = mesh.xr();
const auto node_number = mesh.connectivity().nodeNumber();
@@ -141,7 +139,7 @@ MeshFlatNodeBoundary<const Mesh<Connectivity<3>>>::_getFarestNode(const Mesh<Con
template <>
TinyVector<3, double>
-MeshFlatNodeBoundary<const Mesh<Connectivity<3>>>::_getNormal(const Mesh<Connectivity<3>>& mesh)
+MeshFlatNodeBoundary<const Mesh<3>>::_getNormal(const Mesh<3>& mesh)
{
using R3 = TinyVector<3, double>;
@@ -198,7 +196,7 @@ MeshFlatNodeBoundary<const Mesh<Connectivity<3>>>::_getNormal(const Mesh<Connect
template <>
TinyVector<1, double>
-MeshFlatNodeBoundary<const Mesh<Connectivity<1>>>::_getOutgoingNormal(const Mesh<Connectivity<1>>& mesh)
+MeshFlatNodeBoundary<const Mesh<1>>::_getOutgoingNormal(const Mesh<1>& mesh)
{
using R = TinyVector<1, double>;
@@ -241,7 +239,7 @@ MeshFlatNodeBoundary<const Mesh<Connectivity<1>>>::_getOutgoingNormal(const Mesh
template <>
TinyVector<2, double>
-MeshFlatNodeBoundary<const Mesh<Connectivity<2>>>::_getOutgoingNormal(const Mesh<Connectivity<2>>& mesh)
+MeshFlatNodeBoundary<const Mesh<2>>::_getOutgoingNormal(const Mesh<2>& mesh)
{
using R2 = TinyVector<2, double>;
@@ -284,7 +282,7 @@ MeshFlatNodeBoundary<const Mesh<Connectivity<2>>>::_getOutgoingNormal(const Mesh
template <>
TinyVector<3, double>
-MeshFlatNodeBoundary<const Mesh<Connectivity<3>>>::_getOutgoingNormal(const Mesh<Connectivity<3>>& mesh)
+MeshFlatNodeBoundary<const Mesh<3>>::_getOutgoingNormal(const Mesh<3>& mesh)
{
using R3 = TinyVector<3, double>;
@@ -352,9 +350,6 @@ getMeshFlatNodeBoundary(const MeshType& mesh, const IBoundaryDescriptor& boundar
throw NormalError(ost.str());
}
-template MeshFlatNodeBoundary<const Mesh<Connectivity<1>>> getMeshFlatNodeBoundary(const Mesh<Connectivity<1>>&,
- const IBoundaryDescriptor&);
-template MeshFlatNodeBoundary<const Mesh<Connectivity<2>>> getMeshFlatNodeBoundary(const Mesh<Connectivity<2>>&,
- const IBoundaryDescriptor&);
-template MeshFlatNodeBoundary<const Mesh<Connectivity<3>>> getMeshFlatNodeBoundary(const Mesh<Connectivity<3>>&,
- const IBoundaryDescriptor&);
+template MeshFlatNodeBoundary<const Mesh<1>> getMeshFlatNodeBoundary(const Mesh<1>&, const IBoundaryDescriptor&);
+template MeshFlatNodeBoundary<const Mesh<2>> getMeshFlatNodeBoundary(const Mesh<2>&, const IBoundaryDescriptor&);
+template MeshFlatNodeBoundary<const Mesh<3>> getMeshFlatNodeBoundary(const Mesh<3>&, const IBoundaryDescriptor&);
diff --git a/src/mesh/MeshLineEdgeBoundary.cpp b/src/mesh/MeshLineEdgeBoundary.cpp
index 31c4bb52465870eb4c6def597d1faa89e8df9e73..1d0c3752675689e2b178c70031ca4d0bfe2f8570 100644
--- a/src/mesh/MeshLineEdgeBoundary.cpp
+++ b/src/mesh/MeshLineEdgeBoundary.cpp
@@ -16,7 +16,5 @@ getMeshLineEdgeBoundary(const MeshType& mesh, const IBoundaryDescriptor& boundar
mesh_line_node_boundary.direction()};
}
-template MeshLineEdgeBoundary<const Mesh<Connectivity<2>>> getMeshLineEdgeBoundary(const Mesh<Connectivity<2>>&,
- const IBoundaryDescriptor&);
-template MeshLineEdgeBoundary<const Mesh<Connectivity<3>>> getMeshLineEdgeBoundary(const Mesh<Connectivity<3>>&,
- const IBoundaryDescriptor&);
+template MeshLineEdgeBoundary<const Mesh<2>> getMeshLineEdgeBoundary(const Mesh<2>&, const IBoundaryDescriptor&);
+template MeshLineEdgeBoundary<const Mesh<3>> getMeshLineEdgeBoundary(const Mesh<3>&, const IBoundaryDescriptor&);
diff --git a/src/mesh/MeshLineFaceBoundary.cpp b/src/mesh/MeshLineFaceBoundary.cpp
index e81aef554062fe97509a2ebd770652828743d265..83bfdf8cf069c71db3b2daf0d1f382acb7295d53 100644
--- a/src/mesh/MeshLineFaceBoundary.cpp
+++ b/src/mesh/MeshLineFaceBoundary.cpp
@@ -16,5 +16,4 @@ getMeshLineFaceBoundary(const MeshType& mesh, const IBoundaryDescriptor& boundar
mesh_line_node_boundary.direction()};
}
-template MeshLineFaceBoundary<const Mesh<Connectivity<2>>> getMeshLineFaceBoundary(const Mesh<Connectivity<2>>&,
- const IBoundaryDescriptor&);
+template MeshLineFaceBoundary<const Mesh<2>> getMeshLineFaceBoundary(const Mesh<2>&, const IBoundaryDescriptor&);
diff --git a/src/mesh/MeshLineNodeBoundary.cpp b/src/mesh/MeshLineNodeBoundary.cpp
index a03c6102ba1ceeb37c35b75f3cd02040313b569d..dd145d889a81f1d653add3802089fbbfcdc11262 100644
--- a/src/mesh/MeshLineNodeBoundary.cpp
+++ b/src/mesh/MeshLineNodeBoundary.cpp
@@ -38,7 +38,7 @@ MeshLineNodeBoundary<MeshType>::_checkBoundaryIsLine(const TinyVector<MeshType::
template <>
TinyVector<2>
-MeshLineNodeBoundary<const Mesh<Connectivity<2>>>::_getDirection(const Mesh<Connectivity<2>>& mesh)
+MeshLineNodeBoundary<const Mesh<2>>::_getDirection(const Mesh<2>& mesh)
{
using R2 = TinyVector<2, double>;
@@ -65,7 +65,7 @@ MeshLineNodeBoundary<const Mesh<Connectivity<2>>>::_getDirection(const Mesh<Conn
template <>
TinyVector<3>
-MeshLineNodeBoundary<const Mesh<Connectivity<3>>>::_getDirection(const Mesh<Connectivity<3>>& mesh)
+MeshLineNodeBoundary<const Mesh<3>>::_getDirection(const Mesh<3>& mesh)
{
using R3 = TinyVector<3, double>;
@@ -140,7 +140,5 @@ getMeshLineNodeBoundary(const MeshType& mesh, const IBoundaryDescriptor& boundar
throw NormalError(ost.str());
}
-template MeshLineNodeBoundary<const Mesh<Connectivity<2>>> getMeshLineNodeBoundary(const Mesh<Connectivity<2>>&,
- const IBoundaryDescriptor&);
-template MeshLineNodeBoundary<const Mesh<Connectivity<3>>> getMeshLineNodeBoundary(const Mesh<Connectivity<3>>&,
- const IBoundaryDescriptor&);
+template MeshLineNodeBoundary<const Mesh<2>> getMeshLineNodeBoundary(const Mesh<2>&, const IBoundaryDescriptor&);
+template MeshLineNodeBoundary<const Mesh<3>> getMeshLineNodeBoundary(const Mesh<3>&, const IBoundaryDescriptor&);
diff --git a/src/mesh/MeshNodeBoundary.cpp b/src/mesh/MeshNodeBoundary.cpp
index 97f32e9d433553e4d1826f21be08a46e48a91a78..e750580ce1298e4ea110847d89ddf385ca893c38 100644
--- a/src/mesh/MeshNodeBoundary.cpp
+++ b/src/mesh/MeshNodeBoundary.cpp
@@ -109,17 +109,17 @@ MeshNodeBoundary::MeshNodeBoundary(const MeshType&, const RefNodeList& ref_node_
}
}
-template MeshNodeBoundary::MeshNodeBoundary(const Mesh<Connectivity<1>>&, const RefFaceList&);
-template MeshNodeBoundary::MeshNodeBoundary(const Mesh<Connectivity<2>>&, const RefFaceList&);
-template MeshNodeBoundary::MeshNodeBoundary(const Mesh<Connectivity<3>>&, const RefFaceList&);
+template MeshNodeBoundary::MeshNodeBoundary(const Mesh<1>&, const RefFaceList&);
+template MeshNodeBoundary::MeshNodeBoundary(const Mesh<2>&, const RefFaceList&);
+template MeshNodeBoundary::MeshNodeBoundary(const Mesh<3>&, const RefFaceList&);
-template MeshNodeBoundary::MeshNodeBoundary(const Mesh<Connectivity<1>>&, const RefEdgeList&);
-template MeshNodeBoundary::MeshNodeBoundary(const Mesh<Connectivity<2>>&, const RefEdgeList&);
-template MeshNodeBoundary::MeshNodeBoundary(const Mesh<Connectivity<3>>&, const RefEdgeList&);
+template MeshNodeBoundary::MeshNodeBoundary(const Mesh<1>&, const RefEdgeList&);
+template MeshNodeBoundary::MeshNodeBoundary(const Mesh<2>&, const RefEdgeList&);
+template MeshNodeBoundary::MeshNodeBoundary(const Mesh<3>&, const RefEdgeList&);
-template MeshNodeBoundary::MeshNodeBoundary(const Mesh<Connectivity<1>>&, const RefNodeList&);
-template MeshNodeBoundary::MeshNodeBoundary(const Mesh<Connectivity<2>>&, const RefNodeList&);
-template MeshNodeBoundary::MeshNodeBoundary(const Mesh<Connectivity<3>>&, const RefNodeList&);
+template MeshNodeBoundary::MeshNodeBoundary(const Mesh<1>&, const RefNodeList&);
+template MeshNodeBoundary::MeshNodeBoundary(const Mesh<2>&, const RefNodeList&);
+template MeshNodeBoundary::MeshNodeBoundary(const Mesh<3>&, const RefNodeList&);
template <typename MeshType>
MeshNodeBoundary
@@ -156,6 +156,6 @@ getMeshNodeBoundary(const MeshType& mesh, const IBoundaryDescriptor& boundary_de
throw NormalError(ost.str());
}
-template MeshNodeBoundary getMeshNodeBoundary(const Mesh<Connectivity<1>>&, const IBoundaryDescriptor&);
-template MeshNodeBoundary getMeshNodeBoundary(const Mesh<Connectivity<2>>&, const IBoundaryDescriptor&);
-template MeshNodeBoundary getMeshNodeBoundary(const Mesh<Connectivity<3>>&, const IBoundaryDescriptor&);
+template MeshNodeBoundary getMeshNodeBoundary(const Mesh<1>&, const IBoundaryDescriptor&);
+template MeshNodeBoundary getMeshNodeBoundary(const Mesh<2>&, const IBoundaryDescriptor&);
+template MeshNodeBoundary getMeshNodeBoundary(const Mesh<3>&, const IBoundaryDescriptor&);
diff --git a/src/mesh/MeshNodeBoundaryUtils.cpp b/src/mesh/MeshNodeBoundaryUtils.cpp
index 63df3e36475d2d60029fe5cd27dba9f37eac21cc..d6a0e31cd9a671530216b9dce15747e07f06c9b4 100644
--- a/src/mesh/MeshNodeBoundaryUtils.cpp
+++ b/src/mesh/MeshNodeBoundaryUtils.cpp
@@ -6,7 +6,7 @@
template <>
std::array<TinyVector<2>, 2>
-getBounds(const Mesh<Connectivity<2>>& mesh, const RefNodeList& ref_node_list)
+getBounds(const Mesh<2>& mesh, const RefNodeList& ref_node_list)
{
using R2 = TinyVector<2, double>;
@@ -54,7 +54,7 @@ getBounds(const Mesh<Connectivity<2>>& mesh, const RefNodeList& ref_node_list)
template <>
std::array<TinyVector<3>, 6>
-getBounds(const Mesh<Connectivity<3>>& mesh, const RefNodeList& ref_node_list)
+getBounds(const Mesh<3>& mesh, const RefNodeList& ref_node_list)
{
using R3 = TinyVector<3, double>;
diff --git a/src/mesh/MeshNodeInterface.cpp b/src/mesh/MeshNodeInterface.cpp
index 84a5b2669a27bedbd8e11e8ffbe67d1860051f56..1703f36a024fb81d7df27a6e74204eb087923349 100644
--- a/src/mesh/MeshNodeInterface.cpp
+++ b/src/mesh/MeshNodeInterface.cpp
@@ -1,14 +1,16 @@
#include <mesh/MeshNodeInterface.hpp>
#include <Kokkos_Vector.hpp>
-#include <mesh/Connectivity.hpp>
#include <mesh/Mesh.hpp>
+#include <mesh/MeshTraits.hpp>
#include <utils/Messenger.hpp>
-template <size_t Dimension>
-MeshNodeInterface<Dimension>::MeshNodeInterface(const Mesh<Connectivity<Dimension>>& mesh,
- const RefFaceList& ref_face_list)
+template <typename MeshType>
+MeshNodeInterface::MeshNodeInterface(const MeshType& mesh, const RefFaceList& ref_face_list)
{
+ static_assert(is_polygonal_mesh<MeshType>);
+ constexpr size_t Dimension = MeshType::Dimension;
+
const Array<const FaceId>& face_list = ref_face_list.list();
if (ref_face_list.type() != RefItemListBase::Type::interface) {
std::ostringstream ost;
@@ -51,10 +53,12 @@ MeshNodeInterface<Dimension>::MeshNodeInterface(const Mesh<Connectivity<Dimensio
const_cast<Connectivity<Dimension>&>(mesh.connectivity()).addRefItemList(m_ref_node_list);
}
-template <size_t Dimension>
-MeshNodeInterface<Dimension>::MeshNodeInterface(const Mesh<Connectivity<Dimension>>& mesh,
- const RefEdgeList& ref_edge_list)
+template <typename MeshType>
+MeshNodeInterface::MeshNodeInterface(const MeshType& mesh, const RefEdgeList& ref_edge_list)
{
+ static_assert(is_polygonal_mesh<MeshType>);
+ constexpr size_t Dimension = MeshType::Dimension;
+
const Array<const EdgeId>& edge_list = ref_edge_list.list();
if (ref_edge_list.type() != RefItemListBase::Type::interface) {
std::ostringstream ost;
@@ -96,9 +100,8 @@ MeshNodeInterface<Dimension>::MeshNodeInterface(const Mesh<Connectivity<Dimensio
const_cast<Connectivity<Dimension>&>(mesh.connectivity()).addRefItemList(m_ref_node_list);
}
-template <size_t Dimension>
-MeshNodeInterface<Dimension>::MeshNodeInterface(const Mesh<Connectivity<Dimension>>&, const RefNodeList& ref_node_list)
- : m_ref_node_list(ref_node_list)
+template <typename MeshType>
+MeshNodeInterface::MeshNodeInterface(const MeshType&, const RefNodeList& ref_node_list) : m_ref_node_list(ref_node_list)
{
if (ref_node_list.type() != RefItemListBase::Type::interface) {
std::ostringstream ost;
@@ -108,28 +111,28 @@ MeshNodeInterface<Dimension>::MeshNodeInterface(const Mesh<Connectivity<Dimensio
}
}
-template MeshNodeInterface<1>::MeshNodeInterface(const Mesh<Connectivity<1>>&, const RefFaceList&);
-template MeshNodeInterface<2>::MeshNodeInterface(const Mesh<Connectivity<2>>&, const RefFaceList&);
-template MeshNodeInterface<3>::MeshNodeInterface(const Mesh<Connectivity<3>>&, const RefFaceList&);
+template MeshNodeInterface::MeshNodeInterface(const Mesh<1>&, const RefFaceList&);
+template MeshNodeInterface::MeshNodeInterface(const Mesh<2>&, const RefFaceList&);
+template MeshNodeInterface::MeshNodeInterface(const Mesh<3>&, const RefFaceList&);
-template MeshNodeInterface<1>::MeshNodeInterface(const Mesh<Connectivity<1>>&, const RefEdgeList&);
-template MeshNodeInterface<2>::MeshNodeInterface(const Mesh<Connectivity<2>>&, const RefEdgeList&);
-template MeshNodeInterface<3>::MeshNodeInterface(const Mesh<Connectivity<3>>&, const RefEdgeList&);
+template MeshNodeInterface::MeshNodeInterface(const Mesh<1>&, const RefEdgeList&);
+template MeshNodeInterface::MeshNodeInterface(const Mesh<2>&, const RefEdgeList&);
+template MeshNodeInterface::MeshNodeInterface(const Mesh<3>&, const RefEdgeList&);
-template MeshNodeInterface<1>::MeshNodeInterface(const Mesh<Connectivity<1>>&, const RefNodeList&);
-template MeshNodeInterface<2>::MeshNodeInterface(const Mesh<Connectivity<2>>&, const RefNodeList&);
-template MeshNodeInterface<3>::MeshNodeInterface(const Mesh<Connectivity<3>>&, const RefNodeList&);
+template MeshNodeInterface::MeshNodeInterface(const Mesh<1>&, const RefNodeList&);
+template MeshNodeInterface::MeshNodeInterface(const Mesh<2>&, const RefNodeList&);
+template MeshNodeInterface::MeshNodeInterface(const Mesh<3>&, const RefNodeList&);
-template <size_t Dimension>
-MeshNodeInterface<Dimension>
-getMeshNodeInterface(const Mesh<Connectivity<Dimension>>& mesh, const IInterfaceDescriptor& interface_descriptor)
+template <typename MeshType>
+MeshNodeInterface
+getMeshNodeInterface(const MeshType& mesh, const IInterfaceDescriptor& interface_descriptor)
{
for (size_t i_ref_node_list = 0; i_ref_node_list < mesh.connectivity().template numberOfRefItemList<ItemType::node>();
++i_ref_node_list) {
const auto& ref_node_list = mesh.connectivity().template refItemList<ItemType::node>(i_ref_node_list);
const RefId& ref = ref_node_list.refId();
if (ref == interface_descriptor) {
- return MeshNodeInterface<Dimension>{mesh, ref_node_list};
+ return MeshNodeInterface{mesh, ref_node_list};
}
}
for (size_t i_ref_edge_list = 0; i_ref_edge_list < mesh.connectivity().template numberOfRefItemList<ItemType::edge>();
@@ -137,7 +140,7 @@ getMeshNodeInterface(const Mesh<Connectivity<Dimension>>& mesh, const IInterface
const auto& ref_edge_list = mesh.connectivity().template refItemList<ItemType::edge>(i_ref_edge_list);
const RefId& ref = ref_edge_list.refId();
if (ref == interface_descriptor) {
- return MeshNodeInterface<Dimension>{mesh, ref_edge_list};
+ return MeshNodeInterface{mesh, ref_edge_list};
}
}
for (size_t i_ref_face_list = 0; i_ref_face_list < mesh.connectivity().template numberOfRefItemList<ItemType::face>();
@@ -145,7 +148,7 @@ getMeshNodeInterface(const Mesh<Connectivity<Dimension>>& mesh, const IInterface
const auto& ref_face_list = mesh.connectivity().template refItemList<ItemType::face>(i_ref_face_list);
const RefId& ref = ref_face_list.refId();
if (ref == interface_descriptor) {
- return MeshNodeInterface<Dimension>{mesh, ref_face_list};
+ return MeshNodeInterface{mesh, ref_face_list};
}
}
@@ -155,6 +158,6 @@ getMeshNodeInterface(const Mesh<Connectivity<Dimension>>& mesh, const IInterface
throw NormalError(ost.str());
}
-template MeshNodeInterface<1> getMeshNodeInterface(const Mesh<Connectivity<1>>&, const IInterfaceDescriptor&);
-template MeshNodeInterface<2> getMeshNodeInterface(const Mesh<Connectivity<2>>&, const IInterfaceDescriptor&);
-template MeshNodeInterface<3> getMeshNodeInterface(const Mesh<Connectivity<3>>&, const IInterfaceDescriptor&);
+template MeshNodeInterface getMeshNodeInterface(const Mesh<1>&, const IInterfaceDescriptor&);
+template MeshNodeInterface getMeshNodeInterface(const Mesh<2>&, const IInterfaceDescriptor&);
+template MeshNodeInterface getMeshNodeInterface(const Mesh<3>&, const IInterfaceDescriptor&);
diff --git a/src/mesh/MeshNodeInterface.hpp b/src/mesh/MeshNodeInterface.hpp
index 8fcbb6704df80ffa48697fd20109e4fab36e2b98..4eaafd06b56dfcb07d9e1bc5c0b29d661c0386e1 100644
--- a/src/mesh/MeshNodeInterface.hpp
+++ b/src/mesh/MeshNodeInterface.hpp
@@ -6,56 +6,49 @@
#include <mesh/RefItemList.hpp>
#include <utils/Array.hpp>
-template <size_t Dimension>
-class Connectivity;
-
-template <typename ConnectivityType>
-class Mesh;
-
-template <size_t Dimension>
class [[nodiscard]] MeshNodeInterface // clazy:exclude=copyable-polymorphic
{
protected:
RefNodeList m_ref_node_list;
- std::array<TinyVector<Dimension>, Dimension*(Dimension - 1)> _getBounds(const Mesh<Connectivity<Dimension>>& mesh)
- const;
-
public:
- template <size_t MeshDimension>
- friend MeshNodeInterface<MeshDimension> getMeshNodeInterface(const Mesh<Connectivity<MeshDimension>>& mesh,
- const IInterfaceDescriptor& interface_descriptor);
+ template <typename MeshType>
+ friend MeshNodeInterface getMeshNodeInterface(const MeshType& mesh, const IInterfaceDescriptor& interface_descriptor);
MeshNodeInterface& operator=(const MeshNodeInterface&) = default;
- MeshNodeInterface& operator=(MeshNodeInterface&&) = default;
+ MeshNodeInterface& operator=(MeshNodeInterface&&) = default;
PUGS_INLINE
- const RefNodeList& refNodeList() const
+ const RefNodeList&
+ refNodeList() const
{
return m_ref_node_list;
}
PUGS_INLINE
- const Array<const NodeId>& nodeList() const
+ const Array<const NodeId>&
+ nodeList() const
{
return m_ref_node_list.list();
}
protected:
- MeshNodeInterface(const Mesh<Connectivity<Dimension>>& mesh, const RefFaceList& ref_face_list);
- MeshNodeInterface(const Mesh<Connectivity<Dimension>>& mesh, const RefEdgeList& ref_edge_list);
- MeshNodeInterface(const Mesh<Connectivity<Dimension>>&, const RefNodeList& ref_node_list);
+ template <typename MeshType>
+ MeshNodeInterface(const MeshType& mesh, const RefFaceList& ref_face_list);
+ template <typename MeshType>
+ MeshNodeInterface(const MeshType& mesh, const RefEdgeList& ref_edge_list);
+ template <typename MeshType>
+ MeshNodeInterface(const MeshType&, const RefNodeList& ref_node_list);
public:
MeshNodeInterface(const MeshNodeInterface&) = default;
- MeshNodeInterface(MeshNodeInterface &&) = default;
+ MeshNodeInterface(MeshNodeInterface&&) = default;
MeshNodeInterface() = default;
virtual ~MeshNodeInterface() = default;
};
-template <size_t Dimension>
-MeshNodeInterface<Dimension> getMeshNodeInterface(const Mesh<Connectivity<Dimension>>& mesh,
- const IInterfaceDescriptor& interface_descriptor);
+template <typename MeshType>
+MeshNodeInterface getMeshNodeInterface(const MeshType& mesh, const IInterfaceDescriptor& interface_descriptor);
#endif // MESH_NODE_INTERFACE_HPP
diff --git a/src/mesh/MeshRandomizer.cpp b/src/mesh/MeshRandomizer.cpp
index 9a2a56d5a754d96a195cb88a6e70ac5f0d5c6dab..a5b0e03514ec869230c7ee96d29e18268097b82e 100644
--- a/src/mesh/MeshRandomizer.cpp
+++ b/src/mesh/MeshRandomizer.cpp
@@ -306,7 +306,7 @@ class MeshRandomizerHandler::MeshRandomizer<MeshType>::AxisBoundaryCondition
};
template <>
-class MeshRandomizerHandler::MeshRandomizer<Mesh<Connectivity<1>>>::AxisBoundaryCondition
+class MeshRandomizerHandler::MeshRandomizer<Mesh<1>>::AxisBoundaryCondition
{
public:
AxisBoundaryCondition() = default;
diff --git a/src/mesh/MeshRelaxer.hpp b/src/mesh/MeshRelaxer.hpp
index deb9f8a4afa263d8dc852002d4db5ef748e85596..1e2b082064dadb1f85d1c1e2c93151b331993c1c 100644
--- a/src/mesh/MeshRelaxer.hpp
+++ b/src/mesh/MeshRelaxer.hpp
@@ -3,9 +3,6 @@
class MeshVariant;
-template <typename ConnectivityType>
-class Mesh;
-
#include <memory>
class MeshRelaxer
diff --git a/src/mesh/MeshSmoother.cpp b/src/mesh/MeshSmoother.cpp
index d90b95a4d9356e441155b844bc33eeb754aaee62..5eee712d8c7a3f3b513f24cf48385c239f926a54 100644
--- a/src/mesh/MeshSmoother.cpp
+++ b/src/mesh/MeshSmoother.cpp
@@ -353,7 +353,7 @@ class MeshSmootherHandler::MeshSmoother<MeshType>::AxisBoundaryCondition
};
template <>
-class MeshSmootherHandler::MeshSmoother<Mesh<Connectivity<1>>>::AxisBoundaryCondition
+class MeshSmootherHandler::MeshSmoother<Mesh<1>>::AxisBoundaryCondition
{
public:
AxisBoundaryCondition() = default;
diff --git a/src/mesh/MeshTraits.hpp b/src/mesh/MeshTraits.hpp
index 62ea2282523481325dd01502c93d198a261b6c0f..bf1a8a8f95f508edca7704f7d0d58dd1a5cc18cb 100644
--- a/src/mesh/MeshTraits.hpp
+++ b/src/mesh/MeshTraits.hpp
@@ -3,19 +3,16 @@
#include <cstddef>
-template <size_t>
-class Connectivity;
-
-template <typename ConnectivityType>
+template <size_t Dimension>
class Mesh;
template <typename MeshType>
constexpr inline bool is_polygonal_mesh = false;
template <size_t Dimension>
-constexpr inline bool is_polygonal_mesh<Mesh<Connectivity<Dimension>>> = true;
+constexpr inline bool is_polygonal_mesh<Mesh<Dimension>> = true;
template <size_t Dimension>
-constexpr inline bool is_polygonal_mesh<const Mesh<Connectivity<Dimension>>> = true;
+constexpr inline bool is_polygonal_mesh<const Mesh<Dimension>> = true;
#endif // MESH_TRAITS_HPP
diff --git a/src/mesh/MeshVariant.hpp b/src/mesh/MeshVariant.hpp
index bcc57e4cb601cec9130ce18c43a69b16ad5fafe0..53e4327952334bda21186fb47153722f03e5a924 100644
--- a/src/mesh/MeshVariant.hpp
+++ b/src/mesh/MeshVariant.hpp
@@ -7,18 +7,15 @@
#include <memory>
#include <variant>
-template <size_t>
-class Connectivity;
-
-template <typename ConnectivityType>
+template <size_t Dimension>
class Mesh;
class MeshVariant
{
private:
- using Variant = std::variant<std::shared_ptr<const Mesh<Connectivity<1>>>,
- std::shared_ptr<const Mesh<Connectivity<2>>>,
- std::shared_ptr<const Mesh<Connectivity<3>>>>;
+ using Variant = std::variant<std::shared_ptr<const Mesh<1>>, //
+ std::shared_ptr<const Mesh<2>>, //
+ std::shared_ptr<const Mesh<3>>>;
Variant m_p_mesh;
diff --git a/src/scheme/AcousticSolver.cpp b/src/scheme/AcousticSolver.cpp
index e60372f0c46b1d3eed0852c5fc60486b6f1c611d..3cd71fe5dd4b7bbfbdb959e86851ab8c4259c5b6 100644
--- a/src/scheme/AcousticSolver.cpp
+++ b/src/scheme/AcousticSolver.cpp
@@ -26,7 +26,7 @@ template <size_t Dimension>
double
acoustic_dt(const DiscreteFunctionP0<Dimension, const double>& c)
{
- const Mesh<Connectivity<Dimension>>& mesh = dynamic_cast<const Mesh<Connectivity<Dimension>>&>(*c.mesh());
+ const Mesh<Dimension>& mesh = dynamic_cast<const Mesh<Dimension>&>(*c.mesh());
const auto Vj = MeshDataManager::instance().getMeshData(mesh).Vj();
const auto Sj = MeshDataManager::instance().getMeshData(mesh).sumOverRLjr();
@@ -205,7 +205,7 @@ class AcousticSolverHandler::AcousticSolver final : public AcousticSolverHandler
}
NodeValue<Rd>
- _computeBr(const Mesh<Connectivity<Dimension>>& mesh,
+ _computeBr(const Mesh<Dimension>& mesh,
const NodeValuePerCell<const Rdxd>& Ajr,
const DiscreteVectorFunction& u,
const DiscreteScalarFunction& p) const
@@ -735,7 +735,7 @@ class AcousticSolverHandler::AcousticSolver<MeshType>::PressureBoundaryCondition
};
template <>
-class AcousticSolverHandler::AcousticSolver<const Mesh<Connectivity<1>>>::PressureBoundaryCondition
+class AcousticSolverHandler::AcousticSolver<const Mesh<1>>::PressureBoundaryCondition
{
private:
const MeshNodeBoundary m_mesh_node_boundary;
@@ -794,7 +794,7 @@ class AcousticSolverHandler::AcousticSolver<MeshType>::ExternalFSIVelocityBounda
return m_value_list;
}
- ExternalFSIVelocityBoundaryCondition(const Mesh<Connectivity<Dimension>>& mesh,
+ ExternalFSIVelocityBoundaryCondition(const Mesh<Dimension>& mesh,
const MeshNodeBoundary& mesh_node_boundary,
const std::shared_ptr<const Socket>& socket)
: m_node_to_cell_matrix{mesh.connectivity().nodeToCellMatrix()},
diff --git a/src/scheme/DiscreteFunctionP0.hpp b/src/scheme/DiscreteFunctionP0.hpp
index 0e11b608adf97704bdc345277abde6a49b843fca..9f7fc39c5d9e4af0d66ba276bc6500bc7e467440 100644
--- a/src/scheme/DiscreteFunctionP0.hpp
+++ b/src/scheme/DiscreteFunctionP0.hpp
@@ -3,7 +3,6 @@
#include <language/utils/ASTNodeDataTypeTraits.hpp>
-#include <mesh/Connectivity.hpp>
#include <mesh/ItemValueUtils.hpp>
#include <mesh/Mesh.hpp>
#include <mesh/MeshData.hpp>
@@ -16,7 +15,7 @@ class DiscreteFunctionP0
{
public:
using data_type = DataType;
- using MeshType = Mesh<Connectivity<Dimension>>;
+ using MeshType = Mesh<Dimension>;
friend class DiscreteFunctionP0<Dimension, std::add_const_t<DataType>>;
friend class DiscreteFunctionP0<Dimension, std::remove_const_t<DataType>>;
@@ -714,8 +713,8 @@ class DiscreteFunctionP0
{
Assert(f.m_cell_values.isBuilt());
- const Mesh<Connectivity<Dimension>>& mesh = *f.m_mesh;
- CellValue<const double> cell_volume = MeshDataManager::instance().getMeshData(mesh).Vj();
+ const Mesh<Dimension>& mesh = *f.m_mesh;
+ CellValue<const double> cell_volume = MeshDataManager::instance().getMeshData(mesh).Vj();
CellValue<std::remove_const_t<DataType>> f_v{mesh.connectivity()};
parallel_for(
diff --git a/src/scheme/DiscreteFunctionP0Vector.hpp b/src/scheme/DiscreteFunctionP0Vector.hpp
index 7801626a48a4c9d3961912ad1a7e3f0a28fb37b0..ada6df3adf4a3f899c1d0fdbd87413941a762006 100644
--- a/src/scheme/DiscreteFunctionP0Vector.hpp
+++ b/src/scheme/DiscreteFunctionP0Vector.hpp
@@ -2,7 +2,6 @@
#define DISCRETE_FUNCTION_P0_VECTOR_HPP
#include <algebra/Vector.hpp>
-#include <mesh/Connectivity.hpp>
#include <mesh/ItemArray.hpp>
#include <mesh/Mesh.hpp>
#include <mesh/MeshData.hpp>
@@ -16,7 +15,7 @@ class DiscreteFunctionP0Vector
{
public:
using data_type = DataType;
- using MeshType = Mesh<Connectivity<Dimension>>;
+ using MeshType = Mesh<Dimension>;
friend class DiscreteFunctionP0Vector<Dimension, std::add_const_t<DataType>>;
friend class DiscreteFunctionP0Vector<Dimension, std::remove_const_t<DataType>>;
diff --git a/src/scheme/FluxingAdvectionSolver.cpp b/src/scheme/FluxingAdvectionSolver.cpp
index 25d9cf3de756c2b21f6a85447b4a69359cea4767..0ee4e25f56968b202d2b6a49156555bbbacc7cea 100644
--- a/src/scheme/FluxingAdvectionSolver.cpp
+++ b/src/scheme/FluxingAdvectionSolver.cpp
@@ -96,7 +96,7 @@ class FluxingAdvectionSolver
_storeValueBCList(const size_t i_quantity,
const std::vector<std::shared_ptr<const IBoundaryConditionDescriptor>>& bc_list)
{
- const Mesh<Connectivity<Dimension>>& mesh = *m_old_mesh;
+ const Mesh<Dimension>& mesh = *m_old_mesh;
for (auto& i_bc : bc_list) {
switch (i_bc->type()) {
@@ -139,7 +139,7 @@ class FluxingAdvectionSolver
_storeArrayBCList(const size_t i_quantity,
const std::vector<std::shared_ptr<const IBoundaryConditionDescriptor>>& bc_list)
{
- const Mesh<Connectivity<Dimension>>& mesh = *m_old_mesh;
+ const Mesh<Dimension>& mesh = *m_old_mesh;
for (auto& i_bc : bc_list) {
switch (i_bc->type()) {
@@ -240,8 +240,7 @@ FluxingAdvectionSolver<MeshType>::_computeDonorCells(FaceValue<const double> alg
template <>
void
-FluxingAdvectionSolver<const Mesh<Connectivity<1>>>::_computeDonorCells(
- FaceValue<const double> algebraic_fluxing_volumes)
+FluxingAdvectionSolver<const Mesh<1>>::_computeDonorCells(FaceValue<const double> algebraic_fluxing_volumes)
{
m_donnor_cell = [&] {
const auto face_to_cell_matrix = m_new_mesh->connectivity().faceToCellMatrix();
@@ -275,7 +274,7 @@ FluxingAdvectionSolver<const Mesh<Connectivity<1>>>::_computeDonorCells(
template <>
FaceValue<double>
-FluxingAdvectionSolver<const Mesh<Connectivity<1>>>::_computeAlgebraicFluxingVolume()
+FluxingAdvectionSolver<const Mesh<1>>::_computeAlgebraicFluxingVolume()
{
Array<double> fluxing_volumes{m_new_mesh->numberOfNodes()};
NodeValue<double> nodal_fluxing_volume(m_new_mesh->connectivity(), fluxing_volumes);
@@ -294,7 +293,7 @@ FluxingAdvectionSolver<const Mesh<Connectivity<1>>>::_computeAlgebraicFluxingVol
template <>
FaceValue<double>
-FluxingAdvectionSolver<const Mesh<Connectivity<2>>>::_computeAlgebraicFluxingVolume()
+FluxingAdvectionSolver<const Mesh<2>>::_computeAlgebraicFluxingVolume()
{
const auto face_to_node_matrix = m_old_mesh->connectivity().faceToNodeMatrix();
FaceValue<double> algebraic_fluxing_volume(m_new_mesh->connectivity());
@@ -321,7 +320,7 @@ FluxingAdvectionSolver<const Mesh<Connectivity<2>>>::_computeAlgebraicFluxingVol
template <>
FaceValue<double>
-FluxingAdvectionSolver<const Mesh<Connectivity<3>>>::_computeAlgebraicFluxingVolume()
+FluxingAdvectionSolver<const Mesh<3>>::_computeAlgebraicFluxingVolume()
{
// due to the z component of the jacobian determinant, degree 3
// polynomials must be exactly integrated
diff --git a/src/scheme/HyperelasticSolver.cpp b/src/scheme/HyperelasticSolver.cpp
index 1e0654aebedcdc00108963f7ca05416752422008..edd4e4950e60e0f4e875e3df4ed4a80d49c49a35 100644
--- a/src/scheme/HyperelasticSolver.cpp
+++ b/src/scheme/HyperelasticSolver.cpp
@@ -209,7 +209,7 @@ class HyperelasticSolverHandler::HyperelasticSolver final : public HyperelasticS
}
NodeValue<Rd>
- _computeBr(const Mesh<Connectivity<Dimension>>& mesh,
+ _computeBr(const Mesh<Dimension>& mesh,
const NodeValuePerCell<const Rdxd>& Ajr,
const DiscreteVectorFunction& u,
const DiscreteTensorFunction& sigma) const
@@ -800,7 +800,7 @@ class HyperelasticSolverHandler::HyperelasticSolver<MeshType>::PressureBoundaryC
};
template <>
-class HyperelasticSolverHandler::HyperelasticSolver<const Mesh<Connectivity<1>>>::PressureBoundaryCondition
+class HyperelasticSolverHandler::HyperelasticSolver<const Mesh<1>>::PressureBoundaryCondition
{
private:
const MeshNodeBoundary m_mesh_node_boundary;
@@ -854,7 +854,7 @@ class HyperelasticSolverHandler::HyperelasticSolver<MeshType>::NormalStressBound
};
template <>
-class HyperelasticSolverHandler::HyperelasticSolver<const Mesh<Connectivity<1>>>::NormalStressBoundaryCondition
+class HyperelasticSolverHandler::HyperelasticSolver<const Mesh<1>>::NormalStressBoundaryCondition
{
private:
const MeshNodeBoundary m_mesh_node_boundary;
diff --git a/tests/test_CellIntegrator.cpp b/tests/test_CellIntegrator.cpp
index cce926b6827464a6ef9310586315b5227043c1a1..d5e7dddece7d2e644194ac5c6cb4e62bcd3dea07 100644
--- a/tests/test_CellIntegrator.cpp
+++ b/tests/test_CellIntegrator.cpp
@@ -21,7 +21,7 @@ TEST_CASE("CellIntegrator", "[scheme]")
{
using R1 = TinyVector<1>;
- const auto mesh = MeshDataBaseForTests::get().unordered1DMesh()->get<Mesh<Connectivity<1>>>();
+ const auto mesh = MeshDataBaseForTests::get().unordered1DMesh()->get<Mesh<1>>();
auto f = [](const R1& x) -> double { return x[0] * x[0] + 1; };
Array<const double> int_f_per_cell = [=] {
@@ -236,7 +236,7 @@ TEST_CASE("CellIntegrator", "[scheme]")
{
using R2 = TinyVector<2>;
- const auto mesh = MeshDataBaseForTests::get().hybrid2DMesh()->get<Mesh<Connectivity<2>>>();
+ const auto mesh = MeshDataBaseForTests::get().hybrid2DMesh()->get<Mesh<2>>();
auto f = [](const R2& X) -> double {
const double x = X[0];
@@ -499,7 +499,7 @@ TEST_CASE("CellIntegrator", "[scheme]")
for (const auto& mesh_info : mesh_list) {
auto mesh_name = mesh_info.first;
- auto mesh = mesh_info.second->get<Mesh<Connectivity<3>>>();
+ auto mesh = mesh_info.second->get<Mesh<3>>();
SECTION(mesh_name)
{
@@ -943,7 +943,7 @@ TEST_CASE("CellIntegrator", "[scheme]")
{
using R1 = TinyVector<1>;
- const auto mesh = MeshDataBaseForTests::get().unordered1DMesh()->get<Mesh<Connectivity<1>>>();
+ const auto mesh = MeshDataBaseForTests::get().unordered1DMesh()->get<Mesh<1>>();
auto f = [](const R1& x) -> R2 { return R2{x[0] * x[0] + 1, 2 * x[0]}; };
Array<const R2> int_f_per_cell = [=] {
@@ -1156,7 +1156,7 @@ TEST_CASE("CellIntegrator", "[scheme]")
SECTION("2D")
{
- const auto mesh = MeshDataBaseForTests::get().hybrid2DMesh()->get<Mesh<Connectivity<2>>>();
+ const auto mesh = MeshDataBaseForTests::get().hybrid2DMesh()->get<Mesh<2>>();
auto f = [](const R2& X) -> R2 {
const double x = X[0];
@@ -1418,7 +1418,7 @@ TEST_CASE("CellIntegrator", "[scheme]")
for (const auto& mesh_info : mesh_list) {
auto mesh_name = mesh_info.first;
- auto mesh = mesh_info.second->get<Mesh<Connectivity<3>>>();
+ auto mesh = mesh_info.second->get<Mesh<3>>();
SECTION(mesh_name)
{
@@ -1863,7 +1863,7 @@ TEST_CASE("CellIntegrator", "[scheme]")
{
using R1 = TinyVector<1>;
- const auto mesh = MeshDataBaseForTests::get().unordered1DMesh()->get<Mesh<Connectivity<1>>>();
+ const auto mesh = MeshDataBaseForTests::get().unordered1DMesh()->get<Mesh<1>>();
auto f = [](const R1& x) -> R2x2 { return R2x2{x[0] * x[0] + 1, 2 * x[0], 3 - x[0], 3 * x[0] - 1}; };
Array<const R2x2> int_f_per_cell = [=] {
@@ -2081,7 +2081,7 @@ TEST_CASE("CellIntegrator", "[scheme]")
{
using R2 = TinyVector<2>;
- const auto mesh = MeshDataBaseForTests::get().hybrid2DMesh()->get<Mesh<Connectivity<2>>>();
+ const auto mesh = MeshDataBaseForTests::get().hybrid2DMesh()->get<Mesh<2>>();
auto f = [](const R2& X) -> R2x2 {
const double x = X[0];
@@ -2345,7 +2345,7 @@ TEST_CASE("CellIntegrator", "[scheme]")
for (const auto& mesh_info : mesh_list) {
auto mesh_name = mesh_info.first;
- auto mesh = mesh_info.second->get<Mesh<Connectivity<3>>>();
+ auto mesh = mesh_info.second->get<Mesh<3>>();
SECTION(mesh_name)
{
diff --git a/tests/test_Connectivity.cpp b/tests/test_Connectivity.cpp
index a0277eaf6ac21f12b0025e8197efa829552eaeed..c64988ec2ffb34c273290d65709a0056e392222f 100644
--- a/tests/test_Connectivity.cpp
+++ b/tests/test_Connectivity.cpp
@@ -20,7 +20,7 @@ TEST_CASE("Connectivity", "[mesh]")
{
SECTION("unordered 1D mesh")
{
- const auto& mesh = *MeshDataBaseForTests::get().unordered1DMesh()->get<Mesh<Connectivity<1>>>();
+ const auto& mesh = *MeshDataBaseForTests::get().unordered1DMesh()->get<Mesh<1>>();
if (parallel::size() == 1) {
REQUIRE(mesh.numberOfNodes() == 35);
@@ -37,7 +37,7 @@ TEST_CASE("Connectivity", "[mesh]")
SECTION("cartesian 1D mesh")
{
- const auto& mesh = *MeshDataBaseForTests::get().cartesian1DMesh()->get<Mesh<Connectivity<1>>>();
+ const auto& mesh = *MeshDataBaseForTests::get().cartesian1DMesh()->get<Mesh<1>>();
if (parallel::size() == 1) {
REQUIRE(mesh.numberOfNodes() == 24);
@@ -57,7 +57,7 @@ TEST_CASE("Connectivity", "[mesh]")
{
SECTION("hybrid 2D mesh")
{
- const auto& mesh = *MeshDataBaseForTests::get().hybrid2DMesh()->get<Mesh<Connectivity<2>>>();
+ const auto& mesh = *MeshDataBaseForTests::get().hybrid2DMesh()->get<Mesh<2>>();
if (parallel::size() == 1) {
REQUIRE(mesh.numberOfNodes() == 53);
@@ -74,7 +74,7 @@ TEST_CASE("Connectivity", "[mesh]")
SECTION("cartesian 2D mesh")
{
- const auto& mesh = *MeshDataBaseForTests::get().cartesian2DMesh()->get<Mesh<Connectivity<2>>>();
+ const auto& mesh = *MeshDataBaseForTests::get().cartesian2DMesh()->get<Mesh<2>>();
if (parallel::size() == 1) {
REQUIRE(mesh.numberOfNodes() == 56);
@@ -94,7 +94,7 @@ TEST_CASE("Connectivity", "[mesh]")
{
SECTION("hybrid 3D mesh")
{
- const auto& mesh = *MeshDataBaseForTests::get().hybrid3DMesh()->get<Mesh<Connectivity<3>>>();
+ const auto& mesh = *MeshDataBaseForTests::get().hybrid3DMesh()->get<Mesh<3>>();
if (parallel::size() == 1) {
REQUIRE(mesh.numberOfNodes() == 132);
@@ -111,7 +111,7 @@ TEST_CASE("Connectivity", "[mesh]")
SECTION("cartesian 3D mesh")
{
- const auto& mesh = *MeshDataBaseForTests::get().cartesian3DMesh()->get<Mesh<Connectivity<3>>>();
+ const auto& mesh = *MeshDataBaseForTests::get().cartesian3DMesh()->get<Mesh<3>>();
if (parallel::size() == 1) {
REQUIRE(mesh.numberOfNodes() == 280);
@@ -137,7 +137,7 @@ TEST_CASE("Connectivity", "[mesh]")
for (const auto& named_mesh : mesh_list) {
SECTION(named_mesh.name())
{
- auto mesh_1d = named_mesh.mesh()->get<Mesh<Connectivity<1>>>();
+ auto mesh_1d = named_mesh.mesh()->get<Mesh<1>>();
const Connectivity<1>& connectivity = mesh_1d->connectivity();
SECTION("nodes/edges/faces")
@@ -228,7 +228,7 @@ TEST_CASE("Connectivity", "[mesh]")
for (const auto& named_mesh : mesh_list) {
SECTION(named_mesh.name())
{
- auto mesh_2d = named_mesh.mesh()->get<Mesh<Connectivity<2>>>();
+ auto mesh_2d = named_mesh.mesh()->get<Mesh<2>>();
const Connectivity<2>& connectivity = mesh_2d->connectivity();
SECTION("nodes")
@@ -331,7 +331,7 @@ TEST_CASE("Connectivity", "[mesh]")
for (const auto& named_mesh : mesh_list) {
SECTION(named_mesh.name())
{
- auto mesh_3d = named_mesh.mesh()->get<Mesh<Connectivity<3>>>();
+ auto mesh_3d = named_mesh.mesh()->get<Mesh<3>>();
const Connectivity<3>& connectivity = mesh_3d->connectivity();
@@ -448,7 +448,7 @@ TEST_CASE("Connectivity", "[mesh]")
for (const auto& named_mesh : mesh_list) {
SECTION(named_mesh.name())
{
- auto mesh_1d = named_mesh.mesh()->get<Mesh<Connectivity<1>>>();
+ auto mesh_1d = named_mesh.mesh()->get<Mesh<1>>();
const Connectivity<1>& connectivity = mesh_1d->connectivity();
auto is_boundary_node = connectivity.isBoundaryNode();
@@ -497,7 +497,7 @@ TEST_CASE("Connectivity", "[mesh]")
for (const auto& named_mesh : mesh_list) {
SECTION(named_mesh.name())
{
- auto mesh_2d = named_mesh.mesh()->get<Mesh<Connectivity<2>>>();
+ auto mesh_2d = named_mesh.mesh()->get<Mesh<2>>();
const Connectivity<2>& connectivity = mesh_2d->connectivity();
SECTION("faces/edges")
@@ -575,7 +575,7 @@ TEST_CASE("Connectivity", "[mesh]")
for (const auto& named_mesh : mesh_list) {
SECTION(named_mesh.name())
{
- auto mesh_3d = named_mesh.mesh()->get<Mesh<Connectivity<3>>>();
+ auto mesh_3d = named_mesh.mesh()->get<Mesh<3>>();
const Connectivity<3>& connectivity = mesh_3d->connectivity();
@@ -683,7 +683,7 @@ TEST_CASE("Connectivity", "[mesh]")
{
SECTION("cell -> nodes")
{
- auto mesh = named_mesh.mesh()->get<Mesh<Connectivity<1>>>();
+ auto mesh = named_mesh.mesh()->get<Mesh<1>>();
auto xr = mesh->xr();
const Connectivity<1>& connectivity = mesh->connectivity();
@@ -702,7 +702,7 @@ TEST_CASE("Connectivity", "[mesh]")
SECTION("node -> cells")
{
- auto mesh = named_mesh.mesh()->get<Mesh<Connectivity<1>>>();
+ auto mesh = named_mesh.mesh()->get<Mesh<1>>();
const Connectivity<1>& connectivity = mesh->connectivity();
@@ -731,7 +731,7 @@ TEST_CASE("Connectivity", "[mesh]")
{
SECTION("face -> nodes")
{
- auto mesh = named_mesh.mesh()->get<Mesh<Connectivity<2>>>();
+ auto mesh = named_mesh.mesh()->get<Mesh<2>>();
const Connectivity<2>& connectivity = mesh->connectivity();
@@ -751,7 +751,7 @@ TEST_CASE("Connectivity", "[mesh]")
SECTION("node -> faces")
{
- auto mesh = named_mesh.mesh()->get<Mesh<Connectivity<2>>>();
+ auto mesh = named_mesh.mesh()->get<Mesh<2>>();
const Connectivity<2>& connectivity = mesh->connectivity();
@@ -770,7 +770,7 @@ TEST_CASE("Connectivity", "[mesh]")
SECTION("node -> cells")
{
- auto mesh = named_mesh.mesh()->get<Mesh<Connectivity<2>>>();
+ auto mesh = named_mesh.mesh()->get<Mesh<2>>();
const Connectivity<2>& connectivity = mesh->connectivity();
@@ -789,7 +789,7 @@ TEST_CASE("Connectivity", "[mesh]")
SECTION("face -> cells")
{
- auto mesh = named_mesh.mesh()->get<Mesh<Connectivity<2>>>();
+ auto mesh = named_mesh.mesh()->get<Mesh<2>>();
const Connectivity<2>& connectivity = mesh->connectivity();
@@ -818,7 +818,7 @@ TEST_CASE("Connectivity", "[mesh]")
{
SECTION("edge -> nodes")
{
- auto mesh = named_mesh.mesh()->get<Mesh<Connectivity<3>>>();
+ auto mesh = named_mesh.mesh()->get<Mesh<3>>();
const Connectivity<3>& connectivity = mesh->connectivity();
@@ -838,7 +838,7 @@ TEST_CASE("Connectivity", "[mesh]")
SECTION("face -> nodes")
{
- auto mesh = named_mesh.mesh()->get<Mesh<Connectivity<3>>>();
+ auto mesh = named_mesh.mesh()->get<Mesh<3>>();
const Connectivity<3>& connectivity = mesh->connectivity();
@@ -865,7 +865,7 @@ TEST_CASE("Connectivity", "[mesh]")
SECTION("node -> edges")
{
- auto mesh = named_mesh.mesh()->get<Mesh<Connectivity<3>>>();
+ auto mesh = named_mesh.mesh()->get<Mesh<3>>();
const Connectivity<3>& connectivity = mesh->connectivity();
@@ -887,7 +887,7 @@ TEST_CASE("Connectivity", "[mesh]")
SECTION("node -> faces")
{
- auto mesh = named_mesh.mesh()->get<Mesh<Connectivity<3>>>();
+ auto mesh = named_mesh.mesh()->get<Mesh<3>>();
const Connectivity<3>& connectivity = mesh->connectivity();
@@ -906,7 +906,7 @@ TEST_CASE("Connectivity", "[mesh]")
SECTION("node -> cells")
{
- auto mesh = named_mesh.mesh()->get<Mesh<Connectivity<3>>>();
+ auto mesh = named_mesh.mesh()->get<Mesh<3>>();
const Connectivity<3>& connectivity = mesh->connectivity();
@@ -925,7 +925,7 @@ TEST_CASE("Connectivity", "[mesh]")
SECTION("edge -> faces")
{
- auto mesh = named_mesh.mesh()->get<Mesh<Connectivity<3>>>();
+ auto mesh = named_mesh.mesh()->get<Mesh<3>>();
const Connectivity<3>& connectivity = mesh->connectivity();
@@ -944,7 +944,7 @@ TEST_CASE("Connectivity", "[mesh]")
SECTION("edge -> cells")
{
- auto mesh = named_mesh.mesh()->get<Mesh<Connectivity<3>>>();
+ auto mesh = named_mesh.mesh()->get<Mesh<3>>();
const Connectivity<3>& connectivity = mesh->connectivity();
@@ -963,7 +963,7 @@ TEST_CASE("Connectivity", "[mesh]")
SECTION("face -> cells")
{
- auto mesh = named_mesh.mesh()->get<Mesh<Connectivity<3>>>();
+ auto mesh = named_mesh.mesh()->get<Mesh<3>>();
const Connectivity<3>& connectivity = mesh->connectivity();
@@ -1020,7 +1020,7 @@ TEST_CASE("Connectivity", "[mesh]")
for (const auto& named_mesh : mesh_list) {
SECTION(named_mesh.name())
{
- auto mesh = named_mesh.mesh()->get<Mesh<Connectivity<1>>>();
+ auto mesh = named_mesh.mesh()->get<Mesh<1>>();
const Connectivity<1>& connectivity = mesh->connectivity();
SECTION("node <-> cell")
@@ -1092,7 +1092,7 @@ TEST_CASE("Connectivity", "[mesh]")
for (const auto& named_mesh : mesh_list) {
SECTION(named_mesh.name())
{
- auto mesh = named_mesh.mesh()->get<Mesh<Connectivity<2>>>();
+ auto mesh = named_mesh.mesh()->get<Mesh<2>>();
const Connectivity<2>& connectivity = mesh->connectivity();
SECTION("node <-> cell")
@@ -1194,7 +1194,7 @@ TEST_CASE("Connectivity", "[mesh]")
for (const auto& named_mesh : mesh_list) {
SECTION(named_mesh.name())
{
- auto mesh = named_mesh.mesh()->get<Mesh<Connectivity<3>>>();
+ auto mesh = named_mesh.mesh()->get<Mesh<3>>();
const Connectivity<3>& connectivity = mesh->connectivity();
SECTION("node <-> cell")
diff --git a/tests/test_DiamondDualConnectivityBuilder.cpp b/tests/test_DiamondDualConnectivityBuilder.cpp
index 5c9df535a56dcbfa24741cc662c8bac5c7b78274..2a5d1d4b6f72ff39666536eaa0146c5e80a3029b 100644
--- a/tests/test_DiamondDualConnectivityBuilder.cpp
+++ b/tests/test_DiamondDualConnectivityBuilder.cpp
@@ -37,8 +37,8 @@ TEST_CASE("DiamondDualConnectivityBuilder", "[mesh]")
{
constexpr static size_t Dimension = 2;
- using ConnectivityType = Connectivity<Dimension>;
- using MeshType = Mesh<ConnectivityType>;
+ using MeshType = Mesh<Dimension>;
+ using ConnectivityType = typename MeshType::Connectivity;
std::shared_ptr<const MeshType> mesh = MeshDataBaseForTests::get().hybrid2DMesh()->get<MeshType>();
const ConnectivityType& primal_connectivity = mesh->connectivity();
@@ -147,8 +147,8 @@ TEST_CASE("DiamondDualConnectivityBuilder", "[mesh]")
{
constexpr static size_t Dimension = 3;
- using ConnectivityType = Connectivity<Dimension>;
- using MeshType = Mesh<ConnectivityType>;
+ using MeshType = Mesh<Dimension>;
+ using ConnectivityType = typename MeshType::Connectivity;
std::shared_ptr<const MeshType> mesh = MeshDataBaseForTests::get().hybrid3DMesh()->get<MeshType>();
const ConnectivityType& primal_connectivity = mesh->connectivity();
diff --git a/tests/test_DiamondDualMeshBuilder.cpp b/tests/test_DiamondDualMeshBuilder.cpp
index debdf09f5a4a034745ce2377c9dde2592547152e..c1f95caab35d9ffa2944f8effb5f9ab7fd598770 100644
--- a/tests/test_DiamondDualMeshBuilder.cpp
+++ b/tests/test_DiamondDualMeshBuilder.cpp
@@ -18,18 +18,17 @@ TEST_CASE("DiamondDualMeshBuilder", "[mesh]")
{
constexpr static size_t Dimension = 2;
- using ConnectivityType = Connectivity<Dimension>;
- using MeshType = Mesh<ConnectivityType>;
+ using MeshType = Mesh<Dimension>;
std::shared_ptr primal_mesh_v = MeshDataBaseForTests::get().hybrid2DMesh();
- const MeshType& primal_mesh = *primal_mesh_v->get<Mesh<Connectivity<Dimension>>>();
+ const MeshType& primal_mesh = *primal_mesh_v->get<Mesh<Dimension>>();
REQUIRE(primal_mesh.numberOfNodes() == 53);
REQUIRE(primal_mesh.numberOfFaces() == 110);
REQUIRE(primal_mesh.numberOfCells() == 58);
std::shared_ptr diamond_dual_mesh_v = DualMeshManager::instance().getDiamondDualMesh(primal_mesh_v);
- const MeshType& dual_mesh = *diamond_dual_mesh_v->get<Mesh<Connectivity<Dimension>>>();
+ const MeshType& dual_mesh = *diamond_dual_mesh_v->get<Mesh<Dimension>>();
REQUIRE(dual_mesh.numberOfNodes() == 111);
REQUIRE(dual_mesh.numberOfFaces() == 220);
@@ -85,11 +84,10 @@ TEST_CASE("DiamondDualMeshBuilder", "[mesh]")
{
constexpr static size_t Dimension = 3;
- using ConnectivityType = Connectivity<Dimension>;
- using MeshType = Mesh<ConnectivityType>;
+ using MeshType = Mesh<Dimension>;
std::shared_ptr primal_mesh_v = MeshDataBaseForTests::get().hybrid3DMesh();
- const MeshType& primal_mesh = *primal_mesh_v->get<Mesh<Connectivity<Dimension>>>();
+ const MeshType& primal_mesh = *primal_mesh_v->get<Mesh<Dimension>>();
REQUIRE(primal_mesh.numberOfNodes() == 132);
REQUIRE(primal_mesh.numberOfEdges() == 452);
@@ -97,7 +95,7 @@ TEST_CASE("DiamondDualMeshBuilder", "[mesh]")
REQUIRE(primal_mesh.numberOfCells() == 199);
std::shared_ptr diamond_dual_mesh_v = DualMeshManager::instance().getDiamondDualMesh(primal_mesh_v);
- const MeshType& dual_mesh = *diamond_dual_mesh_v->get<Mesh<Connectivity<Dimension>>>();
+ const MeshType& dual_mesh = *diamond_dual_mesh_v->get<Mesh<Dimension>>();
REQUIRE(dual_mesh.numberOfNodes() == 331);
REQUIRE(dual_mesh.numberOfEdges() == 1461);
diff --git a/tests/test_DiscreteFunctionIntegratorByZone.cpp b/tests/test_DiscreteFunctionIntegratorByZone.cpp
index c556cf6caa623455b238eca7c21463067724936e..2e0d3dad215b2729ae22e0827e3d1b64233be77e 100644
--- a/tests/test_DiscreteFunctionIntegratorByZone.cpp
+++ b/tests/test_DiscreteFunctionIntegratorByZone.cpp
@@ -54,7 +54,7 @@ TEST_CASE("DiscreteFunctionIntegratorByZone", "[scheme]")
std::shared_ptr<const IQuadratureDescriptor> quadrature_descriptor = std::make_shared<GaussQuadratureDescriptor>(3);
auto mesh_1d_v = MeshDataBaseForTests::get().unordered1DMesh();
- auto mesh_1d = mesh_1d_v->get<Mesh<Connectivity<1>>>();
+ auto mesh_1d = mesh_1d_v->get<Mesh<1>>();
std::vector<std::shared_ptr<const IZoneDescriptor>> zone_list;
zone_list.push_back(std::make_shared<NamedZoneDescriptor>("LEFT"));
@@ -384,7 +384,7 @@ let R3x3_non_linear_1d: R^1 -> R^3x3, x -> [[2 * exp(x[0]) * sin(x[0]) + 3, sin(
std::shared_ptr<const IQuadratureDescriptor> quadrature_descriptor = std::make_shared<GaussQuadratureDescriptor>(3);
auto mesh_2d_v = MeshDataBaseForTests::get().hybrid2DMesh();
- auto mesh_2d = mesh_2d_v->get<Mesh<Connectivity<2>>>();
+ auto mesh_2d = mesh_2d_v->get<Mesh<2>>();
std::vector<std::shared_ptr<const IZoneDescriptor>> zone_list;
zone_list.push_back(std::make_shared<NamedZoneDescriptor>("LEFT"));
@@ -714,7 +714,7 @@ let R3x3_non_linear_2d: R^2 -> R^3x3, x -> [[2 * exp(x[0]) * sin(x[1]) + 3, sin(
std::shared_ptr<const IQuadratureDescriptor> quadrature_descriptor = std::make_shared<GaussQuadratureDescriptor>(3);
auto mesh_3d_v = MeshDataBaseForTests::get().hybrid3DMesh();
- auto mesh_3d = mesh_3d_v->get<Mesh<Connectivity<3>>>();
+ auto mesh_3d = mesh_3d_v->get<Mesh<3>>();
std::vector<std::shared_ptr<const IZoneDescriptor>> zone_list;
zone_list.push_back(std::make_shared<NamedZoneDescriptor>("LEFT"));
diff --git a/tests/test_DiscreteFunctionInterpoler.cpp b/tests/test_DiscreteFunctionInterpoler.cpp
index 21bec01b14dc84883e66daa120a2d108bb344d7f..b4c024d3f0f768abb51b34b9772af1e832e7dc5c 100644
--- a/tests/test_DiscreteFunctionInterpoler.cpp
+++ b/tests/test_DiscreteFunctionInterpoler.cpp
@@ -53,7 +53,7 @@ TEST_CASE("DiscreteFunctionInterpoler", "[scheme]")
SECTION(named_mesh.name())
{
auto mesh_1d_v = named_mesh.mesh();
- auto mesh_1d = mesh_1d_v->get<Mesh<Connectivity<Dimension>>>();
+ auto mesh_1d = mesh_1d_v->get<Mesh<Dimension>>();
auto xj = MeshDataManager::instance().getMeshData(*mesh_1d).xj();
@@ -343,7 +343,7 @@ let R3x3_non_linear_1d: R^1 -> R^3x3, x -> [[2 * exp(x[0]) * sin(x[0]) + 3, sin(
SECTION(named_mesh.name())
{
auto mesh_2d_v = named_mesh.mesh();
- auto mesh_2d = mesh_2d_v->get<Mesh<Connectivity<2>>>();
+ auto mesh_2d = mesh_2d_v->get<Mesh<2>>();
auto xj = MeshDataManager::instance().getMeshData(*mesh_2d).xj();
@@ -634,7 +634,7 @@ let R3x3_non_linear_2d: R^2 -> R^3x3, x -> [[2 * exp(x[0]) * sin(x[1]) + 3, sin(
SECTION(named_mesh.name())
{
auto mesh_3d_v = named_mesh.mesh();
- auto mesh_3d = mesh_3d_v->get<Mesh<Connectivity<3>>>();
+ auto mesh_3d = mesh_3d_v->get<Mesh<3>>();
auto xj = MeshDataManager::instance().getMeshData(*mesh_3d).xj();
diff --git a/tests/test_DiscreteFunctionInterpolerByZone.cpp b/tests/test_DiscreteFunctionInterpolerByZone.cpp
index 9a95d7ffd6ff2598bc4c7b9ae3b40287de064cb3..e62ab89fe8951e19d63b1c529237e9dc6ae1cef8 100644
--- a/tests/test_DiscreteFunctionInterpolerByZone.cpp
+++ b/tests/test_DiscreteFunctionInterpolerByZone.cpp
@@ -50,7 +50,7 @@ TEST_CASE("DiscreteFunctionInterpolerByZone", "[scheme]")
constexpr size_t Dimension = 1;
auto mesh_1d_v = MeshDataBaseForTests::get().unordered1DMesh();
- auto mesh_1d = mesh_1d_v->get<Mesh<Connectivity<1>>>();
+ auto mesh_1d = mesh_1d_v->get<Mesh<1>>();
std::vector<std::shared_ptr<const IZoneDescriptor>> zone_list;
zone_list.push_back(std::make_shared<NamedZoneDescriptor>("LEFT"));
@@ -384,7 +384,7 @@ let R3x3_non_linear_1d: R^1 -> R^3x3, x -> [[2 * exp(x[0]) * sin(x[0]) + 3, sin(
constexpr size_t Dimension = 2;
auto mesh_2d_v = MeshDataBaseForTests::get().hybrid2DMesh();
- auto mesh_2d = mesh_2d_v->get<Mesh<Connectivity<2>>>();
+ auto mesh_2d = mesh_2d_v->get<Mesh<2>>();
std::vector<std::shared_ptr<const IZoneDescriptor>> zone_list;
zone_list.push_back(std::make_shared<NamedZoneDescriptor>("LEFT"));
@@ -719,7 +719,7 @@ let R3x3_non_linear_2d: R^2 -> R^3x3, x -> [[2 * exp(x[0]) * sin(x[1]) + 3, sin(
constexpr size_t Dimension = 3;
auto mesh_3d_v = MeshDataBaseForTests::get().hybrid3DMesh();
- auto mesh_3d = mesh_3d_v->get<Mesh<Connectivity<3>>>();
+ auto mesh_3d = mesh_3d_v->get<Mesh<3>>();
std::vector<std::shared_ptr<const IZoneDescriptor>> zone_list;
zone_list.push_back(std::make_shared<NamedZoneDescriptor>("LEFT"));
diff --git a/tests/test_DiscreteFunctionP0.cpp b/tests/test_DiscreteFunctionP0.cpp
index ccf962754e0586d25f995767f73148d3f85a5914..e43e1958de2be4d3ebdc11cecf52d503387ed429 100644
--- a/tests/test_DiscreteFunctionP0.cpp
+++ b/tests/test_DiscreteFunctionP0.cpp
@@ -33,7 +33,7 @@ TEST_CASE("DiscreteFunctionP0", "[scheme]")
for (const auto& named_mesh : mesh_list) {
SECTION(named_mesh.name())
{
- auto mesh = named_mesh.mesh()->get<Mesh<Connectivity<1>>>();
+ auto mesh = named_mesh.mesh()->get<Mesh<1>>();
DiscreteFunctionP0<Dimension, double> f{mesh};
REQUIRE(f.dataType() == ASTNodeDataType::double_t);
REQUIRE(f.descriptor().type() == DiscreteFunctionType::P0);
@@ -81,7 +81,7 @@ TEST_CASE("DiscreteFunctionP0", "[scheme]")
for (const auto& named_mesh : mesh_list) {
SECTION(named_mesh.name())
{
- auto mesh = named_mesh.mesh()->get<Mesh<Connectivity<2>>>();
+ auto mesh = named_mesh.mesh()->get<Mesh<2>>();
DiscreteFunctionP0<Dimension, double> f{mesh};
REQUIRE(f.dataType() == ASTNodeDataType::double_t);
REQUIRE(f.descriptor().type() == DiscreteFunctionType::P0);
@@ -129,7 +129,7 @@ TEST_CASE("DiscreteFunctionP0", "[scheme]")
for (const auto& named_mesh : mesh_list) {
SECTION(named_mesh.name())
{
- auto mesh = named_mesh.mesh()->get<Mesh<Connectivity<3>>>();
+ auto mesh = named_mesh.mesh()->get<Mesh<3>>();
DiscreteFunctionP0<Dimension, double> f{mesh};
REQUIRE(f.dataType() == ASTNodeDataType::double_t);
@@ -191,7 +191,7 @@ TEST_CASE("DiscreteFunctionP0", "[scheme]")
for (const auto& named_mesh : mesh_list) {
SECTION(named_mesh.name())
{
- auto mesh = named_mesh.mesh()->get<Mesh<Connectivity<1>>>();
+ auto mesh = named_mesh.mesh()->get<Mesh<1>>();
DiscreteFunctionP0<Dimension, double> f{mesh};
f.fill(3);
@@ -220,7 +220,7 @@ TEST_CASE("DiscreteFunctionP0", "[scheme]")
for (const auto& named_mesh : mesh_list) {
SECTION(named_mesh.name())
{
- auto mesh = named_mesh.mesh()->get<Mesh<Connectivity<2>>>();
+ auto mesh = named_mesh.mesh()->get<Mesh<2>>();
DiscreteFunctionP0<Dimension, double> f{mesh};
f.fill(3);
@@ -249,7 +249,7 @@ TEST_CASE("DiscreteFunctionP0", "[scheme]")
for (const auto& named_mesh : mesh_list) {
SECTION(named_mesh.name())
{
- auto mesh = named_mesh.mesh()->get<Mesh<Connectivity<3>>>();
+ auto mesh = named_mesh.mesh()->get<Mesh<3>>();
DiscreteFunctionP0<Dimension, double> f{mesh};
f.fill(3);
@@ -291,7 +291,7 @@ TEST_CASE("DiscreteFunctionP0", "[scheme]")
for (const auto& named_mesh : mesh_list) {
SECTION(named_mesh.name())
{
- auto mesh = named_mesh.mesh()->get<Mesh<Connectivity<1>>>();
+ auto mesh = named_mesh.mesh()->get<Mesh<1>>();
SECTION("scalar")
{
@@ -393,7 +393,7 @@ TEST_CASE("DiscreteFunctionP0", "[scheme]")
for (const auto& named_mesh : mesh_list) {
SECTION(named_mesh.name())
{
- auto mesh = named_mesh.mesh()->get<Mesh<Connectivity<2>>>();
+ auto mesh = named_mesh.mesh()->get<Mesh<2>>();
SECTION("scalar")
{
@@ -495,7 +495,7 @@ TEST_CASE("DiscreteFunctionP0", "[scheme]")
for (const auto& named_mesh : mesh_list) {
SECTION(named_mesh.name())
{
- auto mesh = named_mesh.mesh()->get<Mesh<Connectivity<3>>>();
+ auto mesh = named_mesh.mesh()->get<Mesh<3>>();
SECTION("scalar")
{
@@ -599,7 +599,7 @@ TEST_CASE("DiscreteFunctionP0", "[scheme]")
for (const auto& named_mesh : mesh_list) {
SECTION(named_mesh.name())
{
- auto mesh = named_mesh.mesh()->get<Mesh<Connectivity<1>>>();
+ auto mesh = named_mesh.mesh()->get<Mesh<1>>();
auto xj = MeshDataManager::instance().getMeshData(*mesh).xj();
@@ -684,7 +684,7 @@ TEST_CASE("DiscreteFunctionP0", "[scheme]")
for (const auto& named_mesh : mesh_list) {
SECTION(named_mesh.name())
{
- auto mesh = named_mesh.mesh()->get<Mesh<Connectivity<1>>>();
+ auto mesh = named_mesh.mesh()->get<Mesh<1>>();
auto xj = MeshDataManager::instance().getMeshData(*mesh).xj();
@@ -1265,7 +1265,7 @@ TEST_CASE("DiscreteFunctionP0", "[scheme]")
for (const auto& named_mesh : mesh_list) {
SECTION(named_mesh.name())
{
- auto mesh = named_mesh.mesh()->get<Mesh<Connectivity<2>>>();
+ auto mesh = named_mesh.mesh()->get<Mesh<2>>();
auto xj = MeshDataManager::instance().getMeshData(*mesh).xj();
@@ -1851,7 +1851,7 @@ TEST_CASE("DiscreteFunctionP0", "[scheme]")
for (const auto& named_mesh : mesh_list) {
SECTION(named_mesh.name())
{
- auto mesh = named_mesh.mesh()->get<Mesh<Connectivity<3>>>();
+ auto mesh = named_mesh.mesh()->get<Mesh<3>>();
auto xj = MeshDataManager::instance().getMeshData(*mesh).xj();
@@ -2500,7 +2500,7 @@ TEST_CASE("DiscreteFunctionP0", "[scheme]")
for (const auto& named_mesh : mesh_list) {
SECTION(named_mesh.name())
{
- auto mesh = named_mesh.mesh()->get<Mesh<Connectivity<1>>>();
+ auto mesh = named_mesh.mesh()->get<Mesh<1>>();
auto xj = MeshDataManager::instance().getMeshData(*mesh).xj();
@@ -2844,7 +2844,7 @@ TEST_CASE("DiscreteFunctionP0", "[scheme]")
for (const auto& named_mesh : mesh_list) {
SECTION(named_mesh.name())
{
- auto mesh = named_mesh.mesh()->get<Mesh<Connectivity<2>>>();
+ auto mesh = named_mesh.mesh()->get<Mesh<2>>();
auto xj = MeshDataManager::instance().getMeshData(*mesh).xj();
@@ -3188,7 +3188,7 @@ TEST_CASE("DiscreteFunctionP0", "[scheme]")
for (const auto& named_mesh : mesh_list) {
SECTION(named_mesh.name())
{
- auto mesh = named_mesh.mesh()->get<Mesh<Connectivity<3>>>();
+ auto mesh = named_mesh.mesh()->get<Mesh<3>>();
auto xj = MeshDataManager::instance().getMeshData(*mesh).xj();
@@ -3439,7 +3439,7 @@ TEST_CASE("DiscreteFunctionP0", "[scheme]")
mesh->numberOfCells(), PUGS_LAMBDA(const CellId cell_id) {
const double x = xj[cell_id][0];
Ah[cell_id] = TinyMatrix<2>{x + 1, 2 * x - 3, //
- -0.2 * x - 1, 2 + x};
+ -0.2 * x - 1, 2 + x};
});
{
@@ -3461,7 +3461,7 @@ TEST_CASE("DiscreteFunctionP0", "[scheme]")
mesh->numberOfCells(), PUGS_LAMBDA(const CellId cell_id) {
const double x = xj[cell_id][0];
Ah[cell_id] = TinyMatrix<2>{x + 1, 2 * x - 3, //
- -0.2 * x - 1, 2 + x};
+ -0.2 * x - 1, 2 + x};
});
{
@@ -3483,7 +3483,7 @@ TEST_CASE("DiscreteFunctionP0", "[scheme]")
mesh->numberOfCells(), PUGS_LAMBDA(const CellId cell_id) {
const double x = xj[cell_id][0];
Ah[cell_id] = TinyMatrix<2>{x + 1, 2 * x - 3, //
- -0.2 * x - 1, 2 + x};
+ -0.2 * x - 1, 2 + x};
});
{
@@ -3505,7 +3505,7 @@ TEST_CASE("DiscreteFunctionP0", "[scheme]")
mesh->numberOfCells(), PUGS_LAMBDA(const CellId cell_id) {
const double x = xj[cell_id][0];
Ah[cell_id] = TinyMatrix<2>{x + 1, 2 * x - 3, //
- -0.2 * x - 1, 2 + x};
+ -0.2 * x - 1, 2 + x};
});
{
@@ -3626,10 +3626,9 @@ TEST_CASE("DiscreteFunctionP0", "[scheme]")
for (const auto& named_mesh : mesh_list) {
SECTION(named_mesh.name())
{
- auto mesh_1 = named_mesh.mesh()->get<Mesh<Connectivity<1>>>();
+ auto mesh_1 = named_mesh.mesh()->get<Mesh<1>>();
- std::shared_ptr mesh_2 =
- std::make_shared<Mesh<Connectivity<Dimension>>>(mesh_1->shared_connectivity(), mesh_1->xr());
+ std::shared_ptr mesh_2 = std::make_shared<Mesh<Dimension>>(mesh_1->shared_connectivity(), mesh_1->xr());
DiscreteFunctionP0<Dimension, double> f1{mesh_1};
DiscreteFunctionP0<Dimension, double> f2{mesh_2};
@@ -3653,10 +3652,9 @@ TEST_CASE("DiscreteFunctionP0", "[scheme]")
for (const auto& named_mesh : mesh_list) {
SECTION(named_mesh.name())
{
- auto mesh_1 = named_mesh.mesh()->get<Mesh<Connectivity<2>>>();
+ auto mesh_1 = named_mesh.mesh()->get<Mesh<2>>();
- std::shared_ptr mesh_2 =
- std::make_shared<Mesh<Connectivity<Dimension>>>(mesh_1->shared_connectivity(), mesh_1->xr());
+ std::shared_ptr mesh_2 = std::make_shared<Mesh<Dimension>>(mesh_1->shared_connectivity(), mesh_1->xr());
DiscreteFunctionP0<Dimension, double> f1{mesh_1};
DiscreteFunctionP0<Dimension, double> f2{mesh_2};
@@ -3680,10 +3678,9 @@ TEST_CASE("DiscreteFunctionP0", "[scheme]")
for (const auto& named_mesh : mesh_list) {
SECTION(named_mesh.name())
{
- auto mesh_1 = named_mesh.mesh()->get<Mesh<Connectivity<3>>>();
+ auto mesh_1 = named_mesh.mesh()->get<Mesh<3>>();
- std::shared_ptr mesh_2 =
- std::make_shared<Mesh<Connectivity<Dimension>>>(mesh_1->shared_connectivity(), mesh_1->xr());
+ std::shared_ptr mesh_2 = std::make_shared<Mesh<Dimension>>(mesh_1->shared_connectivity(), mesh_1->xr());
DiscreteFunctionP0<Dimension, double> f1{mesh_1};
DiscreteFunctionP0<Dimension, double> f2{mesh_2};
diff --git a/tests/test_DiscreteFunctionP0Vector.cpp b/tests/test_DiscreteFunctionP0Vector.cpp
index 1bf27ddb2648ac3d231e15bae4eca9250109787c..4521e8fa88b9767411c2ff7003f8daecce13eb06 100644
--- a/tests/test_DiscreteFunctionP0Vector.cpp
+++ b/tests/test_DiscreteFunctionP0Vector.cpp
@@ -34,7 +34,7 @@ TEST_CASE("DiscreteFunctionP0Vector", "[scheme]")
for (const auto& named_mesh : mesh_list) {
SECTION(named_mesh.name())
{
- auto mesh = named_mesh.mesh()->get<Mesh<Connectivity<Dimension>>>();
+ auto mesh = named_mesh.mesh()->get<Mesh<Dimension>>();
DiscreteFunctionP0Vector<Dimension, double> f{mesh, size};
REQUIRE(f.dataType() == ASTNodeDataType::double_t);
REQUIRE(f.descriptor().type() == DiscreteFunctionType::P0Vector);
@@ -82,7 +82,7 @@ TEST_CASE("DiscreteFunctionP0Vector", "[scheme]")
for (const auto& named_mesh : mesh_list) {
SECTION(named_mesh.name())
{
- auto mesh = named_mesh.mesh()->get<Mesh<Connectivity<Dimension>>>();
+ auto mesh = named_mesh.mesh()->get<Mesh<Dimension>>();
DiscreteFunctionP0Vector<Dimension, double> f{mesh, size};
REQUIRE(f.dataType() == ASTNodeDataType::double_t);
REQUIRE(f.descriptor().type() == DiscreteFunctionType::P0Vector);
@@ -130,7 +130,7 @@ TEST_CASE("DiscreteFunctionP0Vector", "[scheme]")
for (const auto& named_mesh : mesh_list) {
SECTION(named_mesh.name())
{
- auto mesh = named_mesh.mesh()->get<Mesh<Connectivity<Dimension>>>();
+ auto mesh = named_mesh.mesh()->get<Mesh<Dimension>>();
DiscreteFunctionP0Vector<Dimension, double> f{mesh, size};
REQUIRE(f.dataType() == ASTNodeDataType::double_t);
REQUIRE(f.descriptor().type() == DiscreteFunctionType::P0Vector);
@@ -193,7 +193,7 @@ TEST_CASE("DiscreteFunctionP0Vector", "[scheme]")
SECTION(named_mesh.name())
{
constexpr size_t Dimension = 1;
- auto mesh = named_mesh.mesh()->get<Mesh<Connectivity<Dimension>>>();
+ auto mesh = named_mesh.mesh()->get<Mesh<Dimension>>();
DiscreteFunctionP0Vector<Dimension, double> f{mesh, size};
f.fill(3);
@@ -214,7 +214,7 @@ TEST_CASE("DiscreteFunctionP0Vector", "[scheme]")
for (const auto& named_mesh : mesh_list) {
SECTION(named_mesh.name())
{
- auto mesh = named_mesh.mesh()->get<Mesh<Connectivity<Dimension>>>();
+ auto mesh = named_mesh.mesh()->get<Mesh<Dimension>>();
DiscreteFunctionP0Vector<Dimension, double> f{mesh, size};
f.fill(2.3);
@@ -234,7 +234,7 @@ TEST_CASE("DiscreteFunctionP0Vector", "[scheme]")
for (const auto& named_mesh : mesh_list) {
SECTION(named_mesh.name())
{
- auto mesh = named_mesh.mesh()->get<Mesh<Connectivity<Dimension>>>();
+ auto mesh = named_mesh.mesh()->get<Mesh<Dimension>>();
DiscreteFunctionP0Vector<Dimension, double> f{mesh, size};
f.fill(3.2);
@@ -269,7 +269,7 @@ TEST_CASE("DiscreteFunctionP0Vector", "[scheme]")
for (const auto& named_mesh : mesh_list) {
SECTION(named_mesh.name())
{
- auto mesh = named_mesh.mesh()->get<Mesh<Connectivity<Dimension>>>();
+ auto mesh = named_mesh.mesh()->get<Mesh<Dimension>>();
const size_t size = 3;
const size_t value = parallel::rank() + 1;
@@ -315,7 +315,7 @@ TEST_CASE("DiscreteFunctionP0Vector", "[scheme]")
for (const auto& named_mesh : mesh_list) {
SECTION(named_mesh.name())
{
- auto mesh = named_mesh.mesh()->get<Mesh<Connectivity<Dimension>>>();
+ auto mesh = named_mesh.mesh()->get<Mesh<Dimension>>();
const size_t size = 3;
const size_t value = parallel::rank() + 1;
@@ -362,7 +362,7 @@ TEST_CASE("DiscreteFunctionP0Vector", "[scheme]")
for (const auto& named_mesh : mesh_list) {
SECTION(named_mesh.name())
{
- auto mesh = named_mesh.mesh()->get<Mesh<Connectivity<Dimension>>>();
+ auto mesh = named_mesh.mesh()->get<Mesh<Dimension>>();
const size_t size = 3;
const size_t value = parallel::rank() + 1;
@@ -413,7 +413,7 @@ TEST_CASE("DiscreteFunctionP0Vector", "[scheme]")
for (const auto& named_mesh : mesh_list) {
SECTION(named_mesh.name())
{
- auto mesh = named_mesh.mesh()->get<Mesh<Connectivity<Dimension>>>();
+ auto mesh = named_mesh.mesh()->get<Mesh<Dimension>>();
auto xj = MeshDataManager::instance().getMeshData(*mesh).xj();
@@ -456,7 +456,7 @@ TEST_CASE("DiscreteFunctionP0Vector", "[scheme]")
for (const auto& named_mesh : mesh_list) {
SECTION(named_mesh.name())
{
- auto mesh = named_mesh.mesh()->get<Mesh<Connectivity<Dimension>>>();
+ auto mesh = named_mesh.mesh()->get<Mesh<Dimension>>();
auto xj = MeshDataManager::instance().getMeshData(*mesh).xj();
@@ -500,7 +500,7 @@ TEST_CASE("DiscreteFunctionP0Vector", "[scheme]")
for (const auto& named_mesh : mesh_list) {
SECTION(named_mesh.name())
{
- auto mesh = named_mesh.mesh()->get<Mesh<Connectivity<Dimension>>>();
+ auto mesh = named_mesh.mesh()->get<Mesh<Dimension>>();
auto xj = MeshDataManager::instance().getMeshData(*mesh).xj();
@@ -548,7 +548,7 @@ TEST_CASE("DiscreteFunctionP0Vector", "[scheme]")
for (const auto& named_mesh : mesh_list) {
SECTION(named_mesh.name())
{
- auto mesh = named_mesh.mesh()->get<Mesh<Connectivity<Dimension>>>();
+ auto mesh = named_mesh.mesh()->get<Mesh<Dimension>>();
auto xj = MeshDataManager::instance().getMeshData(*mesh).xj();
@@ -631,7 +631,7 @@ TEST_CASE("DiscreteFunctionP0Vector", "[scheme]")
for (const auto& named_mesh : mesh_list) {
SECTION(named_mesh.name())
{
- auto mesh = named_mesh.mesh()->get<Mesh<Connectivity<Dimension>>>();
+ auto mesh = named_mesh.mesh()->get<Mesh<Dimension>>();
auto xj = MeshDataManager::instance().getMeshData(*mesh).xj();
@@ -714,7 +714,7 @@ TEST_CASE("DiscreteFunctionP0Vector", "[scheme]")
for (const auto& named_mesh : mesh_list) {
SECTION(named_mesh.name())
{
- auto mesh = named_mesh.mesh()->get<Mesh<Connectivity<Dimension>>>();
+ auto mesh = named_mesh.mesh()->get<Mesh<Dimension>>();
auto xj = MeshDataManager::instance().getMeshData(*mesh).xj();
@@ -800,7 +800,7 @@ TEST_CASE("DiscreteFunctionP0Vector", "[scheme]")
for (const auto& named_mesh : mesh_list) {
SECTION(named_mesh.name())
{
- auto mesh = named_mesh.mesh()->get<Mesh<Connectivity<Dimension>>>();
+ auto mesh = named_mesh.mesh()->get<Mesh<Dimension>>();
auto xj = MeshDataManager::instance().getMeshData(*mesh).xj();
@@ -934,7 +934,7 @@ TEST_CASE("DiscreteFunctionP0Vector", "[scheme]")
for (const auto& named_mesh : mesh_list) {
SECTION(named_mesh.name())
{
- auto mesh = named_mesh.mesh()->get<Mesh<Connectivity<Dimension>>>();
+ auto mesh = named_mesh.mesh()->get<Mesh<Dimension>>();
auto xj = MeshDataManager::instance().getMeshData(*mesh).xj();
@@ -1072,7 +1072,7 @@ TEST_CASE("DiscreteFunctionP0Vector", "[scheme]")
for (const auto& named_mesh : mesh_list) {
SECTION(named_mesh.name())
{
- auto mesh = named_mesh.mesh()->get<Mesh<Connectivity<Dimension>>>();
+ auto mesh = named_mesh.mesh()->get<Mesh<Dimension>>();
auto xj = MeshDataManager::instance().getMeshData(*mesh).xj();
diff --git a/tests/test_DiscreteFunctionUtils.cpp b/tests/test_DiscreteFunctionUtils.cpp
index e7b547f40d9de7a7b1c63c7e1a0c9ae2b4a49e88..f7adede0eeea5ec70035c8a1ad92609a9a4235fb 100644
--- a/tests/test_DiscreteFunctionUtils.cpp
+++ b/tests/test_DiscreteFunctionUtils.cpp
@@ -22,7 +22,7 @@ TEST_CASE("DiscreteFunctionUtils", "[scheme]")
SECTION(named_mesh.name())
{
auto mesh_v = named_mesh.mesh();
- auto mesh = mesh_v->get<Mesh<Connectivity<Dimension>>>();
+ auto mesh = mesh_v->get<Mesh<Dimension>>();
std::shared_ptr mesh_copy =
std::make_shared<const std::decay_t<decltype(*mesh)>>(mesh->shared_connectivity(), mesh->xr());
@@ -210,7 +210,7 @@ TEST_CASE("DiscreteFunctionUtils", "[scheme]")
SECTION(named_mesh.name())
{
auto mesh_v = named_mesh.mesh();
- auto mesh = mesh_v->get<Mesh<Connectivity<Dimension>>>();
+ auto mesh = mesh_v->get<Mesh<Dimension>>();
std::shared_ptr mesh_copy =
std::make_shared<const std::decay_t<decltype(*mesh)>>(mesh->shared_connectivity(), mesh->xr());
std::shared_ptr mesh_copy_v = std::make_shared<const MeshVariant>(mesh_copy);
@@ -397,7 +397,7 @@ TEST_CASE("DiscreteFunctionUtils", "[scheme]")
SECTION(named_mesh.name())
{
auto mesh_v = named_mesh.mesh();
- auto mesh = mesh_v->get<Mesh<Connectivity<3>>>();
+ auto mesh = mesh_v->get<Mesh<3>>();
std::shared_ptr mesh_copy =
std::make_shared<const std::decay_t<decltype(*mesh)>>(mesh->shared_connectivity(), mesh->xr());
@@ -586,12 +586,10 @@ TEST_CASE("DiscreteFunctionUtils", "[scheme]")
for (const auto& named_mesh : mesh_list) {
SECTION(named_mesh.name())
{
- auto mesh = named_mesh.mesh()->get<Mesh<Connectivity<Dimension>>>();
+ auto mesh = named_mesh.mesh()->get<Mesh<Dimension>>();
std::shared_ptr other_mesh =
- CartesianMeshBuilder{TinyVector<1>{-1}, TinyVector<1>{3}, TinyVector<1, size_t>{19}}
- .mesh()
- ->get<Mesh<Connectivity<1>>>();
+ CartesianMeshBuilder{TinyVector<1>{-1}, TinyVector<1>{3}, TinyVector<1, size_t>{19}}.mesh()->get<Mesh<1>>();
std::shared_ptr other_mesh_v = std::make_shared<const MeshVariant>(other_mesh);
std::shared_ptr uh = std::make_shared<DiscreteFunctionVariant>(DiscreteFunctionP0<Dimension, double>(mesh));
@@ -605,9 +603,9 @@ TEST_CASE("DiscreteFunctionUtils", "[scheme]")
{
constexpr size_t Dimension = 1;
- std::shared_ptr mesh_1d = MeshDataBaseForTests::get().cartesian1DMesh()->get<Mesh<Connectivity<1>>>();
- std::shared_ptr mesh_2d_v = std::make_shared<const MeshVariant>(
- MeshDataBaseForTests::get().cartesian2DMesh()->get<Mesh<Connectivity<2>>>());
+ std::shared_ptr mesh_1d = MeshDataBaseForTests::get().cartesian1DMesh()->get<Mesh<1>>();
+ std::shared_ptr mesh_2d_v =
+ std::make_shared<const MeshVariant>(MeshDataBaseForTests::get().cartesian2DMesh()->get<Mesh<2>>());
std::shared_ptr uh = std::make_shared<DiscreteFunctionVariant>(DiscreteFunctionP0<Dimension, double>(mesh_1d));
diff --git a/tests/test_DiscreteFunctionVectorIntegrator.cpp b/tests/test_DiscreteFunctionVectorIntegrator.cpp
index 79cd612e5c754752bcf2092813d963b3a3c5cb0a..58ad9cdf21b983138b66566978b754f5b128de5a 100644
--- a/tests/test_DiscreteFunctionVectorIntegrator.cpp
+++ b/tests/test_DiscreteFunctionVectorIntegrator.cpp
@@ -68,7 +68,7 @@ TEST_CASE("DiscreteFunctionVectorIntegrator", "[scheme]")
SECTION(named_mesh.name())
{
auto mesh_1d_v = named_mesh.mesh();
- auto mesh_1d = mesh_1d_v->get<Mesh<Connectivity<Dimension>>>();
+ auto mesh_1d = mesh_1d_v->get<Mesh<Dimension>>();
std::string_view data = R"(
import math;
@@ -158,7 +158,7 @@ let R_scalar_non_linear_1d: R^1 -> R, x -> 2 * exp(x[0]) + 3;
SECTION(named_mesh.name())
{
auto mesh_2d_v = named_mesh.mesh();
- auto mesh_2d = mesh_2d_v->get<Mesh<Connectivity<Dimension>>>();
+ auto mesh_2d = mesh_2d_v->get<Mesh<Dimension>>();
std::string_view data = R"(
import math;
@@ -248,7 +248,7 @@ let R_scalar_non_linear_2d: R^2 -> R, x -> 2 * exp(x[0] + x[1]) + 3;
SECTION(named_mesh.name())
{
auto mesh_3d_v = named_mesh.mesh();
- auto mesh_3d = mesh_3d_v->get<Mesh<Connectivity<3>>>();
+ auto mesh_3d = mesh_3d_v->get<Mesh<3>>();
std::string_view data = R"(
import math;
@@ -336,7 +336,7 @@ let R_scalar_non_linear_3d: R^3 -> R, x -> 2 * exp(x[0] + x[1]) + 3 * x[2];
SECTION(named_mesh.name())
{
auto mesh_3d_v = named_mesh.mesh();
- auto mesh_3d = mesh_3d_v->get<Mesh<Connectivity<3>>>();
+ auto mesh_3d = mesh_3d_v->get<Mesh<3>>();
std::string_view data = R"(
import math;
diff --git a/tests/test_DiscreteFunctionVectorIntegratorByZone.cpp b/tests/test_DiscreteFunctionVectorIntegratorByZone.cpp
index f2672177013526d6070ff03d2ba7aca537f1838a..2522f365166666cd0d1dcbc6f87df1480a1a9825 100644
--- a/tests/test_DiscreteFunctionVectorIntegratorByZone.cpp
+++ b/tests/test_DiscreteFunctionVectorIntegratorByZone.cpp
@@ -65,7 +65,7 @@ TEST_CASE("DiscreteFunctionVectorIntegratorByZone", "[scheme]")
std::shared_ptr<const IQuadratureDescriptor> quadrature_descriptor = std::make_shared<GaussQuadratureDescriptor>(3);
auto mesh_1d_v = MeshDataBaseForTests::get().unordered1DMesh();
- auto mesh_1d = mesh_1d_v->get<Mesh<Connectivity<Dimension>>>();
+ auto mesh_1d = mesh_1d_v->get<Mesh<Dimension>>();
std::vector<std::shared_ptr<const IZoneDescriptor>> zone_list;
zone_list.push_back(std::make_shared<NamedZoneDescriptor>("LEFT"));
@@ -194,7 +194,7 @@ let R_scalar_non_linear_1d: R^1 -> R, x -> 2 * exp(x[0]) + 3;
std::shared_ptr<const IQuadratureDescriptor> quadrature_descriptor = std::make_shared<GaussQuadratureDescriptor>(3);
auto mesh_2d_v = MeshDataBaseForTests::get().hybrid2DMesh();
- auto mesh_2d = mesh_2d_v->get<Mesh<Connectivity<Dimension>>>();
+ auto mesh_2d = mesh_2d_v->get<Mesh<Dimension>>();
std::vector<std::shared_ptr<const IZoneDescriptor>> zone_list;
zone_list.push_back(std::make_shared<NamedZoneDescriptor>("LEFT"));
@@ -323,7 +323,7 @@ let R_scalar_non_linear_2d: R^2 -> R, x -> 2 * exp(x[0] + x[1]) + 3;
std::shared_ptr<const IQuadratureDescriptor> quadrature_descriptor = std::make_shared<GaussQuadratureDescriptor>(3);
auto mesh_3d_v = MeshDataBaseForTests::get().hybrid3DMesh();
- auto mesh_3d = mesh_3d_v->get<Mesh<Connectivity<Dimension>>>();
+ auto mesh_3d = mesh_3d_v->get<Mesh<Dimension>>();
std::vector<std::shared_ptr<const IZoneDescriptor>> zone_list;
zone_list.push_back(std::make_shared<NamedZoneDescriptor>("LEFT"));
diff --git a/tests/test_DiscreteFunctionVectorInterpoler.cpp b/tests/test_DiscreteFunctionVectorInterpoler.cpp
index 0c0d60795834f8518a1306d48e62e072f7e6140b..b86f833a2f19c172d51a5f818db67270271396a2 100644
--- a/tests/test_DiscreteFunctionVectorInterpoler.cpp
+++ b/tests/test_DiscreteFunctionVectorInterpoler.cpp
@@ -63,7 +63,7 @@ TEST_CASE("DiscreteFunctionVectorInterpoler", "[scheme]")
SECTION(named_mesh.name())
{
auto mesh_1d_v = named_mesh.mesh();
- auto mesh_1d = mesh_1d_v->get<Mesh<Connectivity<Dimension>>>();
+ auto mesh_1d = mesh_1d_v->get<Mesh<Dimension>>();
auto xj = MeshDataManager::instance().getMeshData(*mesh_1d).xj();
@@ -168,7 +168,7 @@ let R_scalar_non_linear_1d: R^1 -> R, x -> 2 * exp(x[0]) + 3;
SECTION(named_mesh.name())
{
auto mesh_2d_v = named_mesh.mesh();
- auto mesh_2d = mesh_2d_v->get<Mesh<Connectivity<Dimension>>>();
+ auto mesh_2d = mesh_2d_v->get<Mesh<Dimension>>();
auto xj = MeshDataManager::instance().getMeshData(*mesh_2d).xj();
@@ -273,7 +273,7 @@ let R_scalar_non_linear_2d: R^2 -> R, x -> 2 * exp(x[0] + x[1]) + 3;
SECTION(named_mesh.name())
{
auto mesh_3d_v = named_mesh.mesh();
- auto mesh_3d = mesh_3d_v->get<Mesh<Connectivity<Dimension>>>();
+ auto mesh_3d = mesh_3d_v->get<Mesh<Dimension>>();
auto xj = MeshDataManager::instance().getMeshData(*mesh_3d).xj();
@@ -376,7 +376,7 @@ let R_scalar_non_linear_3d: R^3 -> R, x -> 2 * exp(x[0] + x[1]) + 3 * x[2];
SECTION(named_mesh.name())
{
auto mesh_3d_v = named_mesh.mesh();
- auto mesh_3d = mesh_3d_v->get<Mesh<Connectivity<3>>>();
+ auto mesh_3d = mesh_3d_v->get<Mesh<3>>();
auto xj = MeshDataManager::instance().getMeshData(*mesh_3d).xj();
diff --git a/tests/test_DiscreteFunctionVectorInterpolerByZone.cpp b/tests/test_DiscreteFunctionVectorInterpolerByZone.cpp
index 7468ce99706446d1aaa6f2a8b5647afcb3ab8cb0..c17a9b397bd9570cbbda7e22f9f2bebbbc838b4d 100644
--- a/tests/test_DiscreteFunctionVectorInterpolerByZone.cpp
+++ b/tests/test_DiscreteFunctionVectorInterpolerByZone.cpp
@@ -60,7 +60,7 @@ TEST_CASE("DiscreteFunctionVectorInterpolerByZone", "[scheme]")
constexpr size_t Dimension = 1;
auto mesh_1d_v = MeshDataBaseForTests::get().unordered1DMesh();
- auto mesh_1d = mesh_1d_v->get<Mesh<Connectivity<Dimension>>>();
+ auto mesh_1d = mesh_1d_v->get<Mesh<Dimension>>();
std::vector<std::shared_ptr<const IZoneDescriptor>> zone_list;
zone_list.push_back(std::make_shared<NamedZoneDescriptor>("LEFT"));
@@ -186,7 +186,7 @@ let R_scalar_non_linear_1d: R^1 -> R, x -> 2 * exp(x[0]) + 3;
constexpr size_t Dimension = 2;
auto mesh_2d_v = MeshDataBaseForTests::get().hybrid2DMesh();
- auto mesh_2d = mesh_2d_v->get<Mesh<Connectivity<Dimension>>>();
+ auto mesh_2d = mesh_2d_v->get<Mesh<Dimension>>();
std::vector<std::shared_ptr<const IZoneDescriptor>> zone_list;
zone_list.push_back(std::make_shared<NamedZoneDescriptor>("LEFT"));
@@ -312,7 +312,7 @@ let R_scalar_non_linear_2d: R^2 -> R, x -> 2 * exp(x[0] + x[1]) + 3;
constexpr size_t Dimension = 3;
auto mesh_3d_v = MeshDataBaseForTests::get().hybrid3DMesh();
- auto mesh_3d = mesh_3d_v->get<Mesh<Connectivity<Dimension>>>();
+ auto mesh_3d = mesh_3d_v->get<Mesh<Dimension>>();
std::vector<std::shared_ptr<const IZoneDescriptor>> zone_list;
zone_list.push_back(std::make_shared<NamedZoneDescriptor>("LEFT"));
@@ -441,7 +441,7 @@ let R_scalar_non_linear_3d: R^3 -> R, x -> 2 * exp(x[0] + x[1]) + 3 * x[2];
SECTION(named_mesh.name())
{
auto mesh_3d_v = named_mesh.mesh();
- auto mesh_3d = mesh_3d_v->get<Mesh<Connectivity<3>>>();
+ auto mesh_3d = mesh_3d_v->get<Mesh<3>>();
auto xj = MeshDataManager::instance().getMeshData(*mesh_3d).xj();
diff --git a/tests/test_Dual1DConnectivityBuilder.cpp b/tests/test_Dual1DConnectivityBuilder.cpp
index b8eb91e42fa8743d6c161485dabe9f0be0f10aea..d08d53b2295ec74b793159771132542bebdfef1f 100644
--- a/tests/test_Dual1DConnectivityBuilder.cpp
+++ b/tests/test_Dual1DConnectivityBuilder.cpp
@@ -36,8 +36,8 @@ TEST_CASE("Dual1DConnectivityBuilder", "[mesh]")
{
constexpr static size_t Dimension = 1;
- using ConnectivityType = Connectivity<Dimension>;
- using MeshType = Mesh<ConnectivityType>;
+ using MeshType = Mesh<Dimension>;
+ using ConnectivityType = typename MeshType::Connectivity;
std::shared_ptr<const MeshType> mesh = MeshDataBaseForTests::get().unordered1DMesh()->get<MeshType>();
const ConnectivityType& primal_connectivity = mesh->connectivity();
diff --git a/tests/test_Dual1DMeshBuilder.cpp b/tests/test_Dual1DMeshBuilder.cpp
index 50b130f12dad2b9d23e9d9b18a2e0d2bbdc2b6a9..88ac2de4f0858fbc97ea81b59f4025732e8590fe 100644
--- a/tests/test_Dual1DMeshBuilder.cpp
+++ b/tests/test_Dual1DMeshBuilder.cpp
@@ -18,8 +18,7 @@ TEST_CASE("Dual1DMeshBuilder", "[mesh]")
{
constexpr static size_t Dimension = 1;
- using ConnectivityType = Connectivity<Dimension>;
- using MeshType = Mesh<ConnectivityType>;
+ using MeshType = Mesh<Dimension>;
std::shared_ptr primal_mesh_v = MeshDataBaseForTests::get().unordered1DMesh();
diff --git a/tests/test_DualConnectivityManager.cpp b/tests/test_DualConnectivityManager.cpp
index f3b6703902c2dffae69c12ddc593927150729ef6..8c24eb0798631cbf3b4f365a6e955706fe00bab5 100644
--- a/tests/test_DualConnectivityManager.cpp
+++ b/tests/test_DualConnectivityManager.cpp
@@ -11,8 +11,8 @@
TEST_CASE("DualConnectivityManager", "[mesh]")
{
- using ConnectivityType = Connectivity<2>;
- using MeshType = Mesh<ConnectivityType>;
+ using MeshType = Mesh<2>;
+ using ConnectivityType = typename MeshType::Connectivity;
std::shared_ptr mesh_v = MeshDataBaseForTests::get().hybrid2DMesh();
std::shared_ptr mesh = mesh_v->get<MeshType>();
diff --git a/tests/test_DualMeshManager.cpp b/tests/test_DualMeshManager.cpp
index da9aac55748a14bb54241c99cb158f1334891b90..9c60d8a27eda913ded141a8dc0147c4ccdf0f733 100644
--- a/tests/test_DualMeshManager.cpp
+++ b/tests/test_DualMeshManager.cpp
@@ -26,8 +26,7 @@ TEST_CASE("DualMeshManager", "[mesh]")
SECTION("2D")
{
- using ConnectivityType = Connectivity<2>;
- using MeshType = Mesh<ConnectivityType>;
+ using MeshType = Mesh<2>;
std::shared_ptr mesh_v = MeshDataBaseForTests::get().hybrid2DMesh();
@@ -117,8 +116,7 @@ TEST_CASE("DualMeshManager", "[mesh]")
SECTION("3D")
{
- using ConnectivityType = Connectivity<3>;
- using MeshType = Mesh<ConnectivityType>;
+ using MeshType = Mesh<3>;
std::shared_ptr mesh_v = MeshDataBaseForTests::get().hybrid3DMesh();
std::shared_ptr mesh = mesh_v->get<MeshType>();
diff --git a/tests/test_EdgeIntegrator.cpp b/tests/test_EdgeIntegrator.cpp
index b20cb925f8b386296651fbee23d3787999e66c8d..5415a03f10f16b43574f258233586ec3b0dc3312 100644
--- a/tests/test_EdgeIntegrator.cpp
+++ b/tests/test_EdgeIntegrator.cpp
@@ -38,7 +38,7 @@ TEST_CASE("EdgeIntegrator", "[scheme]")
for (const auto& mesh_info : mesh_list) {
auto mesh_name = mesh_info.first;
- auto mesh = mesh_info.second->get<Mesh<Connectivity<3>>>();
+ auto mesh = mesh_info.second->get<Mesh<3>>();
Array<const double> int_f_per_edge = [=] {
Array<double> int_f(mesh->numberOfEdges());
@@ -288,7 +288,7 @@ TEST_CASE("EdgeIntegrator", "[scheme]")
for (const auto& mesh_info : mesh_list) {
auto mesh_name = mesh_info.first;
- auto mesh = mesh_info.second->get<Mesh<Connectivity<3>>>();
+ auto mesh = mesh_info.second->get<Mesh<3>>();
Array<const R2> int_f_per_edge = [=] {
Array<R2> int_f(mesh->numberOfEdges());
@@ -540,7 +540,7 @@ TEST_CASE("EdgeIntegrator", "[scheme]")
for (const auto& mesh_info : mesh_list) {
auto mesh_name = mesh_info.first;
- auto mesh = mesh_info.second->get<Mesh<Connectivity<3>>>();
+ auto mesh = mesh_info.second->get<Mesh<3>>();
Array<const R2x2> int_f_per_edge = [=] {
Array<R2x2> int_f(mesh->numberOfEdges());
diff --git a/tests/test_EmbeddedDiscreteFunctionMathFunctions1D.cpp b/tests/test_EmbeddedDiscreteFunctionMathFunctions1D.cpp
index aa48d3ffc1143017a6b174efc82242e8f60a7e08..7ebb75f85d0307b2645aead879c529f381b4cd77 100644
--- a/tests/test_EmbeddedDiscreteFunctionMathFunctions1D.cpp
+++ b/tests/test_EmbeddedDiscreteFunctionMathFunctions1D.cpp
@@ -34,10 +34,9 @@ TEST_CASE("EmbeddedDiscreteFunctionVariantMathFunctions1D", "[scheme]")
for (const auto& named_mesh : mesh_list) {
SECTION(named_mesh.name())
{
- auto mesh = named_mesh.mesh()->get<Mesh<Connectivity<Dimension>>>();
+ auto mesh = named_mesh.mesh()->get<Mesh<Dimension>>();
- std::shared_ptr other_mesh =
- std::make_shared<Mesh<Connectivity<Dimension>>>(mesh->shared_connectivity(), mesh->xr());
+ std::shared_ptr other_mesh = std::make_shared<Mesh<Dimension>>(mesh->shared_connectivity(), mesh->xr());
CellValue<const Rd> xj = MeshDataManager::instance().getMeshData(*mesh).xj();
diff --git a/tests/test_EmbeddedDiscreteFunctionMathFunctions2D.cpp b/tests/test_EmbeddedDiscreteFunctionMathFunctions2D.cpp
index 8658aa3da0bc66bc8a06d2e65ac94d67140e4495..989b66d39227ef33e3be242916a5d590c1305cd6 100644
--- a/tests/test_EmbeddedDiscreteFunctionMathFunctions2D.cpp
+++ b/tests/test_EmbeddedDiscreteFunctionMathFunctions2D.cpp
@@ -36,10 +36,9 @@ TEST_CASE("EmbeddedDiscreteFunctionVariantMathFunctions2D", "[scheme]")
for (const auto& named_mesh : mesh_list) {
SECTION(named_mesh.name())
{
- auto mesh = named_mesh.mesh()->get<Mesh<Connectivity<Dimension>>>();
+ auto mesh = named_mesh.mesh()->get<Mesh<Dimension>>();
- std::shared_ptr other_mesh =
- std::make_shared<Mesh<Connectivity<Dimension>>>(mesh->shared_connectivity(), mesh->xr());
+ std::shared_ptr other_mesh = std::make_shared<Mesh<Dimension>>(mesh->shared_connectivity(), mesh->xr());
CellValue<const Rd> xj = MeshDataManager::instance().getMeshData(*mesh).xj();
diff --git a/tests/test_EmbeddedDiscreteFunctionMathFunctions3D.cpp b/tests/test_EmbeddedDiscreteFunctionMathFunctions3D.cpp
index 415173393dfeeb91570bfefb1d98c1535400dde6..25afe5c0a10512fff3c43bc6c78cd28cab7330a7 100644
--- a/tests/test_EmbeddedDiscreteFunctionMathFunctions3D.cpp
+++ b/tests/test_EmbeddedDiscreteFunctionMathFunctions3D.cpp
@@ -34,10 +34,9 @@ TEST_CASE("EmbeddedDiscreteFunctionVariantMathFunctions3D", "[scheme]")
for (const auto& named_mesh : mesh_list) {
SECTION(named_mesh.name())
{
- auto mesh = named_mesh.mesh()->get<Mesh<Connectivity<Dimension>>>();
+ auto mesh = named_mesh.mesh()->get<Mesh<Dimension>>();
- std::shared_ptr other_mesh =
- std::make_shared<Mesh<Connectivity<Dimension>>>(mesh->shared_connectivity(), mesh->xr());
+ std::shared_ptr other_mesh = std::make_shared<Mesh<Dimension>>(mesh->shared_connectivity(), mesh->xr());
CellValue<const Rd> xj = MeshDataManager::instance().getMeshData(*mesh).xj();
diff --git a/tests/test_EmbeddedDiscreteFunctionOperators1D.cpp b/tests/test_EmbeddedDiscreteFunctionOperators1D.cpp
index fbbcc8bc14ada19d5272d8c3c5459e24bb965c56..e790e4608fc58dbb35cfc36ec81d780a88a4513a 100644
--- a/tests/test_EmbeddedDiscreteFunctionOperators1D.cpp
+++ b/tests/test_EmbeddedDiscreteFunctionOperators1D.cpp
@@ -25,10 +25,9 @@ TEST_CASE("EmbeddedDiscreteFunctionOperators1D", "[scheme]")
for (const auto& named_mesh : mesh_list) {
SECTION(named_mesh.name())
{
- auto mesh = named_mesh.mesh()->get<Mesh<Connectivity<Dimension>>>();
+ auto mesh = named_mesh.mesh()->get<Mesh<Dimension>>();
- std::shared_ptr other_mesh =
- std::make_shared<Mesh<Connectivity<Dimension>>>(mesh->shared_connectivity(), mesh->xr());
+ std::shared_ptr other_mesh = std::make_shared<Mesh<Dimension>>(mesh->shared_connectivity(), mesh->xr());
CellValue<const Rd> xj = MeshDataManager::instance().getMeshData(*mesh).xj();
diff --git a/tests/test_EmbeddedDiscreteFunctionOperators2D.cpp b/tests/test_EmbeddedDiscreteFunctionOperators2D.cpp
index d3fece0d8cbad2c949d0e96213f01fab27a2713c..0efe0cfbe5efdac30f11b6ac52d6c4a1030b55b3 100644
--- a/tests/test_EmbeddedDiscreteFunctionOperators2D.cpp
+++ b/tests/test_EmbeddedDiscreteFunctionOperators2D.cpp
@@ -25,10 +25,9 @@ TEST_CASE("EmbeddedDiscreteFunctionOperators2D", "[scheme]")
for (const auto& named_mesh : mesh_list) {
SECTION(named_mesh.name())
{
- auto mesh = named_mesh.mesh()->get<Mesh<Connectivity<Dimension>>>();
+ auto mesh = named_mesh.mesh()->get<Mesh<Dimension>>();
- std::shared_ptr other_mesh =
- std::make_shared<Mesh<Connectivity<Dimension>>>(mesh->shared_connectivity(), mesh->xr());
+ std::shared_ptr other_mesh = std::make_shared<Mesh<Dimension>>(mesh->shared_connectivity(), mesh->xr());
CellValue<const Rd> xj = MeshDataManager::instance().getMeshData(*mesh).xj();
diff --git a/tests/test_EmbeddedDiscreteFunctionOperators3D.cpp b/tests/test_EmbeddedDiscreteFunctionOperators3D.cpp
index eb30d18d5e54f3e1bc50703deaa805675190e4f1..a024c2dde6f69143b74c87324805cc193e40c6a0 100644
--- a/tests/test_EmbeddedDiscreteFunctionOperators3D.cpp
+++ b/tests/test_EmbeddedDiscreteFunctionOperators3D.cpp
@@ -25,10 +25,9 @@ TEST_CASE("EmbeddedDiscreteFunctionOperators3D", "[scheme]")
for (const auto& named_mesh : mesh_list) {
SECTION(named_mesh.name())
{
- auto mesh = named_mesh.mesh()->get<Mesh<Connectivity<Dimension>>>();
+ auto mesh = named_mesh.mesh()->get<Mesh<Dimension>>();
- std::shared_ptr other_mesh =
- std::make_shared<Mesh<Connectivity<Dimension>>>(mesh->shared_connectivity(), mesh->xr());
+ std::shared_ptr other_mesh = std::make_shared<Mesh<Dimension>>(mesh->shared_connectivity(), mesh->xr());
CellValue<const Rd> xj = MeshDataManager::instance().getMeshData(*mesh).xj();
diff --git a/tests/test_EmbeddedDiscreteFunctionUtils.cpp b/tests/test_EmbeddedDiscreteFunctionUtils.cpp
index 3640ea5b5a2dcf3ea692e60d7904b4d4042079ec..464313cb54d55ac4a237533955f115b4288cc5fd 100644
--- a/tests/test_EmbeddedDiscreteFunctionUtils.cpp
+++ b/tests/test_EmbeddedDiscreteFunctionUtils.cpp
@@ -34,7 +34,7 @@ TEST_CASE("EmbeddedDiscreteFunctionUtils", "[language]")
for (const auto& named_mesh : mesh_list) {
SECTION(named_mesh.name())
{
- auto mesh_1d = named_mesh.mesh()->get<Mesh<Connectivity<1>>>();
+ auto mesh_1d = named_mesh.mesh()->get<Mesh<1>>();
REQUIRE(EmbeddedDiscreteFunctionUtils::getOperandTypeName(DiscreteFunctionP0<1, double>{mesh_1d}) ==
"Vh(P0:R)");
@@ -63,7 +63,7 @@ TEST_CASE("EmbeddedDiscreteFunctionUtils", "[language]")
for (const auto& named_mesh : mesh_list) {
SECTION(named_mesh.name())
{
- auto mesh_1d = named_mesh.mesh()->get<Mesh<Connectivity<1>>>();
+ auto mesh_1d = named_mesh.mesh()->get<Mesh<1>>();
REQUIRE(EmbeddedDiscreteFunctionUtils::getOperandTypeName(DiscreteFunctionP0Vector<1, double>{mesh_1d, 2}) ==
"Vh(P0Vector:R)");
diff --git a/tests/test_FaceIntegrator.cpp b/tests/test_FaceIntegrator.cpp
index 092ab55c3b683963004c52724694b594d81fdbaa..73e6ecad94efed88c5c4c7b50ad535a573851beb 100644
--- a/tests/test_FaceIntegrator.cpp
+++ b/tests/test_FaceIntegrator.cpp
@@ -22,7 +22,7 @@ TEST_CASE("FaceIntegrator", "[scheme]")
using R2 = TinyVector<2>;
const auto mesh_v = MeshDataBaseForTests::get().hybrid2DMesh();
- const auto mesh = mesh_v->get<Mesh<Connectivity<2>>>();
+ const auto mesh = mesh_v->get<Mesh<2>>();
auto f = [](const R2& X) -> double {
const double x = X[0];
@@ -265,7 +265,7 @@ TEST_CASE("FaceIntegrator", "[scheme]")
for (const auto& mesh_info : mesh_list) {
auto mesh_name = mesh_info.first;
- auto mesh = mesh_info.second->get<Mesh<Connectivity<3>>>();
+ auto mesh = mesh_info.second->get<Mesh<3>>();
Array<const double> int_f_per_face = [=] {
Array<double> int_f(mesh->numberOfFaces());
@@ -514,7 +514,7 @@ TEST_CASE("FaceIntegrator", "[scheme]")
{
using R2 = TinyVector<2>;
- const auto mesh = MeshDataBaseForTests::get().hybrid2DMesh()->get<Mesh<Connectivity<2>>>();
+ const auto mesh = MeshDataBaseForTests::get().hybrid2DMesh()->get<Mesh<2>>();
auto f = [](const R2& X) -> R2 {
const double x = X[0];
@@ -757,7 +757,7 @@ TEST_CASE("FaceIntegrator", "[scheme]")
for (const auto& mesh_info : mesh_list) {
auto mesh_name = mesh_info.first;
- auto mesh = mesh_info.second->get<Mesh<Connectivity<3>>>();
+ auto mesh = mesh_info.second->get<Mesh<3>>();
Array<const R2> int_f_per_face = [=] {
Array<R2> int_f(mesh->numberOfFaces());
@@ -1009,7 +1009,7 @@ TEST_CASE("FaceIntegrator", "[scheme]")
{
using R2 = TinyVector<2>;
- const auto mesh = MeshDataBaseForTests::get().hybrid2DMesh()->get<Mesh<Connectivity<2>>>();
+ const auto mesh = MeshDataBaseForTests::get().hybrid2DMesh()->get<Mesh<2>>();
auto f = [](const R2& X) -> R2x2 {
const double x = X[0];
@@ -1252,7 +1252,7 @@ TEST_CASE("FaceIntegrator", "[scheme]")
for (const auto& mesh_info : mesh_list) {
auto mesh_name = mesh_info.first;
- auto mesh = mesh_info.second->get<Mesh<Connectivity<3>>>();
+ auto mesh = mesh_info.second->get<Mesh<3>>();
Array<const R2x2> int_f_per_face = [=] {
Array<R2x2> int_f(mesh->numberOfFaces());
diff --git a/tests/test_IntegrateCellArray.cpp b/tests/test_IntegrateCellArray.cpp
index 6e45d92a74e0a8cf02b60216d41111b55cac3c35..738fad05fa74ff45820a5f4071dff9a91abbc42b 100644
--- a/tests/test_IntegrateCellArray.cpp
+++ b/tests/test_IntegrateCellArray.cpp
@@ -55,7 +55,7 @@ TEST_CASE("IntegrateCellArray", "[language]")
SECTION(named_mesh.name())
{
auto mesh_1d_v = named_mesh.mesh();
- auto mesh_1d = mesh_1d_v->get<Mesh<Connectivity<1>>>();
+ auto mesh_1d = mesh_1d_v->get<Mesh<1>>();
std::string_view data = R"(
import math;
@@ -142,7 +142,7 @@ let g: R^1 -> R, x -> 2 * exp(x[0]) + 3;
SECTION(named_mesh.name())
{
auto mesh_2d_v = named_mesh.mesh();
- auto mesh_2d = mesh_2d_v->get<Mesh<Connectivity<2>>>();
+ auto mesh_2d = mesh_2d_v->get<Mesh<2>>();
std::string_view data = R"(
import math;
@@ -240,7 +240,7 @@ let g: R^2 -> R, x -> 2*exp(x[0])*sin(x[1])+3;
SECTION(named_mesh.name())
{
auto mesh_3d_v = named_mesh.mesh();
- auto mesh_3d = mesh_3d_v->get<Mesh<Connectivity<3>>>();
+ auto mesh_3d = mesh_3d_v->get<Mesh<3>>();
std::string_view data = R"(
import math;
@@ -343,7 +343,7 @@ let g: R^3 -> R, x -> 2 * exp(x[0]) * sin(x[1]) * x[2] + 3;
SECTION(named_mesh.name())
{
auto mesh_1d_v = named_mesh.mesh();
- auto mesh_1d = mesh_1d_v->get<Mesh<Connectivity<1>>>();
+ auto mesh_1d = mesh_1d_v->get<Mesh<1>>();
Array<CellId> cell_list{mesh_1d->numberOfCells() / 2 + mesh_1d->numberOfCells() % 2};
@@ -438,7 +438,7 @@ let g: R^1 -> R, x -> 2 * exp(x[0]) + 3;
SECTION(named_mesh.name())
{
auto mesh_2d_v = named_mesh.mesh();
- auto mesh_2d = mesh_2d_v->get<Mesh<Connectivity<2>>>();
+ auto mesh_2d = mesh_2d_v->get<Mesh<2>>();
Array<CellId> cell_list{mesh_2d->numberOfCells() / 2 + mesh_2d->numberOfCells() % 2};
@@ -545,7 +545,7 @@ let g: R^2 -> R, x -> 2*exp(x[0])*sin(x[1])+3;
SECTION(named_mesh.name())
{
auto mesh_3d_v = named_mesh.mesh();
- auto mesh_3d = mesh_3d_v->get<Mesh<Connectivity<3>>>();
+ auto mesh_3d = mesh_3d_v->get<Mesh<3>>();
Array<CellId> cell_list{mesh_3d->numberOfCells() / 2 + mesh_3d->numberOfCells() % 2};
diff --git a/tests/test_IntegrateCellValue.cpp b/tests/test_IntegrateCellValue.cpp
index 88d4e356da9bf364557769e1f851dd22d2a1fa8e..ab2c86f82c7b2b10db0a54555d70a06558460468 100644
--- a/tests/test_IntegrateCellValue.cpp
+++ b/tests/test_IntegrateCellValue.cpp
@@ -53,7 +53,7 @@ TEST_CASE("IntegrateCellValue", "[language]")
SECTION(named_mesh.name())
{
auto mesh_1d_v = named_mesh.mesh();
- auto mesh_1d = mesh_1d_v->get<Mesh<Connectivity<1>>>();
+ auto mesh_1d = mesh_1d_v->get<Mesh<1>>();
std::string_view data = R"(
import math;
@@ -111,7 +111,7 @@ let R2x2_1d: R^1 -> R^2x2, x -> [[2 * exp(x[0]) * sin(x[0]) + 3, sin(x[0] - 2 *
SECTION(named_mesh.name())
{
auto mesh_2d_v = named_mesh.mesh();
- auto mesh_2d = mesh_2d_v->get<Mesh<Connectivity<2>>>();
+ auto mesh_2d = mesh_2d_v->get<Mesh<2>>();
std::string_view data = R"(
import math;
@@ -181,7 +181,7 @@ let R3_2d: R^2 -> R^3, x -> [2*exp(x[0])*sin(x[1])+3, x[0]-2*x[1], 3];
SECTION(named_mesh.name())
{
auto mesh_3d_v = named_mesh.mesh();
- auto mesh_3d = mesh_3d_v->get<Mesh<Connectivity<3>>>();
+ auto mesh_3d = mesh_3d_v->get<Mesh<3>>();
std::string_view data = R"(
import math;
@@ -251,7 +251,7 @@ let scalar_3d: R^3 -> R, x -> 2 * exp(x[0]) * sin(x[1]) * x[2] + 3;
SECTION(named_mesh.name())
{
auto mesh_1d_v = named_mesh.mesh();
- auto mesh_1d = mesh_1d_v->get<Mesh<Connectivity<1>>>();
+ auto mesh_1d = mesh_1d_v->get<Mesh<1>>();
Array<CellId> cell_list{mesh_1d->numberOfCells() / 2 + mesh_1d->numberOfCells() % 2};
@@ -316,7 +316,7 @@ let scalar_1d: R^1 -> R, x -> 2 * exp(x[0]) + 3;
SECTION(named_mesh.name())
{
auto mesh_2d_v = named_mesh.mesh();
- auto mesh_2d = mesh_2d_v->get<Mesh<Connectivity<2>>>();
+ auto mesh_2d = mesh_2d_v->get<Mesh<2>>();
Array<CellId> cell_list{mesh_2d->numberOfCells() / 2 + mesh_2d->numberOfCells() % 2};
@@ -396,7 +396,7 @@ let R3_2d: R^2 -> R^3, x -> [2*exp(x[0])*sin(x[1])+3, x[0]-2*x[1], 3];
SECTION(named_mesh.name())
{
auto mesh_3d_v = named_mesh.mesh();
- auto mesh_3d = mesh_3d_v->get<Mesh<Connectivity<3>>>();
+ auto mesh_3d = mesh_3d_v->get<Mesh<3>>();
Array<CellId> cell_list{mesh_3d->numberOfCells() / 2 + mesh_3d->numberOfCells() % 2};
diff --git a/tests/test_IntegrateOnCells.cpp b/tests/test_IntegrateOnCells.cpp
index f2f46e73dc33c569c7469ce2bf9a5dbaa928e43e..93d1c6806dda456aa17b7f138c6d1f23d8e33315 100644
--- a/tests/test_IntegrateOnCells.cpp
+++ b/tests/test_IntegrateOnCells.cpp
@@ -55,7 +55,7 @@ TEST_CASE("IntegrateOnCells", "[language]")
SECTION(named_mesh.name())
{
auto mesh_1d_v = named_mesh.mesh();
- auto mesh_1d = mesh_1d_v->get<Mesh<Connectivity<1>>>();
+ auto mesh_1d = mesh_1d_v->get<Mesh<1>>();
std::string_view data = R"(
import math;
@@ -156,7 +156,7 @@ let R2x2_1d: R^1 -> R^2x2, x -> [[2 * exp(x[0]) * sin(x[0]) + 3, sin(x[0] - 2 *
SECTION(named_mesh.name())
{
auto mesh_2d_v = named_mesh.mesh();
- auto mesh_2d = mesh_2d_v->get<Mesh<Connectivity<2>>>();
+ auto mesh_2d = mesh_2d_v->get<Mesh<2>>();
std::string_view data = R"(
import math;
@@ -270,7 +270,7 @@ let R2x2_2d: R^2 -> R^2x2, x -> [[2*exp(x[0])*sin(x[1])+3, sin(x[0]-2*x[1])], [3
SECTION(named_mesh.name())
{
auto mesh_3d_v = named_mesh.mesh();
- auto mesh_3d = mesh_3d_v->get<Mesh<Connectivity<3>>>();
+ auto mesh_3d = mesh_3d_v->get<Mesh<3>>();
std::string_view data = R"(
import math;
@@ -389,7 +389,7 @@ let R2x2_3d: R^3 -> R^2x2, x -> [[2 * exp(x[0]) * sin(x[1]) + 3 * cos(x[2]), sin
SECTION(named_mesh.name())
{
auto mesh_1d_v = named_mesh.mesh();
- auto mesh_1d = mesh_1d_v->get<Mesh<Connectivity<1>>>();
+ auto mesh_1d = mesh_1d_v->get<Mesh<1>>();
Array<CellId> cell_list{mesh_1d->numberOfCells() / 2 + mesh_1d->numberOfCells() % 2};
@@ -500,7 +500,7 @@ let R2x2_1d: R^1 -> R^2x2, x -> [[2 * exp(x[0]) * sin(x[0]) + 3, sin(x[0] - 2 *
SECTION(named_mesh.name())
{
auto mesh_2d_v = named_mesh.mesh();
- auto mesh_2d = mesh_2d_v->get<Mesh<Connectivity<2>>>();
+ auto mesh_2d = mesh_2d_v->get<Mesh<2>>();
Array<CellId> cell_list{mesh_2d->numberOfCells() / 2 + mesh_2d->numberOfCells() % 2};
@@ -624,7 +624,7 @@ let R2x2_2d: R^2 -> R^2x2, x -> [[2*exp(x[0])*sin(x[1])+3, sin(x[0]-2*x[1])], [3
SECTION(named_mesh.name())
{
auto mesh_3d_v = named_mesh.mesh();
- auto mesh_3d = mesh_3d_v->get<Mesh<Connectivity<3>>>();
+ auto mesh_3d = mesh_3d_v->get<Mesh<3>>();
Array<CellId> cell_list{mesh_3d->numberOfCells() / 2 + mesh_3d->numberOfCells() % 2};
@@ -758,7 +758,7 @@ let R2x2_3d: R^3 -> R^2x2, x -> [[2 * exp(x[0]) * sin(x[1]) + 3 * cos(x[2]), sin
SECTION(named_mesh.name())
{
auto mesh_1d_v = named_mesh.mesh();
- auto mesh_1d = mesh_1d_v->get<Mesh<Connectivity<1>>>();
+ auto mesh_1d = mesh_1d_v->get<Mesh<1>>();
std::string_view data = R"(
import math;
@@ -859,7 +859,7 @@ let R2x2_1d: R^1 -> R^2x2, x -> [[2 * exp(x[0]) * sin(x[0]) + 3, sin(x[0] - 2 *
SECTION(named_mesh.name())
{
auto mesh_2d_v = named_mesh.mesh();
- auto mesh_2d = mesh_2d_v->get<Mesh<Connectivity<2>>>();
+ auto mesh_2d = mesh_2d_v->get<Mesh<2>>();
std::string_view data = R"(
import math;
@@ -973,7 +973,7 @@ let R2x2_2d: R^2 -> R^2x2, x -> [[2*exp(x[0])*sin(x[1])+3, sin(x[0]-2*x[1])], [3
SECTION(named_mesh.name())
{
auto mesh_3d_v = named_mesh.mesh();
- auto mesh_3d = mesh_3d_v->get<Mesh<Connectivity<3>>>();
+ auto mesh_3d = mesh_3d_v->get<Mesh<3>>();
std::string_view data = R"(
import math;
@@ -1092,7 +1092,7 @@ let R2x2_3d: R^3 -> R^2x2, x -> [[2 * exp(x[0]) * sin(x[1]) + 3 * cos(x[2]), sin
SECTION(named_mesh.name())
{
auto mesh_1d_v = named_mesh.mesh();
- auto mesh_1d = mesh_1d_v->get<Mesh<Connectivity<1>>>();
+ auto mesh_1d = mesh_1d_v->get<Mesh<1>>();
Array<CellId> cell_list{mesh_1d->numberOfCells() / 2 + mesh_1d->numberOfCells() % 2};
@@ -1203,7 +1203,7 @@ let R2x2_1d: R^1 -> R^2x2, x -> [[2 * exp(x[0]) * sin(x[0]) + 3, sin(x[0] - 2 *
SECTION(named_mesh.name())
{
auto mesh_2d_v = named_mesh.mesh();
- auto mesh_2d = mesh_2d_v->get<Mesh<Connectivity<2>>>();
+ auto mesh_2d = mesh_2d_v->get<Mesh<2>>();
Array<CellId> cell_list{mesh_2d->numberOfCells() / 2 + mesh_2d->numberOfCells() % 2};
@@ -1327,7 +1327,7 @@ let R2x2_2d: R^2 -> R^2x2, x -> [[2*exp(x[0])*sin(x[1])+3, sin(x[0]-2*x[1])], [3
SECTION(named_mesh.name())
{
auto mesh_3d_v = named_mesh.mesh();
- auto mesh_3d = mesh_3d_v->get<Mesh<Connectivity<3>>>();
+ auto mesh_3d = mesh_3d_v->get<Mesh<3>>();
Array<CellId> cell_list{mesh_3d->numberOfCells() / 2 + mesh_3d->numberOfCells() % 2};
@@ -1461,7 +1461,7 @@ let R2x2_3d: R^3 -> R^2x2, x -> [[2 * exp(x[0]) * sin(x[1]) + 3 * cos(x[2]), sin
SECTION(named_mesh.name())
{
auto mesh_1d_v = named_mesh.mesh();
- auto mesh_1d = mesh_1d_v->get<Mesh<Connectivity<1>>>();
+ auto mesh_1d = mesh_1d_v->get<Mesh<1>>();
std::string_view data = R"(
import math;
@@ -1562,7 +1562,7 @@ let R2x2_1d: R^1 -> R^2x2, x -> [[2 * exp(x[0]) * sin(x[0]) + 3, sin(x[0] - 2 *
SECTION(named_mesh.name())
{
auto mesh_2d_v = named_mesh.mesh();
- auto mesh_2d = mesh_2d_v->get<Mesh<Connectivity<2>>>();
+ auto mesh_2d = mesh_2d_v->get<Mesh<2>>();
std::string_view data = R"(
import math;
@@ -1676,7 +1676,7 @@ let R2x2_2d: R^2 -> R^2x2, x -> [[2*exp(x[0])*sin(x[1])+3, sin(x[0]-2*x[1])], [3
SECTION(named_mesh.name())
{
auto mesh_3d_v = named_mesh.mesh();
- auto mesh_3d = mesh_3d_v->get<Mesh<Connectivity<3>>>();
+ auto mesh_3d = mesh_3d_v->get<Mesh<3>>();
std::string_view data = R"(
import math;
@@ -1795,7 +1795,7 @@ let R2x2_3d: R^3 -> R^2x2, x -> [[2 * exp(x[0]) * sin(x[1]) + 3 * cos(x[2]), sin
SECTION(named_mesh.name())
{
auto mesh_1d_v = named_mesh.mesh();
- auto mesh_1d = mesh_1d_v->get<Mesh<Connectivity<1>>>();
+ auto mesh_1d = mesh_1d_v->get<Mesh<1>>();
Array<CellId> cell_list{mesh_1d->numberOfCells() / 2 + mesh_1d->numberOfCells() % 2};
@@ -1906,7 +1906,7 @@ let R2x2_1d: R^1 -> R^2x2, x -> [[2 * exp(x[0]) * sin(x[0]) + 3, sin(x[0] - 2 *
SECTION(named_mesh.name())
{
auto mesh_2d_v = named_mesh.mesh();
- auto mesh_2d = mesh_2d_v->get<Mesh<Connectivity<2>>>();
+ auto mesh_2d = mesh_2d_v->get<Mesh<2>>();
Array<CellId> cell_list{mesh_2d->numberOfCells() / 2 + mesh_2d->numberOfCells() % 2};
@@ -2030,7 +2030,7 @@ let R2x2_2d: R^2 -> R^2x2, x -> [[2*exp(x[0])*sin(x[1])+3, sin(x[0]-2*x[1])], [3
SECTION(named_mesh.name())
{
auto mesh_3d_v = named_mesh.mesh();
- auto mesh_3d = mesh_3d_v->get<Mesh<Connectivity<3>>>();
+ auto mesh_3d = mesh_3d_v->get<Mesh<3>>();
Array<CellId> cell_list{mesh_3d->numberOfCells() / 2 + mesh_3d->numberOfCells() % 2};
diff --git a/tests/test_InterpolateItemArray.cpp b/tests/test_InterpolateItemArray.cpp
index c23dedeaa3baff962c1c05a606e85edd13587bdd..cab00ae0a3540baba764495e13ce89be9511923c 100644
--- a/tests/test_InterpolateItemArray.cpp
+++ b/tests/test_InterpolateItemArray.cpp
@@ -57,7 +57,7 @@ TEST_CASE("InterpolateItemArray", "[language]")
SECTION(named_mesh.name())
{
auto mesh_1d_v = named_mesh.mesh();
- auto mesh_1d = mesh_1d_v->get<Mesh<Connectivity<Dimension>>>();
+ auto mesh_1d = mesh_1d_v->get<Mesh<Dimension>>();
auto xj = MeshDataManager::instance().getMeshData(*mesh_1d).xj();
@@ -127,7 +127,7 @@ let scalar_non_linear_1d: R^1 -> R, x -> 2 * exp(x[0]) + 3;
SECTION(named_mesh.name())
{
auto mesh_1d_v = named_mesh.mesh();
- auto mesh_1d = mesh_1d_v->get<Mesh<Connectivity<Dimension>>>();
+ auto mesh_1d = mesh_1d_v->get<Mesh<Dimension>>();
auto xj = MeshDataManager::instance().getMeshData(*mesh_1d).xj();
@@ -187,7 +187,7 @@ let f_1d: R^1 -> (R), x -> (2*x[0] + 2, 2 * exp(x[0]) + 3);
SECTION(named_mesh.name())
{
auto mesh_1d_v = named_mesh.mesh();
- auto mesh_1d = mesh_1d_v->get<Mesh<Connectivity<Dimension>>>();
+ auto mesh_1d = mesh_1d_v->get<Mesh<Dimension>>();
auto xj = MeshDataManager::instance().getMeshData(*mesh_1d).xj();
@@ -256,7 +256,7 @@ let f_1d: R^1 -> (R^1), x -> (2*x[0] + 2, [2 * exp(x[0]) + 3]);
SECTION(named_mesh.name())
{
auto mesh_2d_v = named_mesh.mesh();
- auto mesh_2d = mesh_2d_v->get<Mesh<Connectivity<Dimension>>>();
+ auto mesh_2d = mesh_2d_v->get<Mesh<Dimension>>();
auto xj = MeshDataManager::instance().getMeshData(*mesh_2d).xj();
@@ -326,7 +326,7 @@ let scalar_non_linear_2d: R^2 -> R, x -> 2*exp(x[0])*sin(x[1])+3;
SECTION(named_mesh.name())
{
auto mesh_2d_v = named_mesh.mesh();
- auto mesh_2d = mesh_2d_v->get<Mesh<Connectivity<Dimension>>>();
+ auto mesh_2d = mesh_2d_v->get<Mesh<Dimension>>();
auto xj = MeshDataManager::instance().getMeshData(*mesh_2d).xj();
@@ -393,7 +393,7 @@ let f_2d: R^2 -> (R), x -> (2*x[0] + 3*x[1] + 2, 2*exp(x[0])*sin(x[1])+3);
SECTION(named_mesh.name())
{
auto mesh_3d_v = named_mesh.mesh();
- auto mesh_3d = mesh_3d_v->get<Mesh<Connectivity<Dimension>>>();
+ auto mesh_3d = mesh_3d_v->get<Mesh<Dimension>>();
auto xj = MeshDataManager::instance().getMeshData(*mesh_3d).xj();
@@ -463,7 +463,7 @@ let scalar_non_linear_3d: R^3 -> R, x -> 2 * exp(x[0]) * sin(x[1]) * x[2] + 3;
SECTION(named_mesh.name())
{
auto mesh_3d_v = named_mesh.mesh();
- auto mesh_3d = mesh_3d_v->get<Mesh<Connectivity<Dimension>>>();
+ auto mesh_3d = mesh_3d_v->get<Mesh<Dimension>>();
auto xj = MeshDataManager::instance().getMeshData(*mesh_3d).xj();
@@ -544,7 +544,7 @@ let f_3d: R^3 -> (R), x -> (2 * x[0] + 3 * x[1] + 2 * x[2] - 1, 2 * exp(x[0]) *
SECTION(named_mesh.name())
{
auto mesh_1d_v = named_mesh.mesh();
- auto mesh_1d = mesh_1d_v->get<Mesh<Connectivity<Dimension>>>();
+ auto mesh_1d = mesh_1d_v->get<Mesh<Dimension>>();
auto xj = MeshDataManager::instance().getMeshData(*mesh_1d).xj();
@@ -623,7 +623,7 @@ let scalar_non_linear_1d: R^1 -> R, x -> 2 * exp(x[0]) + 3;
SECTION(named_mesh.name())
{
auto mesh_1d_v = named_mesh.mesh();
- auto mesh_1d = mesh_1d_v->get<Mesh<Connectivity<Dimension>>>();
+ auto mesh_1d = mesh_1d_v->get<Mesh<Dimension>>();
auto xj = MeshDataManager::instance().getMeshData(*mesh_1d).xj();
@@ -699,7 +699,7 @@ let f_1d: R^1 -> (R), x -> (2*x[0] + 2, 2 * exp(x[0]) + 3);
SECTION(named_mesh.name())
{
auto mesh_2d_v = named_mesh.mesh();
- auto mesh_2d = mesh_2d_v->get<Mesh<Connectivity<Dimension>>>();
+ auto mesh_2d = mesh_2d_v->get<Mesh<Dimension>>();
auto xj = MeshDataManager::instance().getMeshData(*mesh_2d).xj();
@@ -775,7 +775,7 @@ let scalar_non_linear_2d: R^2 -> R, x -> 2*exp(x[0])*sin(x[1])+3;
SECTION(named_mesh.name())
{
auto mesh_2d_v = named_mesh.mesh();
- auto mesh_2d = mesh_2d_v->get<Mesh<Connectivity<Dimension>>>();
+ auto mesh_2d = mesh_2d_v->get<Mesh<Dimension>>();
auto xj = MeshDataManager::instance().getMeshData(*mesh_2d).xj();
@@ -841,7 +841,7 @@ let f_2d: R^2 -> (R), x -> (2*x[0] + 3*x[1] + 2, 2*exp(x[0])*sin(x[1])+3);
SECTION(named_mesh.name())
{
auto mesh_2d_v = named_mesh.mesh();
- auto mesh_2d = mesh_2d_v->get<Mesh<Connectivity<Dimension>>>();
+ auto mesh_2d = mesh_2d_v->get<Mesh<Dimension>>();
auto xj = MeshDataManager::instance().getMeshData(*mesh_2d).xj();
@@ -919,7 +919,7 @@ let f_2d: R^2 -> (R^2x2), x -> ([[x[0],0],[2-x[1], x[0]*x[1]]], [[2*x[0], x[1]],
SECTION(named_mesh.name())
{
auto mesh_3d_v = named_mesh.mesh();
- auto mesh_3d = mesh_3d_v->get<Mesh<Connectivity<Dimension>>>();
+ auto mesh_3d = mesh_3d_v->get<Mesh<Dimension>>();
auto xj = MeshDataManager::instance().getMeshData(*mesh_3d).xj();
@@ -995,7 +995,7 @@ let scalar_non_linear_3d: R^3 -> R, x -> 2 * exp(x[0]) * sin(x[1]) * x[2] + 3;
SECTION(named_mesh.name())
{
auto mesh_3d_v = named_mesh.mesh();
- auto mesh_3d = mesh_3d_v->get<Mesh<Connectivity<Dimension>>>();
+ auto mesh_3d = mesh_3d_v->get<Mesh<Dimension>>();
auto xj = MeshDataManager::instance().getMeshData(*mesh_3d).xj();
@@ -1061,7 +1061,7 @@ let f_3d: R^3 -> (R), x -> (2 * x[0] + 3 * x[1] + 2 * x[2] - 1, 2 * exp(x[0]) *
SECTION(named_mesh.name())
{
auto mesh_3d_v = named_mesh.mesh();
- auto mesh_3d = mesh_3d_v->get<Mesh<Connectivity<Dimension>>>();
+ auto mesh_3d = mesh_3d_v->get<Mesh<Dimension>>();
auto xj = MeshDataManager::instance().getMeshData(*mesh_3d).xj();
diff --git a/tests/test_InterpolateItemValue.cpp b/tests/test_InterpolateItemValue.cpp
index 8cb3ac4986459404eb94ff92e2204ecd35a152ba..370b8b4015122ce3227c3ee99da177b1b7e20978 100644
--- a/tests/test_InterpolateItemValue.cpp
+++ b/tests/test_InterpolateItemValue.cpp
@@ -52,7 +52,7 @@ TEST_CASE("InterpolateItemValue", "[language]")
SECTION(named_mesh.name())
{
auto mesh_1d_v = named_mesh.mesh();
- auto mesh_1d = mesh_1d_v->get<Mesh<Connectivity<Dimension>>>();
+ auto mesh_1d = mesh_1d_v->get<Mesh<Dimension>>();
auto xj = MeshDataManager::instance().getMeshData(*mesh_1d).xj();
@@ -224,7 +224,7 @@ let R2x2_non_linear_1d: R^1 -> R^2x2, x -> [[2 * exp(x[0]) * sin(x[0]) + 3, sin(
SECTION(named_mesh.name())
{
auto mesh_2d_v = named_mesh.mesh();
- auto mesh_2d = mesh_2d_v->get<Mesh<Connectivity<Dimension>>>();
+ auto mesh_2d = mesh_2d_v->get<Mesh<Dimension>>();
auto xj = MeshDataManager::instance().getMeshData(*mesh_2d).xj();
@@ -390,7 +390,7 @@ let R2x2_non_linear_2d: R^2 -> R^2x2, x -> [[2*exp(x[0])*sin(x[1])+3, sin(x[0]-2
SECTION(named_mesh.name())
{
auto mesh_3d_v = named_mesh.mesh();
- auto mesh_3d = mesh_3d_v->get<Mesh<Connectivity<Dimension>>>();
+ auto mesh_3d = mesh_3d_v->get<Mesh<Dimension>>();
auto xj = MeshDataManager::instance().getMeshData(*mesh_3d).xj();
@@ -536,7 +536,7 @@ let R2x2_non_linear_3d: R^3 -> R^2x2, x -> [[2 * exp(x[0]) * sin(x[1]) + 3 * cos
cell_value.numberOfItems(), PUGS_LAMBDA(const CellId cell_id) {
const TinyVector<Dimension>& x = xj[cell_id];
cell_value[cell_id] = TinyMatrix<2>{2 * exp(x[0]) * sin(x[1]) + 3 * cos(x[2]),
- sin(x[0] - 2 * x[1] * x[2]), 3, x[0] * x[1] * x[2]};
+ sin(x[0] - 2 * x[1] * x[2]), 3, x[0] * x[1] * x[2]};
});
CellValue<const TinyMatrix<2>> interpolate_value =
InterpolateItemValue<TinyMatrix<2>(TinyVector<Dimension>)>::interpolate(function_symbol_id, xj);
@@ -570,7 +570,7 @@ let R2x2_non_linear_3d: R^3 -> R^2x2, x -> [[2 * exp(x[0]) * sin(x[1]) + 3 * cos
SECTION(named_mesh.name())
{
auto mesh_1d_v = named_mesh.mesh();
- auto mesh_1d = mesh_1d_v->get<Mesh<Connectivity<Dimension>>>();
+ auto mesh_1d = mesh_1d_v->get<Mesh<Dimension>>();
auto xj = MeshDataManager::instance().getMeshData(*mesh_1d).xj();
@@ -753,7 +753,7 @@ let R2x2_non_linear_1d: R^1 -> R^2x2, x -> [[2 * exp(x[0]) * sin(x[0]) + 3, sin(
SECTION(named_mesh.name())
{
auto mesh_2d_v = named_mesh.mesh();
- auto mesh_2d = mesh_2d_v->get<Mesh<Connectivity<Dimension>>>();
+ auto mesh_2d = mesh_2d_v->get<Mesh<Dimension>>();
auto xj = MeshDataManager::instance().getMeshData(*mesh_2d).xj();
@@ -928,7 +928,7 @@ let R2x2_non_linear_2d: R^2 -> R^2x2, x -> [[2*exp(x[0])*sin(x[1])+3, sin(x[0]-2
SECTION(named_mesh.name())
{
auto mesh_3d_v = named_mesh.mesh();
- auto mesh_3d = mesh_3d_v->get<Mesh<Connectivity<Dimension>>>();
+ auto mesh_3d = mesh_3d_v->get<Mesh<Dimension>>();
auto xj = MeshDataManager::instance().getMeshData(*mesh_3d).xj();
diff --git a/tests/test_ItemArray.cpp b/tests/test_ItemArray.cpp
index 313e519727568bf386acbd1826b32e440814839a..212e1214113e929ecf730cfe603f7d1da21963c1 100644
--- a/tests/test_ItemArray.cpp
+++ b/tests/test_ItemArray.cpp
@@ -34,7 +34,7 @@ TEST_CASE("ItemArray", "[mesh]")
SECTION(named_mesh.name())
{
auto mesh_1d_v = named_mesh.mesh();
- auto mesh_1d = mesh_1d_v->get<Mesh<Connectivity<1>>>();
+ auto mesh_1d = mesh_1d_v->get<Mesh<1>>();
const Connectivity<1>& connectivity = mesh_1d->connectivity();
@@ -73,7 +73,7 @@ TEST_CASE("ItemArray", "[mesh]")
SECTION(named_mesh.name())
{
auto mesh_2d_v = named_mesh.mesh();
- auto mesh_2d = mesh_2d_v->get<Mesh<Connectivity<2>>>();
+ auto mesh_2d = mesh_2d_v->get<Mesh<2>>();
const Connectivity<2>& connectivity = mesh_2d->connectivity();
@@ -110,7 +110,7 @@ TEST_CASE("ItemArray", "[mesh]")
SECTION(named_mesh.name())
{
auto mesh_3d_v = named_mesh.mesh();
- auto mesh_3d = mesh_3d_v->get<Mesh<Connectivity<3>>>();
+ auto mesh_3d = mesh_3d_v->get<Mesh<3>>();
const Connectivity<3>& connectivity = mesh_3d->connectivity();
@@ -145,7 +145,7 @@ TEST_CASE("ItemArray", "[mesh]")
SECTION(named_mesh.name())
{
auto mesh_3d_v = named_mesh.mesh();
- auto mesh_3d = mesh_3d_v->get<Mesh<Connectivity<3>>>();
+ auto mesh_3d = mesh_3d_v->get<Mesh<3>>();
const Connectivity<3>& connectivity = mesh_3d->connectivity();
@@ -197,7 +197,7 @@ TEST_CASE("ItemArray", "[mesh]")
SECTION(named_mesh.name())
{
auto mesh_3d_v = named_mesh.mesh();
- auto mesh_3d = mesh_3d_v->get<Mesh<Connectivity<3>>>();
+ auto mesh_3d = mesh_3d_v->get<Mesh<3>>();
const Connectivity<3>& connectivity = mesh_3d->connectivity();
@@ -233,7 +233,7 @@ TEST_CASE("ItemArray", "[mesh]")
SECTION(named_mesh.name())
{
auto mesh_2d_v = named_mesh.mesh();
- auto mesh_2d = mesh_2d_v->get<Mesh<Connectivity<2>>>();
+ auto mesh_2d = mesh_2d_v->get<Mesh<2>>();
const Connectivity<2>& connectivity = mesh_2d->connectivity();
@@ -280,7 +280,7 @@ TEST_CASE("ItemArray", "[mesh]")
SECTION("output")
{
auto mesh_v = MeshDataBaseForTests::get().unordered1DMesh();
- auto mesh = mesh_v->get<Mesh<Connectivity<1>>>();
+ auto mesh = mesh_v->get<Mesh<1>>();
Table<int> table{mesh->numberOfCells(), 3};
for (size_t i = 0; i < table.numberOfRows(); ++i) {
@@ -324,7 +324,7 @@ TEST_CASE("ItemArray", "[mesh]")
SECTION(named_mesh.name())
{
auto mesh_3d_v = named_mesh.mesh();
- auto mesh_3d = mesh_3d_v->get<Mesh<Connectivity<3>>>();
+ auto mesh_3d = mesh_3d_v->get<Mesh<3>>();
const Connectivity<3>& connectivity = mesh_3d->connectivity();
@@ -355,7 +355,7 @@ TEST_CASE("ItemArray", "[mesh]")
SECTION(named_mesh.name())
{
auto mesh_3d_v = named_mesh.mesh();
- auto mesh_3d = mesh_3d_v->get<Mesh<Connectivity<3>>>();
+ auto mesh_3d = mesh_3d_v->get<Mesh<3>>();
const Connectivity<3>& connectivity = mesh_3d->connectivity();
@@ -373,7 +373,7 @@ TEST_CASE("ItemArray", "[mesh]")
SECTION(named_mesh.name())
{
auto mesh_2d_v = named_mesh.mesh();
- auto mesh_2d = mesh_2d_v->get<Mesh<Connectivity<2>>>();
+ auto mesh_2d = mesh_2d_v->get<Mesh<2>>();
const Connectivity<2>& connectivity_2d = mesh_2d->connectivity();
@@ -383,7 +383,7 @@ TEST_CASE("ItemArray", "[mesh]")
SECTION(named_mesh_3d.name())
{
auto mesh_3d_v = named_mesh_3d.mesh();
- auto mesh_3d = mesh_3d_v->get<Mesh<Connectivity<3>>>();
+ auto mesh_3d = mesh_3d_v->get<Mesh<3>>();
const Connectivity<3>& connectivity_3d = mesh_3d->connectivity();
diff --git a/tests/test_ItemArrayUtils.cpp b/tests/test_ItemArrayUtils.cpp
index cd0ea471f1ddef14f2038edf4b8df8e93a740888..8d04842849ed469e1cf0c1e6e217d8e1da4f01f8 100644
--- a/tests/test_ItemArrayUtils.cpp
+++ b/tests/test_ItemArrayUtils.cpp
@@ -26,7 +26,7 @@ TEST_CASE("ItemArrayUtils", "[mesh]")
SECTION(named_mesh.name())
{
auto mesh_1d_v = named_mesh.mesh();
- auto mesh_1d = mesh_1d_v->get<Mesh<Connectivity<1>>>();
+ auto mesh_1d = mesh_1d_v->get<Mesh<1>>();
const Connectivity<1>& connectivity = mesh_1d->connectivity();
@@ -301,7 +301,7 @@ TEST_CASE("ItemArrayUtils", "[mesh]")
SECTION(named_mesh.name())
{
auto mesh_2d_v = named_mesh.mesh();
- auto mesh_2d = mesh_2d_v->get<Mesh<Connectivity<2>>>();
+ auto mesh_2d = mesh_2d_v->get<Mesh<2>>();
const Connectivity<2>& connectivity = mesh_2d->connectivity();
@@ -576,7 +576,7 @@ TEST_CASE("ItemArrayUtils", "[mesh]")
SECTION(named_mesh.name())
{
auto mesh_3d_v = named_mesh.mesh();
- auto mesh_3d = mesh_3d_v->get<Mesh<Connectivity<3>>>();
+ auto mesh_3d = mesh_3d_v->get<Mesh<3>>();
const Connectivity<3>& connectivity = mesh_3d->connectivity();
@@ -854,7 +854,7 @@ TEST_CASE("ItemArrayUtils", "[mesh]")
SECTION(named_mesh.name())
{
auto mesh_1d_v = named_mesh.mesh();
- auto mesh_1d = mesh_1d_v->get<Mesh<Connectivity<1>>>();
+ auto mesh_1d = mesh_1d_v->get<Mesh<1>>();
const Connectivity<1>& connectivity = mesh_1d->connectivity();
@@ -884,7 +884,7 @@ TEST_CASE("ItemArrayUtils", "[mesh]")
SECTION(named_mesh.name())
{
auto mesh_2d_v = named_mesh.mesh();
- auto mesh_2d = mesh_2d_v->get<Mesh<Connectivity<2>>>();
+ auto mesh_2d = mesh_2d_v->get<Mesh<2>>();
const Connectivity<2>& connectivity = mesh_2d->connectivity();
@@ -914,7 +914,7 @@ TEST_CASE("ItemArrayUtils", "[mesh]")
SECTION(named_mesh.name())
{
auto mesh_3d_v = named_mesh.mesh();
- auto mesh_3d = mesh_3d_v->get<Mesh<Connectivity<3>>>();
+ auto mesh_3d = mesh_3d_v->get<Mesh<3>>();
const Connectivity<3>& connectivity = mesh_3d->connectivity();
@@ -947,7 +947,7 @@ TEST_CASE("ItemArrayUtils", "[mesh]")
SECTION(named_mesh.name())
{
auto mesh_1d_v = named_mesh.mesh();
- auto mesh_1d = mesh_1d_v->get<Mesh<Connectivity<1>>>();
+ auto mesh_1d = mesh_1d_v->get<Mesh<1>>();
const Connectivity<1>& connectivity = mesh_1d->connectivity();
@@ -977,7 +977,7 @@ TEST_CASE("ItemArrayUtils", "[mesh]")
SECTION(named_mesh.name())
{
auto mesh_2d_v = named_mesh.mesh();
- auto mesh_2d = mesh_2d_v->get<Mesh<Connectivity<2>>>();
+ auto mesh_2d = mesh_2d_v->get<Mesh<2>>();
const Connectivity<2>& connectivity = mesh_2d->connectivity();
@@ -1007,7 +1007,7 @@ TEST_CASE("ItemArrayUtils", "[mesh]")
SECTION(named_mesh.name())
{
auto mesh_3d_v = named_mesh.mesh();
- auto mesh_3d = mesh_3d_v->get<Mesh<Connectivity<3>>>();
+ auto mesh_3d = mesh_3d_v->get<Mesh<3>>();
const Connectivity<3>& connectivity = mesh_3d->connectivity();
@@ -1040,7 +1040,7 @@ TEST_CASE("ItemArrayUtils", "[mesh]")
SECTION(named_mesh.name())
{
auto mesh_1d_v = named_mesh.mesh();
- auto mesh_1d = mesh_1d_v->get<Mesh<Connectivity<1>>>();
+ auto mesh_1d = mesh_1d_v->get<Mesh<1>>();
const Connectivity<1>& connectivity = mesh_1d->connectivity();
@@ -1070,7 +1070,7 @@ TEST_CASE("ItemArrayUtils", "[mesh]")
SECTION(named_mesh.name() + " for double data")
{
auto mesh_2d_v = named_mesh.mesh();
- auto mesh_2d = mesh_2d_v->get<Mesh<Connectivity<2>>>();
+ auto mesh_2d = mesh_2d_v->get<Mesh<2>>();
const Connectivity<2>& connectivity = mesh_2d->connectivity();
@@ -1098,7 +1098,7 @@ TEST_CASE("ItemArrayUtils", "[mesh]")
SECTION(named_mesh.name() + " for N^2 data")
{
auto mesh_2d_v = named_mesh.mesh();
- auto mesh_2d = mesh_2d_v->get<Mesh<Connectivity<2>>>();
+ auto mesh_2d = mesh_2d_v->get<Mesh<2>>();
const Connectivity<2>& connectivity = mesh_2d->connectivity();
@@ -1131,7 +1131,7 @@ TEST_CASE("ItemArrayUtils", "[mesh]")
SECTION(named_mesh.name() + " for size_t data")
{
auto mesh_3d_v = named_mesh.mesh();
- auto mesh_3d = mesh_3d_v->get<Mesh<Connectivity<3>>>();
+ auto mesh_3d = mesh_3d_v->get<Mesh<3>>();
const Connectivity<3>& connectivity = mesh_3d->connectivity();
@@ -1154,7 +1154,7 @@ TEST_CASE("ItemArrayUtils", "[mesh]")
SECTION(named_mesh.name() + " for N^3x2 data")
{
auto mesh_3d_v = named_mesh.mesh();
- auto mesh_3d = mesh_3d_v->get<Mesh<Connectivity<3>>>();
+ auto mesh_3d = mesh_3d_v->get<Mesh<3>>();
const Connectivity<3>& connectivity = mesh_3d->connectivity();
diff --git a/tests/test_ItemArrayVariant.cpp b/tests/test_ItemArrayVariant.cpp
index 878ee2a7fcf49a1e07af8501edf28c42d58bbbe5..1f0d1d822c594d223cb94075b6be1a87fddb730a 100644
--- a/tests/test_ItemArrayVariant.cpp
+++ b/tests/test_ItemArrayVariant.cpp
@@ -12,7 +12,7 @@
TEST_CASE("ItemArrayVariant", "[mesh]")
{
std::shared_ptr mesh_v = MeshDataBaseForTests::get().hybrid2DMesh();
- std::shared_ptr mesh = mesh_v->get<Mesh<Connectivity<2>>>();
+ std::shared_ptr mesh = mesh_v->get<Mesh<2>>();
const Connectivity<2>& connectivity = *mesh->shared_connectivity();
diff --git a/tests/test_ItemArrayVariantFunctionInterpoler.cpp b/tests/test_ItemArrayVariantFunctionInterpoler.cpp
index 50af9207ea170fac28e74ede2a293bcdcc1d67a8..95f9d9b5c45c0ff66d0c908e5831877f99ef2e29 100644
--- a/tests/test_ItemArrayVariantFunctionInterpoler.cpp
+++ b/tests/test_ItemArrayVariantFunctionInterpoler.cpp
@@ -56,7 +56,7 @@ TEST_CASE("ItemArrayVariantFunctionInterpoler", "[scheme]")
SECTION(named_mesh.name())
{
auto mesh_1d_v = named_mesh.mesh();
- auto mesh_1d = mesh_1d_v->get<Mesh<Connectivity<1>>>();
+ auto mesh_1d = mesh_1d_v->get<Mesh<1>>();
auto xj = MeshDataManager::instance().getMeshData(*mesh_1d).xj();
auto xr = mesh_1d->xr();
@@ -431,7 +431,7 @@ let R3x3_non_linear_1d: R^1 -> R^3x3, x -> [[2 * exp(x[0]) * sin(x[0]) + 3, sin(
SECTION(named_mesh.name())
{
auto mesh_2d_v = named_mesh.mesh();
- auto mesh_2d = mesh_2d_v->get<Mesh<Connectivity<2>>>();
+ auto mesh_2d = mesh_2d_v->get<Mesh<2>>();
auto xl = MeshDataManager::instance().getMeshData(*mesh_2d).xl();
@@ -565,7 +565,7 @@ let R3x3_non_linear_2d: R^2 -> R^3x3, x -> [[2 * exp(x[0]) * sin(x[1]) + 3, sin(
SECTION(named_mesh.name())
{
auto mesh_3d_v = named_mesh.mesh();
- auto mesh_3d = mesh_3d_v->get<Mesh<Connectivity<3>>>();
+ auto mesh_3d = mesh_3d_v->get<Mesh<3>>();
auto xe = MeshDataManager::instance().getMeshData(*mesh_3d).xe();
diff --git a/tests/test_ItemValue.cpp b/tests/test_ItemValue.cpp
index 79b54c22b09dc964acdc4b934c25c4527ad5effe..43a8dd459b59a438b4a3b072ff60b9b543dcf8f2 100644
--- a/tests/test_ItemValue.cpp
+++ b/tests/test_ItemValue.cpp
@@ -32,7 +32,7 @@ TEST_CASE("ItemValue", "[mesh]")
SECTION(named_mesh.name())
{
auto mesh_1d_v = named_mesh.mesh();
- auto mesh_1d = mesh_1d_v->get<Mesh<Connectivity<1>>>();
+ auto mesh_1d = mesh_1d_v->get<Mesh<1>>();
const Connectivity<1>& connectivity = mesh_1d->connectivity();
@@ -66,7 +66,7 @@ TEST_CASE("ItemValue", "[mesh]")
SECTION(named_mesh.name())
{
auto mesh_2d_v = named_mesh.mesh();
- auto mesh_2d = mesh_2d_v->get<Mesh<Connectivity<2>>>();
+ auto mesh_2d = mesh_2d_v->get<Mesh<2>>();
const Connectivity<2>& connectivity = mesh_2d->connectivity();
@@ -96,7 +96,7 @@ TEST_CASE("ItemValue", "[mesh]")
SECTION(named_mesh.name())
{
auto mesh_3d_v = named_mesh.mesh();
- auto mesh_3d = mesh_3d_v->get<Mesh<Connectivity<3>>>();
+ auto mesh_3d = mesh_3d_v->get<Mesh<3>>();
const Connectivity<3>& connectivity = mesh_3d->connectivity();
@@ -121,7 +121,7 @@ TEST_CASE("ItemValue", "[mesh]")
SECTION(named_mesh.name())
{
auto mesh_3d_v = named_mesh.mesh();
- auto mesh_3d = mesh_3d_v->get<Mesh<Connectivity<3>>>();
+ auto mesh_3d = mesh_3d_v->get<Mesh<3>>();
const Connectivity<3>& connectivity = mesh_3d->connectivity();
@@ -152,7 +152,7 @@ TEST_CASE("ItemValue", "[mesh]")
SECTION(named_mesh.name())
{
auto mesh_3d_v = named_mesh.mesh();
- auto mesh_3d = mesh_3d_v->get<Mesh<Connectivity<3>>>();
+ auto mesh_3d = mesh_3d_v->get<Mesh<3>>();
const Connectivity<3>& connectivity = mesh_3d->connectivity();
@@ -191,7 +191,7 @@ TEST_CASE("ItemValue", "[mesh]")
SECTION(named_mesh.name())
{
auto mesh_2d_v = named_mesh.mesh();
- auto mesh_2d = mesh_2d_v->get<Mesh<Connectivity<2>>>();
+ auto mesh_2d = mesh_2d_v->get<Mesh<2>>();
const Connectivity<2>& connectivity = mesh_2d->connectivity();
@@ -226,7 +226,7 @@ TEST_CASE("ItemValue", "[mesh]")
SECTION("output")
{
auto mesh_v = MeshDataBaseForTests::get().unordered1DMesh();
- auto mesh = mesh_v->get<Mesh<Connectivity<1>>>();
+ auto mesh = mesh_v->get<Mesh<1>>();
Array<int> array{mesh->numberOfCells()};
for (size_t i = 0; i < array.size(); ++i) {
@@ -268,7 +268,7 @@ TEST_CASE("ItemValue", "[mesh]")
SECTION(named_mesh.name())
{
auto mesh_3d_v = named_mesh.mesh();
- auto mesh_3d = mesh_3d_v->get<Mesh<Connectivity<3>>>();
+ auto mesh_3d = mesh_3d_v->get<Mesh<3>>();
const Connectivity<3>& connectivity = mesh_3d->connectivity();
@@ -299,7 +299,7 @@ TEST_CASE("ItemValue", "[mesh]")
SECTION(named_mesh.name())
{
auto mesh_3d_v = named_mesh.mesh();
- auto mesh_3d = mesh_3d_v->get<Mesh<Connectivity<3>>>();
+ auto mesh_3d = mesh_3d_v->get<Mesh<3>>();
const Connectivity<3>& connectivity = mesh_3d->connectivity();
@@ -317,7 +317,7 @@ TEST_CASE("ItemValue", "[mesh]")
SECTION(named_mesh.name())
{
auto mesh_2d_v = named_mesh.mesh();
- auto mesh_2d = mesh_2d_v->get<Mesh<Connectivity<2>>>();
+ auto mesh_2d = mesh_2d_v->get<Mesh<2>>();
const Connectivity<2>& connectivity_2d = mesh_2d->connectivity();
@@ -327,7 +327,7 @@ TEST_CASE("ItemValue", "[mesh]")
SECTION(named_mesh_3d.name())
{
auto mesh_3d_v = named_mesh_3d.mesh();
- auto mesh_3d = mesh_3d_v->get<Mesh<Connectivity<3>>>();
+ auto mesh_3d = mesh_3d_v->get<Mesh<3>>();
const Connectivity<3>& connectivity_3d = mesh_3d->connectivity();
diff --git a/tests/test_ItemValueUtils.cpp b/tests/test_ItemValueUtils.cpp
index 6c5d8b0af5a934ba79c76bdc3f695b207b704638..f80819f2c619709584c6574dd185d5a0e2f512fb 100644
--- a/tests/test_ItemValueUtils.cpp
+++ b/tests/test_ItemValueUtils.cpp
@@ -24,7 +24,7 @@ TEST_CASE("ItemValueUtils", "[mesh]")
SECTION(named_mesh.name())
{
auto mesh_2d_v = named_mesh.mesh();
- auto mesh_2d = mesh_2d_v->get<Mesh<Connectivity<2>>>();
+ auto mesh_2d = mesh_2d_v->get<Mesh<2>>();
const Connectivity<2>& connectivity = mesh_2d->connectivity();
@@ -78,7 +78,7 @@ TEST_CASE("ItemValueUtils", "[mesh]")
SECTION(named_mesh.name())
{
auto mesh_1d_v = named_mesh.mesh();
- auto mesh_1d = mesh_1d_v->get<Mesh<Connectivity<1>>>();
+ auto mesh_1d = mesh_1d_v->get<Mesh<1>>();
const Connectivity<1>& connectivity = mesh_1d->connectivity();
@@ -106,7 +106,7 @@ TEST_CASE("ItemValueUtils", "[mesh]")
SECTION(named_mesh.name())
{
auto mesh_2d_v = named_mesh.mesh();
- auto mesh_2d = mesh_2d_v->get<Mesh<Connectivity<2>>>();
+ auto mesh_2d = mesh_2d_v->get<Mesh<2>>();
const Connectivity<2>& connectivity = mesh_2d->connectivity();
@@ -134,7 +134,7 @@ TEST_CASE("ItemValueUtils", "[mesh]")
SECTION(named_mesh.name())
{
auto mesh_3d_v = named_mesh.mesh();
- auto mesh_3d = mesh_3d_v->get<Mesh<Connectivity<3>>>();
+ auto mesh_3d = mesh_3d_v->get<Mesh<3>>();
const Connectivity<3>& connectivity = mesh_3d->connectivity();
@@ -165,7 +165,7 @@ TEST_CASE("ItemValueUtils", "[mesh]")
SECTION(named_mesh.name())
{
auto mesh_1d_v = named_mesh.mesh();
- auto mesh_1d = mesh_1d_v->get<Mesh<Connectivity<1>>>();
+ auto mesh_1d = mesh_1d_v->get<Mesh<1>>();
const Connectivity<1>& connectivity = mesh_1d->connectivity();
@@ -193,7 +193,7 @@ TEST_CASE("ItemValueUtils", "[mesh]")
SECTION(named_mesh.name())
{
auto mesh_2d_v = named_mesh.mesh();
- auto mesh_2d = mesh_2d_v->get<Mesh<Connectivity<2>>>();
+ auto mesh_2d = mesh_2d_v->get<Mesh<2>>();
const Connectivity<2>& connectivity = mesh_2d->connectivity();
@@ -221,7 +221,7 @@ TEST_CASE("ItemValueUtils", "[mesh]")
SECTION(named_mesh.name())
{
auto mesh_3d_v = named_mesh.mesh();
- auto mesh_3d = mesh_3d_v->get<Mesh<Connectivity<3>>>();
+ auto mesh_3d = mesh_3d_v->get<Mesh<3>>();
const Connectivity<3>& connectivity = mesh_3d->connectivity();
@@ -252,7 +252,7 @@ TEST_CASE("ItemValueUtils", "[mesh]")
SECTION(named_mesh.name() + " for size_t data")
{
auto mesh_1d_v = named_mesh.mesh();
- auto mesh_1d = mesh_1d_v->get<Mesh<Connectivity<1>>>();
+ auto mesh_1d = mesh_1d_v->get<Mesh<1>>();
const Connectivity<1>& connectivity = mesh_1d->connectivity();
@@ -280,7 +280,7 @@ TEST_CASE("ItemValueUtils", "[mesh]")
SECTION(named_mesh.name() + " for double data")
{
auto mesh_1d_v = named_mesh.mesh();
- auto mesh_1d = mesh_1d_v->get<Mesh<Connectivity<1>>>();
+ auto mesh_1d = mesh_1d_v->get<Mesh<1>>();
const Connectivity<1>& connectivity = mesh_1d->connectivity();
@@ -317,7 +317,7 @@ TEST_CASE("ItemValueUtils", "[mesh]")
SECTION(named_mesh.name() + "for size_t data")
{
auto mesh_2d_v = named_mesh.mesh();
- auto mesh_2d = mesh_2d_v->get<Mesh<Connectivity<2>>>();
+ auto mesh_2d = mesh_2d_v->get<Mesh<2>>();
const Connectivity<2>& connectivity = mesh_2d->connectivity();
@@ -345,7 +345,7 @@ TEST_CASE("ItemValueUtils", "[mesh]")
SECTION(named_mesh.name() + "for double data")
{
auto mesh_2d_v = named_mesh.mesh();
- auto mesh_2d = mesh_2d_v->get<Mesh<Connectivity<2>>>();
+ auto mesh_2d = mesh_2d_v->get<Mesh<2>>();
const Connectivity<2>& connectivity = mesh_2d->connectivity();
@@ -375,7 +375,7 @@ TEST_CASE("ItemValueUtils", "[mesh]")
SECTION(named_mesh.name() + "for N^3 data")
{
auto mesh_2d_v = named_mesh.mesh();
- auto mesh_2d = mesh_2d_v->get<Mesh<Connectivity<2>>>();
+ auto mesh_2d = mesh_2d_v->get<Mesh<2>>();
const Connectivity<2>& connectivity = mesh_2d->connectivity();
@@ -407,7 +407,7 @@ TEST_CASE("ItemValueUtils", "[mesh]")
SECTION(named_mesh.name() + "for R2 data")
{
auto mesh_2d_v = named_mesh.mesh();
- auto mesh_2d = mesh_2d_v->get<Mesh<Connectivity<2>>>();
+ auto mesh_2d = mesh_2d_v->get<Mesh<2>>();
const Connectivity<2>& connectivity = mesh_2d->connectivity();
@@ -452,7 +452,7 @@ TEST_CASE("ItemValueUtils", "[mesh]")
SECTION(named_mesh.name() + " for size_t data")
{
auto mesh_3d_v = named_mesh.mesh();
- auto mesh_3d = mesh_3d_v->get<Mesh<Connectivity<3>>>();
+ auto mesh_3d = mesh_3d_v->get<Mesh<3>>();
const Connectivity<3>& connectivity = mesh_3d->connectivity();
@@ -475,7 +475,7 @@ TEST_CASE("ItemValueUtils", "[mesh]")
SECTION(named_mesh.name() + " for N^3x2 data")
{
auto mesh_3d_v = named_mesh.mesh();
- auto mesh_3d = mesh_3d_v->get<Mesh<Connectivity<3>>>();
+ auto mesh_3d = mesh_3d_v->get<Mesh<3>>();
const Connectivity<3>& connectivity = mesh_3d->connectivity();
@@ -502,7 +502,7 @@ TEST_CASE("ItemValueUtils", "[mesh]")
SECTION(named_mesh.name() + "for R^2x3 data")
{
auto mesh_3d_v = named_mesh.mesh();
- auto mesh_3d = mesh_3d_v->get<Mesh<Connectivity<3>>>();
+ auto mesh_3d = mesh_3d_v->get<Mesh<3>>();
const Connectivity<3>& connectivity = mesh_3d->connectivity();
diff --git a/tests/test_ItemValueVariant.cpp b/tests/test_ItemValueVariant.cpp
index fb4180f4a5df1961c9ba70a5821357a88f032be1..6eed537b91a0ef1a9009c79653ce89b141e57d27 100644
--- a/tests/test_ItemValueVariant.cpp
+++ b/tests/test_ItemValueVariant.cpp
@@ -12,7 +12,7 @@
TEST_CASE("ItemValueVariant", "[mesh]")
{
auto mesh_2d_v = MeshDataBaseForTests::get().hybrid2DMesh();
- auto mesh_2d = mesh_2d_v->get<Mesh<Connectivity<2>>>();
+ auto mesh_2d = mesh_2d_v->get<Mesh<2>>();
const Connectivity<2>& connectivity = *mesh_2d->shared_connectivity();
diff --git a/tests/test_ItemValueVariantFunctionInterpoler.cpp b/tests/test_ItemValueVariantFunctionInterpoler.cpp
index a55435208c628ffb55d7f6e3daebbc72e4508616..4bc251aba0b6b75d7c65d87474e72010f0cecf1b 100644
--- a/tests/test_ItemValueVariantFunctionInterpoler.cpp
+++ b/tests/test_ItemValueVariantFunctionInterpoler.cpp
@@ -50,7 +50,7 @@ TEST_CASE("ItemValueVariantFunctionInterpoler", "[scheme]")
SECTION(named_mesh.name())
{
auto mesh_1d_v = named_mesh.mesh();
- auto mesh_1d = mesh_1d_v->get<Mesh<Connectivity<1>>>();
+ auto mesh_1d = mesh_1d_v->get<Mesh<1>>();
auto xj = MeshDataManager::instance().getMeshData(*mesh_1d).xj();
auto xr = mesh_1d->xr();
@@ -331,7 +331,7 @@ let R3x3_non_linear_1d: R^1 -> R^3x3, x -> [[2 * exp(x[0]) * sin(x[0]) + 3, sin(
SECTION(named_mesh.name())
{
auto mesh_2d_v = named_mesh.mesh();
- auto mesh_2d = mesh_2d_v->get<Mesh<Connectivity<2>>>();
+ auto mesh_2d = mesh_2d_v->get<Mesh<2>>();
auto xl = MeshDataManager::instance().getMeshData(*mesh_2d).xl();
@@ -451,7 +451,7 @@ let R3x3_non_linear_2d: R^2 -> R^3x3, x -> [[2 * exp(x[0]) * sin(x[1]) + 3, sin(
SECTION(named_mesh.name())
{
auto mesh_3d_v = named_mesh.mesh();
- auto mesh_3d = mesh_3d_v->get<Mesh<Connectivity<3>>>();
+ auto mesh_3d = mesh_3d_v->get<Mesh<3>>();
auto xe = MeshDataManager::instance().getMeshData(*mesh_3d).xe();
diff --git a/tests/test_MedianDualConnectivityBuilder.cpp b/tests/test_MedianDualConnectivityBuilder.cpp
index 34728918dc7d61613ff2d8306d1f4b157d90c89d..ab418f24a8cd59fab535cf5b4997f3e82f1e8972 100644
--- a/tests/test_MedianDualConnectivityBuilder.cpp
+++ b/tests/test_MedianDualConnectivityBuilder.cpp
@@ -37,8 +37,8 @@ TEST_CASE("MedianDualConnectivityBuilder", "[mesh]")
{
constexpr static size_t Dimension = 2;
- using ConnectivityType = Connectivity<Dimension>;
- using MeshType = Mesh<ConnectivityType>;
+ using MeshType = Mesh<Dimension>;
+ using ConnectivityType = typename MeshType::Connectivity;
std::shared_ptr<const MeshType> mesh = MeshDataBaseForTests::get().hybrid2DMesh()->get<MeshType>();
const ConnectivityType& primal_connectivity = mesh->connectivity();
diff --git a/tests/test_MedianDualMeshBuilder.cpp b/tests/test_MedianDualMeshBuilder.cpp
index ac2fe9bb4fb56e1a4f5252d41a1b0fd75cb4ebe0..f4f96d240dcab0f9dbe09cb32c5a25366402574c 100644
--- a/tests/test_MedianDualMeshBuilder.cpp
+++ b/tests/test_MedianDualMeshBuilder.cpp
@@ -18,8 +18,7 @@ TEST_CASE("MedianDualMeshBuilder", "[mesh]")
{
constexpr static size_t Dimension = 2;
- using ConnectivityType = Connectivity<Dimension>;
- using MeshType = Mesh<ConnectivityType>;
+ using MeshType = Mesh<Dimension>;
std::shared_ptr primal_mesh_v = MeshDataBaseForTests::get().hybrid2DMesh();
const MeshType& primal_mesh = *primal_mesh_v->get<MeshType>();
diff --git a/tests/test_MeshEdgeBoundary.cpp b/tests/test_MeshEdgeBoundary.cpp
index 0659a644b7907d5debbfedf3382087728654395e..307d73fdcccdc990a1b8cfa740493dfb21d154d3 100644
--- a/tests/test_MeshEdgeBoundary.cpp
+++ b/tests/test_MeshEdgeBoundary.cpp
@@ -47,8 +47,8 @@ TEST_CASE("MeshEdgeBoundary", "[mesh]")
{
static constexpr size_t Dimension = 1;
- using ConnectivityType = Connectivity<Dimension>;
- using MeshType = Mesh<ConnectivityType>;
+ using MeshType = Mesh<Dimension>;
+ using ConnectivityType = typename MeshType::Connectivity;
SECTION("cartesian 1d")
{
@@ -120,8 +120,8 @@ TEST_CASE("MeshEdgeBoundary", "[mesh]")
{
static constexpr size_t Dimension = 2;
- using ConnectivityType = Connectivity<Dimension>;
- using MeshType = Mesh<ConnectivityType>;
+ using MeshType = Mesh<Dimension>;
+ using ConnectivityType = typename MeshType::Connectivity;
SECTION("cartesian 2d")
{
@@ -192,8 +192,8 @@ TEST_CASE("MeshEdgeBoundary", "[mesh]")
{
static constexpr size_t Dimension = 3;
- using ConnectivityType = Connectivity<Dimension>;
- using MeshType = Mesh<ConnectivityType>;
+ using MeshType = Mesh<Dimension>;
+ using ConnectivityType = typename MeshType::Connectivity;
SECTION("cartesian 3d")
{
@@ -266,8 +266,7 @@ TEST_CASE("MeshEdgeBoundary", "[mesh]")
{
static constexpr size_t Dimension = 3;
- using ConnectivityType = Connectivity<Dimension>;
- using MeshType = Mesh<ConnectivityType>;
+ using MeshType = Mesh<Dimension>;
std::shared_ptr p_mesh = MeshDataBaseForTests::get().hybrid3DMesh();
const MeshType& mesh = *p_mesh->get<MeshType>();
@@ -284,8 +283,7 @@ TEST_CASE("MeshEdgeBoundary", "[mesh]")
{
static constexpr size_t Dimension = 1;
- using ConnectivityType = Connectivity<Dimension>;
- using MeshType = Mesh<ConnectivityType>;
+ using MeshType = Mesh<Dimension>;
std::shared_ptr p_mesh = MeshDataBaseForTests::get().unordered1DMesh();
const MeshType& mesh = *p_mesh->get<MeshType>();
@@ -301,8 +299,7 @@ TEST_CASE("MeshEdgeBoundary", "[mesh]")
{
static constexpr size_t Dimension = 2;
- using ConnectivityType = Connectivity<Dimension>;
- using MeshType = Mesh<ConnectivityType>;
+ using MeshType = Mesh<Dimension>;
std::shared_ptr p_mesh = MeshDataBaseForTests::get().hybrid2DMesh();
const MeshType& mesh = *p_mesh->get<MeshType>();
@@ -318,8 +315,7 @@ TEST_CASE("MeshEdgeBoundary", "[mesh]")
{
static constexpr size_t Dimension = 3;
- using ConnectivityType = Connectivity<Dimension>;
- using MeshType = Mesh<ConnectivityType>;
+ using MeshType = Mesh<Dimension>;
std::shared_ptr p_mesh = MeshDataBaseForTests::get().hybrid3DMesh();
const MeshType& mesh = *p_mesh->get<MeshType>();
diff --git a/tests/test_MeshEdgeInterface.cpp b/tests/test_MeshEdgeInterface.cpp
index b42324c5af4a519da1b7cd17a02c96c09181fca6..328a1ce637f34969994766961f69be5e186ae6a9 100644
--- a/tests/test_MeshEdgeInterface.cpp
+++ b/tests/test_MeshEdgeInterface.cpp
@@ -47,8 +47,8 @@ TEST_CASE("MeshEdgeInterface", "[mesh]")
{
static constexpr size_t Dimension = 1;
- using ConnectivityType = Connectivity<Dimension>;
- using MeshType = Mesh<ConnectivityType>;
+ using MeshType = Mesh<Dimension>;
+ using ConnectivityType = typename MeshType::Connectivity;
SECTION("unordered 1d")
{
@@ -87,8 +87,8 @@ TEST_CASE("MeshEdgeInterface", "[mesh]")
{
static constexpr size_t Dimension = 2;
- using ConnectivityType = Connectivity<Dimension>;
- using MeshType = Mesh<ConnectivityType>;
+ using MeshType = Mesh<Dimension>;
+ using ConnectivityType = typename MeshType::Connectivity;
SECTION("hybrid 2d")
{
@@ -127,8 +127,8 @@ TEST_CASE("MeshEdgeInterface", "[mesh]")
{
static constexpr size_t Dimension = 3;
- using ConnectivityType = Connectivity<Dimension>;
- using MeshType = Mesh<ConnectivityType>;
+ using MeshType = Mesh<Dimension>;
+ using ConnectivityType = typename MeshType::Connectivity;
SECTION("hybrid 3d")
{
@@ -169,8 +169,7 @@ TEST_CASE("MeshEdgeInterface", "[mesh]")
{
static constexpr size_t Dimension = 3;
- using ConnectivityType = Connectivity<Dimension>;
- using MeshType = Mesh<ConnectivityType>;
+ using MeshType = Mesh<Dimension>;
std::shared_ptr p_mesh = MeshDataBaseForTests::get().hybrid3DMesh();
const MeshType& mesh = *p_mesh->get<MeshType>();
diff --git a/tests/test_MeshFaceBoundary.cpp b/tests/test_MeshFaceBoundary.cpp
index 9345506d744921e3330635b8e5e4f3d39b3eab72..40de5b4e1594c454529730cf2fdfb322bf570814 100644
--- a/tests/test_MeshFaceBoundary.cpp
+++ b/tests/test_MeshFaceBoundary.cpp
@@ -47,8 +47,8 @@ TEST_CASE("MeshFaceBoundary", "[mesh]")
{
static constexpr size_t Dimension = 1;
- using ConnectivityType = Connectivity<Dimension>;
- using MeshType = Mesh<ConnectivityType>;
+ using MeshType = Mesh<Dimension>;
+ using ConnectivityType = typename MeshType::Connectivity;
SECTION("cartesian 1d")
{
@@ -120,8 +120,8 @@ TEST_CASE("MeshFaceBoundary", "[mesh]")
{
static constexpr size_t Dimension = 2;
- using ConnectivityType = Connectivity<Dimension>;
- using MeshType = Mesh<ConnectivityType>;
+ using MeshType = Mesh<Dimension>;
+ using ConnectivityType = typename MeshType::Connectivity;
SECTION("cartesian 2d")
{
@@ -192,8 +192,8 @@ TEST_CASE("MeshFaceBoundary", "[mesh]")
{
static constexpr size_t Dimension = 3;
- using ConnectivityType = Connectivity<Dimension>;
- using MeshType = Mesh<ConnectivityType>;
+ using MeshType = Mesh<Dimension>;
+ using ConnectivityType = typename MeshType::Connectivity;
SECTION("cartesian 3d")
{
@@ -266,8 +266,7 @@ TEST_CASE("MeshFaceBoundary", "[mesh]")
{
static constexpr size_t Dimension = 3;
- using ConnectivityType = Connectivity<Dimension>;
- using MeshType = Mesh<ConnectivityType>;
+ using MeshType = Mesh<Dimension>;
std::shared_ptr p_mesh = MeshDataBaseForTests::get().hybrid3DMesh();
const MeshType& mesh = *p_mesh->get<MeshType>();
@@ -284,8 +283,7 @@ TEST_CASE("MeshFaceBoundary", "[mesh]")
{
static constexpr size_t Dimension = 1;
- using ConnectivityType = Connectivity<Dimension>;
- using MeshType = Mesh<ConnectivityType>;
+ using MeshType = Mesh<Dimension>;
std::shared_ptr p_mesh = MeshDataBaseForTests::get().unordered1DMesh();
const MeshType& mesh = *p_mesh->get<MeshType>();
@@ -301,8 +299,7 @@ TEST_CASE("MeshFaceBoundary", "[mesh]")
{
static constexpr size_t Dimension = 2;
- using ConnectivityType = Connectivity<Dimension>;
- using MeshType = Mesh<ConnectivityType>;
+ using MeshType = Mesh<Dimension>;
std::shared_ptr p_mesh = MeshDataBaseForTests::get().hybrid2DMesh();
const MeshType& mesh = *p_mesh->get<MeshType>();
@@ -318,8 +315,7 @@ TEST_CASE("MeshFaceBoundary", "[mesh]")
{
static constexpr size_t Dimension = 3;
- using ConnectivityType = Connectivity<Dimension>;
- using MeshType = Mesh<ConnectivityType>;
+ using MeshType = Mesh<Dimension>;
std::shared_ptr p_mesh = MeshDataBaseForTests::get().hybrid3DMesh();
const MeshType& mesh = *p_mesh->get<MeshType>();
diff --git a/tests/test_MeshFaceInterface.cpp b/tests/test_MeshFaceInterface.cpp
index d762b712c2d5d5c169ba9472972963f630d9468f..b213ffd48d67677700943549b436a4fefbdb9471 100644
--- a/tests/test_MeshFaceInterface.cpp
+++ b/tests/test_MeshFaceInterface.cpp
@@ -47,8 +47,8 @@ TEST_CASE("MeshFaceInterface", "[mesh]")
{
static constexpr size_t Dimension = 1;
- using ConnectivityType = Connectivity<Dimension>;
- using MeshType = Mesh<ConnectivityType>;
+ using MeshType = Mesh<Dimension>;
+ using ConnectivityType = typename MeshType::Connectivity;
SECTION("unordered 1d")
{
@@ -87,8 +87,8 @@ TEST_CASE("MeshFaceInterface", "[mesh]")
{
static constexpr size_t Dimension = 2;
- using ConnectivityType = Connectivity<Dimension>;
- using MeshType = Mesh<ConnectivityType>;
+ using MeshType = Mesh<Dimension>;
+ using ConnectivityType = typename MeshType::Connectivity;
SECTION("hybrid 2d")
{
@@ -127,8 +127,8 @@ TEST_CASE("MeshFaceInterface", "[mesh]")
{
static constexpr size_t Dimension = 3;
- using ConnectivityType = Connectivity<Dimension>;
- using MeshType = Mesh<ConnectivityType>;
+ using MeshType = Mesh<Dimension>;
+ using ConnectivityType = typename MeshType::Connectivity;
SECTION("hybrid 3d")
{
@@ -169,8 +169,7 @@ TEST_CASE("MeshFaceInterface", "[mesh]")
{
static constexpr size_t Dimension = 3;
- using ConnectivityType = Connectivity<Dimension>;
- using MeshType = Mesh<ConnectivityType>;
+ using MeshType = Mesh<Dimension>;
std::shared_ptr p_mesh = MeshDataBaseForTests::get().hybrid3DMesh();
const MeshType& mesh = *p_mesh->get<MeshType>();
@@ -187,8 +186,7 @@ TEST_CASE("MeshFaceInterface", "[mesh]")
{
static constexpr size_t Dimension = 1;
- using ConnectivityType = Connectivity<Dimension>;
- using MeshType = Mesh<ConnectivityType>;
+ using MeshType = Mesh<Dimension>;
std::shared_ptr p_mesh = MeshDataBaseForTests::get().unordered1DMesh();
const MeshType& mesh = *p_mesh->get<MeshType>();
@@ -204,8 +202,7 @@ TEST_CASE("MeshFaceInterface", "[mesh]")
{
static constexpr size_t Dimension = 2;
- using ConnectivityType = Connectivity<Dimension>;
- using MeshType = Mesh<ConnectivityType>;
+ using MeshType = Mesh<Dimension>;
std::shared_ptr p_mesh = MeshDataBaseForTests::get().hybrid2DMesh();
const MeshType& mesh = *p_mesh->get<MeshType>();
@@ -221,8 +218,7 @@ TEST_CASE("MeshFaceInterface", "[mesh]")
{
static constexpr size_t Dimension = 3;
- using ConnectivityType = Connectivity<Dimension>;
- using MeshType = Mesh<ConnectivityType>;
+ using MeshType = Mesh<Dimension>;
std::shared_ptr p_mesh = MeshDataBaseForTests::get().hybrid3DMesh();
const MeshType& mesh = *p_mesh->get<MeshType>();
diff --git a/tests/test_MeshFlatEdgeBoundary.cpp b/tests/test_MeshFlatEdgeBoundary.cpp
index 18070e5db96b51fee1a7f96edd8cf336e9f1b05c..f717565461f9ae77ee8424108013293ec20060f9 100644
--- a/tests/test_MeshFlatEdgeBoundary.cpp
+++ b/tests/test_MeshFlatEdgeBoundary.cpp
@@ -50,8 +50,8 @@ TEST_CASE("MeshFlatEdgeBoundary", "[mesh]")
{
static constexpr size_t Dimension = 1;
- using ConnectivityType = Connectivity<Dimension>;
- using MeshType = Mesh<ConnectivityType>;
+ using MeshType = Mesh<Dimension>;
+ using ConnectivityType = typename MeshType::Connectivity;
using R1 = TinyVector<1>;
@@ -186,8 +186,8 @@ TEST_CASE("MeshFlatEdgeBoundary", "[mesh]")
{
static constexpr size_t Dimension = 2;
- using ConnectivityType = Connectivity<Dimension>;
- using MeshType = Mesh<ConnectivityType>;
+ using MeshType = Mesh<Dimension>;
+ using ConnectivityType = typename MeshType::Connectivity;
using R2 = TinyVector<2>;
@@ -389,8 +389,8 @@ TEST_CASE("MeshFlatEdgeBoundary", "[mesh]")
{
static constexpr size_t Dimension = 3;
- using ConnectivityType = Connectivity<Dimension>;
- using MeshType = Mesh<ConnectivityType>;
+ using MeshType = Mesh<Dimension>;
+ using ConnectivityType = typename MeshType::Connectivity;
using R3 = TinyVector<3>;
@@ -752,8 +752,8 @@ TEST_CASE("MeshFlatEdgeBoundary", "[mesh]")
{
static constexpr size_t Dimension = 2;
- using ConnectivityType = Connectivity<Dimension>;
- using MeshType = Mesh<ConnectivityType>;
+ using MeshType = Mesh<Dimension>;
+ using ConnectivityType = typename MeshType::Connectivity;
using R2 = TinyVector<2>;
@@ -932,8 +932,8 @@ TEST_CASE("MeshFlatEdgeBoundary", "[mesh]")
{
static constexpr size_t Dimension = 3;
- using ConnectivityType = Connectivity<Dimension>;
- using MeshType = Mesh<ConnectivityType>;
+ using MeshType = Mesh<Dimension>;
+ using ConnectivityType = typename MeshType::Connectivity;
using R3 = TinyVector<3>;
@@ -1143,8 +1143,8 @@ TEST_CASE("MeshFlatEdgeBoundary", "[mesh]")
{
static constexpr size_t Dimension = 2;
- using ConnectivityType = Connectivity<Dimension>;
- using MeshType = Mesh<ConnectivityType>;
+ using MeshType = Mesh<Dimension>;
+ using ConnectivityType = typename MeshType::Connectivity;
using R2 = TinyVector<2>;
@@ -1244,8 +1244,8 @@ TEST_CASE("MeshFlatEdgeBoundary", "[mesh]")
{
static constexpr size_t Dimension = 3;
- using ConnectivityType = Connectivity<Dimension>;
- using MeshType = Mesh<ConnectivityType>;
+ using MeshType = Mesh<Dimension>;
+ using ConnectivityType = typename MeshType::Connectivity;
using R3 = TinyVector<3>;
diff --git a/tests/test_MeshFlatFaceBoundary.cpp b/tests/test_MeshFlatFaceBoundary.cpp
index f7a3ea2fabfc3a8e80b649e2527fb31f5784b119..f5630ee1d2b60858525c7b17a60312fcfe55f524 100644
--- a/tests/test_MeshFlatFaceBoundary.cpp
+++ b/tests/test_MeshFlatFaceBoundary.cpp
@@ -50,8 +50,8 @@ TEST_CASE("MeshFlatFaceBoundary", "[mesh]")
{
static constexpr size_t Dimension = 1;
- using ConnectivityType = Connectivity<Dimension>;
- using MeshType = Mesh<ConnectivityType>;
+ using MeshType = Mesh<Dimension>;
+ using ConnectivityType = typename MeshType::Connectivity;
using R1 = TinyVector<1>;
@@ -186,8 +186,8 @@ TEST_CASE("MeshFlatFaceBoundary", "[mesh]")
{
static constexpr size_t Dimension = 2;
- using ConnectivityType = Connectivity<Dimension>;
- using MeshType = Mesh<ConnectivityType>;
+ using MeshType = Mesh<Dimension>;
+ using ConnectivityType = typename MeshType::Connectivity;
using R2 = TinyVector<2>;
@@ -345,8 +345,8 @@ TEST_CASE("MeshFlatFaceBoundary", "[mesh]")
{
static constexpr size_t Dimension = 3;
- using ConnectivityType = Connectivity<Dimension>;
- using MeshType = Mesh<ConnectivityType>;
+ using MeshType = Mesh<Dimension>;
+ using ConnectivityType = typename MeshType::Connectivity;
using R3 = TinyVector<3>;
@@ -532,8 +532,8 @@ TEST_CASE("MeshFlatFaceBoundary", "[mesh]")
{
static constexpr size_t Dimension = 2;
- using ConnectivityType = Connectivity<Dimension>;
- using MeshType = Mesh<ConnectivityType>;
+ using MeshType = Mesh<Dimension>;
+ using ConnectivityType = typename MeshType::Connectivity;
using R2 = TinyVector<2>;
@@ -712,8 +712,8 @@ TEST_CASE("MeshFlatFaceBoundary", "[mesh]")
{
static constexpr size_t Dimension = 3;
- using ConnectivityType = Connectivity<Dimension>;
- using MeshType = Mesh<ConnectivityType>;
+ using MeshType = Mesh<Dimension>;
+ using ConnectivityType = typename MeshType::Connectivity;
using R3 = TinyVector<3>;
@@ -923,8 +923,8 @@ TEST_CASE("MeshFlatFaceBoundary", "[mesh]")
{
static constexpr size_t Dimension = 2;
- using ConnectivityType = Connectivity<Dimension>;
- using MeshType = Mesh<ConnectivityType>;
+ using MeshType = Mesh<Dimension>;
+ using ConnectivityType = typename MeshType::Connectivity;
using R2 = TinyVector<2>;
@@ -1024,8 +1024,8 @@ TEST_CASE("MeshFlatFaceBoundary", "[mesh]")
{
static constexpr size_t Dimension = 3;
- using ConnectivityType = Connectivity<Dimension>;
- using MeshType = Mesh<ConnectivityType>;
+ using MeshType = Mesh<Dimension>;
+ using ConnectivityType = typename MeshType::Connectivity;
using R3 = TinyVector<3>;
diff --git a/tests/test_MeshFlatNodeBoundary.cpp b/tests/test_MeshFlatNodeBoundary.cpp
index bb5cbff8a4903a36e84fdc2996d078cb4643b166..60348c53b204289a0e0a3745d42dc70af132a3a5 100644
--- a/tests/test_MeshFlatNodeBoundary.cpp
+++ b/tests/test_MeshFlatNodeBoundary.cpp
@@ -50,8 +50,8 @@ TEST_CASE("MeshFlatNodeBoundary", "[mesh]")
{
static constexpr size_t Dimension = 1;
- using ConnectivityType = Connectivity<Dimension>;
- using MeshType = Mesh<ConnectivityType>;
+ using MeshType = Mesh<Dimension>;
+ using ConnectivityType = typename MeshType::Connectivity;
using R1 = TinyVector<1>;
@@ -186,8 +186,8 @@ TEST_CASE("MeshFlatNodeBoundary", "[mesh]")
{
static constexpr size_t Dimension = 2;
- using ConnectivityType = Connectivity<Dimension>;
- using MeshType = Mesh<ConnectivityType>;
+ using MeshType = Mesh<Dimension>;
+ using ConnectivityType = typename MeshType::Connectivity;
using R2 = TinyVector<2>;
@@ -389,8 +389,8 @@ TEST_CASE("MeshFlatNodeBoundary", "[mesh]")
{
static constexpr size_t Dimension = 3;
- using ConnectivityType = Connectivity<Dimension>;
- using MeshType = Mesh<ConnectivityType>;
+ using MeshType = Mesh<Dimension>;
+ using ConnectivityType = typename MeshType::Connectivity;
using R3 = TinyVector<3>;
@@ -752,8 +752,8 @@ TEST_CASE("MeshFlatNodeBoundary", "[mesh]")
{
static constexpr size_t Dimension = 2;
- using ConnectivityType = Connectivity<Dimension>;
- using MeshType = Mesh<ConnectivityType>;
+ using MeshType = Mesh<Dimension>;
+ using ConnectivityType = typename MeshType::Connectivity;
using R2 = TinyVector<2>;
@@ -932,8 +932,8 @@ TEST_CASE("MeshFlatNodeBoundary", "[mesh]")
{
static constexpr size_t Dimension = 3;
- using ConnectivityType = Connectivity<Dimension>;
- using MeshType = Mesh<ConnectivityType>;
+ using MeshType = Mesh<Dimension>;
+ using ConnectivityType = typename MeshType::Connectivity;
using R3 = TinyVector<3>;
@@ -1143,8 +1143,8 @@ TEST_CASE("MeshFlatNodeBoundary", "[mesh]")
{
static constexpr size_t Dimension = 2;
- using ConnectivityType = Connectivity<Dimension>;
- using MeshType = Mesh<ConnectivityType>;
+ using MeshType = Mesh<Dimension>;
+ using ConnectivityType = typename MeshType::Connectivity;
using R2 = TinyVector<2>;
@@ -1237,8 +1237,8 @@ TEST_CASE("MeshFlatNodeBoundary", "[mesh]")
{
static constexpr size_t Dimension = 3;
- using ConnectivityType = Connectivity<Dimension>;
- using MeshType = Mesh<ConnectivityType>;
+ using MeshType = Mesh<Dimension>;
+ using ConnectivityType = typename MeshType::Connectivity;
using R3 = TinyVector<3>;
@@ -1350,8 +1350,8 @@ TEST_CASE("MeshFlatNodeBoundary", "[mesh]")
{
static constexpr size_t Dimension = 2;
- using ConnectivityType = Connectivity<Dimension>;
- using MeshType = Mesh<ConnectivityType>;
+ using MeshType = Mesh<Dimension>;
+ using ConnectivityType = typename MeshType::Connectivity;
using R2 = TinyVector<2>;
@@ -1451,8 +1451,8 @@ TEST_CASE("MeshFlatNodeBoundary", "[mesh]")
{
static constexpr size_t Dimension = 3;
- using ConnectivityType = Connectivity<Dimension>;
- using MeshType = Mesh<ConnectivityType>;
+ using MeshType = Mesh<Dimension>;
+ using ConnectivityType = typename MeshType::Connectivity;
using R3 = TinyVector<3>;
diff --git a/tests/test_MeshLineEdgeBoundary.cpp b/tests/test_MeshLineEdgeBoundary.cpp
index 7eb85182b2c8147c694e9dc71486a7567ca63dda..39a93e8880b5f616003a0c5d1c28860af6363124 100644
--- a/tests/test_MeshLineEdgeBoundary.cpp
+++ b/tests/test_MeshLineEdgeBoundary.cpp
@@ -50,8 +50,8 @@ TEST_CASE("MeshLineEdgeBoundary", "[mesh]")
{
static constexpr size_t Dimension = 2;
- using ConnectivityType = Connectivity<Dimension>;
- using MeshType = Mesh<ConnectivityType>;
+ using MeshType = Mesh<Dimension>;
+ using ConnectivityType = typename MeshType::Connectivity;
using R2 = TinyVector<2>;
@@ -253,8 +253,8 @@ TEST_CASE("MeshLineEdgeBoundary", "[mesh]")
{
static constexpr size_t Dimension = 3;
- using ConnectivityType = Connectivity<Dimension>;
- using MeshType = Mesh<ConnectivityType>;
+ using MeshType = Mesh<Dimension>;
+ using ConnectivityType = typename MeshType::Connectivity;
using R3 = TinyVector<3>;
@@ -552,8 +552,8 @@ TEST_CASE("MeshLineEdgeBoundary", "[mesh]")
{
static constexpr size_t Dimension = 2;
- using ConnectivityType = Connectivity<Dimension>;
- using MeshType = Mesh<ConnectivityType>;
+ using MeshType = Mesh<Dimension>;
+ using ConnectivityType = typename MeshType::Connectivity;
using R2 = TinyVector<2>;
@@ -732,8 +732,8 @@ TEST_CASE("MeshLineEdgeBoundary", "[mesh]")
{
static constexpr size_t Dimension = 3;
- using ConnectivityType = Connectivity<Dimension>;
- using MeshType = Mesh<ConnectivityType>;
+ using MeshType = Mesh<Dimension>;
+ using ConnectivityType = typename MeshType::Connectivity;
using R3 = TinyVector<3>;
@@ -929,8 +929,8 @@ TEST_CASE("MeshLineEdgeBoundary", "[mesh]")
{
static constexpr size_t Dimension = 2;
- using ConnectivityType = Connectivity<Dimension>;
- using MeshType = Mesh<ConnectivityType>;
+ using MeshType = Mesh<Dimension>;
+ using ConnectivityType = typename MeshType::Connectivity;
using R2 = TinyVector<2>;
@@ -1030,8 +1030,8 @@ TEST_CASE("MeshLineEdgeBoundary", "[mesh]")
{
static constexpr size_t Dimension = 3;
- using ConnectivityType = Connectivity<Dimension>;
- using MeshType = Mesh<ConnectivityType>;
+ using MeshType = Mesh<Dimension>;
+ using ConnectivityType = typename MeshType::Connectivity;
using R3 = TinyVector<3>;
diff --git a/tests/test_MeshLineFaceBoundary.cpp b/tests/test_MeshLineFaceBoundary.cpp
index dc0a1b8f86fab920de1e12f44d70e8d6cde75221..5599862251d7e9c3e8f8d759e350e2776d0fa013 100644
--- a/tests/test_MeshLineFaceBoundary.cpp
+++ b/tests/test_MeshLineFaceBoundary.cpp
@@ -50,8 +50,8 @@ TEST_CASE("MeshLineFaceBoundary", "[mesh]")
{
static constexpr size_t Dimension = 2;
- using ConnectivityType = Connectivity<Dimension>;
- using MeshType = Mesh<ConnectivityType>;
+ using MeshType = Mesh<Dimension>;
+ using ConnectivityType = typename MeshType::Connectivity;
using R2 = TinyVector<2>;
@@ -256,8 +256,8 @@ TEST_CASE("MeshLineFaceBoundary", "[mesh]")
{
static constexpr size_t Dimension = 2;
- using ConnectivityType = Connectivity<Dimension>;
- using MeshType = Mesh<ConnectivityType>;
+ using MeshType = Mesh<Dimension>;
+ using ConnectivityType = typename MeshType::Connectivity;
using R2 = TinyVector<2>;
@@ -439,8 +439,8 @@ TEST_CASE("MeshLineFaceBoundary", "[mesh]")
{
static constexpr size_t Dimension = 2;
- using ConnectivityType = Connectivity<Dimension>;
- using MeshType = Mesh<ConnectivityType>;
+ using MeshType = Mesh<Dimension>;
+ using ConnectivityType = typename MeshType::Connectivity;
using R2 = TinyVector<2>;
diff --git a/tests/test_MeshLineNodeBoundary.cpp b/tests/test_MeshLineNodeBoundary.cpp
index cccc6568671628382cb44c8be123895535897b11..748c4533528fcbd137a4c88b83601e17ea48793a 100644
--- a/tests/test_MeshLineNodeBoundary.cpp
+++ b/tests/test_MeshLineNodeBoundary.cpp
@@ -50,8 +50,8 @@ TEST_CASE("MeshLineNodeBoundary", "[mesh]")
{
static constexpr size_t Dimension = 2;
- using ConnectivityType = Connectivity<Dimension>;
- using MeshType = Mesh<ConnectivityType>;
+ using MeshType = Mesh<Dimension>;
+ using ConnectivityType = typename MeshType::Connectivity;
using R2 = TinyVector<2>;
@@ -253,8 +253,8 @@ TEST_CASE("MeshLineNodeBoundary", "[mesh]")
{
static constexpr size_t Dimension = 3;
- using ConnectivityType = Connectivity<Dimension>;
- using MeshType = Mesh<ConnectivityType>;
+ using MeshType = Mesh<Dimension>;
+ using ConnectivityType = typename MeshType::Connectivity;
using R3 = TinyVector<3>;
@@ -552,8 +552,8 @@ TEST_CASE("MeshLineNodeBoundary", "[mesh]")
{
static constexpr size_t Dimension = 2;
- using ConnectivityType = Connectivity<Dimension>;
- using MeshType = Mesh<ConnectivityType>;
+ using MeshType = Mesh<Dimension>;
+ using ConnectivityType = typename MeshType::Connectivity;
using R2 = TinyVector<2>;
@@ -732,8 +732,8 @@ TEST_CASE("MeshLineNodeBoundary", "[mesh]")
{
static constexpr size_t Dimension = 3;
- using ConnectivityType = Connectivity<Dimension>;
- using MeshType = Mesh<ConnectivityType>;
+ using MeshType = Mesh<Dimension>;
+ using ConnectivityType = typename MeshType::Connectivity;
using R3 = TinyVector<3>;
@@ -929,8 +929,8 @@ TEST_CASE("MeshLineNodeBoundary", "[mesh]")
{
static constexpr size_t Dimension = 2;
- using ConnectivityType = Connectivity<Dimension>;
- using MeshType = Mesh<ConnectivityType>;
+ using MeshType = Mesh<Dimension>;
+ using ConnectivityType = typename MeshType::Connectivity;
using R2 = TinyVector<2>;
@@ -1030,8 +1030,8 @@ TEST_CASE("MeshLineNodeBoundary", "[mesh]")
{
static constexpr size_t Dimension = 3;
- using ConnectivityType = Connectivity<Dimension>;
- using MeshType = Mesh<ConnectivityType>;
+ using MeshType = Mesh<Dimension>;
+ using ConnectivityType = typename MeshType::Connectivity;
using R3 = TinyVector<3>;
diff --git a/tests/test_MeshNodeBoundary.cpp b/tests/test_MeshNodeBoundary.cpp
index b91de432f324b041bdf5aab4b6eefa726cac661a..68d398095ae3764c30a8fe03c365236453b00dfb 100644
--- a/tests/test_MeshNodeBoundary.cpp
+++ b/tests/test_MeshNodeBoundary.cpp
@@ -47,8 +47,8 @@ TEST_CASE("MeshNodeBoundary", "[mesh]")
{
static constexpr size_t Dimension = 1;
- using ConnectivityType = Connectivity<Dimension>;
- using MeshType = Mesh<ConnectivityType>;
+ using MeshType = Mesh<Dimension>;
+ using ConnectivityType = typename MeshType::Connectivity;
SECTION("cartesian 1d")
{
@@ -120,8 +120,8 @@ TEST_CASE("MeshNodeBoundary", "[mesh]")
{
static constexpr size_t Dimension = 2;
- using ConnectivityType = Connectivity<Dimension>;
- using MeshType = Mesh<ConnectivityType>;
+ using MeshType = Mesh<Dimension>;
+ using ConnectivityType = typename MeshType::Connectivity;
SECTION("cartesian 2d")
{
@@ -194,8 +194,8 @@ TEST_CASE("MeshNodeBoundary", "[mesh]")
{
static constexpr size_t Dimension = 3;
- using ConnectivityType = Connectivity<Dimension>;
- using MeshType = Mesh<ConnectivityType>;
+ using MeshType = Mesh<Dimension>;
+ using ConnectivityType = typename MeshType::Connectivity;
SECTION("cartesian 3d")
{
@@ -284,8 +284,7 @@ TEST_CASE("MeshNodeBoundary", "[mesh]")
{
static constexpr size_t Dimension = 3;
- using ConnectivityType = Connectivity<Dimension>;
- using MeshType = Mesh<ConnectivityType>;
+ using MeshType = Mesh<Dimension>;
std::shared_ptr p_mesh = MeshDataBaseForTests::get().hybrid3DMesh();
const MeshType& mesh = *p_mesh->get<MeshType>();
@@ -302,8 +301,7 @@ TEST_CASE("MeshNodeBoundary", "[mesh]")
{
static constexpr size_t Dimension = 1;
- using ConnectivityType = Connectivity<Dimension>;
- using MeshType = Mesh<ConnectivityType>;
+ using MeshType = Mesh<Dimension>;
std::shared_ptr p_mesh = MeshDataBaseForTests::get().unordered1DMesh();
const MeshType& mesh = *p_mesh->get<MeshType>();
@@ -319,8 +317,7 @@ TEST_CASE("MeshNodeBoundary", "[mesh]")
{
static constexpr size_t Dimension = 2;
- using ConnectivityType = Connectivity<Dimension>;
- using MeshType = Mesh<ConnectivityType>;
+ using MeshType = Mesh<Dimension>;
std::shared_ptr p_mesh = MeshDataBaseForTests::get().hybrid2DMesh();
const MeshType& mesh = *p_mesh->get<MeshType>();
@@ -336,8 +333,7 @@ TEST_CASE("MeshNodeBoundary", "[mesh]")
{
static constexpr size_t Dimension = 3;
- using ConnectivityType = Connectivity<Dimension>;
- using MeshType = Mesh<ConnectivityType>;
+ using MeshType = Mesh<Dimension>;
std::shared_ptr p_mesh = MeshDataBaseForTests::get().hybrid3DMesh();
const MeshType& mesh = *p_mesh->get<MeshType>();
diff --git a/tests/test_MeshNodeInterface.cpp b/tests/test_MeshNodeInterface.cpp
index 88884df79561b6e5fd17de808051bc03ba24a6ec..f8ce8a5fb61c5166db9520868e0d999a0f9e7d0c 100644
--- a/tests/test_MeshNodeInterface.cpp
+++ b/tests/test_MeshNodeInterface.cpp
@@ -47,8 +47,8 @@ TEST_CASE("MeshNodeInterface", "[mesh]")
{
static constexpr size_t Dimension = 1;
- using ConnectivityType = Connectivity<Dimension>;
- using MeshType = Mesh<ConnectivityType>;
+ using MeshType = Mesh<Dimension>;
+ using ConnectivityType = typename MeshType::Connectivity;
SECTION("unordered 1d")
{
@@ -87,8 +87,8 @@ TEST_CASE("MeshNodeInterface", "[mesh]")
{
static constexpr size_t Dimension = 2;
- using ConnectivityType = Connectivity<Dimension>;
- using MeshType = Mesh<ConnectivityType>;
+ using MeshType = Mesh<Dimension>;
+ using ConnectivityType = typename MeshType::Connectivity;
SECTION("hybrid 2d")
{
@@ -127,8 +127,8 @@ TEST_CASE("MeshNodeInterface", "[mesh]")
{
static constexpr size_t Dimension = 3;
- using ConnectivityType = Connectivity<Dimension>;
- using MeshType = Mesh<ConnectivityType>;
+ using MeshType = Mesh<Dimension>;
+ using ConnectivityType = typename MeshType::Connectivity;
SECTION("hybrid 3d")
{
@@ -169,8 +169,7 @@ TEST_CASE("MeshNodeInterface", "[mesh]")
{
static constexpr size_t Dimension = 3;
- using ConnectivityType = Connectivity<Dimension>;
- using MeshType = Mesh<ConnectivityType>;
+ using MeshType = Mesh<Dimension>;
std::shared_ptr p_mesh = MeshDataBaseForTests::get().hybrid3DMesh();
const MeshType& mesh = *p_mesh->get<MeshType>();
@@ -187,8 +186,7 @@ TEST_CASE("MeshNodeInterface", "[mesh]")
{
static constexpr size_t Dimension = 1;
- using ConnectivityType = Connectivity<Dimension>;
- using MeshType = Mesh<ConnectivityType>;
+ using MeshType = Mesh<Dimension>;
std::shared_ptr p_mesh = MeshDataBaseForTests::get().unordered1DMesh();
const MeshType& mesh = *p_mesh->get<MeshType>();
@@ -204,8 +202,7 @@ TEST_CASE("MeshNodeInterface", "[mesh]")
{
static constexpr size_t Dimension = 2;
- using ConnectivityType = Connectivity<Dimension>;
- using MeshType = Mesh<ConnectivityType>;
+ using MeshType = Mesh<Dimension>;
std::shared_ptr p_mesh = MeshDataBaseForTests::get().hybrid2DMesh();
const MeshType& mesh = *p_mesh->get<MeshType>();
@@ -221,8 +218,7 @@ TEST_CASE("MeshNodeInterface", "[mesh]")
{
static constexpr size_t Dimension = 3;
- using ConnectivityType = Connectivity<Dimension>;
- using MeshType = Mesh<ConnectivityType>;
+ using MeshType = Mesh<Dimension>;
std::shared_ptr p_mesh = MeshDataBaseForTests::get().hybrid3DMesh();
const MeshType& mesh = *p_mesh->get<MeshType>();
diff --git a/tests/test_ParallelChecker_read.cpp b/tests/test_ParallelChecker_read.cpp
index b7996db7bc955107c260aa93ae96dffe4356184f..026cd4063fc0837dcec9a83bac312d4d58619561 100644
--- a/tests/test_ParallelChecker_read.cpp
+++ b/tests/test_ParallelChecker_read.cpp
@@ -199,7 +199,7 @@ TEST_CASE("ParallelChecker_read", "[dev]")
{
// ItemValues
{ // 1d
- auto mesh = MeshDataBaseForTests::get().unordered1DMesh()->get<Mesh<Connectivity<1>>>();
+ auto mesh = MeshDataBaseForTests::get().unordered1DMesh()->get<Mesh<1>>();
std::string filename = ParallelChecker::instance().filename();
const Connectivity<1>& connectivity = mesh->connectivity();
@@ -363,7 +363,7 @@ TEST_CASE("ParallelChecker_read", "[dev]")
}
{
- auto other_mesh = MeshDataBaseForTests::get().cartesian1DMesh()->get<Mesh<Connectivity<1>>>();
+ auto other_mesh = MeshDataBaseForTests::get().cartesian1DMesh()->get<Mesh<1>>();
const Connectivity<1>& other_connectivity = other_mesh->connectivity();
CellValue<double> other_shape{other_connectivity};
@@ -377,7 +377,7 @@ TEST_CASE("ParallelChecker_read", "[dev]")
// ItemArray
{ // 1d
- auto mesh = MeshDataBaseForTests::get().unordered1DMesh()->get<Mesh<Connectivity<1>>>();
+ auto mesh = MeshDataBaseForTests::get().unordered1DMesh()->get<Mesh<1>>();
std::string filename = ParallelChecker::instance().filename();
const Connectivity<1>& connectivity = mesh->connectivity();
@@ -503,7 +503,7 @@ TEST_CASE("ParallelChecker_read", "[dev]")
}
{
- auto other_mesh = MeshDataBaseForTests::get().cartesian1DMesh()->get<Mesh<Connectivity<1>>>();
+ auto other_mesh = MeshDataBaseForTests::get().cartesian1DMesh()->get<Mesh<1>>();
const Connectivity<1>& other_connectivity = other_mesh->connectivity();
CellArray<double> other_shape{other_connectivity, 2};
@@ -517,7 +517,7 @@ TEST_CASE("ParallelChecker_read", "[dev]")
// ItemValues
{ // 2d
- auto mesh = MeshDataBaseForTests::get().hybrid2DMesh()->get<Mesh<Connectivity<2>>>();
+ auto mesh = MeshDataBaseForTests::get().hybrid2DMesh()->get<Mesh<2>>();
std::string filename = ParallelChecker::instance().filename();
const Connectivity<2>& connectivity = mesh->connectivity();
@@ -610,7 +610,7 @@ TEST_CASE("ParallelChecker_read", "[dev]")
// ItemArray
{ // 2d
- auto mesh = MeshDataBaseForTests::get().hybrid2DMesh()->get<Mesh<Connectivity<2>>>();
+ auto mesh = MeshDataBaseForTests::get().hybrid2DMesh()->get<Mesh<2>>();
std::string filename = ParallelChecker::instance().filename();
const Connectivity<2>& connectivity = mesh->connectivity();
@@ -747,7 +747,7 @@ TEST_CASE("ParallelChecker_read", "[dev]")
}
{
- auto other_mesh = MeshDataBaseForTests::get().cartesian2DMesh()->get<Mesh<Connectivity<2>>>();
+ auto other_mesh = MeshDataBaseForTests::get().cartesian2DMesh()->get<Mesh<2>>();
const Connectivity<2>& other_connectivity = other_mesh->connectivity();
FaceArray<DataType> other_shape{other_connectivity, 2};
@@ -761,7 +761,7 @@ TEST_CASE("ParallelChecker_read", "[dev]")
// ItemValues
{ // 3d
- auto mesh = MeshDataBaseForTests::get().hybrid3DMesh()->get<Mesh<Connectivity<3>>>();
+ auto mesh = MeshDataBaseForTests::get().hybrid3DMesh()->get<Mesh<3>>();
std::string filename = ParallelChecker::instance().filename();
const Connectivity<3>& connectivity = mesh->connectivity();
@@ -858,7 +858,7 @@ TEST_CASE("ParallelChecker_read", "[dev]")
// ItemArray
{ // 3d
- auto mesh = MeshDataBaseForTests::get().hybrid3DMesh()->get<Mesh<Connectivity<3>>>();
+ auto mesh = MeshDataBaseForTests::get().hybrid3DMesh()->get<Mesh<3>>();
std::string filename = ParallelChecker::instance().filename();
const Connectivity<3>& connectivity = mesh->connectivity();
@@ -991,7 +991,7 @@ TEST_CASE("ParallelChecker_read", "[dev]")
}
{
- auto other_mesh = MeshDataBaseForTests::get().cartesian2DMesh()->get<Mesh<Connectivity<2>>>();
+ auto other_mesh = MeshDataBaseForTests::get().cartesian2DMesh()->get<Mesh<2>>();
const Connectivity<2>& other_connectivity = other_mesh->connectivity();
FaceArray<DataType> other_shape{other_connectivity, 2};
@@ -1007,7 +1007,7 @@ TEST_CASE("ParallelChecker_read", "[dev]")
{
// SubItemValuePerItem
{ // 1d
- auto mesh = MeshDataBaseForTests::get().unordered1DMesh()->get<Mesh<Connectivity<1>>>();
+ auto mesh = MeshDataBaseForTests::get().unordered1DMesh()->get<Mesh<1>>();
std::string filename = ParallelChecker::instance().filename();
const Connectivity<1>& connectivity = mesh->connectivity();
@@ -1198,7 +1198,7 @@ TEST_CASE("ParallelChecker_read", "[dev]")
}
{
- auto other_mesh = MeshDataBaseForTests::get().cartesian1DMesh()->get<Mesh<Connectivity<1>>>();
+ auto other_mesh = MeshDataBaseForTests::get().cartesian1DMesh()->get<Mesh<1>>();
const Connectivity<1>& other_connectivity = other_mesh->connectivity();
CellValuePerNode<double> other_shape{other_connectivity};
@@ -1212,7 +1212,7 @@ TEST_CASE("ParallelChecker_read", "[dev]")
// SubItemArrayPerItem
{ // 1d
- auto mesh = MeshDataBaseForTests::get().unordered1DMesh()->get<Mesh<Connectivity<1>>>();
+ auto mesh = MeshDataBaseForTests::get().unordered1DMesh()->get<Mesh<1>>();
std::string filename = ParallelChecker::instance().filename();
const Connectivity<1>& connectivity = mesh->connectivity();
@@ -1352,7 +1352,7 @@ TEST_CASE("ParallelChecker_read", "[dev]")
}
{
- auto other_mesh = MeshDataBaseForTests::get().cartesian1DMesh()->get<Mesh<Connectivity<1>>>();
+ auto other_mesh = MeshDataBaseForTests::get().cartesian1DMesh()->get<Mesh<1>>();
const Connectivity<1>& other_connectivity = other_mesh->connectivity();
CellArray<double> other_shape{other_connectivity, 2};
@@ -1366,7 +1366,7 @@ TEST_CASE("ParallelChecker_read", "[dev]")
// SubItemValuePerItem
{ // 2d
- auto mesh = MeshDataBaseForTests::get().hybrid2DMesh()->get<Mesh<Connectivity<2>>>();
+ auto mesh = MeshDataBaseForTests::get().hybrid2DMesh()->get<Mesh<2>>();
std::string filename = ParallelChecker::instance().filename();
const Connectivity<2>& connectivity = mesh->connectivity();
@@ -1503,7 +1503,7 @@ TEST_CASE("ParallelChecker_read", "[dev]")
}
{
- auto other_mesh = MeshDataBaseForTests::get().cartesian2DMesh()->get<Mesh<Connectivity<2>>>();
+ auto other_mesh = MeshDataBaseForTests::get().cartesian2DMesh()->get<Mesh<2>>();
const Connectivity<2>& other_connectivity = other_mesh->connectivity();
FaceArray<DataType> other_shape{other_connectivity, 2};
@@ -1517,7 +1517,7 @@ TEST_CASE("ParallelChecker_read", "[dev]")
// SubItemArrayPerItem
{ // 2d
- auto mesh = MeshDataBaseForTests::get().hybrid2DMesh()->get<Mesh<Connectivity<2>>>();
+ auto mesh = MeshDataBaseForTests::get().hybrid2DMesh()->get<Mesh<2>>();
std::string filename = ParallelChecker::instance().filename();
const Connectivity<2>& connectivity = mesh->connectivity();
@@ -1708,7 +1708,7 @@ TEST_CASE("ParallelChecker_read", "[dev]")
}
{
- auto other_mesh = MeshDataBaseForTests::get().cartesian2DMesh()->get<Mesh<Connectivity<2>>>();
+ auto other_mesh = MeshDataBaseForTests::get().cartesian2DMesh()->get<Mesh<2>>();
const Connectivity<2>& other_connectivity = other_mesh->connectivity();
CellValuePerNode<double> other_shape{other_connectivity};
@@ -1722,7 +1722,7 @@ TEST_CASE("ParallelChecker_read", "[dev]")
// SubItemValuePerItem
{ // 3d
- auto mesh = MeshDataBaseForTests::get().hybrid3DMesh()->get<Mesh<Connectivity<3>>>();
+ auto mesh = MeshDataBaseForTests::get().hybrid3DMesh()->get<Mesh<3>>();
std::string filename = ParallelChecker::instance().filename();
const Connectivity<3>& connectivity = mesh->connectivity();
@@ -1855,7 +1855,7 @@ TEST_CASE("ParallelChecker_read", "[dev]")
}
{
- auto other_mesh = MeshDataBaseForTests::get().cartesian2DMesh()->get<Mesh<Connectivity<2>>>();
+ auto other_mesh = MeshDataBaseForTests::get().cartesian2DMesh()->get<Mesh<2>>();
const Connectivity<2>& other_connectivity = other_mesh->connectivity();
FaceArray<DataType> other_shape{other_connectivity, 2};
@@ -1868,7 +1868,7 @@ TEST_CASE("ParallelChecker_read", "[dev]")
// SubItemArrayPerItem
{ // 3d
- auto mesh = MeshDataBaseForTests::get().hybrid3DMesh()->get<Mesh<Connectivity<3>>>();
+ auto mesh = MeshDataBaseForTests::get().hybrid3DMesh()->get<Mesh<3>>();
std::string filename = ParallelChecker::instance().filename();
const Connectivity<3>& connectivity = mesh->connectivity();
@@ -2008,7 +2008,7 @@ TEST_CASE("ParallelChecker_read", "[dev]")
}
{
- auto other_mesh = MeshDataBaseForTests::get().cartesian3DMesh()->get<Mesh<Connectivity<3>>>();
+ auto other_mesh = MeshDataBaseForTests::get().cartesian3DMesh()->get<Mesh<3>>();
const Connectivity<3>& other_connectivity = other_mesh->connectivity();
CellArray<double> other_shape{other_connectivity, 2};
@@ -2046,7 +2046,7 @@ TEST_CASE("ParallelChecker_read", "[dev]")
REQUIRE_NOTHROW(ParallelChecker::instance().setMode(ParallelChecker::Mode::read));
REQUIRE_NOTHROW(not ParallelChecker::instance().isWriting());
- auto mesh = MeshDataBaseForTests::get().unordered1DMesh()->get<Mesh<Connectivity<1>>>();
+ auto mesh = MeshDataBaseForTests::get().unordered1DMesh()->get<Mesh<1>>();
const Connectivity<1>& connectivity = mesh->connectivity();
diff --git a/tests/test_ParallelChecker_write.cpp b/tests/test_ParallelChecker_write.cpp
index d9237c7daa4fb8b96b57498085ab4bdea128b028..e3416599966f65429763e466ac0d51884dfe7f8b 100644
--- a/tests/test_ParallelChecker_write.cpp
+++ b/tests/test_ParallelChecker_write.cpp
@@ -131,7 +131,7 @@ TEST_CASE("ParallelChecker_write", "[dev]")
// ItemValues
{ // 1d
- auto mesh = MeshDataBaseForTests::get().unordered1DMesh()->get<Mesh<Connectivity<1>>>();
+ auto mesh = MeshDataBaseForTests::get().unordered1DMesh()->get<Mesh<1>>();
const Connectivity<1>& connectivity = mesh->connectivity();
@@ -172,7 +172,7 @@ TEST_CASE("ParallelChecker_write", "[dev]")
}
{ // 2d
- auto mesh = MeshDataBaseForTests::get().hybrid2DMesh()->get<Mesh<Connectivity<2>>>();
+ auto mesh = MeshDataBaseForTests::get().hybrid2DMesh()->get<Mesh<2>>();
const Connectivity<2>& connectivity = mesh->connectivity();
@@ -202,7 +202,7 @@ TEST_CASE("ParallelChecker_write", "[dev]")
}
{ // 3d
- auto mesh = MeshDataBaseForTests::get().hybrid3DMesh()->get<Mesh<Connectivity<3>>>();
+ auto mesh = MeshDataBaseForTests::get().hybrid3DMesh()->get<Mesh<3>>();
const Connectivity<3>& connectivity = mesh->connectivity();
@@ -244,7 +244,7 @@ TEST_CASE("ParallelChecker_write", "[dev]")
// ItemArrays
{ // 1d
- auto mesh = MeshDataBaseForTests::get().unordered1DMesh()->get<Mesh<Connectivity<1>>>();
+ auto mesh = MeshDataBaseForTests::get().unordered1DMesh()->get<Mesh<1>>();
const Connectivity<1>& connectivity = mesh->connectivity();
@@ -285,7 +285,7 @@ TEST_CASE("ParallelChecker_write", "[dev]")
}
{ // 2d
- auto mesh = MeshDataBaseForTests::get().hybrid2DMesh()->get<Mesh<Connectivity<2>>>();
+ auto mesh = MeshDataBaseForTests::get().hybrid2DMesh()->get<Mesh<2>>();
const Connectivity<2>& connectivity = mesh->connectivity();
@@ -315,7 +315,7 @@ TEST_CASE("ParallelChecker_write", "[dev]")
}
{ // 3d
- auto mesh = MeshDataBaseForTests::get().hybrid3DMesh()->get<Mesh<Connectivity<3>>>();
+ auto mesh = MeshDataBaseForTests::get().hybrid3DMesh()->get<Mesh<3>>();
const Connectivity<3>& connectivity = mesh->connectivity();
@@ -394,7 +394,7 @@ TEST_CASE("ParallelChecker_write", "[dev]")
// ItemValues
{ // 1d
- auto mesh = MeshDataBaseForTests::get().unordered1DMesh()->get<Mesh<Connectivity<1>>>();
+ auto mesh = MeshDataBaseForTests::get().unordered1DMesh()->get<Mesh<1>>();
const Connectivity<1>& connectivity = mesh->connectivity();
@@ -424,7 +424,7 @@ TEST_CASE("ParallelChecker_write", "[dev]")
}
{ // 2d
- auto mesh = MeshDataBaseForTests::get().hybrid2DMesh()->get<Mesh<Connectivity<2>>>();
+ auto mesh = MeshDataBaseForTests::get().hybrid2DMesh()->get<Mesh<2>>();
const Connectivity<2>& connectivity = mesh->connectivity();
@@ -454,7 +454,7 @@ TEST_CASE("ParallelChecker_write", "[dev]")
}
{ // 3d
- auto mesh = MeshDataBaseForTests::get().hybrid3DMesh()->get<Mesh<Connectivity<3>>>();
+ auto mesh = MeshDataBaseForTests::get().hybrid3DMesh()->get<Mesh<3>>();
const Connectivity<3>& connectivity = mesh->connectivity();
@@ -485,7 +485,7 @@ TEST_CASE("ParallelChecker_write", "[dev]")
// ItemArrays
{ // 1d
- auto mesh = MeshDataBaseForTests::get().unordered1DMesh()->get<Mesh<Connectivity<1>>>();
+ auto mesh = MeshDataBaseForTests::get().unordered1DMesh()->get<Mesh<1>>();
const Connectivity<1>& connectivity = mesh->connectivity();
@@ -515,7 +515,7 @@ TEST_CASE("ParallelChecker_write", "[dev]")
}
{ // 2d
- auto mesh = MeshDataBaseForTests::get().hybrid2DMesh()->get<Mesh<Connectivity<2>>>();
+ auto mesh = MeshDataBaseForTests::get().hybrid2DMesh()->get<Mesh<2>>();
const Connectivity<2>& connectivity = mesh->connectivity();
@@ -545,7 +545,7 @@ TEST_CASE("ParallelChecker_write", "[dev]")
}
{ // 3d
- auto mesh = MeshDataBaseForTests::get().hybrid3DMesh()->get<Mesh<Connectivity<3>>>();
+ auto mesh = MeshDataBaseForTests::get().hybrid3DMesh()->get<Mesh<3>>();
const Connectivity<3>& connectivity = mesh->connectivity();
@@ -620,7 +620,7 @@ TEST_CASE("ParallelChecker_write", "[dev]")
REQUIRE_NOTHROW(ParallelChecker::instance().setMode(ParallelChecker::Mode::automatic));
REQUIRE_NOTHROW(ParallelChecker::instance().isWriting() == (parallel::size() == 1));
- auto mesh = MeshDataBaseForTests::get().unordered1DMesh()->get<Mesh<Connectivity<1>>>();
+ auto mesh = MeshDataBaseForTests::get().unordered1DMesh()->get<Mesh<1>>();
const Connectivity<1>& connectivity = mesh->connectivity();
diff --git a/tests/test_PrimalToDiamondDualConnectivityDataMapper.cpp b/tests/test_PrimalToDiamondDualConnectivityDataMapper.cpp
index 76882047e56b119d5f9363d3d1643d8f814cb82c..68a2301df4dd28e3cc847f065c014bfc44ded6cc 100644
--- a/tests/test_PrimalToDiamondDualConnectivityDataMapper.cpp
+++ b/tests/test_PrimalToDiamondDualConnectivityDataMapper.cpp
@@ -17,8 +17,8 @@ TEST_CASE("PrimalToDiamondDualConnectivityDataMapper", "[mesh]")
{
constexpr static size_t Dimension = 1;
- using ConnectivityType = Connectivity<Dimension>;
- using MeshType = Mesh<ConnectivityType>;
+ using MeshType = Mesh<Dimension>;
+ using ConnectivityType = typename MeshType::Connectivity;
std::shared_ptr<const MeshType> mesh = MeshDataBaseForTests::get().unordered1DMesh()->get<MeshType>();
const ConnectivityType& primal_connectivity = mesh->connectivity();
@@ -44,8 +44,8 @@ TEST_CASE("PrimalToDiamondDualConnectivityDataMapper", "[mesh]")
{
constexpr static size_t Dimension = 2;
- using ConnectivityType = Connectivity<Dimension>;
- using MeshType = Mesh<ConnectivityType>;
+ using MeshType = Mesh<Dimension>;
+ using ConnectivityType = typename MeshType::Connectivity;
std::shared_ptr<const MeshType> mesh = MeshDataBaseForTests::get().hybrid2DMesh()->get<MeshType>();
const ConnectivityType& primal_connectivity = mesh->connectivity();
@@ -223,8 +223,8 @@ TEST_CASE("PrimalToDiamondDualConnectivityDataMapper", "[mesh]")
{
constexpr static size_t Dimension = 3;
- using ConnectivityType = Connectivity<Dimension>;
- using MeshType = Mesh<ConnectivityType>;
+ using MeshType = Mesh<Dimension>;
+ using ConnectivityType = typename MeshType::Connectivity;
std::shared_ptr<const MeshType> mesh = MeshDataBaseForTests::get().hybrid3DMesh()->get<MeshType>();
const ConnectivityType& primal_connectivity = mesh->connectivity();
diff --git a/tests/test_PrimalToDual1DConnectivityDataMapper.cpp b/tests/test_PrimalToDual1DConnectivityDataMapper.cpp
index df4e580464d401d04a16578db5bc6f7c29226a93..10a6fc4b7a2691e356de6e9b146447296c1e7e6f 100644
--- a/tests/test_PrimalToDual1DConnectivityDataMapper.cpp
+++ b/tests/test_PrimalToDual1DConnectivityDataMapper.cpp
@@ -14,8 +14,8 @@ TEST_CASE("PrimalToDual1DConnectivityDataMapper", "[mesh]")
{
constexpr static size_t Dimension = 1;
- using ConnectivityType = Connectivity<Dimension>;
- using MeshType = Mesh<ConnectivityType>;
+ using MeshType = Mesh<Dimension>;
+ using ConnectivityType = typename MeshType::Connectivity;
std::shared_ptr<const MeshType> mesh = MeshDataBaseForTests::get().unordered1DMesh()->get<MeshType>();
const ConnectivityType& primal_connectivity = mesh->connectivity();
diff --git a/tests/test_PrimalToMedianDualConnectivityDataMapper.cpp b/tests/test_PrimalToMedianDualConnectivityDataMapper.cpp
index 7c9a25d68827c5fea6c2da4fc5ccbdaf2103444a..7fb1e02a06c57b8657570a6f51e4312b40ca9fe2 100644
--- a/tests/test_PrimalToMedianDualConnectivityDataMapper.cpp
+++ b/tests/test_PrimalToMedianDualConnectivityDataMapper.cpp
@@ -17,8 +17,8 @@ TEST_CASE("PrimalToMedianDualConnectivityDataMapper", "[mesh]")
{
constexpr static size_t Dimension = 1;
- using ConnectivityType = Connectivity<Dimension>;
- using MeshType = Mesh<ConnectivityType>;
+ using MeshType = Mesh<Dimension>;
+ using ConnectivityType = typename MeshType::Connectivity;
std::shared_ptr<const MeshType> mesh = MeshDataBaseForTests::get().unordered1DMesh()->get<MeshType>();
const ConnectivityType& primal_connectivity = mesh->connectivity();
@@ -44,8 +44,8 @@ TEST_CASE("PrimalToMedianDualConnectivityDataMapper", "[mesh]")
{
constexpr static size_t Dimension = 2;
- using ConnectivityType = Connectivity<Dimension>;
- using MeshType = Mesh<ConnectivityType>;
+ using MeshType = Mesh<Dimension>;
+ using ConnectivityType = typename MeshType::Connectivity;
std::shared_ptr<const MeshType> mesh = MeshDataBaseForTests::get().hybrid2DMesh()->get<MeshType>();
const ConnectivityType& primal_connectivity = mesh->connectivity();
diff --git a/tests/test_SubItemArrayPerItem.cpp b/tests/test_SubItemArrayPerItem.cpp
index f027c56be1899e730c93fe73421c376c1d52d746..23f32a41bc9cfd6a503a7e78e726106a8c79fdf9 100644
--- a/tests/test_SubItemArrayPerItem.cpp
+++ b/tests/test_SubItemArrayPerItem.cpp
@@ -69,7 +69,7 @@ TEST_CASE("SubItemArrayPerItem", "[mesh]")
for (const auto& named_mesh : mesh_list) {
SECTION(named_mesh.name())
{
- auto mesh_1d = named_mesh.mesh()->get<Mesh<Connectivity<1>>>();
+ auto mesh_1d = named_mesh.mesh()->get<Mesh<1>>();
const Connectivity<1>& connectivity = mesh_1d->connectivity();
@@ -192,7 +192,7 @@ TEST_CASE("SubItemArrayPerItem", "[mesh]")
for (const auto& named_mesh : mesh_list) {
SECTION(named_mesh.name())
{
- auto mesh_2d = named_mesh.mesh()->get<Mesh<Connectivity<2>>>();
+ auto mesh_2d = named_mesh.mesh()->get<Mesh<2>>();
const Connectivity<2>& connectivity = mesh_2d->connectivity();
@@ -358,7 +358,7 @@ TEST_CASE("SubItemArrayPerItem", "[mesh]")
for (const auto& named_mesh : mesh_list) {
SECTION(named_mesh.name())
{
- auto mesh_3d = named_mesh.mesh()->get<Mesh<Connectivity<3>>>();
+ auto mesh_3d = named_mesh.mesh()->get<Mesh<3>>();
const Connectivity<3>& connectivity = mesh_3d->connectivity();
@@ -555,7 +555,7 @@ TEST_CASE("SubItemArrayPerItem", "[mesh]")
for (const auto& named_mesh : mesh_list) {
SECTION(named_mesh.name())
{
- auto mesh_1d = named_mesh.mesh()->get<Mesh<Connectivity<1>>>();
+ auto mesh_1d = named_mesh.mesh()->get<Mesh<1>>();
const Connectivity<1>& connectivity = mesh_1d->connectivity();
@@ -628,7 +628,7 @@ TEST_CASE("SubItemArrayPerItem", "[mesh]")
for (const auto& named_mesh : mesh_list) {
SECTION(named_mesh.name())
{
- auto mesh_2d = named_mesh.mesh()->get<Mesh<Connectivity<2>>>();
+ auto mesh_2d = named_mesh.mesh()->get<Mesh<2>>();
const Connectivity<2>& connectivity = mesh_2d->connectivity();
@@ -700,7 +700,7 @@ TEST_CASE("SubItemArrayPerItem", "[mesh]")
for (const auto& named_mesh : mesh_list) {
SECTION(named_mesh.name())
{
- auto mesh_3d = named_mesh.mesh()->get<Mesh<Connectivity<3>>>();
+ auto mesh_3d = named_mesh.mesh()->get<Mesh<3>>();
const Connectivity<3>& connectivity = mesh_3d->connectivity();
@@ -772,7 +772,7 @@ TEST_CASE("SubItemArrayPerItem", "[mesh]")
for (const auto& named_mesh : mesh_list) {
SECTION(named_mesh.name())
{
- auto mesh_3d = named_mesh.mesh()->get<Mesh<Connectivity<3>>>();
+ auto mesh_3d = named_mesh.mesh()->get<Mesh<3>>();
const Connectivity<3>& connectivity = mesh_3d->connectivity();
@@ -883,7 +883,7 @@ TEST_CASE("SubItemArrayPerItem", "[mesh]")
for (const auto& named_mesh : mesh_list) {
SECTION(named_mesh.name())
{
- auto mesh_3d = named_mesh.mesh()->get<Mesh<Connectivity<3>>>();
+ auto mesh_3d = named_mesh.mesh()->get<Mesh<3>>();
const Connectivity<3>& connectivity = mesh_3d->connectivity();
@@ -994,7 +994,7 @@ TEST_CASE("SubItemArrayPerItem", "[mesh]")
for (const auto& named_mesh : mesh_list) {
SECTION(named_mesh.name())
{
- auto mesh_2d = named_mesh.mesh()->get<Mesh<Connectivity<2>>>();
+ auto mesh_2d = named_mesh.mesh()->get<Mesh<2>>();
const Connectivity<2>& connectivity = mesh_2d->connectivity();
@@ -1068,7 +1068,7 @@ TEST_CASE("SubItemArrayPerItem", "[mesh]")
SECTION("checking invalid table size in constructor")
{
- auto mesh_3d = MeshDataBaseForTests::get().hybrid3DMesh()->get<Mesh<Connectivity<3>>>();
+ auto mesh_3d = MeshDataBaseForTests::get().hybrid3DMesh()->get<Mesh<3>>();
const Connectivity<3>& connectivity = mesh_3d->connectivity();
@@ -1090,7 +1090,7 @@ TEST_CASE("SubItemArrayPerItem", "[mesh]")
for (const auto& named_mesh : mesh_list) {
SECTION(named_mesh.name())
{
- auto mesh_3d = named_mesh.mesh()->get<Mesh<Connectivity<3>>>();
+ auto mesh_3d = named_mesh.mesh()->get<Mesh<3>>();
const Connectivity<3>& connectivity = mesh_3d->connectivity();
diff --git a/tests/test_SubItemArrayPerItemUtils.cpp b/tests/test_SubItemArrayPerItemUtils.cpp
index 8d4c7516762f8efb3ec95654c32cef4a728ffa5a..1ff8bb17f7bf170c748cfa15aa3fc82f25f4a879 100644
--- a/tests/test_SubItemArrayPerItemUtils.cpp
+++ b/tests/test_SubItemArrayPerItemUtils.cpp
@@ -23,7 +23,7 @@ TEST_CASE("SubItemArrayPerItemUtils", "[mesh]")
for (const auto& named_mesh : mesh_list) {
SECTION(named_mesh.name())
{
- auto mesh_2d = named_mesh.mesh()->get<Mesh<Connectivity<2>>>();
+ auto mesh_2d = named_mesh.mesh()->get<Mesh<2>>();
const Connectivity<2>& connectivity = mesh_2d->connectivity();
@@ -95,7 +95,7 @@ TEST_CASE("SubItemArrayPerItemUtils", "[mesh]")
for (const auto& named_mesh : mesh_list) {
SECTION(named_mesh.name())
{
- auto mesh_1d = named_mesh.mesh()->get<Mesh<Connectivity<1>>>();
+ auto mesh_1d = named_mesh.mesh()->get<Mesh<1>>();
const Connectivity<1>& connectivity = mesh_1d->connectivity();
@@ -127,7 +127,7 @@ TEST_CASE("SubItemArrayPerItemUtils", "[mesh]")
for (const auto& named_mesh : mesh_list) {
SECTION(named_mesh.name())
{
- auto mesh_2d = named_mesh.mesh()->get<Mesh<Connectivity<2>>>();
+ auto mesh_2d = named_mesh.mesh()->get<Mesh<2>>();
const Connectivity<2>& connectivity = mesh_2d->connectivity();
@@ -159,7 +159,7 @@ TEST_CASE("SubItemArrayPerItemUtils", "[mesh]")
for (const auto& named_mesh : mesh_list) {
SECTION(named_mesh.name())
{
- auto mesh_3d = named_mesh.mesh()->get<Mesh<Connectivity<3>>>();
+ auto mesh_3d = named_mesh.mesh()->get<Mesh<3>>();
const Connectivity<3>& connectivity = mesh_3d->connectivity();
@@ -194,7 +194,7 @@ TEST_CASE("SubItemArrayPerItemUtils", "[mesh]")
for (const auto& named_mesh : mesh_list) {
SECTION(named_mesh.name())
{
- auto mesh_1d = named_mesh.mesh()->get<Mesh<Connectivity<1>>>();
+ auto mesh_1d = named_mesh.mesh()->get<Mesh<1>>();
const Connectivity<1>& connectivity = mesh_1d->connectivity();
@@ -226,7 +226,7 @@ TEST_CASE("SubItemArrayPerItemUtils", "[mesh]")
for (const auto& named_mesh : mesh_list) {
SECTION(named_mesh.name())
{
- auto mesh_2d = named_mesh.mesh()->get<Mesh<Connectivity<2>>>();
+ auto mesh_2d = named_mesh.mesh()->get<Mesh<2>>();
const Connectivity<2>& connectivity = mesh_2d->connectivity();
@@ -258,7 +258,7 @@ TEST_CASE("SubItemArrayPerItemUtils", "[mesh]")
for (const auto& named_mesh : mesh_list) {
SECTION(named_mesh.name())
{
- auto mesh_3d = named_mesh.mesh()->get<Mesh<Connectivity<3>>>();
+ auto mesh_3d = named_mesh.mesh()->get<Mesh<3>>();
const Connectivity<3>& connectivity = mesh_3d->connectivity();
@@ -294,7 +294,7 @@ TEST_CASE("SubItemArrayPerItemUtils", "[mesh]")
for (const auto& named_mesh : mesh_list) {
SECTION(named_mesh.name())
{
- auto mesh_1d = named_mesh.mesh()->get<Mesh<Connectivity<1>>>();
+ auto mesh_1d = named_mesh.mesh()->get<Mesh<1>>();
const Connectivity<1>& connectivity = mesh_1d->connectivity();
@@ -324,7 +324,7 @@ TEST_CASE("SubItemArrayPerItemUtils", "[mesh]")
for (const auto& named_mesh : mesh_list) {
SECTION(named_mesh.name() + " for size_t data")
{
- auto mesh_2d = named_mesh.mesh()->get<Mesh<Connectivity<2>>>();
+ auto mesh_2d = named_mesh.mesh()->get<Mesh<2>>();
const Connectivity<2>& connectivity = mesh_2d->connectivity();
@@ -347,7 +347,7 @@ TEST_CASE("SubItemArrayPerItemUtils", "[mesh]")
SECTION(named_mesh.name() + " for N^2 data")
{
- auto mesh_2d = named_mesh.mesh()->get<Mesh<Connectivity<2>>>();
+ auto mesh_2d = named_mesh.mesh()->get<Mesh<2>>();
const Connectivity<2>& connectivity = mesh_2d->connectivity();
@@ -373,7 +373,7 @@ TEST_CASE("SubItemArrayPerItemUtils", "[mesh]")
SECTION(named_mesh.name() + " for float data")
{
- auto mesh_2d = named_mesh.mesh()->get<Mesh<Connectivity<2>>>();
+ auto mesh_2d = named_mesh.mesh()->get<Mesh<2>>();
const Connectivity<2>& connectivity = mesh_2d->connectivity();
@@ -417,7 +417,7 @@ TEST_CASE("SubItemArrayPerItemUtils", "[mesh]")
for (const auto& named_mesh : mesh_list) {
SECTION(named_mesh.name() + " for size_t data")
{
- auto mesh_3d = named_mesh.mesh()->get<Mesh<Connectivity<3>>>();
+ auto mesh_3d = named_mesh.mesh()->get<Mesh<3>>();
const Connectivity<3>& connectivity = mesh_3d->connectivity();
@@ -440,7 +440,7 @@ TEST_CASE("SubItemArrayPerItemUtils", "[mesh]")
SECTION(named_mesh.name() + " for N^2x3 data")
{
- auto mesh_3d = named_mesh.mesh()->get<Mesh<Connectivity<3>>>();
+ auto mesh_3d = named_mesh.mesh()->get<Mesh<3>>();
const Connectivity<3>& connectivity = mesh_3d->connectivity();
@@ -466,7 +466,7 @@ TEST_CASE("SubItemArrayPerItemUtils", "[mesh]")
SECTION(named_mesh.name() + " for R^2x3 data")
{
- auto mesh_3d = named_mesh.mesh()->get<Mesh<Connectivity<3>>>();
+ auto mesh_3d = named_mesh.mesh()->get<Mesh<3>>();
const Connectivity<3>& connectivity = mesh_3d->connectivity();
diff --git a/tests/test_SubItemArrayPerItemVariant.cpp b/tests/test_SubItemArrayPerItemVariant.cpp
index e4038670133b167a33b56a3ba71c4794e45eb8ec..7a224e10aa03cdd58ed184b1b92d13d470151259 100644
--- a/tests/test_SubItemArrayPerItemVariant.cpp
+++ b/tests/test_SubItemArrayPerItemVariant.cpp
@@ -11,9 +11,9 @@
TEST_CASE("SubItemArrayPerItemVariant", "[mesh]")
{
- std::shared_ptr mesh = MeshDataBaseForTests::get().hybrid3DMesh()->get<Mesh<Connectivity<3>>>();
+ std::shared_ptr mesh = MeshDataBaseForTests::get().hybrid3DMesh()->get<Mesh<3>>();
- const Connectivity<3>& connectivity = *mesh->shared_connectivity();
+ const Connectivity<3>& connectivity = mesh->connectivity();
using R1 = TinyVector<1>;
using R2 = TinyVector<2>;
diff --git a/tests/test_SubItemValuePerItem.cpp b/tests/test_SubItemValuePerItem.cpp
index 92569d9823ec90205f0ff9ff908b9c094a036d42..72b794bb1e8b5e0fdc3b07bdb6b63da9a4782683 100644
--- a/tests/test_SubItemValuePerItem.cpp
+++ b/tests/test_SubItemValuePerItem.cpp
@@ -79,7 +79,7 @@ TEST_CASE("SubItemValuePerItem", "[mesh]")
for (const auto& named_mesh : mesh_list) {
SECTION(named_mesh.name())
{
- auto mesh_1d = named_mesh.mesh()->get<Mesh<Connectivity<1>>>();
+ auto mesh_1d = named_mesh.mesh()->get<Mesh<1>>();
const Connectivity<1>& connectivity = mesh_1d->connectivity();
@@ -215,7 +215,7 @@ TEST_CASE("SubItemValuePerItem", "[mesh]")
for (const auto& named_mesh : mesh_list) {
SECTION(named_mesh.name())
{
- auto mesh_2d = named_mesh.mesh()->get<Mesh<Connectivity<2>>>();
+ auto mesh_2d = named_mesh.mesh()->get<Mesh<2>>();
const Connectivity<2>& connectivity = mesh_2d->connectivity();
@@ -391,7 +391,7 @@ TEST_CASE("SubItemValuePerItem", "[mesh]")
for (const auto& named_mesh : mesh_list) {
SECTION(named_mesh.name())
{
- auto mesh_3d = named_mesh.mesh()->get<Mesh<Connectivity<3>>>();
+ auto mesh_3d = named_mesh.mesh()->get<Mesh<3>>();
const Connectivity<3>& connectivity = mesh_3d->connectivity();
@@ -596,7 +596,7 @@ TEST_CASE("SubItemValuePerItem", "[mesh]")
for (const auto& named_mesh : mesh_list) {
SECTION(named_mesh.name())
{
- auto mesh_1d = named_mesh.mesh()->get<Mesh<Connectivity<1>>>();
+ auto mesh_1d = named_mesh.mesh()->get<Mesh<1>>();
const Connectivity<1>& connectivity = mesh_1d->connectivity();
@@ -712,7 +712,7 @@ TEST_CASE("SubItemValuePerItem", "[mesh]")
for (const auto& named_mesh : mesh_list) {
SECTION(named_mesh.name())
{
- auto mesh_2d = named_mesh.mesh()->get<Mesh<Connectivity<2>>>();
+ auto mesh_2d = named_mesh.mesh()->get<Mesh<2>>();
const Connectivity<2>& connectivity = mesh_2d->connectivity();
@@ -898,7 +898,7 @@ TEST_CASE("SubItemValuePerItem", "[mesh]")
for (const auto& named_mesh : mesh_list) {
SECTION(named_mesh.name())
{
- auto mesh_3d = named_mesh.mesh()->get<Mesh<Connectivity<3>>>();
+ auto mesh_3d = named_mesh.mesh()->get<Mesh<3>>();
const Connectivity<3>& connectivity = mesh_3d->connectivity();
@@ -1119,7 +1119,7 @@ TEST_CASE("SubItemValuePerItem", "[mesh]")
for (const auto& named_mesh : mesh_list) {
SECTION(named_mesh.name())
{
- auto mesh_1d = named_mesh.mesh()->get<Mesh<Connectivity<1>>>();
+ auto mesh_1d = named_mesh.mesh()->get<Mesh<1>>();
const Connectivity<1>& connectivity = mesh_1d->connectivity();
@@ -1164,7 +1164,7 @@ TEST_CASE("SubItemValuePerItem", "[mesh]")
for (const auto& named_mesh : mesh_list) {
SECTION(named_mesh.name())
{
- auto mesh_2d = named_mesh.mesh()->get<Mesh<Connectivity<2>>>();
+ auto mesh_2d = named_mesh.mesh()->get<Mesh<2>>();
const Connectivity<2>& connectivity = mesh_2d->connectivity();
@@ -1208,7 +1208,7 @@ TEST_CASE("SubItemValuePerItem", "[mesh]")
for (const auto& named_mesh : mesh_list) {
SECTION(named_mesh.name())
{
- auto mesh_3d = named_mesh.mesh()->get<Mesh<Connectivity<3>>>();
+ auto mesh_3d = named_mesh.mesh()->get<Mesh<3>>();
const Connectivity<3>& connectivity = mesh_3d->connectivity();
@@ -1254,7 +1254,7 @@ TEST_CASE("SubItemValuePerItem", "[mesh]")
for (const auto& named_mesh : mesh_list) {
SECTION(named_mesh.name())
{
- auto mesh_3d = named_mesh.mesh()->get<Mesh<Connectivity<3>>>();
+ auto mesh_3d = named_mesh.mesh()->get<Mesh<3>>();
const Connectivity<3>& connectivity = mesh_3d->connectivity();
@@ -1353,7 +1353,7 @@ TEST_CASE("SubItemValuePerItem", "[mesh]")
for (const auto& named_mesh : mesh_list) {
SECTION(named_mesh.name())
{
- auto mesh_2d = named_mesh.mesh()->get<Mesh<Connectivity<2>>>();
+ auto mesh_2d = named_mesh.mesh()->get<Mesh<2>>();
const Connectivity<2>& connectivity = mesh_2d->connectivity();
@@ -1420,7 +1420,7 @@ TEST_CASE("SubItemValuePerItem", "[mesh]")
SECTION("checking invalid array size in constructor")
{
- auto mesh_3d = MeshDataBaseForTests::get().hybrid3DMesh()->get<Mesh<Connectivity<3>>>();
+ auto mesh_3d = MeshDataBaseForTests::get().hybrid3DMesh()->get<Mesh<3>>();
const Connectivity<3>& connectivity = mesh_3d->connectivity();
@@ -1442,7 +1442,7 @@ TEST_CASE("SubItemValuePerItem", "[mesh]")
for (const auto& named_mesh : mesh_list) {
SECTION(named_mesh.name())
{
- auto mesh_3d = named_mesh.mesh()->get<Mesh<Connectivity<3>>>();
+ auto mesh_3d = named_mesh.mesh()->get<Mesh<3>>();
const Connectivity<3>& connectivity = mesh_3d->connectivity();
diff --git a/tests/test_SubItemValuePerItemUtils.cpp b/tests/test_SubItemValuePerItemUtils.cpp
index ee73b60721bd1d3a1bee958c25481b8149057603..69214c1c933ba6abb0c93e094317901263e70eed 100644
--- a/tests/test_SubItemValuePerItemUtils.cpp
+++ b/tests/test_SubItemValuePerItemUtils.cpp
@@ -23,7 +23,7 @@ TEST_CASE("SubItemValuePerItemUtils", "[mesh]")
for (const auto& named_mesh : mesh_list) {
SECTION(named_mesh.name())
{
- auto mesh_2d = named_mesh.mesh()->get<Mesh<Connectivity<2>>>();
+ auto mesh_2d = named_mesh.mesh()->get<Mesh<2>>();
const Connectivity<2>& connectivity = mesh_2d->connectivity();
@@ -90,7 +90,7 @@ TEST_CASE("SubItemValuePerItemUtils", "[mesh]")
for (const auto& named_mesh : mesh_list) {
SECTION(named_mesh.name())
{
- auto mesh_1d = named_mesh.mesh()->get<Mesh<Connectivity<1>>>();
+ auto mesh_1d = named_mesh.mesh()->get<Mesh<1>>();
const Connectivity<1>& connectivity = mesh_1d->connectivity();
@@ -120,7 +120,7 @@ TEST_CASE("SubItemValuePerItemUtils", "[mesh]")
for (const auto& named_mesh : mesh_list) {
SECTION(named_mesh.name())
{
- auto mesh_2d = named_mesh.mesh()->get<Mesh<Connectivity<2>>>();
+ auto mesh_2d = named_mesh.mesh()->get<Mesh<2>>();
const Connectivity<2>& connectivity = mesh_2d->connectivity();
@@ -150,7 +150,7 @@ TEST_CASE("SubItemValuePerItemUtils", "[mesh]")
for (const auto& named_mesh : mesh_list) {
SECTION(named_mesh.name())
{
- auto mesh_3d = named_mesh.mesh()->get<Mesh<Connectivity<3>>>();
+ auto mesh_3d = named_mesh.mesh()->get<Mesh<3>>();
const Connectivity<3>& connectivity = mesh_3d->connectivity();
@@ -183,7 +183,7 @@ TEST_CASE("SubItemValuePerItemUtils", "[mesh]")
for (const auto& named_mesh : mesh_list) {
SECTION(named_mesh.name())
{
- auto mesh_1d = named_mesh.mesh()->get<Mesh<Connectivity<1>>>();
+ auto mesh_1d = named_mesh.mesh()->get<Mesh<1>>();
const Connectivity<1>& connectivity = mesh_1d->connectivity();
@@ -213,7 +213,7 @@ TEST_CASE("SubItemValuePerItemUtils", "[mesh]")
for (const auto& named_mesh : mesh_list) {
SECTION(named_mesh.name())
{
- auto mesh_2d = named_mesh.mesh()->get<Mesh<Connectivity<2>>>();
+ auto mesh_2d = named_mesh.mesh()->get<Mesh<2>>();
const Connectivity<2>& connectivity = mesh_2d->connectivity();
@@ -243,7 +243,7 @@ TEST_CASE("SubItemValuePerItemUtils", "[mesh]")
for (const auto& named_mesh : mesh_list) {
SECTION(named_mesh.name())
{
- auto mesh_3d = named_mesh.mesh()->get<Mesh<Connectivity<3>>>();
+ auto mesh_3d = named_mesh.mesh()->get<Mesh<3>>();
const Connectivity<3>& connectivity = mesh_3d->connectivity();
@@ -278,7 +278,7 @@ TEST_CASE("SubItemValuePerItemUtils", "[mesh]")
for (const auto& named_mesh : mesh_list) {
SECTION(named_mesh.name())
{
- auto mesh_1d = named_mesh.mesh()->get<Mesh<Connectivity<1>>>();
+ auto mesh_1d = named_mesh.mesh()->get<Mesh<1>>();
const Connectivity<1>& connectivity = mesh_1d->connectivity();
@@ -307,7 +307,7 @@ TEST_CASE("SubItemValuePerItemUtils", "[mesh]")
for (const auto& named_mesh : mesh_list) {
SECTION(named_mesh.name() + " for size_t data")
{
- auto mesh_2d = named_mesh.mesh()->get<Mesh<Connectivity<2>>>();
+ auto mesh_2d = named_mesh.mesh()->get<Mesh<2>>();
const Connectivity<2>& connectivity = mesh_2d->connectivity();
@@ -329,7 +329,7 @@ TEST_CASE("SubItemValuePerItemUtils", "[mesh]")
SECTION(named_mesh.name() + " for R^2 data")
{
- auto mesh_2d = named_mesh.mesh()->get<Mesh<Connectivity<2>>>();
+ auto mesh_2d = named_mesh.mesh()->get<Mesh<2>>();
const Connectivity<2>& connectivity = mesh_2d->connectivity();
@@ -366,7 +366,7 @@ TEST_CASE("SubItemValuePerItemUtils", "[mesh]")
for (const auto& named_mesh : mesh_list) {
SECTION(named_mesh.name() + " for size_t data")
{
- auto mesh_3d = named_mesh.mesh()->get<Mesh<Connectivity<3>>>();
+ auto mesh_3d = named_mesh.mesh()->get<Mesh<3>>();
const Connectivity<3>& connectivity = mesh_3d->connectivity();
@@ -388,7 +388,7 @@ TEST_CASE("SubItemValuePerItemUtils", "[mesh]")
SECTION(named_mesh.name() + " for N^2x3 data")
{
- auto mesh_3d = named_mesh.mesh()->get<Mesh<Connectivity<3>>>();
+ auto mesh_3d = named_mesh.mesh()->get<Mesh<3>>();
const Connectivity<3>& connectivity = mesh_3d->connectivity();
@@ -413,7 +413,7 @@ TEST_CASE("SubItemValuePerItemUtils", "[mesh]")
SECTION(named_mesh.name() + " for R^3x2 float")
{
- auto mesh_3d = named_mesh.mesh()->get<Mesh<Connectivity<3>>>();
+ auto mesh_3d = named_mesh.mesh()->get<Mesh<3>>();
const Connectivity<3>& connectivity = mesh_3d->connectivity();
diff --git a/tests/test_SubItemValuePerItemVariant.cpp b/tests/test_SubItemValuePerItemVariant.cpp
index 01ff40651e5e63fbfd00ebf840c79ffad863a67d..28f62ace44001ac1a95c0f81e59b9bd1bf3e1215 100644
--- a/tests/test_SubItemValuePerItemVariant.cpp
+++ b/tests/test_SubItemValuePerItemVariant.cpp
@@ -11,9 +11,9 @@
TEST_CASE("SubItemValuePerItemVariant", "[mesh]")
{
- std::shared_ptr mesh = MeshDataBaseForTests::get().hybrid3DMesh()->get<Mesh<Connectivity<3>>>();
+ std::shared_ptr mesh = MeshDataBaseForTests::get().hybrid3DMesh()->get<Mesh<3>>();
- const Connectivity<3>& connectivity = *mesh->shared_connectivity();
+ const Connectivity<3>& connectivity = mesh->connectivity();
using R1 = TinyVector<1>;
using R2 = TinyVector<2>;
diff --git a/tests/test_Synchronizer.cpp b/tests/test_Synchronizer.cpp
index e14aff66a0ab8d82ee46368c527e918a9a18ca81..b54f9321aa3a01a9ce6114ef9d0194e01d4ee367 100644
--- a/tests/test_Synchronizer.cpp
+++ b/tests/test_Synchronizer.cpp
@@ -29,7 +29,7 @@ TEST_CASE("Synchronizer", "[mesh]")
using ConnectivityType = Connectivity<Dimension>;
const ConnectivityType& connectivity =
- MeshDataBaseForTests::get().unordered1DMesh()->get<Mesh<Connectivity<1>>>()->connectivity();
+ MeshDataBaseForTests::get().unordered1DMesh()->get<Mesh<1>>()->connectivity();
SECTION("synchonize NodeValue")
{
@@ -137,8 +137,7 @@ TEST_CASE("Synchronizer", "[mesh]")
constexpr size_t Dimension = 2;
using ConnectivityType = Connectivity<Dimension>;
- const ConnectivityType& connectivity =
- MeshDataBaseForTests::get().hybrid2DMesh()->get<Mesh<Connectivity<2>>>()->connectivity();
+ const ConnectivityType& connectivity = MeshDataBaseForTests::get().hybrid2DMesh()->get<Mesh<2>>()->connectivity();
SECTION("synchonize NodeValue")
{
@@ -246,8 +245,7 @@ TEST_CASE("Synchronizer", "[mesh]")
constexpr size_t Dimension = 3;
using ConnectivityType = Connectivity<Dimension>;
- const ConnectivityType& connectivity =
- MeshDataBaseForTests::get().hybrid3DMesh()->get<Mesh<Connectivity<3>>>()->connectivity();
+ const ConnectivityType& connectivity = MeshDataBaseForTests::get().hybrid3DMesh()->get<Mesh<3>>()->connectivity();
SECTION("synchonize NodeValue")
{
@@ -370,7 +368,7 @@ TEST_CASE("Synchronizer", "[mesh]")
using ConnectivityType = Connectivity<Dimension>;
const ConnectivityType& connectivity =
- MeshDataBaseForTests::get().unordered1DMesh()->get<Mesh<Connectivity<1>>>()->connectivity();
+ MeshDataBaseForTests::get().unordered1DMesh()->get<Mesh<1>>()->connectivity();
SECTION("synchonize NodeArray")
{
@@ -502,8 +500,7 @@ TEST_CASE("Synchronizer", "[mesh]")
constexpr size_t Dimension = 2;
using ConnectivityType = Connectivity<Dimension>;
- const ConnectivityType& connectivity =
- MeshDataBaseForTests::get().hybrid2DMesh()->get<Mesh<Connectivity<2>>>()->connectivity();
+ const ConnectivityType& connectivity = MeshDataBaseForTests::get().hybrid2DMesh()->get<Mesh<2>>()->connectivity();
SECTION("synchonize NodeArray")
{
@@ -635,8 +632,7 @@ TEST_CASE("Synchronizer", "[mesh]")
constexpr size_t Dimension = 3;
using ConnectivityType = Connectivity<Dimension>;
- const ConnectivityType& connectivity =
- MeshDataBaseForTests::get().hybrid3DMesh()->get<Mesh<Connectivity<3>>>()->connectivity();
+ const ConnectivityType& connectivity = MeshDataBaseForTests::get().hybrid3DMesh()->get<Mesh<3>>()->connectivity();
SECTION("synchonize NodeArray")
{
@@ -796,7 +792,7 @@ TEST_CASE("Synchronizer", "[mesh]")
using ConnectivityType = Connectivity<Dimension>;
const ConnectivityType& connectivity =
- MeshDataBaseForTests::get().unordered1DMesh()->get<Mesh<Connectivity<1>>>()->connectivity();
+ MeshDataBaseForTests::get().unordered1DMesh()->get<Mesh<1>>()->connectivity();
SECTION("synchonize NodeValuePerCell")
{
@@ -1071,8 +1067,7 @@ TEST_CASE("Synchronizer", "[mesh]")
constexpr size_t Dimension = 2;
using ConnectivityType = Connectivity<Dimension>;
- const ConnectivityType& connectivity =
- MeshDataBaseForTests::get().hybrid2DMesh()->get<Mesh<Connectivity<2>>>()->connectivity();
+ const ConnectivityType& connectivity = MeshDataBaseForTests::get().hybrid2DMesh()->get<Mesh<2>>()->connectivity();
SECTION("synchonize NodeValuePerCell")
{
@@ -1525,8 +1520,7 @@ TEST_CASE("Synchronizer", "[mesh]")
constexpr size_t Dimension = 3;
using ConnectivityType = Connectivity<Dimension>;
- const ConnectivityType& connectivity =
- MeshDataBaseForTests::get().hybrid3DMesh()->get<Mesh<Connectivity<3>>>()->connectivity();
+ const ConnectivityType& connectivity = MeshDataBaseForTests::get().hybrid3DMesh()->get<Mesh<3>>()->connectivity();
SECTION("synchonize NodeValuePerCell")
{
@@ -2103,7 +2097,7 @@ TEST_CASE("Synchronizer", "[mesh]")
using ConnectivityType = Connectivity<Dimension>;
const ConnectivityType& connectivity =
- MeshDataBaseForTests::get().unordered1DMesh()->get<Mesh<Connectivity<1>>>()->connectivity();
+ MeshDataBaseForTests::get().unordered1DMesh()->get<Mesh<1>>()->connectivity();
SECTION("synchonize NodeArrayPerCell")
{
@@ -2408,8 +2402,7 @@ TEST_CASE("Synchronizer", "[mesh]")
constexpr size_t Dimension = 2;
using ConnectivityType = Connectivity<Dimension>;
- const ConnectivityType& connectivity =
- MeshDataBaseForTests::get().hybrid2DMesh()->get<Mesh<Connectivity<2>>>()->connectivity();
+ const ConnectivityType& connectivity = MeshDataBaseForTests::get().hybrid2DMesh()->get<Mesh<2>>()->connectivity();
SECTION("synchonize NodeArrayPerCell")
{
@@ -2912,8 +2905,7 @@ TEST_CASE("Synchronizer", "[mesh]")
constexpr size_t Dimension = 3;
using ConnectivityType = Connectivity<Dimension>;
- const ConnectivityType& connectivity =
- MeshDataBaseForTests::get().hybrid3DMesh()->get<Mesh<Connectivity<3>>>()->connectivity();
+ const ConnectivityType& connectivity = MeshDataBaseForTests::get().hybrid3DMesh()->get<Mesh<3>>()->connectivity();
SECTION("synchonize NodeArrayPerCell")
{