63 double locValue = Sacado::scalarValue(value);
64 double glbValue = 0.0;
68 Teuchos::reduceAll(*this->getComm(), Teuchos::REDUCE_MAX,
static_cast<Thyra::Ordinal
>(1), &locValue,&glbValue);
70 Teuchos::reduceAll(*this->getComm(), Teuchos::REDUCE_MIN,
static_cast<Thyra::Ordinal
>(1), &locValue,&glbValue);
74#ifdef PANZER_HAVE_EPETRA_STACK
76 if(this->useEpetra()) {
78 this->getEpetraVector()[0] = glbValue;
84 TEUCHOS_ASSERT(this->useThyra());
86 this->getThyraVector()[0] = glbValue;
94 using Teuchos::rcp_dynamic_cast;
96 Teuchos::RCP<Thyra::MultiVectorBase<double> > dgdx_unique = getDerivative();
98 uniqueContainer_ = linObjFactory_->buildLinearObjContainer();
99 Teuchos::rcp_dynamic_cast<ThyraObjContainer<double> >(uniqueContainer_)->set_x_th(dgdx_unique->col(0));
103 uniqueContainer_ = Teuchos::null;
111 using Teuchos::rcp_dynamic_cast;
113 Teuchos::RCP<Thyra::MultiVectorBase<double> > dgdx_unique = getDerivative();
115 uniqueContainer_ = linObjFactory_->buildLinearObjContainer();
116 Teuchos::rcp_dynamic_cast<ThyraObjContainer<double> >(uniqueContainer_)->set_x_th(dgdx_unique->col(0));
120 uniqueContainer_ = Teuchos::null;
128 const int n = value.size();
129 const int num_deriv = this->numDeriv();
130 TEUCHOS_ASSERT(n == 0 || n == num_deriv);
135 Teuchos::reduceAll(*this->getComm(), Teuchos::REDUCE_MAX, Thyra::Ordinal(1), &value.val(), &glbValue.val());
137 Teuchos::reduceAll(*this->getComm(), Teuchos::REDUCE_MIN, Thyra::Ordinal(1), &value.val(), &glbValue.val());
141 int locProc = value.val() == glbValue.val() ? this->getComm()->getRank() : this->getComm()->getSize();
143 Teuchos::reduceAll(*this->getComm(), Teuchos::REDUCE_MIN, Thyra::Ordinal(1), &locProc, &glbProc);
146 Teuchos::broadcast(*this->getComm(), glbProc, Thyra::Ordinal(n), &glbValue.fastAccessDx(0));
152#ifdef PANZER_HAVE_EPETRA_STACK
153 if(this->useEpetra()) {
156 for (
int i=0; i<num_deriv; ++i)
157 deriv[i][0] = glbValue.dx(i);
163 TEUCHOS_ASSERT(this->useThyra());
164 Thyra::ArrayRCP< Thyra::ArrayRCP<double> > deriv = this->getThyraMultiVector();
165 for (
int i=0; i<num_deriv; ++i)
166 deriv[i][0] = glbValue.dx(i);