/*

 * Copyright (c) 2004 David A. Schleef <ds@schleef.org>

 * Copyright (c) 2005 Eric Anholt <anholt@FreeBSD.org>

 * All rights reserved.

 *

 * Redistribution and use in source and binary forms, with or without

 * modification, are permitted provided that the following conditions

 * are met:

 * 1. Redistributions of source code must retain the above copyright

 *    notice, this list of conditions and the following disclaimer.

 * 2. Redistributions in binary form must reproduce the above copyright

 *    notice, this list of conditions and the following disclaimer in the

 *    documentation and/or other materials provided with the distribution.

 * 

 * THIS SOFTWARE IS PROVIDED BY THE AUTHOR ``AS IS'' AND ANY EXPRESS OR

 * IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED

 * WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE

 * ARE DISCLAIMED.  IN NO EVENT SHALL THE AUTHOR BE LIABLE FOR ANY DIRECT,

 * INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES

 * (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR

 * SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION)

 * HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT,

 * STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING

 * IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE

 * POSSIBILITY OF SUCH DAMAGE.

 */

/* Tests math functions against the reference, failing if they differ by more

 * than some epsilon, and printing the difference.

 */

#ifdef HAVE_CONFIG_H

#include "config.h"

#endif

#include <stdio.h>

#include <liboil/liboil.h>

#include <ctype.h>

#include <stdlib.h>

#include <string.h>

#include <math.h>

#ifdef HAVE_INTTYPES_H

#include <inttypes.h>

#endif

#include <liboil/liboilprototype.h>

#include <liboil/liboiltest.h>

#include <liboil/liboilcpu.h>

#include <liboil/liboilparameter.h>

#define LOG_FILE "c:\\logs\\testsuite_align_log.txt"

#include "std_log_result.h"

#define LOG_FILENAME_LINE __FILE__, __LINE__

void create_xml(int result)

{

    if(result)

        assert_failed = 1;

    testResultXml("testsuite_align");

    close_log_file();

}

#ifndef PRIx8

#define PRIx8  "x"

#define PRIx16 "x"

#define PRIx32 "x"

#define PRIx64 "llx"

#define PRId8  "d"

#define PRId16 "d"

#define PRId32 "d"

#define PRId64 "lld"

#define PRIu8  "u"

#define PRIu16 "u"

#define PRIu32 "u"

#define PRIu64 "llu"

#endif

/* Amount by which results of different types are allowed to deviate from the

 * reference.

 */

#define INT_EPSILON 1

#define FLOAT_EPSILON 0.0001

void

dump_array (void *data, void *ref_data, OilType type, int pre_n, int stride,

    int post_n)

{

  int i, j;

  int s2 = oil_type_sizeof (type);

  double x;

#define DUMP(type, format, int) do { \

  for(i=0;i<post_n;i++){ \

    float epsilon = (int) ? INT_EPSILON : FLOAT_EPSILON; \

    std_log(LOG_FILENAME_LINE,"    "); \

    for(j=0;j<pre_n;j++){ \

      x = fabs(OIL_GET(data, i*stride + j*s2, type) - \

          OIL_GET(ref_data, i*stride + j*s2, type)); \

      if (x > epsilon) { \

        std_log(LOG_FILENAME_LINE,"*" format "* (" format ") ", \

        OIL_GET(data, i*stride + j*s2, type), \

        OIL_GET(ref_data, i*stride + j*s2, type)); \

      } else { \

        std_log(LOG_FILENAME_LINE," " format " ", OIL_GET(data, i*stride + j*s2, type)); \

      } \

    } \

    std_log(LOG_FILENAME_LINE,"\n"); \

  } \

} while(0)

  switch(type) {

    case OIL_TYPE_s8p:

    case OIL_TYPE_u8p:

      DUMP(int8_t, "0x%02" PRIx8, 1);

      break;

    case OIL_TYPE_s16p:

    case OIL_TYPE_u16p:

      DUMP(uint16_t, "0x%04" PRIx16, 1);

      break;

    case OIL_TYPE_s32p:

    case OIL_TYPE_u32p:

      DUMP(uint32_t, "0x%08" PRIx32, 1);

      break;

    case OIL_TYPE_f32p:

      DUMP(float, "%g", 0);

      break;

    case OIL_TYPE_s64p:

    case OIL_TYPE_u64p:

      DUMP(uint64_t, "0x%016" PRIx64, 1);

      break;

    case OIL_TYPE_f64p:

      DUMP(double, "%g", 0);

      break;

    default:

      break;

  }

}

void

dump_source (OilTest *test)

{

  int i;

  for(i=0;i<OIL_ARG_LAST;i++){

    OilParameter *p = &test->params[i];

    if (p->is_pointer) {

      if (p->direction == 'i' || p->direction == 's') {

        printf ("  %s:\n", p->parameter_name);

        dump_array (p->src_data + p->test_header,

            p->src_data + p->test_header,

            p->type, p->pre_n, p->stride, p->post_n);

      }

    }

  }

}

void

dump_dest_ref (OilTest *test)

{

  int i;

  for(i=0;i<OIL_ARG_LAST;i++){

    OilParameter *p = &test->params[i];

    if (p->is_pointer) {

      if (p->direction == 'd') {

        printf ("  %s:\n", p->parameter_name);

        dump_array (p->ref_data + p->test_header,

            p->ref_data + p->test_header,

            p->type, p->pre_n, p->stride, p->post_n);

      }

    }

  }

}

