main.cc

Go to the documentation of this file.

// Calls the unit tests for the PRISMS-PF code
// Orignal author: Stephen DeWitt (stvdwtt@umich.edu)

#include <iostream>
#include "unitTest.h"
#include "../../include/userInputParameters.h"
#include "../../src/userInputParameters/userInputParameters.cc"
#include "../../include/vectorBCFunction.h"
#include "../../src/utilities/vectorBCFunction.cc"
#include "../../include/initialConditions.h"
#include "initialConditions.cc"
#include "boundaryConditions.cc"
#include "unit_test_inputs.cc"

int main(int argc, char **argv)
{
  dealii::deallog.depth_console(3);
  //init MPI
  Utilities::MPI::MPI_InitFinalize mpi_initialization(argc, argv, numbers::invalid_unsigned_int);

  // Load input
  variableAttributeLoader variable_attributes;
  inputFileReader input_file_reader("parameters_test.in",variable_attributes);
  userInputParameters<2> userInputs(input_file_reader,input_file_reader.parameter_handler,variable_attributes);
  load_unit_test_inputs<2>(userInputs);

  std::cout << "Beginning unit tests..." << std::endl;
  bool pass = false;
  int tests_passed = 0, total_tests = 0;

  // Unit tests for the method "parse_line" in "inputFileReader"
  total_tests++;
  unitTest<2,double> parse_line_tester;
  pass = parse_line_tester.test_parse_line();
  tests_passed += pass;

  // Unit tests for the method "get_subsection_entry_list" in "inputFileReader"
  total_tests++;
  unitTest<2,double> get_subsection_entry_list_tester;
  pass = get_subsection_entry_list_tester.test_get_subsection_entry_list();
  tests_passed += pass;

  // Unit tests for the method "get_entry_name_ending_list" in "inputFileReader"
  total_tests++;
  unitTest<2,double> get_entry_name_ending_list_tester;
  pass = get_entry_name_ending_list_tester.test_get_entry_name_ending_list();
  tests_passed += pass;

  // Unit tests for the method "load_BC_list" in "userInputParameters"
  total_tests++;
  unitTest<2,double> load_BC_list_tester;
  pass = load_BC_list_tester.test_load_BC_list();
  tests_passed += pass;

  // Unit tests for the method "setOutputTimeSteps" for all four types of spacing, indirectly through "loadInputParameters"
  total_tests++;
  unitTest<2,double> setOutputTimeSteps_tester_eq;
  pass = setOutputTimeSteps_tester_eq.test_setOutputTimeSteps();
  tests_passed += pass;

  // Unit tests for the method "computeInvM"
  total_tests++;
  unitTest<2,double> computeInvM_tester_2D;
  pass = computeInvM_tester_2D.test_computeInvM(argc, argv, userInputs);
  tests_passed += pass;

  // Unit tests for the method "outputResults" (this test doesn't actually do anything)
  // total_tests++;
  // unitTest<2,double> outputResults_tester_2D;
  // pass = outputResults_tester_2D.test_outputResults(argc, argv, userInputs);
  // tests_passed += pass;

  // Unit tests for the method "computeStress"
  total_tests++;
  unitTest<1,dealii::VectorizedArray<double>[1][1]> computeStress_tester_1D;
  pass = computeStress_tester_1D.test_computeStress();
  tests_passed += pass;

  total_tests++;
  unitTest<1,dealii::Table<2, double> > computeStress_tester_1DT;
  pass = computeStress_tester_1DT.test_computeStress();
  tests_passed += pass;

  total_tests++;
  unitTest<2,dealii::VectorizedArray<double>[3][3]> computeStress_tester_2D;
  pass = computeStress_tester_2D.test_computeStress();
  tests_passed += pass;

  total_tests++;
  unitTest<2,dealii::Table<2, double> > computeStress_tester_2DT;
  pass = computeStress_tester_2DT.test_computeStress();
  tests_passed += pass;

  total_tests++;
  unitTest<3,dealii::VectorizedArray<double>[6][6]> computeStress_tester_3D;
  pass = computeStress_tester_3D.test_computeStress();
  tests_passed += pass;

  total_tests++;
  unitTest<3,dealii::Table<2, double> > computeStress_tester_3DT;
  pass = computeStress_tester_3DT.test_computeStress();
  tests_passed += pass;

  // Unit tests for the method "setRigidBodyModeConstraints"
  total_tests++;
  unitTest<2,double> setRigidBodyModeConstraints_tester_null;
  std::vector<int> rigidBodyModeComponents;
  pass = setRigidBodyModeConstraints_tester_null.test_setRigidBodyModeConstraints(rigidBodyModeComponents, userInputs);
  tests_passed += pass;

  total_tests++;
  unitTest<2,double> setRigidBodyModeConstraints_tester_one;
  rigidBodyModeComponents.clear();
  rigidBodyModeComponents.push_back(0);
  pass = setRigidBodyModeConstraints_tester_one.test_setRigidBodyModeConstraints(rigidBodyModeComponents, userInputs);
  tests_passed += pass;

  // In debug mode this test has an exception because it is trying to access components higher than 0 on the mesh for a scalar.
  // To get this test working again, I'll need a DoFHandler for a vector field.
  // total_tests++;
  // unitTest<2,double> setRigidBodyModeConstraints_tester_three;
  // rigidBodyModeComponents.clear();
  // rigidBodyModeComponents.push_back(0);
  // rigidBodyModeComponents.push_back(1);
  // rigidBodyModeComponents.push_back(2);
  // pass = setRigidBodyModeConstraints_tester_three.test_setRigidBodyModeConstraints(rigidBodyModeComponents, userInputs);
  // tests_passed += pass;

  // Unit tests for the "LinearSolverParameters" class
  total_tests++;
  unitTest<2,double> LinearSolverParameters_tester;
  pass = LinearSolverParameters_tester.test_LinearSolverParameters();
  tests_passed += pass;

  // Unit tests for the "NonlinearSolverParameters" class
  total_tests++;
  unitTest<2,double> NonlinearSolverParameters_tester;
  pass = NonlinearSolverParameters_tester.test_NonlinearSolverParameters();
  tests_passed += pass;

  // Unit tests for the "parseDependences" method in the "UserInputParameters" class
  total_tests++;
  unitTest<2,double> EquationDependencyParser_tester1;
  pass = EquationDependencyParser_tester1.test_EquationDependencyParser_variables_and_residuals_needed();
  tests_passed += pass;

  // Unit tests for the "parseDependences" method in the "UserInputParameters" class
  total_tests++;
  unitTest<2,double> EquationDependencyParser_tester2;
  pass = EquationDependencyParser_tester2.test_EquationDependencyParser_nonlinear();
  tests_passed += pass;

  // Unit tests for the "parseDependences" method in the "UserInputParameters" class
  total_tests++;
  unitTest<2,double> EquationDependencyParser_tester3;
  pass = EquationDependencyParser_tester3.test_EquationDependencyParser_postprocessing();
  tests_passed += pass;

  // Unit tests for the FloodFiller class
  total_tests++;
  unitTest<2,double> FloodFiller_tester;
  pass = FloodFiller_tester.test_FloodFiller();
  tests_passed += pass;

  // Unit tests for the SimplifiedGrainRepresentation class
  total_tests++;
  unitTest<2,double> SimplifiedGrainRepresentation_tester;
  pass = SimplifiedGrainRepresentation_tester.test_SimplifiedGrainRepresentation();
  tests_passed += pass;

  // Unit tests for the "transferGrainIds" method in the "SimplifiedGrainManipulator" class
  total_tests++;
  unitTest<2,double> SimplifiedGrainManipulator_tester1;
  pass = SimplifiedGrainManipulator_tester1.test_SimplifiedGrainManipulator_transferGrainIds();
  tests_passed += pass;

  // Unit tests for the "reassignGrains" method in the "SimplifiedGrainManipulator" class
  total_tests++;
  unitTest<2,double> SimplifiedGrainManipulator_tester2;
  pass = SimplifiedGrainManipulator_tester2.test_SimplifiedGrainManipulator_reassignGrains();
  tests_passed += pass;

  // Unit tests for the OrderParameterRemapper class
  total_tests++;
  unitTest<2,double> OrderParameterRemapper_tester;
  pass = OrderParameterRemapper_tester.test_OrderParameterRemapper();
  tests_passed += pass;

  // Print out results
  char buffer[100];
  sprintf(buffer, "\n\nNumber of tests passed: %u/%u \n\n", tests_passed, total_tests);
  std::cout << buffer;

  // Write results to a file
  std::ofstream result_file ("unit_test_results.txt");
  result_file << tests_passed << std::endl;
  result_file << total_tests << std::endl;
  result_file.close();


  return 0;
}