57 #include "Teuchos_oblackholestream.hpp" 58 #include "Teuchos_XMLParameterListHelpers.hpp" 59 #include "Teuchos_GlobalMPISession.hpp" 60 #include "Teuchos_Comm.hpp" 61 #include "Teuchos_DefaultComm.hpp" 62 #include "Teuchos_CommHelpers.hpp" 77 int main(
int argc,
char *argv[]) {
79 Teuchos::GlobalMPISession mpiSession(&argc, &argv);
80 Teuchos::RCP<const Teuchos::Comm<int> > comm
81 = Teuchos::DefaultComm<int>::getComm();
84 int iprint = argc - 1;
85 bool print = (iprint>0);
86 Teuchos::RCP<std::ostream> outStream;
87 Teuchos::oblackholestream bhs;
89 outStream = Teuchos::rcp(&std::cout,
false);
91 outStream = Teuchos::rcp(&bhs,
false);
93 bool print0 = print && !comm->getRank();
94 Teuchos::RCP<std::ostream> outStream0;
96 outStream0 = Teuchos::rcp(&std::cout,
false);
98 outStream0 = Teuchos::rcp(&bhs,
false);
113 Teuchos::RCP<BurgersFEM<RealT> > fem
115 fem->test_inverse_mass(*outStream0);
116 fem->test_inverse_H1(*outStream0);
120 Teuchos::RCP<std::vector<RealT> > ud_rcp
121 = Teuchos::rcp(
new std::vector<RealT> (nx, 1.0) );
122 Teuchos::RCP<ROL::Vector<RealT> > ud
124 Teuchos::RCP<ROL::Objective_SimOpt<RealT> > pobj
130 Teuchos::RCP<ROL::EqualityConstraint_SimOpt<RealT> > pcon
136 Teuchos::RCP<std::vector<RealT> > z_rcp
137 = Teuchos::rcp(
new std::vector<RealT> (nx+2, 1.0) );
138 Teuchos::RCP<std::vector<RealT> > gz_rcp
139 = Teuchos::rcp(
new std::vector<RealT> (nx+2, 1.0) );
140 Teuchos::RCP<std::vector<RealT> > yz_rcp
141 = Teuchos::rcp(
new std::vector<RealT> (nx+2, 1.0) );
142 for (
int i=0; i<nx+2; i++) {
143 (*yz_rcp)[i] = 2.0*random<RealT>(comm)-1.0;
145 Teuchos::RCP<ROL::Vector<RealT> > zp
147 Teuchos::RCP<ROL::Vector<RealT> > gzp
149 Teuchos::RCP<ROL::Vector<RealT> > yzp
151 std::vector<RealT> zvar(1,0.0*random<RealT>(comm));
152 std::vector<RealT> gvar(1,random<RealT>(comm));
153 std::vector<RealT> yvar(1,random<RealT>(comm));
156 Teuchos::RCP<std::vector<RealT> > u_rcp
157 = Teuchos::rcp(
new std::vector<RealT> (nx, 1.0) );
158 Teuchos::RCP<std::vector<RealT> > gu_rcp
159 = Teuchos::rcp(
new std::vector<RealT> (nx, 1.0) );
160 Teuchos::RCP<ROL::Vector<RealT> > up
162 Teuchos::RCP<ROL::Vector<RealT> > gup
165 Teuchos::RCP<std::vector<RealT> > c_rcp
166 = Teuchos::rcp(
new std::vector<RealT> (nx, 1.0) );
167 Teuchos::RCP<std::vector<RealT> > l_rcp
168 = Teuchos::rcp(
new std::vector<RealT> (nx, 1.0) );
169 for (
int i=0; i<nx; i++) {
170 (*l_rcp)[i] = random<RealT>(comm);
172 Teuchos::RCP<ROL::Vector<RealT> > cp
174 Teuchos::RCP<ROL::Vector<RealT> > lp
179 int dim = 4, nSamp = 1000;
180 std::vector<RealT> tmp(2,0.0); tmp[0] = -1.0; tmp[1] = 1.0;
181 std::vector<std::vector<RealT> > bounds(dim,tmp);
182 Teuchos::RCP<ROL::BatchManager<RealT> > bman
184 Teuchos::RCP<ROL::SampleGenerator<RealT> > sampler
186 nSamp,bounds,bman,
false,
false,100));
190 bool storage =
true, fdhess =
false;
191 Teuchos::RCP<ROL::Objective<RealT> > robj
193 pobj,pcon,up,lp,gup,cp,storage,fdhess));
194 RealT order = 2.0, prob = 0.95;
195 Teuchos::RCP<ROL::Objective<RealT> > obj
197 robj,order,prob,sampler,storage));
202 bool derivcheck =
false;
204 int nranks = sampler->numBatches();
205 for (
int pid = 0; pid < nranks; pid++) {
206 if ( pid == sampler->batchID() ) {
207 for (
int i = sampler->start(); i < sampler->numMySamples(); i++) {
208 *outStream <<
"Sample " << i <<
" Rank " << sampler->batchID() <<
"\n";
209 *outStream <<
"(" << sampler->getMyPoint(i)[0] <<
", " 210 << sampler->getMyPoint(i)[1] <<
", " 211 << sampler->getMyPoint(i)[2] <<
", " 212 << sampler->getMyPoint(i)[3] <<
")\n";
213 pcon->setParameter(sampler->getMyPoint(i));
214 pcon->checkSolve(*up,*zp,*cp,print,*outStream);
215 robj->setParameter(sampler->getMyPoint(i));
217 robj->checkGradient(*zp,*gzp,*yzp,print,*outStream);
218 robj->checkHessVec(*zp,*gzp,*yzp,print,*outStream);
219 *outStream <<
"\n\n";
225 obj->checkGradient(z,g,y,print0,*outStream0);
226 obj->checkHessVec(z,g,y,print0,*outStream0);
231 std::string filename =
"input.xml";
232 Teuchos::RCP<Teuchos::ParameterList> parlist
233 = Teuchos::rcp(
new Teuchos::ParameterList() );
234 Teuchos::updateParametersFromXmlFile( filename, parlist.ptr() );
239 algo.run(z, g, *obj, print0, *outStream0);
244 for (
int i = 0; i < nx+2; i++ ) {
245 *outStream0 << std::scientific << std::setprecision(10);
246 *outStream0 << std::setw(20) << std::left << (
RealT)i/((
RealT)nx+1.0);
247 *outStream0 << std::setw(20) << std::left << (*z_rcp)[i];
251 *outStream0 <<
"Scalar Parameter: " << z.getStatistic(0) <<
"\n";
253 catch (std::logic_error err) {
254 *outStream << err.what() <<
"\n";
260 std::cout <<
"End Result: TEST FAILED\n";
262 std::cout <<
"End Result: TEST PASSED\n";
H1VectorDual< RealT > PrimalConstraintVector
int main(int argc, char *argv[])
H1VectorPrimal< RealT > DualConstraintVector
L2VectorDual< RealT > DualControlVector
Provides an interface to run optimization algorithms.
H1VectorDual< RealT > DualStateVector
L2VectorPrimal< RealT > PrimalControlVector
H1VectorPrimal< RealT > PrimalStateVector