TEST_F(NumLibMeshComponentMapTest, DISABLED_CheckOrderByComponent)
#endif
{
// - Entries in the vector are arranged in the order of a component type and then node ID
// - For example, x=[(node 0, comp 0) (node 1, comp 0) ... (node n, comp0), (node 0, comp1) ... ]
cmap = new MeshComponentMap(components,
NumLib::ComponentOrder::BY_COMPONENT);
ASSERT_EQ(2 * mesh->getNumberOfNodes(), cmap->dofSizeWithGhosts());
for (std::size_t i = 0; i < mesh_size; i++)
{
// Test global indices for the different components of the node.
ASSERT_EQ(i , giAtNodeForComponent(i, comp0_id));
ASSERT_EQ(mesh_size + 1 + i, giAtNodeForComponent(i, comp1_id));
// Test component ids of the node.
std::vector<std::size_t> const vecCompIDs = cmap->getComponentIDs(
Location(mesh->getID(), MeshItemType::Node, i));
ASSERT_EQ(2u, vecCompIDs.size());
ASSERT_EQ(0u, vecCompIDs[0]);
ASSERT_EQ(1u, vecCompIDs[1]);
}
}
TEST_F(AssemblerLibMeshComponentMapTest, DISABLED_CheckOrderByLocation)
#endif
{
// - Entries in the vector are arranged in the order of node ID and then a component type
// - For example, x=[(node 0, comp 0) (node 0, comp 1) ... (node n, comp0), (node n, comp1) ]
cmap = new MeshComponentMap(components,
AssemblerLib::ComponentOrder::BY_LOCATION);
ASSERT_EQ(2 * mesh->getNNodes(), cmap->size());
for (std::size_t i = 0; i < mesh_size; i++)
{
// Test global indices for the different components of the node.
ASSERT_EQ(2 * i , giAtNodeForComponent(i, comp0_id));
ASSERT_EQ(2 * i + 1, giAtNodeForComponent(i, comp1_id));
// Test component ids of the node.
std::vector<std::size_t> const vecCompIDs = cmap->getComponentIDs(
Location(mesh->getID(), MeshItemType::Node, i));
ASSERT_EQ(2u, vecCompIDs.size());
ASSERT_EQ(0u, vecCompIDs[0]);
ASSERT_EQ(1u, vecCompIDs[1]);
}
}
