82 const std::vector<Teuchos::RCP<const GlobalIndexer>>&
84 const Teuchos::ParameterList& p)
87 useTimeDerivativeSolutionVector_(false),
88 globalDataKey_(
"Tangent Gather Container")
99 const vector<string>& names = *(p.get<RCP<vector<string>>>(
"DOF Names"));
101 RCP<PureBasis> basis = p.get<RCP<PureBasis>>(
"Basis");
102 if (p.isType<
bool>(
"Use Time Derivative Solution Vector"))
104 p.get<
bool>(
"Use Time Derivative Solution Vector");
105 if (p.isType<
string>(
"Global Data Key"))
114 MDField<ScalarT, Cell, NODE>(names[fd], basis->functional);
119 string firstName(
"<none>");
121 firstName = names[0];
122 string n(
"GatherTangent (Blocked Epetra): " + firstName +
" (" +
123 print<EvalT>() +
")");
191 typename TRAITS::EvalData workset)
197 using Teuchos::ArrayRCP;
198 using Teuchos::ptrFromRef;
200 using Teuchos::rcp_dynamic_cast;
202 using Thyra::SpmdVectorBase;
206 if (xBvRoGed_.is_null())
210 string blockId(this->wda(workset).block_id);
211 const vector<size_t>& localCellIds = this->wda(workset).cell_local_ids;
212 int numFields(gatherFields_.size()), numCells(localCellIds.size());
215 for (
int fieldIndex(0); fieldIndex <
numFields; ++fieldIndex)
217 MDField<ScalarT, Cell, NODE>&
field = gatherFields_[fieldIndex];
218 int indexerId(indexerIds_[fieldIndex]),
219 subFieldNum(subFieldIds_[fieldIndex]);
222 auto xEvRoGed = xBvRoGed_->getGEDBlock(indexerId);
223 auto subRowIndexer = indexers_[indexerId];
224 const vector<int>& elmtOffset =
225 subRowIndexer->getGIDFieldOffsets(blockId, subFieldNum);
226 int numBases(elmtOffset.size());
229 for (
int cell(0); cell < numCells; ++cell)
231 LO cellLocalId = localCellIds[cell];
232 auto LIDs = subRowIndexer->getElementLIDs(cellLocalId);
235 for (
int basis(0); basis < numBases; ++basis)
237 int offset(elmtOffset[basis]), lid(LIDs[offset]);
238 field(cell, basis) = (*xEvRoGed)[lid];
PHX::MDField< ScalarT, panzer::Cell, panzer::BASIS > field
A field to which we'll contribute, or in which we'll store, the result of computing this integral.