mirror of
https://github.com/OPM/opm-simulators.git
synced 2024-11-25 18:50:19 -06:00
a626ab869e
This works by having a "focus degree of freedom" during linearization. When evaluating the local residual, all derivatives of the residual/fluxes are with regard to the primary variables of that DOF. The two main offenders were the Forchheimer velocity model and the model for the Stokes equations. To ensure that they continue to work, the "powerinjection" and the "stokestest2c" problems are now both compiled and tested with both, automatic differentiation and finite differences, and the results of these tests is compared against the same reference solution. The majority of the time required to develop this patch was actually required for testing: All tests compile and pass with debugging and aggressive optimization flags with at least GCC 5, GCC 7 and clang 3.8, as well as Dune 2.3 and 2.4. Also, the results of flow_ebos stay identical for Norne whilst the performance difference is below the measurement noise on my machine. (the version with this patch applied was actually about 1% faster.)
56 lines
2.0 KiB
C++
56 lines
2.0 KiB
C++
// -*- mode: C++; tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 4 -*-
|
|
// vi: set et ts=4 sw=4 sts=4:
|
|
/*
|
|
This file is part of the Open Porous Media project (OPM).
|
|
|
|
OPM is free software: you can redistribute it and/or modify
|
|
it under the terms of the GNU General Public License as published by
|
|
the Free Software Foundation, either version 2 of the License, or
|
|
(at your option) any later version.
|
|
|
|
OPM is distributed in the hope that it will be useful,
|
|
but WITHOUT ANY WARRANTY; without even the implied warranty of
|
|
MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
|
|
GNU General Public License for more details.
|
|
|
|
You should have received a copy of the GNU General Public License
|
|
along with OPM. If not, see <http://www.gnu.org/licenses/>.
|
|
|
|
Consult the COPYING file in the top-level source directory of this
|
|
module for the precise wording of the license and the list of
|
|
copyright holders.
|
|
*/
|
|
/*!
|
|
* \file
|
|
*
|
|
* \brief Two-phase test for the immiscible model which uses the element-centered finite
|
|
* volume discretization in conjunction with automatic differentiation
|
|
*/
|
|
#include "config.h"
|
|
|
|
#include <ewoms/common/start.hh>
|
|
#include <ewoms/models/immiscible/immisciblemodel.hh>
|
|
#include <ewoms/disc/ecfv/ecfvdiscretization.hh>
|
|
#include "problems/lensproblem.hh"
|
|
|
|
namespace Ewoms {
|
|
namespace Properties {
|
|
NEW_TYPE_TAG(LensProblemEcfvAd, INHERITS_FROM(ImmiscibleTwoPhaseModel, LensBaseProblem));
|
|
|
|
// use the element centered finite volume spatial discretization
|
|
SET_TAG_PROP(LensProblemEcfvAd, SpatialDiscretizationSplice, EcfvDiscretization);
|
|
|
|
// use automatic differentiation for this simulator
|
|
SET_TAG_PROP(LensProblemEcfvAd, LocalLinearizerSplice, AutoDiffLocalLinearizer);
|
|
|
|
// this problem works fine if the linear solver uses single precision scalars
|
|
SET_TYPE_PROP(LensProblemEcfvAd, LinearSolverScalar, float);
|
|
|
|
}}
|
|
|
|
int main(int argc, char **argv)
|
|
{
|
|
typedef TTAG(LensProblemEcfvAd) ProblemTypeTag;
|
|
return Ewoms::start<ProblemTypeTag>(argc, argv);
|
|
}
|