#ifndef GSL_MATRIX_USHORT_H
#define GSL_MATRIX_USHORT_H

#include <stdlib.h>
#include <gsl_errno.h>
#include <gsl_vector_ushort.h>

typedef struct gsl_matrix_ushort_struct gsl_matrix_ushort;

struct gsl_matrix_ushort_struct
{
  size_t size1;
  size_t size2;
  size_t dim2;
  unsigned short * data;
  gsl_block_ushort * block;
} ;

gsl_matrix_ushort * 
gsl_matrix_ushort_alloc (size_t n1, size_t n2);

gsl_matrix_ushort * 
gsl_matrix_ushort_calloc (size_t n1, size_t n2);

gsl_matrix_ushort * 
gsl_matrix_ushort_alloc_from_block (gsl_block_ushort * b, size_t offset, 
                                   size_t n1, size_t n2, size_t d2);

gsl_matrix_ushort * 
gsl_matrix_ushort_alloc_from_matrix (gsl_matrix_ushort * m,
                                    size_t k1, size_t k2,
                                    size_t n1, size_t n2);

gsl_vector_ushort * 
gsl_vector_ushort_alloc_row_from_matrix (gsl_matrix_ushort * m,
                                        size_t i);

gsl_vector_ushort * 
gsl_vector_ushort_alloc_col_from_matrix (gsl_matrix_ushort * m,
                                        size_t j);

void gsl_matrix_ushort_free (gsl_matrix_ushort * m);

unsigned short * gsl_matrix_ushort_ptr(const gsl_matrix_ushort * m, size_t i, size_t j);
unsigned short   gsl_matrix_ushort_get(const gsl_matrix_ushort * m, size_t i, size_t j);
void    gsl_matrix_ushort_set(gsl_matrix_ushort * m, size_t i, size_t j, unsigned short x);

int gsl_matrix_ushort_fread (FILE * stream, gsl_matrix_ushort * m) ;
int gsl_matrix_ushort_fwrite (FILE * stream, const gsl_matrix_ushort * m) ;
int gsl_matrix_ushort_fscanf (FILE * stream, gsl_matrix_ushort * m);
int gsl_matrix_ushort_fprintf (FILE * stream, const gsl_matrix_ushort * m, const char * format);
 
int gsl_matrix_ushort_copy_row(gsl_vector_ushort * v, const gsl_matrix_ushort * m, size_t i);
int gsl_matrix_ushort_copy_col(gsl_vector_ushort * v, const gsl_matrix_ushort * m, size_t j);
int gsl_matrix_ushort_set_row(gsl_matrix_ushort * m, size_t i, const gsl_vector_ushort * v);
int gsl_matrix_ushort_set_col(gsl_matrix_ushort * m, size_t j, const gsl_vector_ushort * v);

int gsl_matrix_ushort_swap_rows(gsl_matrix_ushort * m, size_t i, size_t j);
int gsl_matrix_ushort_swap_cols(gsl_matrix_ushort * m, size_t i, size_t j);
int gsl_matrix_ushort_swap_rowcol(gsl_matrix_ushort * m, size_t i, size_t j);

int gsl_matrix_ushort_copy(gsl_matrix_ushort * dest, const gsl_matrix_ushort * src);

int gsl_vector_ushort_view_row_from_matrix (gsl_vector_ushort * v, gsl_matrix_ushort * m, size_t i);
int gsl_vector_ushort_view_col_from_matrix (gsl_vector_ushort * v, gsl_matrix_ushort * m, size_t j);

int gsl_matrix_ushort_view_from_vector (gsl_matrix_ushort * m, 
                                       gsl_vector_ushort * base,
                                       size_t offset, 
                                       size_t n1, size_t n2, size_t d2);

extern int gsl_check_range ;

/* inline functions if you are using GCC */

#ifdef HAVE_INLINE
extern inline 
unsigned short
gsl_matrix_ushort_get(const gsl_matrix_ushort * m, 
		     const size_t i, const size_t j)
{
#ifndef GSL_RANGE_CHECK_OFF
  if (i >= m->size1)
    {
      GSL_ERROR_RETURN("first index out of range", GSL_EINVAL, 0) ;
    }
  else if (j >= m->size2)
    {
      GSL_ERROR_RETURN("second index out of range", GSL_EINVAL, 0) ;
    }
#endif
  return m->data[i * m->dim2 + j] ;
} 

extern inline 
void
gsl_matrix_ushort_set(gsl_matrix_ushort * m, 
		     const size_t i, const size_t j, const unsigned short x)
{
#ifndef GSL_RANGE_CHECK_OFF
  if (i >= m->size1)
    {
      GSL_ERROR_RETURN_NOTHING("first index out of range", GSL_EINVAL) ;
    }
  else if (j >= m->size2)
    {
      GSL_ERROR_RETURN_NOTHING("second index out of range", GSL_EINVAL) ;
    }
#endif
  m->data[i * m->dim2 + j] = x ;
}
#endif

#endif /* GSL_MATRIX_USHORT_H */