scaling_method_with_program_args.cpp File Reference

#include <MoFEM.hpp>

Go to the source code of this file.

Functions
int	main (int argc, char *argv[])

Variables
static char	help [] = "...\n\n"

Function Documentation

main()

int main	(	int	argc,
		char *	argv[]
	)

Definition at line 13 of file scaling_method_with_program_args.cpp.

                                 {
 
  MoFEM::Core::Initialize(&argc, &argv, (char *)0, help);
  std::vector<double> scalar_values = {1.1, 2.4, 3.6,  4.1,  3.1,
                                       5.1, 9.1, 10.5, 11.2, 15.3};
  try {
    auto time_scale = std::make_shared<TimeScale>();
    for (int i = 1; i <= scalar_values.size(); i++) {
      if (std::fabs(time_scale->getScale(double(i)) - scalar_values[i - 1]) >
          std::numeric_limits<double>::epsilon()) {
        SETERRQ2(PETSC_COMM_SELF, MOFEM_ATOM_TEST_INVALID,
                 "Validation for data scaling from csv "
                 "failed for time: %d value: %d",
                 double(i), time_scale->getScale(i));
      }
    }
    double time1 = 3.0;
    double time0 = 2.0;
    double scale1 = scalar_values[2];
    double scale0 = scalar_values[1];
    double input_time = 2.5;
    double interp_t = (input_time - time0) / (time1 - time0);
    double expected_scale = scale0 + (scale1 - scale0) * interp_t;
    double actual_scale = time_scale->getScale(2.5);
    if (std::fabs(expected_scale - actual_scale) >
        std::numeric_limits<double>::epsilon()) {
      SETERRQ2(
          PETSC_COMM_SELF, MOFEM_ATOM_TEST_INVALID,
          "Validation for data scaling from csv failed for time: %f value: %f",
          2.5, time_scale->getScale(2.5));
    }
  }
  CATCH_ERRORS;
 
  MoFEM::Core::Finalize();
 
  return 0;
}

Variable Documentation

help

char help[] = "...\n\n"

static

Definition at line 11 of file scaling_method_with_program_args.cpp.