Commit 464c9580 authored by Markus Wolff's avatar Markus Wolff
Browse files

fixed bugs to due changes in impetproblem.hh

   - update vector is now resized in the update() function of the
     transport model



git-svn-id: svn://svn.iws.uni-stuttgart.de/DUMUX/dumux/trunk@7413 2fb0f335-1f38-0410-981e-8018bf24f1b0
parent cba112d9
......@@ -159,11 +159,16 @@ public:
/*! For an immiscible IMPES scheme, this is the saturation. For Miscible simulations, however,
* the total concentration of all components is transported.
*/
TransportSolutionType& transportedQuantity()
TransportSolutionType& transportedQuantity() DUNE_DEPRECATED
{
return totalConcentration_;
}
void getTransportedQuantity(TransportSolutionType& transportedQuantity)
{
transportedQuantity = totalConcentration_;
}
const TransportSolutionType& totalConcentration() const
DUNE_DEPRECATED
{
......@@ -222,7 +227,10 @@ void FVTransport2P2C<TypeTag>::update(const Scalar t, Scalar& dt, TransportSolut
// initialize dt very large
dt = 1E100;
// set update vector to zero
// resize update vector and set to zero
updateVec.resize(GET_PROP_VALUE(TypeTag, NumComponents));
updateVec[wCompIdx].resize(problem_.gridView().size(0));
updateVec[nCompIdx].resize(problem_.gridView().size(0));
updateVec[wCompIdx] = 0;
updateVec[nCompIdx] = 0;
......
......@@ -139,7 +139,10 @@ void FVTransport2P2CMultiPhysics<TypeTag>::update(const Scalar t, Scalar& dt, Tr
// initialize dt very large
dt = 1E100;
// set update vector to zero
// resize update vector and set to zero
updateVec.resize(GET_PROP_VALUE(TypeTag, NumComponents));
updateVec[wCompIdx].resize(problem().gridView().size(0));
updateVec[nCompIdx].resize(problem().gridView().size(0));
updateVec[wCompIdx] = 0;
updateVec[nCompIdx] = 0;
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment