matdb-suite.c

Go to the documentation of this file.

/*
** matdb-suite.c
** 
** Made by (Johnny Q. Hacker)
** Login   <solarion@borkborkbork>
** 

    Copyright (C) 2009 Joseph Pingenot

    This program is free software: you can redistribute it and/or modify
    it under the terms of the GNU Affero General Public License as published by
    the Free Software Foundation, either version 3 of the License, or
    (at your option) any later version.

    This program is distributed in the hope that it will be useful,
    but WITHOUT ANY WARRANTY; without even the implied warranty of
    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
    GNU Affero General Public License for more details.

    You should have received a copy of the GNU Affero General Public License
    along with this program.  If not, see <http://www.gnu.org/licenses/>.

** Started on  Thu Mar 19 12:52:41 2009 Johnny Q. Hacker
** Last update Sun May 12 01:17:25 2002 Speed Blue
*/

#include <math.h>
#include <stdio.h>
#include <CUnit/CUnit.h>
#include "argdefaults.h"
#include "matdb.h"
#include "tests/matdb-suite.h"

struct matdb *mdb;

void create_matdb(void) {
  mdb = matdb_new();
  CU_ASSERT(mdb != NULL);
  CU_ASSERT(mdb->materials != NULL);
  CU_ASSERT(mdb->bowings != NULL);
}

void destroy_matdb(void) {
  matdb_destroy(mdb);
}

/*Test to set and then get material 1 parameter 1*/
void set_get_material_1_1(void) {
  //fprintf(stderr, "set_get_material_1_1:");
  double d = 0.5;
  int err;
  set_material_param(mdb, "material_1", "param_1", d, &err);
  CU_ASSERT(err == 0);
  double e = get_material_param(mdb, "material_1", "param_1", &err);
  //fprintf(stderr, "d=%g e=%g err=%d\n", d, e, err);
  CU_ASSERT(d == e);
  CU_ASSERT(err == 0);
}

/*Test to set and then get material 1 parameter 2*/
void set_get_material_1_2(void) {
  //fprintf(stderr, "set_get_material_1_2\n");
  double d = 73;
  int err;
  set_material_param(mdb, "material_1", "param_2", d, &err);
  CU_ASSERT(err == 0);
  double e = get_material_param(mdb, "material_1", "param_2", &err);
  CU_ASSERT(d == e);
  CU_ASSERT(err == 0);
}

/*Test to set and then get material 2 parameter 1*/
void set_get_material_2_1(void) {
  //fprintf(stderr, "set_get_material_2_2\n");
  double d = 1.3;
  int err;
  set_material_param(mdb, "material_2", "param_1", d, &err);
  CU_ASSERT(err == 0);
  double e = get_material_param(mdb, "material_2", "param_1", &err);
  CU_ASSERT(d == e);
  CU_ASSERT(err == 0);
}
/*Test to set and then get material 2 parameter 2*/
void set_get_material_2_2(void) {
  //fprintf(stderr, "set_get_material_2_2\n");
  double d = 0.17;
  int err;
  set_material_param(mdb, "material_2", "param_2", d, &err);
  CU_ASSERT(err == 0);
  double e = get_material_param(mdb, "material_2", "param_2", &err);
  CU_ASSERT(d == e);
  CU_ASSERT(err == 0);
}

/*Test to set and then get alloy of 1 and 2, with bowing
 *Also checks the other way (alloy of 2 and 1 with same param)
 */
