← Back to branch summary

~ubuntu-branches/debian/sid/gsl/sid

« back to all changes in this revision

Viewing changes to matrix/matrix_source.c

  • Committer: Bazaar Package Importer
  • Author(s): Dirk Eddelbuettel
  • Date: 2008-12-16 06:17:55 UTC
  • mfrom: (1.3.2 upstream) (3.1.15 jaunty)
  • Revision ID: james.westby@ubuntu.com-20081216061755-9la7p0qwrhopk8pl
Tags: 1.12+dfsg-1
* New upstream version released today

* doc/*: As before, removed the 'non-free' documentation to create a 
  source package that complies with Debian's interpretation of what is free. 

Show diffs side-by-side

added added

removed removed

Lines of Context:
matrix/matrix_source.c
1
 
/* matrix/matrix_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_matrix, get) (const TYPE (gsl_matrix) * m,
23
 
                            const size_t i, const size_t j)
24
 
{
25
 
  BASE zero = ZERO;
26
 
 
27
 
  if (gsl_check_range)
28
 
    {
29
 
      if (i >= m->size1)        /* size_t is unsigned, can't be negative */
30
 
        {
31
 
          GSL_ERROR_VAL ("first index out of range", GSL_EINVAL, zero);
32
 
        }
33
 
      else if (j >= m->size2)   /* size_t is unsigned, can't be negative */
34
 
        {
35
 
          GSL_ERROR_VAL ("second index out of range", GSL_EINVAL, zero);
36
 
        }
37
 
    }
38
 
  return *(BASE *) (m->data + MULTIPLICITY * (i * m->tda + j));
39
 
}
40
 
 
41
 
void
42
 
FUNCTION (gsl_matrix, set) (TYPE (gsl_matrix) * m,
43
 
                            const size_t i, const size_t j,
44
 
                            const BASE x)
45
 
{
46
 
  if (gsl_check_range)
47
 
    {
48
 
      if (i >= m->size1)        /* size_t is unsigned, can't be negative */
49
 
        {
50
 
          GSL_ERROR_VOID ("first index out of range", GSL_EINVAL);
51
 
        }
52
 
      else if (j >= m->size2)   /* size_t is unsigned, can't be negative */
53
 
        {
54
 
          GSL_ERROR_VOID ("second index out of range", GSL_EINVAL);
55
 
        }
56
 
    }
57
 
  *(BASE *) (m->data + MULTIPLICITY * (i * m->tda + j)) = x;
58
 
}
59
 
 
60
 
 
61
 
BASE *
62
 
FUNCTION (gsl_matrix, ptr) (TYPE (gsl_matrix) * m,
63
 
                            const size_t i, const size_t j)
64
 
{
65
 
  if (gsl_check_range)
66
 
    {
67
 
      if (i >= m->size1)        /* size_t is unsigned, can't be negative */
68
 
        {
69
 
          GSL_ERROR_NULL ("first index out of range", GSL_EINVAL);
70
 
        }
71
 
      else if (j >= m->size2)   /* size_t is unsigned, can't be negative */
72
 
        {
73
 
          GSL_ERROR_NULL ("second index out of range", GSL_EINVAL);
74
 
        }
75
 
    }
76
 
  return (BASE *) (m->data + MULTIPLICITY * (i * m->tda + j));
77
 
}
78
 
 
79
 
 
80
 
const BASE *
81
 
FUNCTION (gsl_matrix, const_ptr) (const TYPE (gsl_matrix) * m,
82
 
                                  const size_t i, const size_t j)
83
 
{
84
 
  if (gsl_check_range)
85
 
    {
86
 
      if (i >= m->size1)        /* size_t is unsigned, can't be negative */
87
 
        {
88
 
          GSL_ERROR_NULL ("first index out of range", GSL_EINVAL);
89
 
        }
90
 
      else if (j >= m->size2)   /* size_t is unsigned, can't be negative */
91
 
        {
92
 
          GSL_ERROR_NULL ("second index out of range", GSL_EINVAL);
93
 
        }
94
 
    }
95
 
  return (const BASE *) (m->data + MULTIPLICITY * (i * m->tda + j));
96
 
}
97
 
#endif