#include "aubio_priv.h"
#include "fmat.h"
#include "tensor.h"

aubio_tensor_t *new_aubio_tensor(uint_t ndim, uint_t *shape)
{
  aubio_tensor_t *c = AUBIO_NEW(aubio_tensor_t);
  uint_t i;

  if ((sint_t)ndim <= 0) goto failure;
  for (i = 0; i < ndim; i++) {
    if ((sint_t)shape[i] <= 0) goto failure;
  }

  c->ndim = ndim;
  c->items_per_row = 1;
  //c->shape = AUBIO_ARRAY(uint_t, ndim);
  c->shape[0] = shape[0];
  for (i = 1; i < ndim; i++) {
    c->shape[i] = shape[i];
    c->items_per_row *= shape[i];
  }
  c->n_items = c->items_per_row * shape[0];
  c->data = AUBIO_ARRAY(smpl_t*, shape[0]);
  c->data[0] = AUBIO_ARRAY(smpl_t, c->n_items);
  for (i = 1; i < c->shape[0]; i++) {
    c->data[i] = c->data[0] + i * c->items_per_row;
  }

  return c;

failure:
  del_aubio_tensor(c);
  return NULL;
}

void del_aubio_tensor(aubio_tensor_t *c)
{
  AUBIO_ASSERT(c);
  if (c->data) {
    if (c->data[0]) {
      AUBIO_FREE(c->data[0]);
    }
    AUBIO_FREE(c->data);
  }
  //if (c->shape)
  //  AUBIO_FREE(c->shape);
  AUBIO_FREE(c);
}

uint_t aubio_tensor_as_fvec(aubio_tensor_t *c, fvec_t *o) {
  if (c->ndim  != 1) return AUBIO_FAIL;
  if (c->shape[0] <= 0) return AUBIO_FAIL;
  o->length = c->shape[0];
  o->data = c->data[0];
  return AUBIO_OK;
}

uint_t aubio_fvec_as_tensor(fvec_t *o, aubio_tensor_t *c) {
  if (o == NULL) return AUBIO_FAIL;
  c->ndim = 1;
  c->shape[0] = o->length;
  c->data = &o->data;
  return AUBIO_OK;
}

uint_t aubio_tensor_as_fmat(aubio_tensor_t *c, fmat_t *o) {
  if (c->ndim  != 2) return AUBIO_FAIL;
  if (c->shape[0] <= 0) return AUBIO_FAIL;
  if (c->shape[1] <= 0) return AUBIO_FAIL;
  o->height = c->shape[0];
  o->length = c->shape[1];
  o->data = c->data;
  return AUBIO_OK;
}

uint_t aubio_fmat_as_tensor(fmat_t *o, aubio_tensor_t *c) {
  if (o == NULL) return AUBIO_FAIL;
  if (c == NULL) return AUBIO_FAIL;
  c->ndim = 2;
  c->shape[0] = o->height;
  c->shape[1] = o->length;
  c->data = o->data;
  return AUBIO_OK;
}

smpl_t aubio_tensor_max(aubio_tensor_t *t)
{
  uint_t i;
  smpl_t max = -1000000;
  for (i = 0; i < t->n_items; i++) {
    max = MAX(t->data[0][i], max);
  }
  return max;
}