void set_get_alloy_1_2_1_b(void) {
  //fprintf(stderr, "set_get_alloy_1_2_1_b: ");
  double b = 7;
  /*25% mat1 and 75% mat2*/
  double x=0.25;
  int err;
  set_bowing_param(mdb, "material_1", "material_2", "param_1", b, &err);
  //fprintf(stderr, "b=%g x=%g err=%d ", b, x, err);
  CU_ASSERT(err == 0);
  double d1 = get_material_param(mdb, "material_1", "param_1", &err);
  CU_ASSERT(err == 0);
  double d2 = get_material_param(mdb, "material_2", "param_1", &err);
  CU_ASSERT(err == 0);
  double v12 = get_alloy_param(mdb, "material_1", "material_2", x, "param_1", &err);
  CU_ASSERT(err == 0);
  double v21 = get_alloy_param(mdb, "material_2", "material_1", (1.0-x), "param_1", &err);
  CU_ASSERT(err == 0);
  CU_ASSERT(fabs(v12 - v21) <= fabs(CLOSE_ENOUGH*(v12)));
  double v_12 = alloy_bowed(d1, d2, b, x);
  double v_21 = alloy_bowed(d2, d1, b, (1-x));
  CU_ASSERT(fabs(v_12 - v_21) <= fabs(CLOSE_ENOUGH*(v_12)));
  CU_ASSERT(fabs(v_12 - v21) <= fabs(CLOSE_ENOUGH*(v_12)));
  CU_ASSERT(fabs(v12 - v_21) <= fabs(CLOSE_ENOUGH*(v_12)));
  CU_ASSERT(fabs(v21 - v_12) <= fabs(CLOSE_ENOUGH*(v_12)));
  CU_ASSERT(fabs(v_21 - v12) <= fabs(CLOSE_ENOUGH*(v_12)));
  //fprintf(stderr, "v12=%g v21=%g v_12=%g v_21=%g\n", v12, v21, v_12, v_21);
}

/*Test to set and then get alloy of 1 and 2, with bowing
 *Also checks the other way (alloy of 2 and 1 with same param)
 */
void set_get_alloy_1_2_2_l(void) {
  //fprintf(stderr, "set_get_alloy_1_2_2_b\n");
  /*25% mat1 and 75% mat2*/
  double x=0.25;
  int err;
  double d1 = get_material_param(mdb, "material_1", "param_2", &err);
  CU_ASSERT(err == 0);
  double d2 = get_material_param(mdb, "material_2", "param_2", &err);
  CU_ASSERT(err == 0);
  double v12 = get_alloy_param(mdb, "material_1", "material_2", x, "param_2", &err);
  CU_ASSERT(err == 0);
  double v21 = get_alloy_param(mdb, "material_2", "material_1", (1.0-x), "param_2", &err);
  CU_ASSERT(err == 0);
  CU_ASSERT(fabs(v12 - v21) <= fabs(CLOSE_ENOUGH*(v12)));
  double v_12 = alloy_linear(d1, d2, x);
  double v_21 = alloy_linear(d2, d1, (1-x));
  CU_ASSERT(fabs(v_12 - v_21) <= fabs(CLOSE_ENOUGH*(v_12)));
  CU_ASSERT(fabs(v_12 - v21) <= fabs(CLOSE_ENOUGH*(v_12)));
  CU_ASSERT(fabs(v12 - v_21) <= fabs(CLOSE_ENOUGH*(v_12)));
  CU_ASSERT(fabs(v21 - v_12) <= fabs(CLOSE_ENOUGH*(v_12)));
  CU_ASSERT(fabs(v_21 - v12) <= fabs(CLOSE_ENOUGH*(v_12)));
}

/*Like the above, but shouldn't be called until after 1 is set for
  both params 1 and 2, but 2 only has param1..*/
