00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023
00024
00025
00026
00027
00028
00029
00030
00031
00032
00033
00034
00035
00036
00037
00038
00039
00040
00041
00042 #ifndef NOX_ABSTRACT_MULTIVECTOR_H
00043 #define NOX_ABSTRACT_MULTIVECTOR_H
00044
00045 #include "NOX_Common.H"
00046
00047 #include "Teuchos_SerialDenseMatrix.hpp"
00048
00049 #include "NOX_Abstract_Vector.H"
00050
00051 namespace NOX {
00052
00053 namespace Abstract {
00054
00058 class MultiVector {
00059
00060 public:
00061
00063 typedef Teuchos::SerialDenseMatrix<int,double> DenseMatrix;
00064
00066 MultiVector() {};
00067
00069 virtual ~MultiVector() {};
00070
00072
00074 virtual NOX::Abstract::MultiVector&
00075 init(double gamma) = 0;
00076
00078 virtual NOX::Abstract::MultiVector&
00079 random(bool useSeed = false, int seed = 1) = 0;
00080
00082 virtual NOX::Abstract::MultiVector&
00083 operator=(const NOX::Abstract::MultiVector& source) = 0;
00084
00090 virtual NOX::Abstract::MultiVector&
00091 setBlock(const NOX::Abstract::MultiVector& source,
00092 const vector<int>& index) = 0;
00093
00095 virtual NOX::Abstract::MultiVector&
00096 augment(const NOX::Abstract::MultiVector& source) = 0;
00097
00099
00101
00110 virtual NOX::Abstract::Vector& operator [] (int i) = 0;
00111
00120 virtual const NOX::Abstract::Vector& operator [] (int i) const = 0;
00121
00123
00125
00127 virtual NOX::Abstract::MultiVector&
00128 scale(double gamma) = 0;
00129
00134 virtual NOX::Abstract::MultiVector&
00135 update(double alpha, const NOX::Abstract::MultiVector& a,
00136 double gamma = 0.0) = 0;
00137
00142 virtual NOX::Abstract::MultiVector&
00143 update(double alpha, const NOX::Abstract::MultiVector& a,
00144 double beta, const NOX::Abstract::MultiVector& b,
00145 double gamma = 0.0) = 0;
00146
00153 virtual NOX::Abstract::MultiVector&
00154 update(Teuchos::ETransp transb, double alpha,
00155 const NOX::Abstract::MultiVector& a,
00156 const DenseMatrix& b, double gamma = 0.0) = 0;
00157
00159
00161
00175 virtual Teuchos::RCP<NOX::Abstract::MultiVector>
00176 clone(NOX::CopyType type = NOX::DeepCopy) const = 0;
00177
00179 virtual Teuchos::RCP<NOX::Abstract::MultiVector>
00180 clone(int numvecs) const = 0;
00181
00186 virtual Teuchos::RCP<NOX::Abstract::MultiVector>
00187 subCopy(const vector<int>& index) const = 0;
00188
00193 virtual Teuchos::RCP<NOX::Abstract::MultiVector>
00194 subView(const vector<int>& index) const = 0;
00195
00197
00199
00201 virtual void
00202 norm(vector<double>& result,
00203 NOX::Abstract::Vector::NormType type =
00204 NOX::Abstract::Vector::TwoNorm) const = 0;
00205
00207
00209
00213 virtual void multiply(double alpha, const NOX::Abstract::MultiVector& y,
00214 DenseMatrix& b) const = 0;
00215
00217
00219 virtual int length() const = 0;
00220
00222 virtual int numVectors() const = 0;
00223
00225 virtual void print(std::ostream& stream) const = 0;
00226
00227 };
00228 }
00229 }
00230
00231 #endif // NOX_ABSTRACT_MULTIVECTOR_H