AquiferConstantFlux: add restart serialization support

This commit is contained in:
Arne Morten Kvarving
2023-02-22 14:13:10 +01:00
parent 9418110d3e
commit da4784b314
3 changed files with 40 additions and 0 deletions

View File

@@ -52,6 +52,14 @@ public:
connection_flux_.resize(this->connections_.size(), {0});
}
static AquiferConstantFlux serializationTestObject(const Simulator& ebos_simulator)
{
AquiferConstantFlux<TypeTag> result({}, {}, ebos_simulator);
result.cumulative_flux_ = 1.0;
return result;
}
virtual ~AquiferConstantFlux() = default;
void updateAquifer(const SingleAquiferFlux& aquifer) {
@@ -113,6 +121,17 @@ public:
+= this->connection_flux_[idx] / model.dofTotalVolume(cellIdx);
}
template<class Serializer>
void serializeOp(Serializer& serializer)
{
serializer(cumulative_flux_);
}
bool operator==(const AquiferConstantFlux& rhs) const
{
return this->cumulative_flux_ == rhs.cumulative_flux_;
}
private:
const std::vector<Aquancon::AquancCell>& connections_;
SingleAquiferFlux aquifer_data_;

View File

@@ -244,12 +244,15 @@ serializeOp(Serializer& serializer)
auto* ct = dynamic_cast<AquiferCarterTracy<TypeTag>*>(aiPtr.get());
auto* fetp = dynamic_cast<AquiferFetkovich<TypeTag>*>(aiPtr.get());
auto* num = dynamic_cast<AquiferNumerical<TypeTag>*>(aiPtr.get());
auto* flux = dynamic_cast<AquiferConstantFlux<TypeTag>*>(aiPtr.get());
if (ct) {
serializer(*ct);
} else if (fetp) {
serializer(*fetp);
} else if (num) {
serializer(*num);
} else if (flux) {
serializer(*flux);
} else {
OPM_THROW(std::logic_error, "Error serializing BlackoilAquiferModel: unknown aquifer type");
}