39 using CellCenterSolutionVector =
typename std::decay_t<decltype(std::declval<SolutionVector>()[GridGeometry::cellCenterIdx()])>;
59 const std::string& name =
"newton_convergence")
60 : gridGeometry_(gridGeometry)
61 , ccWriter_(gridGeometry.gridView(), name,
"",
"")
63 , faceSequenceWriter_(faceWriter_, name +
"-face",
"",
"",
64 gridGeometry.gridView().comm().rank(),
65 gridGeometry.gridView().comm().size())
69 for (
int eqIdx = 0; eqIdx < numEqCellCenter; ++eqIdx)
71 ccWriter_.addCellData(xCellCenter_[eqIdx],
"x_" + std::to_string(eqIdx));
72 ccWriter_.addCellData(deltaCellCenter_[eqIdx],
"delta_" + std::to_string(eqIdx));
73 ccWriter_.addCellData(defCellCenter_[eqIdx],
"defect_" + std::to_string(eqIdx));
80 const auto numCellCenterDofs = gridGeometry_.numCellCenterDofs();
81 const auto numFaceDofs = gridGeometry_.numFaceDofs();
84 for (
int eqIdx = 0; eqIdx < numEqCellCenter; ++eqIdx)
86 defCellCenter_[eqIdx].resize(numCellCenterDofs);
87 deltaCellCenter_[eqIdx].resize(numCellCenterDofs);
88 xCellCenter_[eqIdx].resize(numCellCenterDofs);
92 for (
int eqIdx = 0; eqIdx < numEqFace; ++eqIdx)
94 defFace_[eqIdx].resize(numFaceDofs);
95 deltaFace_[eqIdx].resize(numFaceDofs);
96 xFace_[eqIdx].resize(numFaceDofs);
99 coordinates_.resize(numFaceDofs);
100 for (
auto&& facet : facets(gridGeometry_.gridView()))
102 const auto dofIdxGlobal = gridGeometry_.gridView().indexSet().index(facet);
103 coordinates_[dofIdxGlobal] = facet.geometry().center();
112 void write(
const SolutionVector& uLastIter,
113 const ResidualVector& deltaU,
114 const ResidualVector& residual)
override
116 assert(uLastIter.size() == deltaU.size() && uLastIter.size() == residual.size());
118 for (std::size_t dofIdxGlobal = 0; dofIdxGlobal < deltaU[GridGeometry::cellCenterIdx()].size(); ++dofIdxGlobal)
120 for (
int eqIdx = 0; eqIdx < numEqCellCenter; ++eqIdx)
122 xCellCenter_[eqIdx][dofIdxGlobal] = uLastIter[GridGeometry::cellCenterIdx()][dofIdxGlobal][eqIdx];
123 deltaCellCenter_[eqIdx][dofIdxGlobal] = - deltaU[GridGeometry::cellCenterIdx()][dofIdxGlobal][eqIdx];
124 defCellCenter_[eqIdx][dofIdxGlobal] = residual[GridGeometry::cellCenterIdx()][dofIdxGlobal][eqIdx];
128 for (
int eqIdx = 0; eqIdx < numEqFace; ++eqIdx)
130 faceWriter_->addPointData(xFace_[eqIdx],
"x_" + std::to_string(eqIdx));
131 faceWriter_->addPointData(deltaFace_[eqIdx],
"delta_" + std::to_string(eqIdx));
132 faceWriter_->addPointData(defFace_[eqIdx],
"defect_" + std::to_string(eqIdx));
135 for (std::size_t dofIdxGlobal = 0; dofIdxGlobal < deltaU[GridGeometry::faceIdx()].size(); ++dofIdxGlobal)
137 for (
int eqIdx = 0; eqIdx < numEqFace; ++eqIdx)
139 xFace_[eqIdx][dofIdxGlobal] = uLastIter[GridGeometry::faceIdx()][dofIdxGlobal][eqIdx];
140 deltaFace_[eqIdx][dofIdxGlobal] = - deltaU[GridGeometry::faceIdx()][dofIdxGlobal][eqIdx];
141 defFace_[eqIdx][dofIdxGlobal] = residual[GridGeometry::faceIdx()][dofIdxGlobal][eqIdx];
145 ccWriter_.write(
static_cast<double>(id_) +
static_cast<double>(iteration_)/1000);
146 faceSequenceWriter_.write(
static_cast<double>(id_) +
static_cast<double>(iteration_)/1000);