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);
360 for(
int i = vec.
size() - 1; i >= 0; --i )
397 for(
int i = vec.
size() - 1; i >= 0; --i )
419 for(
int i =
size() - 1; i >= 0; --i )
435 int j = w.
size() - 1;
438 if( i >= 0 && j >= 0 )
443 while( i != 0 && j != 0 )
459 while( i != 0 && vi != wj )
462 while( j != 0 && vi != wj )
474 template <
class S,
class T >
478 template <
class S,
class T >
493 template <
class S,
class T >
497 template <
class S,
class T >
501 template <
class S,
class T >
507 template <
class S,
class T >
517 for(
int i = 0; i <
num; ++i )
538 for(
int i = 0; i <
num; ++i )
570 if(
index(i) >= newdim )
604 for(
int i = 0; i <
num; ++i )
619 #ifdef ENABLE_CONSISTENCY_CHECKS 628 for(
int i = 0; i < VectorBase<R>::dim(); ++i )
634 MSG_ERROR( std::cerr <<
"ESSVEC01 i = " << i
663 len = (p_dim < 1) ? 1 : p_dim;
729 for(
int i =
size() - 1; i >= 0; --i )
740 for(
int i = 0; i < d; ++i )
742 if( rhs.
val[i] != 0 )
787 for(
int i =
size() - 1; i >= 0; --i )
798 for(
int i = 0; i < d; ++i )
833 for(
int i =
size() - 1; i >= 0; --i )
843 for(
int i = 0; i < rhs.
dim(); ++i )
894 template <
class S,
class T >
898 template <
class S,
class T >
902 template <
class S,
class T >
909 #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.