int

test_difference (void *data, void *ref_data, OilType type, int pre_n, int stride,

    int post_n)

{

  int i, j;

  int s2 = oil_type_sizeof (type);

  double x;

#define CHECK(type, is_int) do { \

  float epsilon = (is_int) ? INT_EPSILON : FLOAT_EPSILON; \

  for(i=0;i<post_n;i++){ \

    for(j=0;j<pre_n;j++){ \

      x = fabs(OIL_GET(data, i*stride + j*s2, type) - \

          OIL_GET(ref_data, i*stride + j*s2, type)); \

      if (x > epsilon) { \

        return 1; \

      } \

    } \

  } \

  return 0; \

} while(0)

  switch(type) {

    case OIL_TYPE_s8p:

      CHECK(int8_t, 1);

      break;

    case OIL_TYPE_u8p:

      CHECK(uint8_t, 1);

      break;

    case OIL_TYPE_s16p:

      CHECK(int16_t, 1);

      break;

    case OIL_TYPE_u16p:

      CHECK(uint16_t, 1);

      break;

    case OIL_TYPE_s32p:

      CHECK(int32_t, 1);

      break;

    case OIL_TYPE_u32p:

      CHECK(uint32_t, 1);

      break;

    case OIL_TYPE_s64p:

      CHECK(int64_t, 1);

      break;

    case OIL_TYPE_u64p:

      CHECK(uint64_t, 1);

      break;

    case OIL_TYPE_f32p:

      CHECK(float, 0);

      break;

    case OIL_TYPE_f64p:

      CHECK(double, 0);

      break;

    default:

      return 1;

  }

}

int

check_test (OilTest *test)

{

  int i, failed = 0;

  for(i=0;i<OIL_ARG_LAST;i++){

    OilParameter *p = &test->params[i];

    if (p->is_pointer) {

      if (p->direction == 'i' || p->direction == 'd') {

	if (!test_difference(p->test_data + p->test_header,

            p->ref_data + p->test_header,

            p->type, p->pre_n, p->stride, p->post_n))

	  continue;

        printf (" Failure in %s (marked by *, ref in ()):\n",

	    p->parameter_name);

        dump_array (p->test_data + p->test_header,

            p->ref_data + p->test_header,

            p->type, p->pre_n, p->stride, p->post_n);

	failed = 1;

      }

    }

  }

  return failed;

}

int check_class_with_alignment (OilFunctionClass *klass,

    OilArgType arg, int n, int align)

{

  OilParameter *p;

  int align_offset;

  int test_failed = 0;

  OilTest *test;

  OilFunctionImpl *impl;

  test = oil_test_new(klass);

  p = &test->params[arg];

  align_offset = align * oil_type_sizeof(p->type);

  oil_test_set_test_header(test, p, OIL_TEST_HEADER + align_offset);

  oil_test_set_iterations(test, 1);

  test->n = n;

  test->m = n;

  impl = klass->reference_impl;

  oil_test_check_impl (test, impl);

  for (impl = klass->first_impl; impl; impl = impl->next) {

    if (impl == klass->reference_impl)

      continue;

    if (oil_impl_is_runnable (impl)) {

      if (!oil_test_check_impl (test, impl)) {

        printf ("impl %s with arg %d offset %d, n=%d\n", impl->name, arg,

            align_offset, n);

        std_log(LOG_FILENAME_LINE,"dests for %s:\n", klass->name);

        dump_dest_ref(test);

        std_log(LOG_FILENAME_LINE,"sources for %s:\n", klass->name);

        dump_source(test);

      }

    }

  }

  oil_test_free(test);

  if (test_failed)

      std_log(LOG_FILENAME_LINE,"Failed in class %s", klass->name);

  return test_failed;

}

/* Check a function class for all implementations matching the reference when

 * each parameter is varied in its offset from malloc's alignment by 0 - 3 units

 * times size of the type, and with the number of elements varying between 8 and

 * 11.

 */

int check_class(OilFunctionClass *klass)

{

  OilTest *test;

  int failed = 0;

  int i, n;

  oil_class_optimize (klass);

  std_log(LOG_FILENAME_LINE,"checking class %s\n", klass->name);

  test = oil_test_new(klass);

  for (i=0; i < OIL_ARG_LAST; i++) {

    OilParameter *p;

    int align;

    p = &test->params[i];

    if (!p->is_pointer) {

      continue;

    }

    for (n = 8; n <= 11; n++) {

      for (align = 0; align <= 3; align++) {

        failed |= check_class_with_alignment (klass, i, n, align);

      }

    }

  }

  oil_test_free (test);

  return failed;

}

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

{

  int failed = 0;

  int i, n;

  std_log(LOG_FILENAME_LINE,"Test started testsuite_align");

  oil_init ();

  n = oil_class_get_n_classes ();

  for (i = 0; i < n; i++) {

    OilFunctionClass *klass = oil_class_get_by_index(i);

    failed |= check_class(klass);

  }

  if(failed)

      assert_failed = failed;

  if(assert_failed)

            std_log(LOG_FILENAME_LINE,"Test Fail");

  else

            std_log(LOG_FILENAME_LINE,"Test Successful");

  create_xml(0);

  return 0;

}