Skip to content
Snippets Groups Projects
Commit 73189fd2 authored by Andreas Lauser's avatar Andreas Lauser
Browse files

1p box model: make it compile with Scalar != double

git-svn-id: svn://svn.iws.uni-stuttgart.de/DUMUX/dumux/trunk@6570 2fb0f335-1f38-0410-981e-8018bf24f1b0
parent 598a77ec
No related branches found
No related tags found
No related merge requests found
...@@ -135,7 +135,10 @@ public: ...@@ -135,7 +135,10 @@ public:
Vector tmpVec; Vector tmpVec;
fluxVars.intrinsicPermeability().mv(fluxVars.potentialGrad(), fluxVars.intrinsicPermeability().mv(fluxVars.potentialGrad(),
tmpVec); tmpVec);
Scalar normalFlux = - (tmpVec*fluxVars.face().normal); Scalar normalFlux = 0;
for (int i = 0; i < Vector::size; ++i)
normalFlux += tmpVec[i]*fluxVars.face().normal[i];
normalFlux *= -1;
const VolumeVariables &up = this->curVolVars_(fluxVars.upstreamIdx(normalFlux)); const VolumeVariables &up = this->curVolVars_(fluxVars.upstreamIdx(normalFlux));
const VolumeVariables &dn = this->curVolVars_(fluxVars.downstreamIdx(normalFlux)); const VolumeVariables &dn = this->curVolVars_(fluxVars.downstreamIdx(normalFlux));
......
...@@ -93,14 +93,14 @@ public: ...@@ -93,14 +93,14 @@ public:
void addOutputVtkFields(const SolutionVector &sol, void addOutputVtkFields(const SolutionVector &sol,
MultiWriter &writer) MultiWriter &writer)
{ {
typedef Dune::BlockVector<Dune::FieldVector<Scalar, 1> > ScalarField; typedef Dune::BlockVector<Dune::FieldVector<double, 1> > ScalarField;
// create the required scalar fields // create the required scalar fields
unsigned numVertices = this->problem_().gridView().size(dim); unsigned numVertices = this->problem_().gridView().size(dim);
ScalarField *p = writer.allocateManagedBuffer (numVertices); ScalarField *p = writer.allocateManagedBuffer(numVertices);
unsigned numElements = this->gridView_().size(0); unsigned numElements = this->gridView_().size(0);
ScalarField *rank = writer.allocateManagedBuffer (numElements); ScalarField *rank = writer.allocateManagedBuffer(numElements);
FVElementGeometry fvElemGeom; FVElementGeometry fvElemGeom;
VolumeVariables volVars; VolumeVariables volVars;
......
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