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
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
|
/*
* author: Bruno Levy, INRIA, project ALICE
* website: http://www.loria.fr/~levy/software
*
* This library is free software; you can redistribute it and/or
* modify it under the terms of the GNU Lesser General Public
* License version 2.1 as published by the Free Software Foundation
*
* This library is distributed in the hope that it will be useful,
* but WITHOUT ANY WARRANTY; without even the implied warranty of
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
* Lesser General Public License for more details.
*
* You should have received a copy of the GNU Lesser General Public
* License along with this library; if not, write to the Free Software
* Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
*
* Scientific work that use this software can reference the website and
* the following publication:
*
* @INPROCEEDINGS {levy:NMDGP:05,
* AUTHOR = Bruno Levy,
* TITLE = Numerical Methods for Digital Geometry Processing,
* BOOKTITLE =Israel Korea Bi-National Conference,
* YEAR=November 2005,
* URL=http://www.loria.fr/~levy/php/article.php?pub=../publications/papers/2005/Numerics
* }
*
* Laurent Saboret 01/2005: Change for CGAL:
* - Added OpenNL namespace
* - FullVector is now a model of the SparseLinearAlgebraTraits_d::Vector concept
* - Coefficients are initialized with zeros
*/
#ifndef __OPENNL_FULL_VECTOR__
#define __OPENNL_FULL_VECTOR__
#include <OpenNL/blas.h>
#include <cstdlib>
#include <cassert>
namespace OpenNL {
// Class FullVector
// Model of the SparseLinearAlgebraTraits_d::Vector concept
template <class T> class FullVector
{
// Public types
public:
typedef T CoeffType ;
// for SparseLinearAlgebraTraits_d::Vector concept
typedef T NT;
// Create a vector initialized with zeros
FullVector(unsigned int dim) {
dimension_ = dim ; coeff_ = new T[dim] ;
// init with zeros (for SparseLinearAlgebraTraits_d::Vector concept)
for (unsigned int i=0; i < dimension_; i++)
coeff_[i] = 0;
}
// Copy constructor
FullVector(const FullVector& toCopy) {
dimension_ = toCopy.dimension_ ;
coeff_ = new T[dimension_] ;
for (unsigned int i=0; i < dimension_; i++)
coeff_[i] = toCopy.coeff_[i];
}
// operator =()
FullVector& operator =(const FullVector& toCopy) {
delete[] coeff_ ;
dimension_ = toCopy.dimension_ ;
coeff_ = new T[dimension_] ;
for (unsigned int i=0; i < dimension_; i++)
coeff_[i] = toCopy.coeff_[i];
return *this;
}
~FullVector() {
delete[] coeff_ ;
coeff_ = NULL ;
}
// Return the vector's number of coefficients
unsigned int dimension() const {
return dimension_ ;
}
// Read/write access to 1 vector coefficient.
//
// Preconditions:
// * 0 <= i < dimension()
const T& operator[](unsigned int i) const {
assert(i < dimension_) ;
return coeff_[i] ;
}
T& operator[](unsigned int i) {
assert(i < dimension_) ;
return coeff_[i] ;
}
private:
unsigned int dimension_ ;
T* coeff_ ;
} ;
template <class T> class BLAS< FullVector<T> > {
public:
typedef FullVector<T> VectorType ;
typedef T CoeffType ;
/** y <- y + a*x */
static void axpy(CoeffType a, const VectorType& x, VectorType& y) {
assert(x.dimension() == y.dimension()) ;
for(unsigned int i=0; i<x.dimension(); i++) {
y[i] += a * x[i] ;
}
}
/** x <- a*x */
static void scal(CoeffType a, VectorType& x) {
for(unsigned int i=0; i<x.dimension(); i++) {
x[i] *= a ;
}
}
/** y <- x */
static void copy(const VectorType& x, VectorType& y) {
assert(x.dimension() == y.dimension()) ;
for(unsigned int i=0; i<x.dimension(); i++) {
y[i] = x[i] ;
}
}
/** returns x^t * y */
static CoeffType dot(const VectorType& x, const VectorType& y) {
CoeffType result = 0 ;
assert(x.dimension() == y.dimension()) ;
for(unsigned int i=0; i<x.dimension(); i++) {
result += y[i] * x[i] ;
}
return result ;
}
} ;
}; // namespace OpenNL
#endif // __OPENNL_FULL_VECTOR__
|
