https://gitlab.opengeosys.org/ogs/ogs.git
Tip revision: 6d59cc40a45c71576076de5921ee4a6e06c49bc1 authored by Tobias Meisel on 06 April 2021, 12:13:45 UTC
[MeL/IO] Add missing xdmf parameter documentation.
[MeL/IO] Add missing xdmf parameter documentation.
Tip revision: 6d59cc4
TestDuplicateGeometry.cpp
/**
* \copyright
* Copyright (c) 2012-2021, OpenGeoSys Community (http://www.opengeosys.org)
* Distributed under a Modified BSD License.
* See accompanying file LICENSE.txt or
* http://www.opengeosys.org/project/license
*
*
*/
#include <gtest/gtest.h>
#include <cstdlib>
#include <ctime>
#include <random>
#include <string>
#include "GeoLib/DuplicateGeometry.h"
#include "GeoLib/GEOObjects.h"
#include "GeoLib/Point.h"
#include "GeoLib/Polyline.h"
#include "GeoLib/Surface.h"
#include "GeoLib/Triangle.h"
TEST(GeoLib, DuplicateGeometry)
{
GeoLib::GEOObjects geo;
std::string input_name("input");
// generate points
std::srand(static_cast<unsigned>(std::time(nullptr)));
std::size_t n_pnts(rand() % 1000 + 100);
int box_size(std::rand());
double half_box_size(box_size / 2);
auto pnts = std::make_unique<std::vector<GeoLib::Point*>>();
pnts->reserve(n_pnts);
for (int k(0); k < static_cast<int>(n_pnts); k++)
{
pnts->push_back(
new GeoLib::Point(std::rand() % box_size - half_box_size,
std::rand() % box_size - half_box_size,
std::rand() % box_size - half_box_size));
}
geo.addPointVec(std::move(pnts), input_name);
// update number of points in case because possibly non-unique points have
// been removed
n_pnts = geo.getPointVec(input_name)->size();
std::string output("output_geometry");
// duplicate points
{
GeoLib::DuplicateGeometry dup(geo, input_name, output);
std::vector<GeoLib::Point*> const* const pnts(
geo.getPointVec(input_name));
std::vector<GeoLib::Point*> const* const new_pnts(
geo.getPointVec(output));
std::vector<GeoLib::Point*>& mod_pnts(dup.getPointVectorCopy());
ASSERT_EQ(n_pnts, new_pnts->size());
ASSERT_EQ(n_pnts, mod_pnts.size());
for (std::size_t i = 0; i < n_pnts; ++i)
{
ASSERT_EQ((*(*pnts)[i])[0], (*(*new_pnts)[i])[0]);
ASSERT_EQ((*(*pnts)[i])[1], (*(*new_pnts)[i])[1]);
ASSERT_EQ((*(*pnts)[i])[2], (*(*new_pnts)[i])[2]);
}
mod_pnts.push_back(new GeoLib::Point(0, 0, 0));
ASSERT_EQ(mod_pnts.size(), pnts->size() + 1);
ASSERT_EQ(mod_pnts.size(), new_pnts->size());
}
// add polylines
{
int const n_plys = std::rand() % 10 + 1;
auto plys = std::make_unique<std::vector<GeoLib::Polyline*>>();
for (int i = 0; i < n_plys; ++i)
{
int const n_ply_pnts = std::rand() % 100 + 2;
auto* line = new GeoLib::Polyline(*geo.getPointVec(input_name));
for (int j = 0; j < n_ply_pnts; ++j)
{
line->addPoint(std::rand() % n_pnts);
}
plys->push_back(line);
}
geo.addPolylineVec(std::move(plys), input_name);
}
// duplicate polylines
{
GeoLib::DuplicateGeometry dup(geo, input_name, output);
std::string const& output2 = dup.getFinalizedOutputName();
ASSERT_FALSE(output == output2);
std::vector<GeoLib::Polyline*> const* const plys(
geo.getPolylineVec(input_name));
std::vector<GeoLib::Polyline*> const* const new_plys(
geo.getPolylineVec(output2));
ASSERT_EQ(plys->size(), new_plys->size());
for (std::size_t i = 0; i < plys->size(); ++i)
{
std::size_t const n_ply_pnts((*new_plys)[i]->getNumberOfPoints());
ASSERT_EQ(n_ply_pnts, (*plys)[i]->getNumberOfPoints());
ASSERT_EQ((*new_plys)[i]->getNumberOfSegments(),
(*plys)[i]->getNumberOfSegments());
for (std::size_t j = 0; j < n_ply_pnts; ++j)
{
ASSERT_EQ((*plys)[i]->getPointID(j),
(*new_plys)[i]->getPointID(j));
}
}
std::vector<GeoLib::Polyline*>& mod_plys(dup.getPolylineVectorCopy());
mod_plys.push_back(new GeoLib::Polyline(*(*new_plys)[0]));
ASSERT_EQ(mod_plys.size(), plys->size() + 1);
ASSERT_EQ(mod_plys.size(), new_plys->size());
}
// add surfaces
{
int const n_sfcs = std::rand() % 10 + 1;
auto sfcs = std::make_unique<std::vector<GeoLib::Surface*>>();
for (int i = 0; i < n_sfcs; ++i)
{
int const n_tris = std::rand() % 10 + 1;
auto* sfc = new GeoLib::Surface(*geo.getPointVec(input_name));
while (sfc->getNumberOfTriangles() <= n_tris)
{
sfc->addTriangle(std::rand() % n_pnts,
std::rand() % n_pnts,
std::rand() % n_pnts);
}
ASSERT_GT(sfc->getNumberOfTriangles(), 0);
sfcs->push_back(sfc);
}
geo.addSurfaceVec(std::move(sfcs), input_name);
}
// duplicate surfaces
{
GeoLib::DuplicateGeometry dup(geo, input_name, output);
std::string const& output2 = dup.getFinalizedOutputName();
ASSERT_FALSE(output == output2);
std::vector<GeoLib::Surface*> const* sfcs(
geo.getSurfaceVec(input_name));
std::vector<GeoLib::Surface*> const* new_sfcs(
geo.getSurfaceVec(output2));
ASSERT_EQ(sfcs->size(), new_sfcs->size());
for (std::size_t i = 0; i < sfcs->size(); ++i)
{
std::size_t const n_tris((*new_sfcs)[i]->getNumberOfTriangles());
ASSERT_EQ(n_tris, (*sfcs)[i]->getNumberOfTriangles());
for (std::size_t j = 0; j < n_tris; ++j)
{
for (std::size_t k = 0; k < 3; ++k)
{
ASSERT_EQ((*(*(*sfcs)[i])[j])[k],
(*(*(*new_sfcs)[i])[j])[k]);
}
}
}
std::vector<GeoLib::Point*>& mod_pnts(dup.getPointVectorCopy());
std::vector<GeoLib::Surface*>& mod_sfcs(dup.getSurfaceVectorCopy());
std::size_t n_pnts = mod_pnts.size();
mod_pnts.push_back(new GeoLib::Point(1, 0, 0, n_pnts));
mod_pnts.push_back(new GeoLib::Point(0, 1, 0, n_pnts + 1));
mod_pnts.push_back(new GeoLib::Point(0, 0, 1, n_pnts + 2));
auto* sfc = new GeoLib::Surface(mod_pnts);
sfc->addTriangle(n_pnts, n_pnts + 1, n_pnts + 2);
mod_sfcs.push_back(sfc);
ASSERT_EQ(mod_sfcs.size(), sfcs->size() + 1);
ASSERT_EQ(mod_sfcs.size(), new_sfcs->size());
}
}