tests_lukas.cpp

Go to the documentation of this file.

#include "problem.h"
#include "gelib.h"

#include "meso2d.h"
#include "meso3d.h"
#include "mesoface.h"

// #include "lukas\synthesis.h"

using namespace mumech;


double deg2rad (double deg) { return deg * (PI / 180.0); }
void rotate_point_about_X_axis (const double alpha, meso3d::point_rec & point, rotation r);
void rotate_point_about_Y_axis (const double alpha, meso3d::point_rec & point, rotation r);
void rotate_point_about_Z_axis (const double alpha, meso3d::point_rec & point, rotation r);
void rotate_point_by_euler_angles_in_313_notation (const double * eAngles, meso3d::point_rec & point);
void test_point_rotation ();

void fill_ellipsoid_rec (Inclusion* inc, meso3d::ellipsoid_rec & L);
void check_ellipsoid_rec_consistency (meso3d::ellipsoid_rec & L);
void transfom_ellipsoid_rec (meso3d::ellipsoid_rec & L);
bool ellipsoids_overlap (Inclusion* inc1, Inclusion* inc2);
void test_prolinani_elipsoidu ();

void fill_ellipse_rec (Inclusion* inc, meso2d::ellipse_rec & L);
void check_ellipse_rec_consistency (meso2d::ellipse_rec & L);
void transfom_ellipse_rec (meso2d::ellipse_rec & L);
bool ellipses_overlap (Inclusion* inc1, Inclusion* inc2);
void test_prolinani_elips ();

//
void ovlivneni(void);

// XXXXXXXXXX XXXXXXXXXX XXXXXXXXXX XXXXXXXXXX XXXXXXXXXX XXXXXXXXXX XXXXXXXXXX

void tests_lukas (void)
{
  fprintf (stdout, "lukasovo test - start\n\n");

  // test_point_rotation ();
  // test_prolinani_elipsoidu ();
  //test_prolinani_elips ();

  //if (1) {
  //    Synthesis test ("debugInput.txt");
  //}

  //ovlivneni();

  fprintf (stdout, "\nlukasovo test - end\n");
  return;
}

void test_point_rotation ()
{
  meso3d::point_rec bod;

  double uhly[3];
    
  uhly[0] = deg2rad (30.0); uhly[1] = deg2rad (45.0); uhly[2] = deg2rad (60.0);
  bod.x = 3.0; bod.y = 2.0; bod.z = 1.0;
  //decasne zakom.
  //mesoface::rotate_point_by_euler_angles_in_313_notation (uhly, bod);
  printf ("%.4lf %.4lf %.4lf\n", bod.x, bod.y, bod.z);
  printf ("-1.1195 1.4750 3.2513\n");

  printf ("\n");

  uhly[0] = deg2rad (10.0); uhly[1] = deg2rad (15.0); uhly[2] = deg2rad (130.0);
  bod.x = 1.0; bod.y = 2.0; bod.z = 3.0;
  //decasne zakom.
  MesoFace::rotate_point_by_euler_angles_in_313_notation (uhly, bod);
  printf ("%.4lf %.4lf %.4lf\n", bod.x, bod.y, bod.z);
  printf ("-1.9199 -1.6365 2.7633\n");

  return;
}








/*end of file*/