TEST(FloatBoxTest, PositionTest)
{
FloatBox box(0, 0, 0, 4, 4, 4);
box.move(FloatPoint3D(1, 2, 3));
EXPECT_PRED_FORMAT2(FloatBoxTest::AssertAlmostEqual, FloatBox(1, 2, 3, 4, 4, 4), box);
box.setOrigin(FloatPoint3D(-1, -2, -3));
box.move(FloatPoint3D(-1, -2, -3));
EXPECT_PRED_FORMAT2(FloatBoxTest::AssertAlmostEqual, FloatBox(-2, -4, -6, 4, 4, 4), box);
}
TEST(FloatBoxTest, UnionTest)
{
FloatBox box;
EXPECT_TRUE(box.isEmpty());
FloatBox unionedBox(3, 5, 6, 5, 3, 9);
box.unionBounds(unionedBox);
EXPECT_PRED_FORMAT2(FloatBoxTest::AssertAlmostEqual, unionedBox, box);
box.unionBounds(FloatBox());
EXPECT_PRED_FORMAT2(FloatBoxTest::AssertAlmostEqual, unionedBox, box);
box.unionBounds(FloatBox(0, 0, 0, 1, 1, 1));
EXPECT_PRED_FORMAT2(FloatBoxTest::AssertAlmostEqual, FloatBox(0, 0, 0, 8, 8, 15), box);
}
inline void createBoundingBoxesForElementsInElementBlocks(const stk::mesh::BulkData &bulk, FlaotBoxVector& domainBoxes)
{
stk::mesh::EntityVector elements;
stk::mesh::get_selected_entities(bulk.mesh_meta_data().locally_owned_part(), bulk.buckets(stk::topology::ELEM_RANK), elements);
size_t numberBoundingBoxes = elements.size();
domainBoxes.resize(numberBoundingBoxes);
stk::mesh::FieldBase const * coords = bulk.mesh_meta_data().coordinate_field();
std::vector<double> boxCoordinates(6);
for(size_t i=0;i<elements.size();++i)
{
unsigned num_nodes = bulk.num_nodes(elements[i]);
std::vector<double> coordinates(3*num_nodes,0);
const stk::mesh::Entity* nodes = bulk.begin_nodes(elements[i]);
for(unsigned j=0;j<num_nodes;++j)
{
double* data = static_cast<double*>(stk::mesh::field_data(*coords, nodes[j]));
coordinates[3*j] = data[0];
coordinates[3*j+1] = data[1];
coordinates[3*j+2] = data[2];
}
findBoundingBoxCoordinates(coordinates, boxCoordinates);
unsigned id = bulk.identifier(elements[i]);
Ident domainBoxId(id, bulk.parallel_rank());
domainBoxes[i] = std::make_pair(FloatBox(boxCoordinates[0], boxCoordinates[1], boxCoordinates[2],
boxCoordinates[3], boxCoordinates[4], boxCoordinates[5]),
domainBoxId);
}
}
TEST(FloatBoxTest, CopyTest)
{
FloatBox box(1, 2, 3, 4, 4, 4);
FloatBox box2(box);
EXPECT_EQ(box, box2);
box.setSize(FloatPoint3D(3, 3, 3));
EXPECT_NE(box, box2);
EXPECT_PRED_FORMAT2(FloatBoxTest::AssertAlmostEqual, FloatBox(1, 2, 3, 3, 3, 3), box);
}
TEST(FloatBoxTest, ExpandTests)
{
FloatBox box;
box.expandTo(FloatPoint3D(10, -3, 2));
EXPECT_PRED_FORMAT2(FloatBoxTest::AssertAlmostEqual, FloatBox(0, -3, 0, 10, 3, 2), box);
box.expandTo(FloatPoint3D(-15, 6, 8));
EXPECT_PRED_FORMAT2(FloatBoxTest::AssertAlmostEqual, FloatBox(-15, -3, 0, 25, 9, 8), box);
box = FloatBox();
box.expandTo(FloatPoint3D(-3, 6, 9), FloatPoint3D(-2, 10, 11));
EXPECT_PRED_FORMAT2(FloatBoxTest::AssertAlmostEqual, FloatBox(-3, 0, 0, 3, 10, 11), box);
box = FloatBox();
box.expandTo(FloatBox(-10, -10, -10, 3, 30, 40));
box.expandTo(FloatBox(-11, 3, 50, 10, 15, 1));
EXPECT_PRED_FORMAT2(FloatBoxTest::AssertAlmostEqual, FloatBox(-11, -10, -10, 11, 30, 61), box);
}
TEST(FloatBoxTest, FlattenTest)
{
FloatBox box(1, 2, 3, 4, 4, 4);
box.flatten();
EXPECT_PRED_FORMAT2(FloatBoxTest::AssertAlmostEqual, FloatBox(1, 2, 0, 4, 4, 0), box);
}
