htswanalysis/MACS/lib/gsl/gsl-1.11/vector/vector_source.c

   1 /* vector/vector_source.c
   2  *
   3  * Copyright (C) 1996, 1997, 1998, 1999, 2000, 2007 Gerard Jungman, Brian Gough
   4  *
   5  * This program is free software; you can redistribute it and/or modify
   6  * it under the terms of the GNU General Public License as published by
   7  * the Free Software Foundation; either version 3 of the License, or (at
   8  * your option) any later version.
   9  *
  10  * This program is distributed in the hope that it will be useful, but
  11  * WITHOUT ANY WARRANTY; without even the implied warranty of
  12  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  13  * General Public License for more details.
  14  *
  15  * You should have received a copy of the GNU General Public License
  16  * along with this program; if not, write to the Free Software
  17  * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301, USA.
  18  */
  19
  20 #ifndef HIDE_INLINE_STATIC
  21 BASE
  22 FUNCTION (gsl_vector, get) (const TYPE (gsl_vector) * v, const size_t i)
  23 {
  24   if (gsl_check_range)
  25     {
  26       if (i >= v->size)         /* size_t is unsigned, can't be negative */
  27         {
  28           BASE zero = ZERO;
  29           GSL_ERROR_VAL ("index out of range", GSL_EINVAL, zero);
  30         }
  31     }
  32
  33   /* The following line is a generalization of return v->data[i] */
  34
  35   return *(BASE *) (v->data + MULTIPLICITY * i * v->stride);
  36 }
  37
  38 void
  39 FUNCTION (gsl_vector, set) (TYPE (gsl_vector) * v, const size_t i, BASE x)
  40 {
  41   if (gsl_check_range)
  42     {
  43       if (i >= v->size)         /* size_t is unsigned, can't be negative */
  44         {
  45           GSL_ERROR_VOID ("index out of range", GSL_EINVAL);
  46         }
  47     }
  48
  49   /* The following line is a generalization of v->data[i] = x */
  50
  51   *(BASE *) (v->data + MULTIPLICITY * i * v->stride) = x;
  52 }
  53
  54 BASE *
  55 FUNCTION (gsl_vector, ptr) (TYPE (gsl_vector) * v, const size_t i)
  56 {
  57   if (gsl_check_range)
  58     {
  59       if (i >= v->size)         /* size_t is unsigned, can't be negative */
  60         {
  61           GSL_ERROR_NULL ("index out of range", GSL_EINVAL);
  62         }
  63     }
  64
  65   return (BASE *) (v->data + MULTIPLICITY * i * v->stride);
  66 }
  67
  68 const BASE *
  69 FUNCTION (gsl_vector, const_ptr) (const TYPE (gsl_vector) * v, const size_t i)
  70 {
  71   if (gsl_check_range)
  72     {
  73       if (i >= v->size)         /* size_t is unsigned, can't be negative */
  74         {
  75           GSL_ERROR_NULL ("index out of range", GSL_EINVAL);
  76         }
  77     }
  78
  79   /* The following line is a generalization of return v->data[i] */
  80
  81   return (const BASE *) (v->data + MULTIPLICITY * i * v->stride);
  82 }
  83 #endif