Skip to content
Snippets Groups Projects
Commit 6d045c54 authored by Bernd Flemisch's avatar Bernd Flemisch
Browse files

fix results from DimVector/GlobalPosition name changes

git-svn-id: svn://svn.iws.uni-stuttgart.de/DUMUX/dumux/trunk@11956 2fb0f335-1f38-0410-981e-8018bf24f1b0
parent 0634d30a
No related branches found
No related tags found
No related merge requests found
...@@ -178,11 +178,11 @@ protected: ...@@ -178,11 +178,11 @@ protected:
{ {
// initial guess of velocity: Darcy relation // initial guess of velocity: Darcy relation
// first taken from base class, later overwritten in base class // first taken from base class, later overwritten in base class
DimVector velocity = this-> velocity(phaseIdx); GlobalPosition velocity = this-> velocity(phaseIdx);
DimVector deltaV; // the change in velocity between Newton iterations GlobalPosition deltaV; // the change in velocity between Newton iterations
DimVector residual(10e10); // the residual (function value that is to be minimized) GlobalPosition residual(10e10); // the residual (function value that is to be minimized)
DimVector tmp; // temporary variable for numerical differentiation GlobalPosition tmp; // temporary variable for numerical differentiation
Tensor gradF; // slope of equation that is to be solved Tensor gradF; // slope of equation that is to be solved
// search by means of the Newton method for a root of Forchheimer equation // search by means of the Newton method for a root of Forchheimer equation
...@@ -259,13 +259,13 @@ protected: ...@@ -259,13 +259,13 @@ protected:
* \param gradPotential The gradient in potential * \param gradPotential The gradient in potential
* \param phaseIdx The index of the currently considered phase * \param phaseIdx The index of the currently considered phase
*/ */
void forchheimerResidual_(DimVector & residual, void forchheimerResidual_(GlobalPosition & residual,
const Scalar forchCoeff, const Scalar forchCoeff,
const DimMatrix & sqrtK, const DimMatrix & sqrtK,
const DimMatrix & K, const DimMatrix & K,
const DimVector & velocity, const GlobalPosition & velocity,
const ElementVolumeVariables & elemVolVars, const ElementVolumeVariables & elemVolVars,
const DimVector & gradPotential , const GlobalPosition & gradPotential ,
const unsigned int phaseIdx) const const unsigned int phaseIdx) const
{ {
const VolumeVariables upVolVars = elemVolVars[this->upstreamIdx(phaseIdx)]; const VolumeVariables upVolVars = elemVolVars[this->upstreamIdx(phaseIdx)];
...@@ -332,7 +332,7 @@ protected: ...@@ -332,7 +332,7 @@ protected:
void forchheimerDerivative_(Tensor & derivative, void forchheimerDerivative_(Tensor & derivative,
const Scalar forchCoeff, const Scalar forchCoeff,
const DimMatrix & sqrtK, const DimMatrix & sqrtK,
const DimVector & velocity, const GlobalPosition & velocity,
const ElementVolumeVariables & elemVolVars, const ElementVolumeVariables & elemVolVars,
const unsigned int phaseIdx) const const unsigned int phaseIdx) const
{ {
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment