Skip to content

Commit

Permalink
Merge pull request #12074 from KratosMultiphysics/core/clean-up-cpp-t…
Browse files Browse the repository at this point in the history 
…ests

[Core] Clean up C++ tests for processes
  • Loading branch information
@loumalouomega
loumalouomega authored Feb 19, 2024
2 parents 042a62d + a011f78 commit 0ddd7c9
Show file tree
Hide file tree
Showing 19 changed files with 5,633 additions and 5,714 deletions.
378 changes: 187 additions & 191 deletions kratos/tests/cpp_tests/processes/test_apply_ray_casting_process.cpp
View file

Large diffs are not rendered by default.

3,674 changes: 1,836 additions & 1,838 deletions kratos/tests/cpp_tests/processes/test_calculate_discontinuous_distance_to_skin_process.cpp
View file

Large diffs are not rendered by default.

1,049 changes: 519 additions & 530 deletions kratos/tests/cpp_tests/processes/test_calculate_distance_to_skin_process.cpp
View file

Large diffs are not rendered by default.

270 changes: 134 additions & 136 deletions ...os/tests/cpp_tests/processes/test_calculate_embedded_nodal_variable_from_skin_process.cpp
View file
Original file line number Diff line number Diff line change
Expand Up @@ -23,144 +23,142 @@
#include "spaces/ublas_space.h"
#include "testing/testing.h"