void get_alloy_early_1_n2_1_b(void) {
  //fprintf(stderr, "get_alloy_early_1_n2_1_b: ");
  double b = 7;
  /*25% mat1 and 75% mat2*/
  double x=0.25;
  int err;
  /*Should fail at setting the bowing; 4 ways of failing*/
  set_bowing_param(mdb, "material_1", "material_2", "param_2", b, &err);
  //fprintf(stderr, "err=%d ", err);
  CU_ASSERT(err == 4);
  set_bowing_param(mdb, "material_2", "material_1", "param_2", b, &err);
  //fprintf(stderr, "err=%d ", err);
  CU_ASSERT(err == 2);
  set_bowing_param(mdb, "material_1", "material_172", "param_2", b, &err);
  //fprintf(stderr, "err=%d ", err);
  CU_ASSERT(err == 3);
  set_bowing_param(mdb, "material_172", "material_1", "param_2", b, &err);
  //fprintf(stderr, "err=%d ", err);
  CU_ASSERT(err == 1);

  /*Should get material1, but not material2.*/
  get_material_param(mdb, "material_1", "param_2", &err);
  //fprintf(stderr, "err=%d ", err);
  CU_ASSERT(err == 0);
  get_material_param(mdb, "material_2", "param_2", &err);
  //fprintf(stderr, "err=%d ", err);
  CU_ASSERT(err == 2);
  /*Should fail to get "material_172", since it has no entry whatsoever*/
  get_material_param(mdb, "material_172", "param_2", &err);
  //fprintf(stderr, "err=%d ", err);
  CU_ASSERT(err == 1);

  /*Should fail to get the bowing; 4 ways of failing.*/
  get_alloy_param(mdb, "material_1", "material_2", x, "param_2", &err);
  //fprintf(stderr, "err=%d ", err);
  CU_ASSERT(err == 5);
  get_alloy_param(mdb, "material_2", "material_1", x, "param_2", &err);
  //fprintf(stderr, "err=%d ", err);
  CU_ASSERT(err == 3);
  get_alloy_param(mdb, "material_1", "material_172", x, "param_2", &err);
  //fprintf(stderr, "err=%d ", err);
  CU_ASSERT(err == 4);
  get_alloy_param(mdb, "material_172", "material_1", x, "param_2", &err);
  //fprintf(stderr, "err=%d ", err);
  CU_ASSERT(err == 2);
  //fprintf(stderr, "\n");
}

/*Like the above, but shouldn't be called until after 1 is set, but
  not 2.*/
void get_alloy_early_1_n2_2_l(void) {
  ////fprintf(stderr, "get_alloy_early_1_n2_2_l\n");
  double b = 7;
  /*25% mat1 and 75% mat2*/
  double x=0.25;
  int err;
  /*Should get material1, but not material2.*/
  get_material_param(mdb, "material_1", "param_2", &err);
  CU_ASSERT(err == 0);
  get_material_param(mdb, "material_2", "param_2", &err);
  CU_ASSERT(err == 2);
  /*Should fail to get "material_172", since it has no entry whatsoever*/
  get_material_param(mdb, "material_172", "param_2", &err);
  CU_ASSERT(err == 1);

  /*Should fail to get the bowing; 4 ways of failing.*/
  get_alloy_param(mdb, "material_1", "material_2", x, "param_2", &err);
  CU_ASSERT(err == 5);
  get_alloy_param(mdb, "material_2", "material_1", x, "param_2", &err);
  CU_ASSERT(err == 3);
  get_alloy_param(mdb, "material_1", "material_172", x, "param_2", &err);
  CU_ASSERT(err == 4);
  get_alloy_param(mdb, "material_172", "material_1", x, "param_2", &err);
  CU_ASSERT(err == 2);
}


void matdb_suite_add_tests(CU_pSuite suite) {
  CU_add_test(suite, "create", create_matdb);
  CU_add_test(suite, "set_get_material_1_1", set_get_material_1_1);
  CU_add_test(suite, "set_get_material_1_2", set_get_material_1_2);
  CU_add_test(suite, "set_get_material_2_1", set_get_material_2_1);
  CU_add_test(suite, "get_alloy_early_1_n2_1_b", get_alloy_early_1_n2_1_b);
  CU_add_test(suite, "get_alloy_early_1_n2_2_l", get_alloy_early_1_n2_2_l);
  CU_add_test(suite, "set_get_material_2_2", set_get_material_2_2);
  CU_add_test(suite, "set_get_alloy_1_2_1_b", set_get_alloy_1_2_1_b);
  CU_add_test(suite, "set_get_alloy_1_2_2_l", set_get_alloy_1_2_2_l);
  CU_add_test(suite, "destroy", destroy_matdb);
}