ModelEvaluatorDefaultBaseTypes::DefaultDerivMvAdjointSupport
ModelEvaluatorDefaultBase<Scalar>::determineDefaultDerivMvAdjointSupport(
const ModelEvaluatorBase::DerivativeSupport &derivSupportImpl,
const VectorSpaceBase<Scalar> &fnc_space,
const VectorSpaceBase<Scalar> &var_space
)
{
typedef ModelEvaluatorBase MEB;
// Here we will support the adjoint copy for of a multi-vector if both
// spaces give rise to in-core vectors.
const bool implSupportsMv =
( derivSupportImpl.supports(MEB::DERIV_MV_BY_COL)
|| derivSupportImpl.supports(MEB::DERIV_TRANS_MV_BY_ROW) );
const bool bothSpacesHaveInCoreViews =
( fnc_space.hasInCoreView() && var_space.hasInCoreView() );
if ( implSupportsMv && bothSpacesHaveInCoreViews ) {
return DefaultDerivMvAdjointSupport(
derivSupportImpl.supports(MEB::DERIV_MV_BY_COL)
? MEB::DERIV_TRANS_MV_BY_ROW
: MEB::DERIV_MV_BY_COL
);
}
// We can't provide an adjoint copy or such a copy is not needed!
return DefaultDerivMvAdjointSupport();
}
bool VectorSpaceTester<Scalar>::check(
const VectorSpaceBase<Scalar> &vs
,Teuchos::FancyOStream *out_arg
) const
{
using std::endl;
using Teuchos::describe;
using Teuchos::FancyOStream;
using Teuchos::OSTab;
typedef Teuchos::ScalarTraits<Scalar> ST;
//typedef typename ST::magnitudeType ScalarMag;
Teuchos::RCP<FancyOStream> out = Teuchos::rcp(out_arg,false);
const Teuchos::EVerbosityLevel verbLevel = (dump_all()?Teuchos::VERB_EXTREME:Teuchos::VERB_MEDIUM);
OSTab tab(out,1,"THYRA");
bool result, success = true;
if(out.get()) *out <<endl<< "*** Entering Thyra::VectorSpaceTester<"<<ST::name()<<">::check(vs,...) ...\n";
if(out.get()) *out <<endl<< "Testing a vector space vs described as:\n" << describe(vs,verbLevel);
if(out.get())
*out
<<endl<< "A) Calling basic query functions ...\n"
<<endl<< "vs.dim() = " << vs.dim()
<<endl<< "vs.hasInCoreView() = " << vs.hasInCoreView() << std::endl;
if(out.get()) *out <<endl<< "B) Checking that vs is compatible with itself ...\n";
if(out.get()) *out <<endl<< "vs.isCompatible(vs)=";
result = vs.isCompatible(vs);
if(!result) success = false;
if(out.get()) *out << result << " == true : " << passfail(result) << std::endl;
if(out.get()) *out <<endl<< "C) Creating a randomized vector member v ...\n";
Teuchos::RCP<Thyra::VectorBase<Scalar> >
v = createMember(vs);
randomize(Scalar(-ST::one()),Scalar(+ST::one()),v.ptr());
if(out.get()) *out <<endl<< "D) Testing the VectorBase interface of v ...\n";
result = vectorTester_.check(*v,out.get());
if(!result) success = false;
if(out.get()) *out <<endl<< "C) Creating a randomized MultiVector member mv ...\n";
Teuchos::RCP<Thyra::MultiVectorBase<Scalar> >
mv = createMembers(vs,num_mv_cols());
randomize(Scalar(-ST::one()),Scalar(+ST::one()),mv.ptr());
if(out.get()) *out <<endl<< "D) Testing the MultiVectorBase interface of mv ...\n";
result = vectorTester_.multiVectorTester().check(*mv, out.ptr());
if(!result) success = false;
if(out.get()) {
if(success)
*out << endl <<"Congratulations, this VectorSpaceBase object seems to check out!\n";
else
*out << endl <<"Oh no, at least one of the tests performed with this VectorSpaceBase object failed (see above failures)!\n";
*out << endl << "*** Leaving VectorSpaceTester<"<<ST::name()<<">::check(vs,...)\n";
}
return success;
}
