array.c

basic.c */ /**//* Project:NeuroBasic, basic package*//**/ /* Survey:This is a simple Basic b-code

/****************************************************************/
/*								*/
/* Name:	array.c 					*/
/*								*/
/* Project:	NeuroBasic, basic package	          	*/
/*								*/
/* Survey:	Contains functions for neuro-object "array".	*/
/*		(allocate, load/save, nput/nget, etc.)		*/
/*								*/
/* Author:	Beat Flepp					*/
/*		Electronics Laboratory, ETH Zuerich		*/
/*		Switzerland					*/
/*								*/
/* Created:	November 13, 1997       			*/
/* Modified:	                 				*/
/*								*/
/****************************************************************/


/* system header files */
#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include <ctype.h>
#include <math.h>
#include <musiclib.h>

/* NeuroBasic header files */
#include "host.h"		
#include "allnfcts.h"
#include "allnobjs.h"
#include "neurolib.h"
#include "array.h"



static void release_array(MINT nobj)
/*================================*/
{
  if (neuro_objs[nobj].array.pdata != NULL)
    dfree(neuro_objs[nobj].array.pdata);
  free_nobj(nobj);
} /* end of release_array() */




static MINT save_array(MINT nobj)
/*=============================*/
{
  MINT		ndim, xdim, ydim, zdim, block_size;
  MFLOAT	*pcons;
  comm_def_t	comm_def;

  ndim = 2;
  xdim = neuro_objs[nobj].array.height;
  ydim = neuro_objs[nobj].array.width;
  zdim = 1;
  block_size = xdim * ydim;
  if (!setup_save(ndim, xdim, ydim, zdim, block_size))
    return 0;					/* abort from host */

  comm_def.dim.x = xdim;
  comm_def.dim.y = ydim;
  comm_def.dim.z = 1;
  comm_def.elem_size = sizeof(MFLOAT);
  comm_def.prod.part.x = comm_def.dim.x;
  comm_def.prod.offset.x = 0;
  comm_def.prod.part.y = comm_def.dim.y;
  comm_def.prod.offset.y = 0;
  comm_def.prod.part.z = 1;
  comm_def.prod.offset.z = 0;
  comm_def.prod.nelements = comm_def.prod.part.x
    * comm_def.prod.part.y * comm_def.prod.part.z;
  if (comm_def.prod.nelements == 0)
    comm_def.prod.offset = comm_def.dim;
  Complete_cons_window(&comm_def, HOST_ONLY, 0, 0, 0);
  pcons = neuro_objs[nobj].array.pdata;

  Init_comm(&comm_def, pcons, HOST_IO, COMM_NORM);
  Data_ready(ALL_DATA);

  return 0;
} /* end of save_array() */




static MINT nput_array(MINT nobj, MINT ni, MINT nj, MFLOAT value)
/*=============================================================*/
{
  MINT		index;
  
  if (ni < 0 || ni >= neuro_objs[nobj].array.height
      || nj < 0 || nj >= neuro_objs[nobj].array.width)
  {
    fn_error = ERR_PARCONF;
    return 0;
  }

  /*===== change data =====*/
  index = ni *  neuro_objs[nobj].array.width + nj;
  neuro_objs[nobj].array.pdata[index] = value;

  return 0;
} /* end of nput_array() */




MINT new_array(MINT width, MINT height, MINT format)
/*================================================*/
/* only FLOAT format suppported */
{
  MFLOAT	*pd;
  MINT		failure, nobj;
  MINT	        allocsize, i;
  

  /*===== test parameters for plausibility =====*/
  if ((format != FF_FLOAT32) || 
      (width < 1) || (height < 1))
  {
    fn_error = ERR_PARCONF;
    return NO_OBJECT;
  }

  nobj = alloc_nobj(0);
  if (nobj < 0)
  {
    fn_error = ERR_OBJTABFULL;
    return NO_OBJECT;
  }

  /*===== try to allocate space for matrix =====*/
  failure = FALSE;
  pd = NULL;
  allocsize  = 0;
  switch ((int)format)
  {
    case FF_FLOAT32:
      allocsize = width * height * sizeof(MFLOAT);
      break;
    default: 
      break;
  } /* end switch () */

  if (allocsize > 0)
  {
    pd = (MFLOAT*)dmalloc(allocsize, MT_SLOW);
    if (pd == NULL) failure = TRUE;
  }
  
  if (!failure)
  {						/* success */
    neuro_objs[nobj].public.obtype = N_ARRAY;
    neuro_objs[nobj].public.fn_release = release_array;
    neuro_objs[nobj].public.fn_load = NULL;
    neuro_objs[nobj].public.fn_save = save_array;
    neuro_objs[nobj].public.fn_nput = nput_array;
    neuro_objs[nobj].public.fn_nget = NULL;
    neuro_objs[nobj].public.fn_rand = NULL;

    /*===== put in some data for debugging =====*/
    for (i=0;i<width * height;i++){
      pd[i] = 0.0;
    }

    neuro_objs[nobj].array.format = format;
    neuro_objs[nobj].array.width = width;
    neuro_objs[nobj].array.height = height;
    neuro_objs[nobj].array.pdata = pd;
    neuro_objs[nobj].array.nvalues = width * height;

    return nobj;
  }
  else
  {						/* failure */
    release_array(nobj);
    fn_error = ERR_NOMEM;
    return NO_OBJECT;
  }
} /* end of new_array() */