1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
|
/**
*/
/*
*
* Template Numerical Toolkit (TNT): Linear Algebra Module
*
* Mathematical and Computational Sciences Division
* National Institute of Technology,
* Gaithersburg, MD USA
*
*
* This software was developed at the National Institute of Standards and
* Technology (NIST) by employees of the Federal Government in the course
* of their official duties. Pursuant to title 17 Section 105 of the
* United States Code, this software is not subject to copyright protection
* and is in the public domain. The Template Numerical Toolkit (TNT) is
* an experimental system. NIST assumes no responsibility whatsoever for
* its use by other parties, and makes no guarantees, expressed or implied,
* about its quality, reliability, or any other characteristic.
*
* BETA VERSION INCOMPLETE AND SUBJECT TO CHANGE
* see http://math.nist.gov/tnt for latest updates.
*
*/
// Vector/Matrix/Array Index Module
#ifndef INDEX_H
#define INDEX_H
#include "subscript.h"
namespace TNT
{
class Index1D
{
Subscript lbound_;
Subscript ubound_;
public:
Subscript lbound() const { return lbound_; }
Subscript ubound() const { return ubound_; }
Index1D(const Index1D &D) : lbound_(D.lbound_), ubound_(D.ubound_) {}
Index1D(Subscript i1, Subscript i2) : lbound_(i1), ubound_(i2) {}
Index1D & operator=(const Index1D &D)
{
lbound_ = D.lbound_;
ubound_ = D.ubound_;
return *this;
}
};
inline Index1D operator+(const Index1D &D, Subscript i)
{
return Index1D(i+D.lbound(), i+D.ubound());
}
inline Index1D operator+(Subscript i, const Index1D &D)
{
return Index1D(i+D.lbound(), i+D.ubound());
}
inline Index1D operator-(Index1D &D, Subscript i)
{
return Index1D(D.lbound()-i, D.ubound()-i);
}
inline Index1D operator-(Subscript i, Index1D &D)
{
return Index1D(i-D.lbound(), i-D.ubound());
}
} // namespace TNT
#endif
|
