Tempus  Version of the Day
Time Integration
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
Thyra_MultiVectorLinearOpWithSolveFactory.hpp
Go to the documentation of this file.
1 //@HEADER
2 // *****************************************************************************
3 // Tempus: Time Integration and Sensitivity Analysis Package
4 //
5 // Copyright 2017 NTESS and the Tempus contributors.
6 // SPDX-License-Identifier: BSD-3-Clause
7 // *****************************************************************************
8 //@HEADER
9 
10 #ifndef Thyra_MultiVectorLinearOpWithSolveFactory_hpp
11 #define Thyra_MultiVectorLinearOpWithSolveFactory_hpp
12 
13 #include "Thyra_LinearOpWithSolveFactoryBase.hpp"
14 #include "Thyra_MultiVectorLinearOp.hpp"
15 #include "Thyra_MultiVectorPreconditioner.hpp"
16 #include "Thyra_MultiVectorPreconditionerFactory.hpp"
17 #include "Thyra_DefaultMultiVectorLinearOpWithSolve.hpp"
18 #include "Thyra_DefaultMultiVectorProductVectorSpace.hpp"
19 #include "Thyra_DefaultLinearOpSource.hpp"
20 
21 namespace Thyra {
22 
25 template <class Scalar>
26 class MultiVectorLinearOpWithSolveFactory
27  : virtual public LinearOpWithSolveFactoryBase<Scalar> {
28  public:
31 
33  MultiVectorLinearOpWithSolveFactory() {}
34 
47  void nonconstInitialize(
48  const RCP<LinearOpWithSolveFactoryBase<Scalar> > &lowsf,
49  const RCP<const DefaultMultiVectorProductVectorSpace<Scalar> >
50  &multiVecRange,
51  const RCP<const DefaultMultiVectorProductVectorSpace<Scalar> >
52  &multiVecDomain);
53 
66  void initialize(const RCP<const LinearOpWithSolveFactoryBase<Scalar> > &lowsf,
67  const RCP<const DefaultMultiVectorProductVectorSpace<Scalar> >
68  &multiVecRange,
69  const RCP<const DefaultMultiVectorProductVectorSpace<Scalar> >
70  &multiVecDomain);
71 
72  RCP<LinearOpWithSolveFactoryBase<Scalar> > getUnderlyingLOWSF();
73 
74  RCP<const LinearOpWithSolveFactoryBase<Scalar> > getUnderlyingLOWSF() const;
75 
77 
80 
81  std::string description() const;
82 
84 
88 
89  void setParameterList(RCP<ParameterList> const &paramList);
90  RCP<ParameterList> getNonconstParameterList();
91  RCP<ParameterList> unsetParameterList();
92  RCP<const ParameterList> getParameterList() const;
93  RCP<const ParameterList> getValidParameters() const;
94 
96 
99 
101  virtual bool acceptsPreconditionerFactory() const;
102 
104  virtual void setPreconditionerFactory(
105  const RCP<PreconditionerFactoryBase<Scalar> > &precFactory,
106  const std::string &precFactoryName);
107 
109  virtual RCP<PreconditionerFactoryBase<Scalar> > getPreconditionerFactory()
110  const;
111 
113  virtual void unsetPreconditionerFactory(
114  RCP<PreconditionerFactoryBase<Scalar> > *precFactory,
115  std::string *precFactoryName);
116 
117  virtual bool isCompatible(const LinearOpSourceBase<Scalar> &fwdOpSrc) const;
118 
119  virtual RCP<LinearOpWithSolveBase<Scalar> > createOp() const;
120 
121  virtual void initializeOp(
122  const RCP<const LinearOpSourceBase<Scalar> > &fwdOpSrc,
123  LinearOpWithSolveBase<Scalar> *Op,
124  const ESupportSolveUse supportSolveUse) const;
125 
126  virtual void initializeAndReuseOp(
127  const RCP<const LinearOpSourceBase<Scalar> > &fwdOpSrc,
128  LinearOpWithSolveBase<Scalar> *Op) const;
129 
130  virtual void uninitializeOp(
131  LinearOpWithSolveBase<Scalar> *Op,
132  RCP<const LinearOpSourceBase<Scalar> > *fwdOpSrc,
133  RCP<const PreconditionerBase<Scalar> > *prec,
134  RCP<const LinearOpSourceBase<Scalar> > *approxFwdOpSrc,
135  ESupportSolveUse *supportSolveUse) const;
136 
137  virtual bool supportsPreconditionerInputType(
138  const EPreconditionerInputType precOpType) const;
139 
140  virtual void initializePreconditionedOp(
141  const RCP<const LinearOpSourceBase<Scalar> > &fwdOpSrc,
142  const RCP<const PreconditionerBase<Scalar> > &prec,
143  LinearOpWithSolveBase<Scalar> *Op,
144  const ESupportSolveUse supportSolveUse) const;
145 
146  virtual void initializeApproxPreconditionedOp(
147  const RCP<const LinearOpSourceBase<Scalar> > &fwdOpSrc,
148  const RCP<const LinearOpSourceBase<Scalar> > &approxFwdOpSrc,
149  LinearOpWithSolveBase<Scalar> *Op,
150  const ESupportSolveUse supportSolveUse) const;
151 
153 
154  protected:
157 
158  void informUpdatedVerbosityState() const;
159 
161 
162  private:
163  typedef Teuchos::ConstNonconstObjectContainer<
164  LinearOpWithSolveFactoryBase<Scalar> >
165  LOWSF_t;
166 
167  LOWSF_t lowsf_;
168  RCP<const DefaultMultiVectorProductVectorSpace<Scalar> > multiVecRange_;
169  RCP<const DefaultMultiVectorProductVectorSpace<Scalar> > multiVecDomain_;
170 };
171 
176 template <class Scalar>
177 RCP<MultiVectorLinearOpWithSolveFactory<Scalar> >
178 nonconstMultiVectorLinearOpWithSolveFactory(
179  const RCP<LinearOpWithSolveFactoryBase<Scalar> > &lowsf,
180  const RCP<const DefaultMultiVectorProductVectorSpace<Scalar> >
181  &multiVecRange,
182  const RCP<const DefaultMultiVectorProductVectorSpace<Scalar> >
183  &multiVecDomain)
184 {
185  RCP<MultiVectorLinearOpWithSolveFactory<Scalar> > mvlowsf =
186  Teuchos::rcp(new MultiVectorLinearOpWithSolveFactory<Scalar>);
187  mvlowsf->nonconstInitialize(lowsf, multiVecRange, multiVecDomain);
188  return mvlowsf;
189 }
190 
195 template <class Scalar>
196 RCP<MultiVectorLinearOpWithSolveFactory<Scalar> >
197 nonconstMultiVectorLinearOpWithSolveFactory(
198  const RCP<LinearOpWithSolveFactoryBase<Scalar> > &lowsf,
199  const int num_blocks)
200 {
201  RCP<LinearOpWithSolveBase<Scalar> > op = lowsf->createOp();
202  RCP<const Thyra::DefaultMultiVectorProductVectorSpace<Scalar> > mv_domain =
203  Thyra::multiVectorProductVectorSpace(op->domain(), num_blocks);
204  RCP<const Thyra::DefaultMultiVectorProductVectorSpace<Scalar> > mv_range =
205  Thyra::multiVectorProductVectorSpace(op->range(), num_blocks);
206  return nonconstMultiVectorLinearOpWithSolveFactory(lowsf, mv_range,
207  mv_domain);
208 }
209 
214 template <class Scalar>
215 RCP<MultiVectorLinearOpWithSolveFactory<Scalar> >
216 multiVectorLinearOpWithSolveFactory(
217  const RCP<const LinearOpWithSolveFactoryBase<Scalar> > &lowsf,
218  const RCP<const DefaultMultiVectorProductVectorSpace<Scalar> >
219  &multiVecRange,
220  const RCP<const DefaultMultiVectorProductVectorSpace<Scalar> >
221  &multiVecDomain)
222 {
223  RCP<MultiVectorLinearOpWithSolveFactory<Scalar> > mvlowsf =
224  Teuchos::rcp(new MultiVectorLinearOpWithSolveFactory<Scalar>);
225  mvlowsf->initialize(lowsf, multiVecRange, multiVecDomain);
226  return mvlowsf;
227 }
228 
233 template <class Scalar>
234 RCP<MultiVectorLinearOpWithSolveFactory<Scalar> >
235 multiVectorLinearOpWithSolveFactory(
236  const RCP<const LinearOpWithSolveFactoryBase<Scalar> > &lowsf,
237  const int num_blocks)
238 {
239  RCP<LinearOpWithSolveBase<Scalar> > op = lowsf->createOp();
240  RCP<const Thyra::DefaultMultiVectorProductVectorSpace<Scalar> > mv_domain =
241  Thyra::multiVectorProductVectorSpace(op->domain(), num_blocks);
242  RCP<const Thyra::DefaultMultiVectorProductVectorSpace<Scalar> > mv_range =
243  Thyra::multiVectorProductVectorSpace(op->range(), num_blocks);
244  return multiVectorLinearOpWithSolveFactory(lowsf, mv_range, mv_domain);
245 }
246 
247 // Overridden from Constructors/Initializers/Accessors
248 
249 template <class Scalar>
250 void MultiVectorLinearOpWithSolveFactory<Scalar>::nonconstInitialize(
251  const RCP<LinearOpWithSolveFactoryBase<Scalar> > &lowsf,
252  const RCP<const DefaultMultiVectorProductVectorSpace<Scalar> >
253  &multiVecRange,
254  const RCP<const DefaultMultiVectorProductVectorSpace<Scalar> >
255  &multiVecDomain)
256 {
257 #ifdef TEUCHOS_DEBUG
258  TEUCHOS_TEST_FOR_EXCEPT(is_null(lowsf));
259 #endif
260  lowsf_.initialize(lowsf);
261  multiVecRange_ = multiVecRange;
262  multiVecDomain_ = multiVecDomain;
263 }
264 
265 template <class Scalar>
266 void MultiVectorLinearOpWithSolveFactory<Scalar>::initialize(
267  const RCP<const LinearOpWithSolveFactoryBase<Scalar> > &lowsf,
268  const RCP<const DefaultMultiVectorProductVectorSpace<Scalar> >
269  &multiVecRange,
270  const RCP<const DefaultMultiVectorProductVectorSpace<Scalar> >
271  &multiVecDomain)
272 {
273 #ifdef TEUCHOS_DEBUG
274  TEUCHOS_TEST_FOR_EXCEPT(is_null(lowsf));
275 #endif
276  lowsf_.initialize(lowsf);
277  multiVecRange_ = multiVecRange;
278  multiVecDomain_ = multiVecDomain;
279 }
280 
281 template <class Scalar>
282 RCP<LinearOpWithSolveFactoryBase<Scalar> >
283 MultiVectorLinearOpWithSolveFactory<Scalar>::getUnderlyingLOWSF()
284 {
285  return lowsf_.getNonconstObj();
286 }
287 
288 template <class Scalar>
289 RCP<const LinearOpWithSolveFactoryBase<Scalar> >
290 MultiVectorLinearOpWithSolveFactory<Scalar>::getUnderlyingLOWSF() const
291 {
292  return lowsf_.getConstObj();
293 }
294 
295 // Overridden from Teuchos::Describable
296 
297 template <class Scalar>
298 std::string MultiVectorLinearOpWithSolveFactory<Scalar>::description() const
299 {
300  std::ostringstream oss;
301  oss << this->Teuchos::Describable::description() << "{"
302  << "lowsf=";
303  if (!is_null(lowsf_.getConstObj()))
304  oss << lowsf_.getConstObj()->description();
305  else
306  oss << "NULL";
307  oss << "}";
308  return oss.str();
309 }
310 
311 // Overridden from ParameterListAcceptor
312 
313 template <class Scalar>
314 void MultiVectorLinearOpWithSolveFactory<Scalar>::setParameterList(
315  RCP<ParameterList> const &paramList)
316 {
317  lowsf_.getNonconstObj()->setParameterList(paramList);
318 }
319 
320 template <class Scalar>
321 RCP<ParameterList>
322 MultiVectorLinearOpWithSolveFactory<Scalar>::getNonconstParameterList()
323 {
324  return lowsf_.getNonconstObj()->getNonconstParameterList();
325 }
326 
327 template <class Scalar>
328 RCP<ParameterList>
329 MultiVectorLinearOpWithSolveFactory<Scalar>::unsetParameterList()
330 {
331  return lowsf_.getNonconstObj()->unsetParameterList();
332 }
333 
334 template <class Scalar>
335 RCP<const ParameterList>
336 MultiVectorLinearOpWithSolveFactory<Scalar>::getParameterList() const
337 {
338  return lowsf_.getConstObj()->getParameterList();
339 }
340 
341 template <class Scalar>
342 RCP<const ParameterList>
343 MultiVectorLinearOpWithSolveFactory<Scalar>::getValidParameters() const
344 {
345  return lowsf_.getConstObj()->getValidParameters();
346 }
347 
348 // Overridden from LinearOpWithSolveFactoyBase
349 
350 template <class Scalar>
351 bool MultiVectorLinearOpWithSolveFactory<Scalar>::acceptsPreconditionerFactory()
352  const
353 {
354  return lowsf_.getConstObj()->acceptsPreconditionerFactory();
355 }
356 
357 template <class Scalar>
358 void MultiVectorLinearOpWithSolveFactory<Scalar>::setPreconditionerFactory(
359  const RCP<PreconditionerFactoryBase<Scalar> > &precFactory,
360  const std::string &precFactoryName)
361 {
362  typedef MultiVectorPreconditionerFactory<Scalar> MVPF;
363  RCP<MVPF> mvpf = Teuchos::rcp_dynamic_cast<MVPF>(precFactory);
364  lowsf_.getNonconstObj()->setPreconditionerFactory(
365  mvpf->getNonconstPreconditionerFactory(), precFactoryName);
366 }
367 
368 template <class Scalar>
369 RCP<PreconditionerFactoryBase<Scalar> >
370 MultiVectorLinearOpWithSolveFactory<Scalar>::getPreconditionerFactory() const
371 {
372  RCP<PreconditionerFactoryBase<Scalar> > prec_fac =
373  lowsf_.getConstObj()->getPreconditionerFactory();
374  if (prec_fac == Teuchos::null)
375  return Teuchos::null;
376  else
377  return nonconstMultiVectorPreconditionerFactory(prec_fac, multiVecRange_,
378  multiVecDomain_);
379 }
380 
381 template <class Scalar>
382 void MultiVectorLinearOpWithSolveFactory<Scalar>::unsetPreconditionerFactory(
383  RCP<PreconditionerFactoryBase<Scalar> > *precFactory,
384  std::string *precFactoryName)
385 {
386  RCP<PreconditionerFactoryBase<Scalar> > inner_precFactory;
387  lowsf_.getNonconstObj()->unsetPreconditionerFactory(
388  precFactory ? &inner_precFactory : NULL, precFactoryName);
389  if (precFactory)
390  *precFactory = nonconstMultiVectorPreconditionerFactory(
391  inner_precFactory, multiVecRange_, multiVecDomain_);
392 }
393 
394 template <class Scalar>
395 bool MultiVectorLinearOpWithSolveFactory<Scalar>::isCompatible(
396  const LinearOpSourceBase<Scalar> &fwdOpSrc) const
397 {
398  typedef MultiVectorLinearOp<Scalar> MVLO;
399  RCP<const MVLO> mvlo =
400  Teuchos::rcp_dynamic_cast<const MVLO>(fwdOpSrc.getOp().assert_not_null());
401  RCP<const LinearOpSourceBase<Scalar> > inner_fwdOpSrc =
402  defaultLinearOpSource<Scalar>(mvlo->getLinearOp());
403  return lowsf_.getConstObj()->isCompatible(*inner_fwdOpSrc);
404 }
405 
406 template <class Scalar>
407 RCP<LinearOpWithSolveBase<Scalar> >
408 MultiVectorLinearOpWithSolveFactory<Scalar>::createOp() const
409 {
410  return nonconstMultiVectorLinearOpWithSolve<Scalar>(
411  lowsf_.getConstObj()->createOp(), multiVecRange_, multiVecDomain_);
412 }
413 
414 template <class Scalar>
415 void MultiVectorLinearOpWithSolveFactory<Scalar>::initializeOp(
416  const RCP<const LinearOpSourceBase<Scalar> > &fwdOpSrc,
417  LinearOpWithSolveBase<Scalar> *Op,
418  const ESupportSolveUse supportSolveUse) const
419 {
420  using Teuchos::dyn_cast;
421  using Teuchos::rcp_dynamic_cast;
422 
423 #ifdef TEUCHOS_DEBUG
424  TEUCHOS_TEST_FOR_EXCEPT(0 == Op);
425 #endif
426 
427  // Set the verbosity settings for the wrapped LOWSF object!
428  lowsf_.getConstObj()->setOStream(this->getOStream());
429  lowsf_.getConstObj()->setVerbLevel(this->getVerbLevel());
430 
431  typedef MultiVectorLinearOp<Scalar> MVLO;
432  typedef DefaultMultiVectorLinearOpWithSolve<Scalar> MVLOWS;
433  const RCP<const MVLO> mvlo =
434  rcp_dynamic_cast<const MVLO>(fwdOpSrc->getOp().assert_not_null());
435  MVLOWS &mvlows = dyn_cast<MVLOWS>(*Op);
436 
437  lowsf_.getConstObj()->initializeOp(
438  defaultLinearOpSource<Scalar>(mvlo->getLinearOp()),
439  mvlows.getNonconstLinearOpWithSolve().get(), supportSolveUse);
440 }
441 
442 template <class Scalar>
443 void MultiVectorLinearOpWithSolveFactory<Scalar>::initializeAndReuseOp(
444  const RCP<const LinearOpSourceBase<Scalar> > &fwdOpSrc,
445  LinearOpWithSolveBase<Scalar> *Op) const
446 {
447  using Teuchos::dyn_cast;
448  using Teuchos::rcp_dynamic_cast;
449 
450 #ifdef TEUCHOS_DEBUG
451  TEUCHOS_TEST_FOR_EXCEPT(0 == Op);
452 #endif
453 
454  // Set the verbosity settings for the wrapped LOWSF object!
455  lowsf_.getConstObj()->setOStream(this->getOStream());
456  lowsf_.getConstObj()->setVerbLevel(this->getVerbLevel());
457 
458  typedef MultiVectorLinearOp<Scalar> MVLO;
459  typedef DefaultMultiVectorLinearOpWithSolve<Scalar> MVLOWS;
460  const RCP<const MVLO> mvlo =
461  rcp_dynamic_cast<const MVLO>(fwdOpSrc->getOp().assert_not_null());
462  MVLOWS &mvlows = dyn_cast<MVLOWS>(*Op);
463 
464  lowsf_.getConstObj()->initializeAndReuseOp(
465  defaultLinearOpSource<Scalar>(mvlo->getLinearOp()),
466  mvlows.getNonconstLinearOpWithSolve().get());
467 }
468 
469 template <class Scalar>
470 void MultiVectorLinearOpWithSolveFactory<Scalar>::uninitializeOp(
471  LinearOpWithSolveBase<Scalar> *Op,
472  RCP<const LinearOpSourceBase<Scalar> > *fwdOpSrc,
473  RCP<const PreconditionerBase<Scalar> > *prec,
474  RCP<const LinearOpSourceBase<Scalar> > *approxFwdOpSrc,
475  ESupportSolveUse *supportSolveUse) const
476 {
477  using Teuchos::dyn_cast;
478 
479 #ifdef TEUCHOS_DEBUG
480  TEUCHOS_TEST_FOR_EXCEPT(0 == Op);
481 #endif
482  typedef DefaultMultiVectorLinearOpWithSolve<Scalar> MVLOWS;
483  MVLOWS &mvlowsOp = dyn_cast<MVLOWS>(*Op);
484  RCP<const LinearOpSourceBase<Scalar> > inner_fwdOpSrc;
485  RCP<const PreconditionerBase<Scalar> > inner_prec;
486  RCP<const LinearOpSourceBase<Scalar> > inner_approxFwdOpSrc;
487  lowsf_.getConstObj()->uninitializeOp(
488  mvlowsOp.getNonconstLinearOpWithSolve().get(),
489  fwdOpSrc ? &inner_fwdOpSrc : NULL, prec ? &inner_prec : NULL,
490  approxFwdOpSrc ? &inner_approxFwdOpSrc : NULL, supportSolveUse);
491  if (fwdOpSrc)
492  *fwdOpSrc = defaultLinearOpSource<Scalar>(multiVectorLinearOp(
493  inner_fwdOpSrc->getOp(), multiVecRange_, multiVecDomain_));
494  if (prec)
495  *prec =
496  multiVectorPreconditioner(inner_prec, multiVecRange_, multiVecDomain_);
497  if (fwdOpSrc)
498  *approxFwdOpSrc = defaultLinearOpSource<Scalar>(multiVectorLinearOp(
499  inner_approxFwdOpSrc->getOp(), multiVecRange_, multiVecDomain_));
500 }
501 
502 template <class Scalar>
503 bool MultiVectorLinearOpWithSolveFactory<
504  Scalar>::supportsPreconditionerInputType(const EPreconditionerInputType
505  precOpType) const
506 {
507  return lowsf_.getConstObj()->supportsPreconditionerInputType(precOpType);
508 }
509 
510 template <class Scalar>
511 void MultiVectorLinearOpWithSolveFactory<Scalar>::initializePreconditionedOp(
512  const RCP<const LinearOpSourceBase<Scalar> > &fwdOpSrc,
513  const RCP<const PreconditionerBase<Scalar> > &prec,
514  LinearOpWithSolveBase<Scalar> *Op,
515  const ESupportSolveUse supportSolveUse) const
516 {
517  using Teuchos::dyn_cast;
518  using Teuchos::rcp_dynamic_cast;
519 
520 #ifdef TEUCHOS_DEBUG
521  TEUCHOS_TEST_FOR_EXCEPT(0 == Op);
522 #endif
523 
524  // Set the verbosity settings for the wrapped LOWSF object!
525  lowsf_.getConstObj()->setOStream(this->getOStream());
526  lowsf_.getConstObj()->setVerbLevel(this->getVerbLevel());
527 
528  typedef MultiVectorLinearOp<Scalar> MVLO;
529  typedef MultiVectorPreconditioner<Scalar> MVP;
530  typedef DefaultMultiVectorLinearOpWithSolve<Scalar> MVLOWS;
531  const RCP<const MVLO> mvlo =
532  rcp_dynamic_cast<const MVLO>(fwdOpSrc->getOp().assert_not_null());
533  const RCP<const MVP> mvp = rcp_dynamic_cast<const MVP>(prec);
534  MVLOWS &mvlows = dyn_cast<MVLOWS>(*Op);
535 
536  lowsf_.getConstObj()->initializePreconditionedOp(
537  defaultLinearOpSource<Scalar>(mvlo->getLinearOp()),
538  mvp->getPreconditioner(), mvlows.getNonconstLinearOpWithSolve().get(),
539  supportSolveUse);
540 }
541 
542 template <class Scalar>
543 void MultiVectorLinearOpWithSolveFactory<Scalar>::
544  initializeApproxPreconditionedOp(
545  const RCP<const LinearOpSourceBase<Scalar> > &fwdOpSrc,
546  const RCP<const LinearOpSourceBase<Scalar> > &approxFwdOpSrc,
547  LinearOpWithSolveBase<Scalar> *Op,
548  const ESupportSolveUse supportSolveUse) const
549 {
550  using Teuchos::dyn_cast;
551  using Teuchos::rcp_dynamic_cast;
552 
553 #ifdef TEUCHOS_DEBUG
554  TEUCHOS_TEST_FOR_EXCEPT(0 == Op);
555 #endif
556 
557  // Set the verbosity settings for the wrapped LOWSF object!
558  lowsf_.getConstObj()->setOStream(this->getOStream());
559  lowsf_.getConstObj()->setVerbLevel(this->getVerbLevel());
560 
561  typedef MultiVectorLinearOp<Scalar> MVLO;
562  typedef DefaultMultiVectorLinearOpWithSolve<Scalar> MVLOWS;
563  const RCP<const MVLO> mvlo =
564  rcp_dynamic_cast<const MVLO>(fwdOpSrc->getOp().assert_not_null());
565  const RCP<const MVLO> amvlo =
566  rcp_dynamic_cast<const MVLO>(approxFwdOpSrc->getOp().assert_not_null());
567  MVLOWS &mvlows = dyn_cast<MVLOWS>(*Op);
568 
569  lowsf_.getConstObj()->initializeApproxPreconditionedOp(
570  defaultLinearOpSource<Scalar>(mvlo->getLinearOp()),
571  defaultLinearOpSource<Scalar>(amvlo->getLinearOp()),
572  mvlows.getNonconstLinearOpWithSolve().get(), supportSolveUse);
573 }
574 
575 // protected
576 
577 template <class Scalar>
578 void MultiVectorLinearOpWithSolveFactory<Scalar>::informUpdatedVerbosityState()
579  const
580 {
581  lowsf_.getConstObj()->setVerbLevel(this->getVerbLevel());
582  lowsf_.getConstObj()->setOStream(this->getOStream());
583 }
584 
585 } // namespace Thyra
586 
587 #endif
Thyra::MultiVectorLinearOpWithSolveFactory::initializePreconditionedOp
virtual void initializePreconditionedOp(const RCP< const LinearOpSourceBase< Scalar > > &fwdOpSrc, const RCP< const PreconditionerBase< Scalar > > &prec, LinearOpWithSolveBase< Scalar > *Op, const ESupportSolveUse supportSolveUse) const
Definition: Thyra_MultiVectorLinearOpWithSolveFactory.hpp:511
Thyra::MultiVectorPreconditionerFactory
Concrete PreconditionerFactoryBase subclass that wraps a preconditioner in MultiVectorPreconditioner...
Definition: Thyra_MultiVectorPreconditionerFactory.hpp:25
Thyra::MultiVectorLinearOpWithSolveFactory::multiVecRange_
RCP< const DefaultMultiVectorProductVectorSpace< Scalar > > multiVecRange_
Definition: Thyra_MultiVectorLinearOpWithSolveFactory.hpp:168
is_null
bool is_null(const boost::shared_ptr< T > &p)
Thyra::MultiVectorLinearOpWithSolveFactory::multiVecDomain_
RCP< const DefaultMultiVectorProductVectorSpace< Scalar > > multiVecDomain_
Definition: Thyra_MultiVectorLinearOpWithSolveFactory.hpp:169
Thyra::MultiVectorLinearOpWithSolveFactory::supportsPreconditionerInputType
virtual bool supportsPreconditionerInputType(const EPreconditionerInputType precOpType) const
Definition: Thyra_MultiVectorLinearOpWithSolveFactory.hpp:504
Thyra::MultiVectorLinearOpWithSolveFactory::initializeOp
virtual void initializeOp(const RCP< const LinearOpSourceBase< Scalar > > &fwdOpSrc, LinearOpWithSolveBase< Scalar > *Op, const ESupportSolveUse supportSolveUse) const
Definition: Thyra_MultiVectorLinearOpWithSolveFactory.hpp:415
Thyra::MultiVectorLinearOpWithSolveFactory::multiVectorLinearOpWithSolveFactory
RCP< MultiVectorLinearOpWithSolveFactory< Scalar > > multiVectorLinearOpWithSolveFactory(const RCP< const LinearOpWithSolveFactoryBase< Scalar > > &lowsf, const RCP< const DefaultMultiVectorProductVectorSpace< Scalar > > &multiVecRange, const RCP< const DefaultMultiVectorProductVectorSpace< Scalar > > &multiVecDomain)
Nonmember constructor.
Definition: Thyra_MultiVectorLinearOpWithSolveFactory.hpp:216
Teuchos::dyn_cast
T_To & dyn_cast(T_From &from)
Thyra::MultiVectorLinearOpWithSolveFactory::createOp
virtual RCP< LinearOpWithSolveBase< Scalar > > createOp() const
Definition: Thyra_MultiVectorLinearOpWithSolveFactory.hpp:408
Thyra::MultiVectorLinearOpWithSolveFactory
Create a LinearOpWithSolveFactory for a flattened-out multi-vector.
Definition: Thyra_MultiVectorLinearOpWithSolveFactory.hpp:26
Thyra_MultiVectorLinearOp.hpp
Thyra::MultiVectorLinearOpWithSolveFactory::getUnderlyingLOWSF
RCP< LinearOpWithSolveFactoryBase< Scalar > > getUnderlyingLOWSF()
Definition: Thyra_MultiVectorLinearOpWithSolveFactory.hpp:283
Thyra::MultiVectorLinearOpWithSolveFactory::setPreconditionerFactory
virtual void setPreconditionerFactory(const RCP< PreconditionerFactoryBase< Scalar > > &precFactory, const std::string &precFactoryName)
Throws exception.
Definition: Thyra_MultiVectorLinearOpWithSolveFactory.hpp:358
Thyra::MultiVectorLinearOpWithSolveFactory::nonconstMultiVectorLinearOpWithSolveFactory
RCP< MultiVectorLinearOpWithSolveFactory< Scalar > > nonconstMultiVectorLinearOpWithSolveFactory(const RCP< LinearOpWithSolveFactoryBase< Scalar > > &lowsf, const RCP< const DefaultMultiVectorProductVectorSpace< Scalar > > &multiVecRange, const RCP< const DefaultMultiVectorProductVectorSpace< Scalar > > &multiVecDomain)
Nonmember constructor.
Definition: Thyra_MultiVectorLinearOpWithSolveFactory.hpp:178
Thyra::MultiVectorLinearOpWithSolveFactory::getPreconditionerFactory
virtual RCP< PreconditionerFactoryBase< Scalar > > getPreconditionerFactory() const
Returns null .
Definition: Thyra_MultiVectorLinearOpWithSolveFactory.hpp:370
Teuchos::rcp
TEUCHOS_DEPRECATED RCP< T > rcp(T *p, Dealloc_T dealloc, bool owns_mem)
Thyra::MultiVectorLinearOpWithSolveFactory::nonconstMultiVectorLinearOpWithSolveFactory
RCP< MultiVectorLinearOpWithSolveFactory< Scalar > > nonconstMultiVectorLinearOpWithSolveFactory(const RCP< LinearOpWithSolveFactoryBase< Scalar > > &lowsf, const int num_blocks)
Nonmember constructor.
Definition: Thyra_MultiVectorLinearOpWithSolveFactory.hpp:197
Thyra::MultiVectorLinearOpWithSolveFactory::nonconstInitialize
void nonconstInitialize(const RCP< LinearOpWithSolveFactoryBase< Scalar > > &lowsf, const RCP< const DefaultMultiVectorProductVectorSpace< Scalar > > &multiVecRange, const RCP< const DefaultMultiVectorProductVectorSpace< Scalar > > &multiVecDomain)
Initialize given a single non-const LOWSFB object.
Definition: Thyra_MultiVectorLinearOpWithSolveFactory.hpp:250
Teuchos::Describable::description
virtual std::string description() const
Thyra::MultiVectorLinearOpWithSolveFactory::initializeAndReuseOp
virtual void initializeAndReuseOp(const RCP< const LinearOpSourceBase< Scalar > > &fwdOpSrc, LinearOpWithSolveBase< Scalar > *Op) const
Definition: Thyra_MultiVectorLinearOpWithSolveFactory.hpp:443
ESupportSolveUse
ESupportSolveUse
Thyra::MultiVectorLinearOpWithSolveFactory::initialize
void initialize(const RCP< const LinearOpWithSolveFactoryBase< Scalar > > &lowsf, const RCP< const DefaultMultiVectorProductVectorSpace< Scalar > > &multiVecRange, const RCP< const DefaultMultiVectorProductVectorSpace< Scalar > > &multiVecDomain)
Initialize given a single const LOWSFB object.
Definition: Thyra_MultiVectorLinearOpWithSolveFactory.hpp:266
Thyra::MultiVectorLinearOpWithSolveFactory::unsetPreconditionerFactory
virtual void unsetPreconditionerFactory(RCP< PreconditionerFactoryBase< Scalar > > *precFactory, std::string *precFactoryName)
Throws exception.
Definition: Thyra_MultiVectorLinearOpWithSolveFactory.hpp:382
Thyra_MultiVectorPreconditioner.hpp
Thyra::MultiVectorLinearOp
Implicit concrete LinearOpBase subclass that takes a flattended out multi-vector and performs a multi...
Definition: Thyra_MultiVectorLinearOp.hpp:30
Thyra::MultiVectorPreconditioner
Concrete PreconditionerBase subclass that wraps a preconditioner operator in MultiVectorLinearOp.
Definition: Thyra_MultiVectorPreconditioner.hpp:24
Thyra::MultiVectorLinearOpWithSolveFactory::isCompatible
virtual bool isCompatible(const LinearOpSourceBase< Scalar > &fwdOpSrc) const
Definition: Thyra_MultiVectorLinearOpWithSolveFactory.hpp:395
Thyra_MultiVectorPreconditionerFactory.hpp
Thyra::MultiVectorLinearOpWithSolveFactory::uninitializeOp
virtual void uninitializeOp(LinearOpWithSolveBase< Scalar > *Op, RCP< const LinearOpSourceBase< Scalar > > *fwdOpSrc, RCP< const PreconditionerBase< Scalar > > *prec, RCP< const LinearOpSourceBase< Scalar > > *approxFwdOpSrc, ESupportSolveUse *supportSolveUse) const
Definition: Thyra_MultiVectorLinearOpWithSolveFactory.hpp:470
EPreconditionerInputType
EPreconditionerInputType
Thyra::MultiVectorLinearOpWithSolveFactory::initializeApproxPreconditionedOp
virtual void initializeApproxPreconditionedOp(const RCP< const LinearOpSourceBase< Scalar > > &fwdOpSrc, const RCP< const LinearOpSourceBase< Scalar > > &approxFwdOpSrc, LinearOpWithSolveBase< Scalar > *Op, const ESupportSolveUse supportSolveUse) const
Definition: Thyra_MultiVectorLinearOpWithSolveFactory.hpp:544
Thyra::MultiVectorLinearOpWithSolveFactory::multiVectorLinearOpWithSolveFactory
RCP< MultiVectorLinearOpWithSolveFactory< Scalar > > multiVectorLinearOpWithSolveFactory(const RCP< const LinearOpWithSolveFactoryBase< Scalar > > &lowsf, const int num_blocks)
Nonmember constructor.
Definition: Thyra_MultiVectorLinearOpWithSolveFactory.hpp:235
Thyra::LinearOpSourceBase::getOp
virtual Teuchos::RCP< const LinearOpBase< Scalar > > getOp() const =0
Thyra::MultiVectorLinearOpWithSolveFactory::LOWSF_t
Teuchos::ConstNonconstObjectContainer< LinearOpWithSolveFactoryBase< Scalar > > LOWSF_t
Definition: Thyra_MultiVectorLinearOpWithSolveFactory.hpp:165
TEUCHOS_TEST_FOR_EXCEPT
#define TEUCHOS_TEST_FOR_EXCEPT(throw_exception_test)
Thyra::MultiVectorLinearOpWithSolveFactory::informUpdatedVerbosityState
void informUpdatedVerbosityState() const
Overridden from Teuchos::VerboseObjectBase.
Definition: Thyra_MultiVectorLinearOpWithSolveFactory.hpp:578
Generated on Fri Nov 22 2024 09:23:27 for Tempus by   doxygen 1.8.5