diff --git a/dumux/boxmodels/1p/1plocalresidual.hh b/dumux/boxmodels/1p/1plocalresidual.hh index cbac28aed2ad9c4d697706f8adcc9fc9b3c1c9a2..502f0d5123a15720a5b693de5bc17cb39d5ced02 100644 --- a/dumux/boxmodels/1p/1plocalresidual.hh +++ b/dumux/boxmodels/1p/1plocalresidual.hh @@ -135,7 +135,10 @@ public: Vector tmpVec; fluxVars.intrinsicPermeability().mv(fluxVars.potentialGrad(), 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 &dn = this->curVolVars_(fluxVars.downstreamIdx(normalFlux)); diff --git a/dumux/boxmodels/1p/1pmodel.hh b/dumux/boxmodels/1p/1pmodel.hh index 21a8cabcd33591ce91e529d067b63f2cef084bda..f49799f193ea2edfdb46c005804df8014391d89a 100644 --- a/dumux/boxmodels/1p/1pmodel.hh +++ b/dumux/boxmodels/1p/1pmodel.hh @@ -93,14 +93,14 @@ public: void addOutputVtkFields(const SolutionVector &sol, MultiWriter &writer) { - typedef Dune::BlockVector<Dune::FieldVector<Scalar, 1> > ScalarField; + typedef Dune::BlockVector<Dune::FieldVector<double, 1> > ScalarField; // create the required scalar fields unsigned numVertices = this->problem_().gridView().size(dim); - ScalarField *p = writer.allocateManagedBuffer (numVertices); + ScalarField *p = writer.allocateManagedBuffer(numVertices); unsigned numElements = this->gridView_().size(0); - ScalarField *rank = writer.allocateManagedBuffer (numElements); + ScalarField *rank = writer.allocateManagedBuffer(numElements); FVElementGeometry fvElemGeom; VolumeVariables volVars;