OilfieldToolsNet/opm-simulators
4
0
Fork 0
mirror of https://github.com/OPM/opm-simulators.git synced 2025-02-25 18:55:30 -06:00
Code Issues Packages Projects Releases Wiki Activity

RegionAttributeHelpers: template Scalar type

Browse Source
This commit is contained in:
Arne Morten Kvarving 2024-02-20 13:40:33 +01:00
parent e3891843df
commit 4e7ab26841
1 changed files with 15 additions and 16 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

31
opm/simulators/wells/RegionAttributeHelpers.hpp
View File

@ -28,7 +28,6 @@
#include <dune/grid/common/gridenums.hh>
#include <algorithm>
#include <cmath>
#include <memory>
#include <stdexcept>
#include <type_traits>
@ -63,7 +62,7 @@ namespace Opm {
* In a parallel run only cells owned contribute to the cell average.
* \tparam is_parallel Whether this is a parallel run.
*/
template<bool is_parallel>
template<class Scalar, bool is_parallel>
struct AverageIncrementCalculator
{
/**
@ -77,12 +76,12 @@ namespace Opm {
* by this process (value 1), or not (value 0).
* \param cell The cell index.
*/
std::tuple<double, double, double, double, int>
operator()(const std::vector<double>& pressure,
const std::vector<double>& temperature,
const std::vector<double>& rs,
const std::vector<double>& rv,
const std::vector<double>& ownership,
std::tuple<Scalar, Scalar, Scalar, Scalar, int>
operator()(const std::vector<Scalar>& pressure,
const std::vector<Scalar>& temperature,
const std::vector<Scalar>& rs,
const std::vector<Scalar>& rv,
const std::vector<Scalar>& ownership,
std::size_t cell){
if ( ownership[cell] )
{
@ -98,15 +97,15 @@ namespace Opm {
}
}
};
template<>
struct AverageIncrementCalculator<false>
template<class Scalar>
struct AverageIncrementCalculator<Scalar, false>
{
std::tuple<double, double, double, double, int>
operator()(const std::vector<double>& pressure,
const std::vector<double>& temperature,
const std::vector<double>& rs,
const std::vector<double>& rv,
const std::vector<double>&,
std::tuple<Scalar, Scalar, Scalar, Scalar, int>
operator()(const std::vector<Scalar>& pressure,
const std::vector<Scalar>& temperature,
const std::vector<Scalar>& rs,
const std::vector<Scalar>& rv,
const std::vector<Scalar>&,
std::size_t cell){
return std::make_tuple(pressure[cell],
temperature[cell],