20 #ifndef _SSVECTORBASE_H_ 21 #define _SSVECTORBASE_H_ 32 template <
class R >
class SVectorBase;
141 for(
int i = 0; i < d; ++i )
185 assert(n >= 0 && n <
size());
272 assert(
index(n) >= 0);
359 for(
int i = vec.
size() - 1; i >= 0; --i )
396 for(
int i = vec.
size() - 1; i >= 0; --i )
418 for(
int i =
size() - 1; i >= 0; --i )
434 int j = w.
size() - 1;
437 if( i >= 0 && j >= 0 )
442 while( i != 0 && j != 0 )
458 while( i != 0 && vi != wj )
461 while( j != 0 && vi != wj )
473 template <
class S,
class T >
477 template <
class S,
class T >
492 template <
class S,
class T >
496 template <
class S,
class T >
500 template <
class S,
class T >
506 template <
class S,
class T >
516 for(
int i = 0; i <
num; ++i )
537 for(
int i = 0; i <
num; ++i )
569 if(
index(i) >= newdim )
603 for(
int i = 0; i <
num; ++i )
618 #ifdef ENABLE_CONSISTENCY_CHECKS 627 for(
int i = 0; i < VectorBase<R>::dim(); ++i )
633 MSG_ERROR( std::cerr <<
"ESSVEC01 i = " << i
662 len = (p_dim < 1) ? 1 : p_dim;
728 for(
int i =
size() - 1; i >= 0; --i )
739 for(
int i = 0; i < d; ++i )
741 if( rhs.
val[i] != 0 )
786 for(
int i =
size() - 1; i >= 0; --i )
797 for(
int i = 0; i < d; ++i )
832 for(
int i =
size() - 1; i >= 0; --i )
842 for(
int i = 0; i < rhs.
dim(); ++i )
893 template <
class S,
class T >
897 template <
class S,
class T >
901 template <
class S,
class T >
908 #endif // _SSVECTORBASE_H_ Rational spxAbs(const Rational &r)
Absolute.
R maxAbs() const
Maximum absolute value, i.e., infinity norm.
void add(int n)
appends n uninitialized indices.
SSVectorBase< R > & assign2productFull(const SVSetBase< S > &A, const SSVectorBase< T > &x)
Assignment helper.
int * idx
array of indices
bool isSetup() const
Returns setup status.
void reDim(int newdim, const bool setZero=true)
Resets DVectorBase's dimension to newdim.
bool isConsistent() const
Consistency check.
VectorBase< R > & operator=(const VectorBase< S > &vec)
Assignment operator.
Memory allocation routines.
int max() const
returns the maximal number of indices which can be stored in IdxSet.
SSVectorBase< R > & assign(const SVectorBase< S > &rhs)
Assigns only the elements of rhs.
R length2() const
Squared norm.
Dynamic dense vectors.Class DVectorBase is a derived class of VectorBase adding automatic memory mana...
Dense vector.Class VectorBase provides dense linear algebra vectors. It does not provide memory manag...
int pos(int i) const
Finds the position of index i in the IdxSet, or -1 if i doesn't exist.
SSVectorBase< R > & operator-=(const VectorBase< S > &vec)
Subtraction.
SSVectorBase< R > & operator+=(const SSVectorBase< S > &vec)
Addition.
SSVectorBase< R > & assign2product1(const SVSetBase< S > &A, const SSVectorBase< T > &x)
Assignment helper.
Dynamic sparse vectors.Class DSVectorBase implements dynamic sparse vectors, i.e. SVectorBases with a...
void unSetup()
Makes SSVectorBase not setup.
void setValue(int i, R x)
Sets i 'th element to x.
void clear()
removes all indices.
SSVectorBase< R > & operator=(const SSVectorBase< R > &rhs)
Assignment operator.
void reMem(int newsize)
Resets memory consumption to newsize.
SSVectorBase< R > & operator*=(S x)
Scaling.
IdxSet()
default constructor.
SSVectorBase< R > & multAdd(S x, const VectorBase< T > &vec)
Addition of a scaled vector.
R operator*(const SSVectorBase< S > &w)
int * altIndexMem()
Returns array indices.
bool setupStatus
Is the SSVectorBase set up?
VectorBase< R > & operator+=(const VectorBase< S > &vec)
Addition.
void addIdx(int i)
appends index i.
Semi sparse vector.This class implements semi-sparse vectors. Such are DVectorBases where the indices...
bool isConsistent() const
consistency check.
int dim() const
Dimension of VectorBase.
SSVectorBase< R > & assign2productAndSetup(const SVSetBase< S > &A, SSVectorBase< T > &x)
Assigns SSVectorBase to thereby setting up x.
void spx_alloc(T &p, int n=1)
Allocate memory.
R * altValues()
Returns array values.
R * get_ptr()
Conversion to C-style pointer.
void reSize(int newsize)
Resets DVectorBase's memory size to newsize.
Real spxLdexp(Real x, int exp)
returns x * 2^exp
const R * values() const
Returns array values.
SSVectorBase< R > & assign2product(const SSVectorBase< S > &x, const SVSetBase< T > &A)
Assigns to SSVectorBase.
int pos(int i) const
returns the position of index i.
int size() const
returns the number of used indices.
#define MSG_ERROR(x)
Prints out message x if the verbosity level is at least SPxOut::ERROR.
int len
length of array idx
const IdxSet & indices() const
Returns indices.
void clearNum(int n)
Sets n 'th nonzero element to 0 (index n must exist).
IdxSet & operator=(const IdxSet &set)
assignment operator.
R operator[](int i) const
Returns i 'th value.
Real spxSqrt(Real a)
returns square root
void clear()
Clears vector.
IdxSet & altIndices()
Returns indices.
void setEpsilon(R eps)
Changes the non-zero epsilon, invalidating the setup. */.
bool isConsistent() const
consistency check.
SSVectorBase< R > & multAdd(S xx, const SVectorBase< T > &vec)
Addition of a scaled vector.
int index(int n) const
Returns index of the n 'th nonzero element.
SSVectorBase< R > & operator-=(const SSVectorBase< S > &vec)
Subtraction.
const int * indexMem() const
Returns array indices.
void reDim(int newdim)
Resets dimension to newdim.
Debugging, floating point type and parameter definitions.
void setSize(int n)
Sets number of nonzeros (thereby unSetup SSVectorBase).
void clearIdx(int i)
Clears element i.
R * get_ptr()
Only used in slufactor.cpp.
void spx_realloc(T &p, int n)
Change amount of allocated memory.
R length2() const
Squared euclidian norm.
VectorBase< R > & multAdd(const S &x, const VectorBase< T > &vec)
Addition of scaled vector.
void setup_and_assign(SSVectorBase< S > &rhs)
Sets up rhs vector, and assigns it.
int size() const
Returns the number of nonzeros.
Everything should be within this namespace.
int dim() const
Dimension of the vector defined as maximal index + 1.
void setMax(int newmax)
Allocates enough space to accommodate newmax values.
R getEpsilon() const
Returns the non-zero epsilon used.
SSVectorBase< R > & operator=(const SSVectorBase< S > &rhs)
Assignment operator.
SSVectorBase< R > & operator+=(const VectorBase< S > &vec)
Addition.
VectorBase< R > & operator-=(const VectorBase< S > &vec)
Subtraction.
void setup()
Initializes nonzero indices for elements with absolute values above epsilon and sets all other elemen...
void clear()
Set vector to 0.
void forceSetup()
Forces setup status.
SSVectorBase< R > & operator=(const VectorBase< S > &rhs)
Assignment operator.
Real length() const
Floating point approximation of euclidian norm (without any approximation guarantee).
int dim() const
returns the maximal index.
Sparse vectors.Class SVectorBase provides packed sparse vectors. Such are a sparse vectors...
SSVectorBase< R > & assignPWproduct4setup(const SSVectorBase< S > &x, const SSVectorBase< T > &y)
Assigns pair wise vector product to SSVectorBase.
#define MSGinconsistent(name)
R maxAbs() const
Maximum absolute value, i.e., infinity norm.
R epsilon
A value x with |x| < epsilon is considered zero.
int index(int n) const
access n 'th index.
int num
number of used indices
void add(int i, R x)
Adds nonzero (i, x) to SSVectorBase.
SSVectorBase< R > & assign2productShort(const SVSetBase< S > &A, const SSVectorBase< T > &x)
Assignment helper.
Set of indices.Class IdxSet provides a set of indices. At construction it must be given an array of i...
void spx_free(T &p)
Release memory.
Sparse vector set.Class SVSetBase provides a set of sparse vectors SVectorBase. All SVectorBases in a...
void scaleValue(int i, int scaleExp)
Scale i 'th element by a.
R value(int n) const
Returns value of the n 'th nonzero element.
SSVectorBase< R > & assign2product4setup(const SVSetBase< S > &A, const SSVectorBase< T > &x)
Assigns SSVectorBase to for a setup x.