Commit 8659721c authored by Benjamin Faigle's avatar Benjamin Faigle
Browse files

removed potentials from dec2p2c because it is not used and access to the...

removed potentials from dec2p2c because it is not used and access to the vectro incredibly slow for large domains

git-svn-id: svn://svn.iws.uni-stuttgart.de/DUMUX/dumux/trunk@5125 2fb0f335-1f38-0410-981e-8018bf24f1b0
parent 73ee6550
......@@ -111,16 +111,6 @@ public:
template<class MultiWriter>
void addOutputVtkFields(MultiWriter &writer)
{
#if DUNE_MINIMAL_DEBUG_LEVEL <= 3
// add debug stuff
Dune::BlockVector<Dune::FieldVector<double,1> > *potentialDifferenceW = writer.template createField<double, 1> (problem_.variables().gridSize());
*potentialDifferenceW = potential_[wPhaseIdx];
writer.addCellData(potentialDifferenceW, "potentialDifferenceW pre/post pressure");
Dune::BlockVector<Dune::FieldVector<double,1> > *potentialDifferenceNW = writer.template createField<double, 1> (problem_.variables().gridSize());
*potentialDifferenceNW = potential_[nPhaseIdx];
writer.addCellData(potentialDifferenceNW, "potentialDifferenceNW pre/post pressure");
#endif
return;
}
......@@ -154,10 +144,6 @@ public:
{
DUNE_THROW(Dune::NotImplemented, "Velocity type not supported!");
}
potential_.resize(2);
potential_[0].resize(problem_.variables().gridSize());
potential_[1].resize(problem_.variables().gridSize());
}
~FVTransport2P2C()
......@@ -167,8 +153,6 @@ public:
private:
Problem& problem_;
TransportSolutionType potential_;
static const int pressureType = GET_PROP_VALUE(TypeTag, PTAG(PressureFormulation)); //!< gives kind of pressure used (\f$ 0 = p_w \f$, \f$ 1 = p_n \f$, \f$ 2 = p_{global} \f$)
bool switchNormals;
};
......
......@@ -329,31 +329,6 @@ public:
return saturation_[Idx][0];
}
//! Return vector of wetting phase potential gradients
/*! \param Idx1 Element index
* \param isIdx Intersection index */
Scalar& potentialWetting(int Idx1, int isIdx)
{
return this->potential(Idx1, isIdx)[wPhaseIdx];
}
//! \copydoc Dumux::VariableClass2P2C::potentialWetting()
const Scalar& potentialWetting(int Idx1, int isIdx) const
{
return this->potential(Idx1, isIdx)[wPhaseIdx];
}
//! Return vector of nonwetting phase potential gradients
/*! \param Idx1 Element index
* \param isIdx Intersection index */
Scalar& potentialNonwetting(int Idx1, int isIdx)
{
return this->potential(Idx1, isIdx)[nPhaseIdx];
}
//! \copydoc Dumux::VariableClass2P2C::potentialNonwetting()
const Scalar& potentialNonwetting(int Idx1, int isIdx) const
{
return this->potential(Idx1, isIdx)[nPhaseIdx];
}
//! Return vector of wetting phase mobilities
/*! \param Idx Element index*/
Scalar& mobilityWetting(int Idx)
......
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