namespace Kratos {
namespace Testing {

typedef UblasSpace<double, Matrix, Vector> LocalSpaceType;
typedef UblasSpace<double, CompressedMatrix, Vector> SparseSpaceType;
typedef LinearSolver<SparseSpaceType, LocalSpaceType> LinearSolverType;
typedef CalculateEmbeddedNodalVariableFromSkinProcess<double, SparseSpaceType, LocalSpaceType, LinearSolverType> EmbeddedNodalVariableProcessDouble;
typedef CalculateEmbeddedNodalVariableFromSkinProcess<array_1d<double, 3>, SparseSpaceType, LocalSpaceType, LinearSolverType> EmbeddedNodalVariableProcessArray;

KRATOS_TEST_CASE_IN_SUITE(CalculateEmbeddedNodalVariableFromSkinProcessDouble, KratosCoreFastSuite)
namespace Kratos::Testing {

using LocalSpaceType = UblasSpace<double, Matrix, Vector>;
using SparseSpaceType = UblasSpace<double, CompressedMatrix, Vector>;
using LinearSolverType = LinearSolver<SparseSpaceType, LocalSpaceType>;
using EmbeddedNodalVariableProcessDouble = CalculateEmbeddedNodalVariableFromSkinProcess<double, SparseSpaceType, LocalSpaceType, LinearSolverType>;
using EmbeddedNodalVariableProcessArray = CalculateEmbeddedNodalVariableFromSkinProcess<array_1d<double, 3>, SparseSpaceType, LocalSpaceType, LinearSolverType>;

KRATOS_TEST_CASE_IN_SUITE(CalculateEmbeddedNodalVariableFromSkinProcessDouble, KratosCoreFastSuite)
{
// Generate a volume mesh (done with the StructuredMeshGeneratorProcess)
Node::Pointer p_point_1 = Kratos::make_intrusive<Node>(1, 0.00, 0.00, 0.00);
Node::Pointer p_point_2 = Kratos::make_intrusive<Node>(2, 0.00, 1.00, 0.00);
Node::Pointer p_point_3 = Kratos::make_intrusive<Node>(3, 1.00, 1.00, 0.00);
Node::Pointer p_point_4 = Kratos::make_intrusive<Node>(4, 1.00, 0.00, 0.00);

Quadrilateral2D4<Node> geometry(p_point_1, p_point_2, p_point_3, p_point_4);

Parameters mesher_parameters(R"({
"number_of_divisions": 7,
"element_name": "Element2D3N"
})");

Model current_model;
ModelPart &surface_part = current_model.CreateModelPart("Volume");
surface_part.AddNodalSolutionStepVariable(DISTANCE);
surface_part.AddNodalSolutionStepVariable(TEMPERATURE);
StructuredMeshGeneratorProcess(geometry, surface_part, mesher_parameters).Execute();

// Generate the skin
ModelPart &skin_part = current_model.CreateModelPart("Skin");
skin_part.AddNodalSolutionStepVariable(TEMPERATURE);
skin_part.CreateNewNode(901, 0.24, 0.34, 0.0);
skin_part.CreateNewNode(902, 0.76, 0.34, 0.0);
skin_part.CreateNewNode(903, 0.76, 0.66, 0.0);
skin_part.CreateNewNode(904, 0.24, 0.66, 0.0);
Properties::Pointer p_properties(new Properties(0));
skin_part.CreateNewElement("Element2D2N", 901, {{901, 902}}, p_properties);
skin_part.CreateNewElement("Element2D2N", 902, {{902, 903}}, p_properties);
skin_part.CreateNewElement("Element2D2N", 903, {{903, 904}}, p_properties);
skin_part.CreateNewElement("Element2D2N", 904, {{904, 901}}, p_properties);
skin_part.GetNode(901).FastGetSolutionStepValue(TEMPERATURE) = 1.0;
skin_part.GetNode(902).FastGetSolutionStepValue(TEMPERATURE) = 2.0;
skin_part.GetNode(903).FastGetSolutionStepValue(TEMPERATURE) = 3.0;
skin_part.GetNode(904).FastGetSolutionStepValue(TEMPERATURE) = 2.0;

// Compute distance
CalculateDistanceToSkinProcess<2>(surface_part, skin_part).Execute();

// Compute the embedded nodal variable values
Parameters settings(R"(
{
// Generate a volume mesh (done with the StructuredMeshGeneratorProcess)
Node::Pointer p_point_1 = Kratos::make_intrusive<Node>(1, 0.00, 0.00, 0.00);
Node::Pointer p_point_2 = Kratos::make_intrusive<Node>(2, 0.00, 1.00, 0.00);
Node::Pointer p_point_3 = Kratos::make_intrusive<Node>(3, 1.00, 1.00, 0.00);
Node::Pointer p_point_4 = Kratos::make_intrusive<Node>(4, 1.00, 0.00, 0.00);

Quadrilateral2D4<Node> geometry(p_point_1, p_point_2, p_point_3, p_point_4);

Parameters mesher_parameters(R"({
"number_of_divisions": 7,
"element_name": "Element2D3N"
})");

Model current_model;
ModelPart &surface_part = current_model.CreateModelPart("Volume");
surface_part.AddNodalSolutionStepVariable(DISTANCE);
surface_part.AddNodalSolutionStepVariable(TEMPERATURE);
StructuredMeshGeneratorProcess(geometry, surface_part, mesher_parameters).Execute();

// Generate the skin
ModelPart &skin_part = current_model.CreateModelPart("Skin");
skin_part.AddNodalSolutionStepVariable(TEMPERATURE);
skin_part.CreateNewNode(901, 0.24, 0.34, 0.0);
skin_part.CreateNewNode(902, 0.76, 0.34, 0.0);
skin_part.CreateNewNode(903, 0.76, 0.66, 0.0);
skin_part.CreateNewNode(904, 0.24, 0.66, 0.0);
Properties::Pointer p_properties(new Properties(0));
skin_part.CreateNewElement("Element2D2N", 901, {{901, 902}}, p_properties);
skin_part.CreateNewElement("Element2D2N", 902, {{902, 903}}, p_properties);
skin_part.CreateNewElement("Element2D2N", 903, {{903, 904}}, p_properties);
skin_part.CreateNewElement("Element2D2N", 904, {{904, 901}}, p_properties);
skin_part.GetNode(901).FastGetSolutionStepValue(TEMPERATURE) = 1.0;
skin_part.GetNode(902).FastGetSolutionStepValue(TEMPERATURE) = 2.0;
skin_part.GetNode(903).FastGetSolutionStepValue(TEMPERATURE) = 3.0;
skin_part.GetNode(904).FastGetSolutionStepValue(TEMPERATURE) = 2.0;

// Compute distance
CalculateDistanceToSkinProcess<2>(surface_part, skin_part).Execute();

// Compute the embedded nodal variable values
Parameters settings(R"(
{
"base_model_part_name": "Volume",
"skin_model_part_name": "Skin",
"skin_variable_name": "TEMPERATURE",
"embedded_nodal_variable_name": "TEMPERATURE",
"gradient_penalty_coefficient": 1.0e-5
}
)");
EmbeddedNodalVariableProcessDouble emb_nod_var_from_skin_proc(current_model, settings);
emb_nod_var_from_skin_proc.Execute();
emb_nod_var_from_skin_proc.Clear();

// Check values
std::vector<std::size_t> check_nodes_ids = {19, 20, 28, 45, 46, 54};
std::vector<double> expected_values = {0.883925, 1.42073, 1.69546, 2.57927, 3.11607, 3.47553};
for (std::size_t i_node = 0; i_node < check_nodes_ids.size(); ++i_node) {
const auto p_node = surface_part.pGetNode(check_nodes_ids[i_node]);
KRATOS_EXPECT_NEAR(p_node->FastGetSolutionStepValue(TEMPERATURE), expected_values[i_node], 1e-5);
}
"base_model_part_name": "Volume",
"skin_model_part_name": "Skin",
"skin_variable_name": "TEMPERATURE",
"embedded_nodal_variable_name": "TEMPERATURE",
"gradient_penalty_coefficient": 1.0e-5
}

KRATOS_TEST_CASE_IN_SUITE(CalculateEmbeddedNodalVariableFromSkinProcessArray, KratosCoreFastSuite)
)");
EmbeddedNodalVariableProcessDouble emb_nod_var_from_skin_proc(current_model, settings);
emb_nod_var_from_skin_proc.Execute();
emb_nod_var_from_skin_proc.Clear();

// Check values
std::vector<std::size_t> check_nodes_ids = {19, 20, 28, 45, 46, 54};
std::vector<double> expected_values = {0.883925, 1.42073, 1.69546, 2.57927, 3.11607, 3.47553};
for (std::size_t i_node = 0; i_node < check_nodes_ids.size(); ++i_node) {
const auto p_node = surface_part.pGetNode(check_nodes_ids[i_node]);
KRATOS_EXPECT_NEAR(p_node->FastGetSolutionStepValue(TEMPERATURE), expected_values[i_node], 1e-5);
}
}

KRATOS_TEST_CASE_IN_SUITE(CalculateEmbeddedNodalVariableFromSkinProcessArray, KratosCoreFastSuite)
{
// Generate a volume mesh (done with the StructuredMeshGeneratorProcess)
Node::Pointer p_point_1 = Kratos::make_intrusive<Node>(1, 0.00, 0.00, 0.00);
Node::Pointer p_point_2 = Kratos::make_intrusive<Node>(2, 0.00, 1.00, 0.00);
Node::Pointer p_point_3 = Kratos::make_intrusive<Node>(3, 1.00, 1.00, 0.00);
Node::Pointer p_point_4 = Kratos::make_intrusive<Node>(4, 1.00, 0.00, 0.00);

Quadrilateral2D4<Node> geometry(p_point_1, p_point_2, p_point_3, p_point_4);

Parameters mesher_parameters(R"({
"number_of_divisions": 7,
"element_name": "Element2D3N"
})");

Model current_model;
ModelPart &surface_part = current_model.CreateModelPart("Volume");
surface_part.AddNodalSolutionStepVariable(DISTANCE);
surface_part.AddNodalSolutionStepVariable(DISPLACEMENT);
StructuredMeshGeneratorProcess(geometry, surface_part, mesher_parameters).Execute();

// Generate the skin
ModelPart &skin_part = current_model.CreateModelPart("Skin");
skin_part.AddNodalSolutionStepVariable(DISPLACEMENT);
skin_part.CreateNewNode(901, 0.24, 0.34, 0.0);
skin_part.CreateNewNode(902, 0.76, 0.34, 0.0);
skin_part.CreateNewNode(903, 0.76, 0.66, 0.0);
skin_part.CreateNewNode(904, 0.24, 0.66, 0.0);
Properties::Pointer p_properties(new Properties(0));
skin_part.CreateNewElement("Element2D2N", 901, {{901, 902}}, p_properties);
skin_part.CreateNewElement("Element2D2N", 902, {{902, 903}}, p_properties);
skin_part.CreateNewElement("Element2D2N", 903, {{903, 904}}, p_properties);
skin_part.CreateNewElement("Element2D2N", 904, {{904, 901}}, p_properties);
skin_part.GetNode(901).FastGetSolutionStepValue(DISPLACEMENT_X) = 1.0;
skin_part.GetNode(902).FastGetSolutionStepValue(DISPLACEMENT_Y) = 2.0;
skin_part.GetNode(903).FastGetSolutionStepValue(DISPLACEMENT_X) = 3.0;
skin_part.GetNode(904).FastGetSolutionStepValue(DISPLACEMENT_Y) = 2.0;

// Compute distance
CalculateDiscontinuousDistanceToSkinProcess<2>(surface_part, skin_part).Execute();

// Compute the embedded nodal variable values
Parameters settings(R"(
{
// Generate a volume mesh (done with the StructuredMeshGeneratorProcess)
Node::Pointer p_point_1 = Kratos::make_intrusive<Node>(1, 0.00, 0.00, 0.00);
Node::Pointer p_point_2 = Kratos::make_intrusive<Node>(2, 0.00, 1.00, 0.00);
Node::Pointer p_point_3 = Kratos::make_intrusive<Node>(3, 1.00, 1.00, 0.00);
Node::Pointer p_point_4 = Kratos::make_intrusive<Node>(4, 1.00, 0.00, 0.00);

Quadrilateral2D4<Node> geometry(p_point_1, p_point_2, p_point_3, p_point_4);

Parameters mesher_parameters(R"({
"number_of_divisions": 7,
"element_name": "Element2D3N"
})");

Model current_model;
ModelPart &surface_part = current_model.CreateModelPart("Volume");
surface_part.AddNodalSolutionStepVariable(DISTANCE);
surface_part.AddNodalSolutionStepVariable(DISPLACEMENT);
StructuredMeshGeneratorProcess(geometry, surface_part, mesher_parameters).Execute();

// Generate the skin
ModelPart &skin_part = current_model.CreateModelPart("Skin");
skin_part.AddNodalSolutionStepVariable(DISPLACEMENT);
skin_part.CreateNewNode(901, 0.24, 0.34, 0.0);
skin_part.CreateNewNode(902, 0.76, 0.34, 0.0);
skin_part.CreateNewNode(903, 0.76, 0.66, 0.0);
skin_part.CreateNewNode(904, 0.24, 0.66, 0.0);
Properties::Pointer p_properties(new Properties(0));
skin_part.CreateNewElement("Element2D2N", 901, {{901, 902}}, p_properties);
skin_part.CreateNewElement("Element2D2N", 902, {{902, 903}}, p_properties);
skin_part.CreateNewElement("Element2D2N", 903, {{903, 904}}, p_properties);
skin_part.CreateNewElement("Element2D2N", 904, {{904, 901}}, p_properties);
skin_part.GetNode(901).FastGetSolutionStepValue(DISPLACEMENT_X) = 1.0;
skin_part.GetNode(902).FastGetSolutionStepValue(DISPLACEMENT_Y) = 2.0;
skin_part.GetNode(903).FastGetSolutionStepValue(DISPLACEMENT_X) = 3.0;
skin_part.GetNode(904).FastGetSolutionStepValue(DISPLACEMENT_Y) = 2.0;

// Compute distance
CalculateDiscontinuousDistanceToSkinProcess<2>(surface_part, skin_part).Execute();

// Compute the embedded nodal variable values
Parameters settings(R"(
{
"base_model_part_name": "Volume",
"skin_model_part_name": "Skin",
"skin_variable_name": "DISPLACEMENT",
"embedded_nodal_variable_name": "DISPLACEMENT",
"gradient_penalty_coefficient": 1.0e-5
}
)");
EmbeddedNodalVariableProcessArray emb_nod_var_from_skin_proc(current_model, settings);
emb_nod_var_from_skin_proc.Execute();
emb_nod_var_from_skin_proc.Clear();

// Check values
const std::vector<std::size_t> check_nodes_ids = {19, 20, 28, 45, 46, 54};
const std::vector<double> expected_values_x = {0.344868, 1.83755, 1.56283, 2.99609, 2.57702, 1.75274};
const std::vector<double> expected_values_y = {0.539057, -0.416819, 0.132631, -0.416819, 0.539057, 1.72279};
for (std::size_t i_node = 0; i_node < check_nodes_ids.size(); ++i_node) {
const auto p_node = surface_part.pGetNode(check_nodes_ids[i_node]);
KRATOS_EXPECT_NEAR(p_node->FastGetSolutionStepValue(DISPLACEMENT_X), expected_values_x[i_node], 1e-5);
KRATOS_EXPECT_NEAR(p_node->FastGetSolutionStepValue(DISPLACEMENT_Y), expected_values_y[i_node], 1e-5);
}
"base_model_part_name": "Volume",
"skin_model_part_name": "Skin",
"skin_variable_name": "DISPLACEMENT",
"embedded_nodal_variable_name": "DISPLACEMENT",
"gradient_penalty_coefficient": 1.0e-5
}
)");
EmbeddedNodalVariableProcessArray emb_nod_var_from_skin_proc(current_model, settings);
emb_nod_var_from_skin_proc.Execute();
emb_nod_var_from_skin_proc.Clear();

// Check values
const std::vector<std::size_t> check_nodes_ids = {19, 20, 28, 45, 46, 54};
const std::vector<double> expected_values_x = {0.344868, 1.83755, 1.56283, 2.99609, 2.57702, 1.75274};
const std::vector<double> expected_values_y = {0.539057, -0.416819, 0.132631, -0.416819, 0.539057, 1.72279};
for (std::size_t i_node = 0; i_node < check_nodes_ids.size(); ++i_node) {
const auto p_node = surface_part.pGetNode(check_nodes_ids[i_node]);
KRATOS_EXPECT_NEAR(p_node->FastGetSolutionStepValue(DISPLACEMENT_X), expected_values_x[i_node], 1e-5);
KRATOS_EXPECT_NEAR(p_node->FastGetSolutionStepValue(DISPLACEMENT_Y), expected_values_y[i_node], 1e-5);
}
}

} // namespace Testing.
} // namespace Kratos.
} // namespace Kratos::Testing.
Loading

0 comments on commit 0ddd7c9

Please sign in to comment.