19 #ifndef _SVECTORBASE_H_ 20 #define _SVECTORBASE_H_ 30 template <
class R >
class VectorBase;
31 template <
class R >
class SSVectorBase;
157 assert(m_elem != 0 || memused == 0);
164 assert(m_elem != 0 || memused == 0);
177 d = (d > e->
idx) ? d : e->
idx;
194 for(
int p = 0; p < n; ++p)
196 if(m_elem[p].
idx == i)
198 assert(index(p) == i);
213 return m_elem[n].
val;
242 return m_elem[n].
idx;
251 return m_elem[n].
idx;
260 return m_elem[n].
val;
269 return m_elem[n].
val;
273 void add(
int i,
const R& v)
276 assert(size() < max());
286 assert(size() <= max());
294 assert(size() < max());
301 assert(size() <= max());
311 void add(
int n,
const int i[],
const R v[])
313 assert(n + size() <= max());
336 set_size(size() + newnnz);
341 void add(
int n,
const int i[],
const S v[])
343 assert(n + size() <= max());
366 set_size(size() + newnnz);
372 assert(n + size() <= max());
392 set_size(size() + newnnz);
396 void remove(
int n,
int m)
405 cpy = (size() - m >= cpy) ? cpy : size() - m;
410 set_size(size() - cpy);
425 int newsize = size() - 1;
429 m_elem[n] = m_elem[newsize];
449 for(l = s, w = s + 1; w < e; l = w, ++w)
462 while(l->
idx > dummy.
idx);
481 for(
int i = size() - 1; i >= 0; --i)
484 maxi =
spxAbs(m_elem[i].val);
497 for(
int i = size() - 1; i >= 0; --i)
500 mini =
spxAbs(m_elem[i].val);
511 return std::sqrt(R(length2()));
629 assert(max() >= sv.
size());
661 assert(max() >= sv.
size());
692 assert(max() >= sv.
size());
694 for(
int i = 0; i < sv.
size(); ++i)
700 assert(isConsistent());
708 bool negateExp =
false)
712 assert(max() >= sv.
size());
716 for(
int i = 0; i < sv.
size(); ++i)
724 for(
int i = 0; i < sv.
size(); ++i)
731 assert(isConsistent());
742 assert(max() >= rowSize);
746 for(i = 0; i < rowSize && i < max(); i++)
748 m_elem[i].
val = rowValues[i];
749 m_elem[i].
idx = rowIndices[i];
776 assert(m_elem != 0 || s == 0);
783 assert(m_elem != 0 || m == 0);
791 assert(n == 0 || elmem != 0);
807 #ifdef ENABLE_CONSISTENCY_CHECKS 811 const int my_size = size();
812 const int my_max = max();
814 if(my_size < 0 || my_max < 0 || my_size > my_max)
817 for(
int i = 1; i < my_size; ++i)
819 for(
int j = 0; j < i; ++j)
822 if(m_elem[i].
idx == m_elem[j].
idx && m_elem[i].
val != 0)
893 #endif // _SVECTORBASE_H_ Rational spxAbs(const Rational &r)
Absolute.
const Nonzero< R > & element(int n) const
The n 'th nonzero element.
Sparse vector nonzero element.
R operator*(const VectorBase< R > &w) const
Inner product.
THREADLOCAL const Real infinity
R maxAbs() const
Maximum absolute value, i.e., infinity norm.
Dense vector.Class VectorBase provides dense linear algebra vectors. Internally, VectorBase wraps std...
void setMem(int n, Nonzero< R > *elmem)
Set the memory area where the nonzeros will be stored.
void set_max(int m)
Set the maximum number of nonzeros in the vector.
SVectorBase< R > & operator=(const SVectorBase< R > &sv)
Assignment operator.
int size() const
Number of used indices.
Nonzero< R > & element(int n)
Reference to the n 'th nonzero element.
R minAbs() const
Minimum absolute value.
void add(int n, const int i[], const R v[])
Append n nonzeros.
void set_size(int s)
Set size of the vector.
void add(int n, const Nonzero< R > e[])
Append n nonzeros.
R & value(int n)
Reference to value of n 'th nonzero.
SVectorBase< R > & operator=(const SVectorBase< S > &sv)
Assignment operator.
R length() const
Floating point approximation of euclidian norm (without any approximation guarantee).
int idx
Index of nonzero element.
Nonzero< R > * mem() const
get pointer to internal memory.
Semi sparse vector.This class implements semi-sparse vectors. Such are VectorBases where the indices ...
Rational spxLdexp(Rational x, int exp)
const R & value(int n) const
Value of n 'th nonzero.
SVectorBase< R > & operator*=(const R &x)
Scaling.
DSVectorBase< R > operator*(const SVectorBase< R > &v, R x)
Scaling.
R operator*(const SVectorBase< S > &w) const
inner product for sparse vectors
int & index(int n)
Reference to index of n 'th nonzero.
R operator[](int i) const
Value to index i.
void add(int i, const R &v)
Append one nonzero (i,v).
int index(int n) const
Index of n 'th nonzero.
bool isConsistent() const
Consistency check.
Everything should be within this namespace.
int dim() const
Dimension of the vector defined as maximal index + 1.
int max() const
Maximal number of indices.
void clear()
Remove all indices.
SVectorBase< Real > & scaleAssign(int scaleExp, const SVectorBase< Real > &sv)
scale and assign
Nonzero< R > & operator=(const Nonzero< S > &vec)
void add(const SVectorBase &sv)
Append nonzeros of sv.
SVectorBase< R > & assignArray(const S *rowValues, const int *rowIndices, int rowSize)
Assignment operator.
void sort()
Sort nonzeros to increasing indices.
int pos(int i) const
Position of index i.
R val
Value of nonzero element.
void add(int i)
Append one uninitialized nonzero.
Sparse vectors.Class SVectorBase provides packed sparse vectors. Such are a sparse vectors...
#define MSGinconsistent(name)
SVectorBase< Real > & scaleAssign(const int *scaleExp, const SVectorBase< Real > &sv, bool negateExp=false)
scale and assign
void add(int n, const int i[], const S v[])
Append n nonzeros.
R length2() const
Squared norm.