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 #include "SundanceEvalVectorArray.hpp"
00033 #include "SundanceSparsitySuperset.hpp"
00034 #include "PlayaTabs.hpp"
00035
00036 using namespace Sundance;
00037 using namespace Sundance;
00038
00039 using namespace Teuchos;
00040
00041 EvalVectorArray::EvalVectorArray()
00042 : Array<RCP<EvalVector> >()
00043 {;}
00044
00045 void EvalVectorArray::copy(const RCP<EvalVectorArray>& other)
00046 {
00047 resize(other->size());
00048
00049 for (int i=0; i<other->size(); i++)
00050 {
00051 (*this)[i]->copy((*other)[i]);
00052 }
00053 }
00054
00055 void EvalVectorArray::steal(const RCP<EvalVectorArray>& other)
00056 {
00057 resize(other->size());
00058
00059 for (int i=0; i<other->size(); i++)
00060 {
00061 (*this)[i] = (*other)[i];
00062 }
00063 }
00064
00065 ostream& EvalVectorArray::print(std::ostream& os,
00066 const SparsitySuperset* derivs) const
00067 {
00068 Tabs tab;
00069 TEUCHOS_TEST_FOR_EXCEPTION(derivs->numDerivs() != this->size(),
00070 std::logic_error,
00071 "mismatch between deriv set size=" << derivs->numDerivs()
00072 << "and result vector size " << this->size()
00073 << "in EvalVectorArray::print");
00074
00075 int maxlen = 25;
00076 for (int i=0; i<derivs->numDerivs(); i++)
00077 {
00078 int s = derivs->deriv(i).toString().length();
00079 if (s > maxlen) maxlen = s;
00080 }
00081
00082 for (int i=0; i<derivs->numDerivs(); i++)
00083 {
00084 os << tab;
00085 os.width(maxlen);
00086 os.setf(ios_base::left, ios_base::adjustfield);
00087 os << derivs->deriv(i).toString() << "\t\t";
00088 (*this)[i]->print(os);
00089 os << std::endl;
00090 }
00091 return os;
00092 }