51 #include "Teuchos_oblackholestream.hpp" 52 #include "Teuchos_GlobalMPISession.hpp" 63 typedef typename PV::size_type size_type;
65 const PV eb = Teuchos::dyn_cast<
const PV>(x);
68 for(size_type k=0; k<n; ++k) {
69 std::cout <<
"[subvector " << k <<
"]" << std::endl;
70 Teuchos::RCP<const V> vec = eb.get(k);
71 Teuchos::RCP<const std::vector<Real> > vp =
72 Teuchos::dyn_cast<
const SV>(*vec).getVector();
73 for(size_type i=0;i<vp->size();++i) {
74 std::cout << (*vp)[i] << std::endl;
80 int main(
int argc,
char *argv[]) {
82 using namespace Teuchos;
88 GlobalMPISession mpiSession(&argc, &argv);
90 int iprint = argc - 1;
92 RCP<std::ostream> outStream;
96 outStream = rcp(&std::cout,
false);
98 outStream = rcp(&bhs,
false);
102 RealT errtol = ROL::ROL_THRESHOLD<RealT>();
106 PV::size_type nvec = 3;
107 std::vector<int> dim;
115 RealT left = -1e0, right = 1e0;
117 std::vector<RCP<V> > x_rcp;
118 std::vector<RCP<V> > y_rcp;
119 std::vector<RCP<V> > z_rcp;
121 for( PV::size_type k=0; k<nvec; ++k ) {
122 RCP<std::vector<RealT> > xk_rcp = rcp(
new std::vector<RealT>(dim[k]) );
123 RCP<std::vector<RealT> > yk_rcp = rcp(
new std::vector<RealT>(dim[k]) );
124 RCP<std::vector<RealT> > zk_rcp = rcp(
new std::vector<RealT>(dim[k]) );
126 for(
int i=0; i<dim[k]; ++i ) {
127 (*xk_rcp)[i] = ( (
RealT)rand() / (
RealT)RAND_MAX ) * (right - left) + left;
128 (*yk_rcp)[i] = ( (
RealT)rand() / (
RealT)RAND_MAX ) * (right - left) + left;
129 (*zk_rcp)[i] = ( (
RealT)rand() / (
RealT)RAND_MAX ) * (right - left) + left;
132 RCP<V> xk = rcp(
new SV( xk_rcp ) );
133 RCP<V> yk = rcp(
new SV( yk_rcp ) );
134 RCP<V> zk = rcp(
new SV( zk_rcp ) );
144 RCP<V> y = ROL::CreatePartitionedVector<RealT>(y_rcp[0],y_rcp[1],y_rcp[2]);
148 std::vector<RealT> consistency = x.checkVector(*y, z,
true, *outStream);
150 if (checkvec.
norm() > std::sqrt(errtol)) {
156 Teuchos::RCP<ROL::Vector<RealT> > zp = x.clone();
158 RealT znorm = zp->norm();
159 *outStream <<
"Norm of ROL::Vector z (first basis vector): " << znorm <<
"\n";
160 if ( std::abs(znorm-1.0) > errtol ) {
161 *outStream <<
"---> POSSIBLE ERROR ABOVE!\n";
165 zp = x.basis(total_dim/2);
167 *outStream <<
"\nNorm of ROL::Vector z ('middle' basis vector): " << znorm <<
"\n";
168 if ( std::abs(znorm-1.0) > errtol ) {
169 *outStream <<
"---> POSSIBLE ERROR ABOVE!\n";
173 zp = x.basis(total_dim-1);
175 *outStream <<
"\nNorm of ROL::Vector z (last basis vector): " << znorm <<
"\n";
176 if ( std::abs(znorm-1.0) > errtol ) {
177 *outStream <<
"---> POSSIBLE ERROR ABOVE!\n";
183 consistency = x.checkVector(x, x,
true, *outStream);
184 if (checkvec.
norm() > 0.0) {
189 int m = atoi(argv[1]);
198 catch (std::logic_error err) {
199 *outStream << err.what() <<
"\n";
204 std::cout <<
"End Result: TEST FAILED\n";
206 std::cout <<
"End Result: TEST PASSED\n";
Defines the linear algebra of vector space on a generic partitioned vector.
Contains definitions of custom data types in ROL.
Defines the linear algebra or vector space interface.
void print_vector(const ROL::Vector< Real > &x)
int main(int argc, char *argv[])
size_type numVectors() const
Real norm() const
Returns where .