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 #include "NOX_Abstract_Vector.H"
00043
00044
00045 #include "NOX_MultiVector.H"
00046
00047 NOX::Abstract::Vector& NOX::Abstract::Vector::random(bool useSeed, int seed)
00048 {
00049 cerr << "NOX::Abstract::Vector::random() function not implemented" << endl;
00050 throw "NOX Error";
00051 return *this;
00052 }
00053
00054 void NOX::Abstract::Vector::print(std::ostream& stream) const
00055 {
00056 return;
00057 }
00058
00059 Teuchos::RCP<NOX::Abstract::MultiVector>
00060 NOX::Abstract::Vector::createMultiVector(
00061 const NOX::Abstract::Vector* const* vecs,
00062 int numVecs, NOX::CopyType type) const
00063 {
00064 if (numVecs < 0) {
00065 cerr << "NOX::Abstract::Vector::createMultiVector: Error! Multivector"
00066 << " must have postive number of columns!" << endl;
00067 throw "NOX Error";
00068 }
00069
00070 const NOX::Abstract::Vector** tmp =
00071 new const NOX::Abstract::Vector*[numVecs+1];
00072
00073 tmp[0] = this;
00074 for (int i=0; i<numVecs; i++)
00075 tmp[i+1] = vecs[i];
00076
00077 Teuchos::RCP<NOX::MultiVector> mv =
00078 Teuchos::rcp(new NOX::MultiVector(tmp, numVecs+1, type));
00079
00080 delete [] tmp;
00081
00082 return mv;
00083 }
00084
00085 Teuchos::RCP<NOX::Abstract::MultiVector>
00086 NOX::Abstract::Vector::createMultiVector(int numVecs, NOX::CopyType type) const
00087 {
00088 if (numVecs <= 0) {
00089 cerr << "NOX::Abstract::Vector::createMultiVector: Error! Multivector"
00090 << " must have postive number of columns!" << endl;
00091 throw "NOX Error";
00092 }
00093
00094 return Teuchos::rcp(new NOX::MultiVector(*this, numVecs, type));
00095 }