// -*- 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 This is the unit test for the black oil fluid system
 *
 * This test requires the presence of opm-parser.
 */
#include "config.h"

#if !HAVE_ECL_INPUT
#error "The test for the black oil fluid system classes requires ecl input support in opm-common"
#endif

#include <boost/mpl/list.hpp>

#define BOOST_TEST_MODULE EclBlackOilFluidSystem
#include <boost/test/unit_test.hpp>

#include <opm/material/fluidsystems/BlackOilFluidSystem.hpp>
#include <opm/material/fluidstates/BlackOilFluidState.hpp>
#include <opm/material/densead/Evaluation.hpp>

#include <opm/material/common/Valgrind.hpp>

#include <opm/input/eclipse/Parser/Parser.hpp>
#include <opm/input/eclipse/Deck/Deck.hpp>
#include <opm/input/eclipse/EclipseState/EclipseState.hpp>
#include <opm/input/eclipse/Python/Python.hpp>
#include <opm/input/eclipse/Schedule/Schedule.hpp>

#include <type_traits>
#include <cmath>

// values of strings based on the SPE1 and NORNE cases of opm-data.
static constexpr const char* deckString1 =
    "RUNSPEC\n"
    "\n"
    "DIMENS\n"
    "   10 10 3 /\n"
    "\n"
    "METRIC\n"
    "\n"
    "TABDIMS\n"
    " * 2 /\n"
    "\n"
    "OIL\n"
    "GAS\n"
    "WATER\n"
    "\n"
    "DISGAS\n"
    "VAPOIL\n"
    "\n"
    "GRID\n"
    "\n"
    "DX\n"
    "   	300*1000 /\n"
    "DY\n"
    "	300*1000 /\n"
    "DZ\n"
    "	100*20 100*30 100*50 /\n"
    "\n"
    "TOPS\n"
    "	100*1234 /\n"
    "\n"
    "PORO\n"
    "  300*0.15 /\n"
    "PROPS\n"
    "\n"
    "-- Copyright (C) 2015 Statoil\n"
    "\n"
    "PVTG\n"
    "\n"
    "-- PRESSURE       RSG        B-GAS     VISCOSITY\n"
    "--   BAR                                 (CP)\n"
    "\n"
    "     50.00    0.00000497   0.024958     0.01441\n"
    "              0.00000248   0.024958     0.01440\n"
    "              0.00000000   0.024958     0.01440 /\n"
    "\n"
    "     70.00    0.00000521   0.017639     0.01491\n"
    "              0.00000261   0.017641     0.01490\n"
    "              0.00000000   0.017643     0.01490 /\n"
    "\n"
    "     90.00    0.00000627   0.013608     0.01547\n"
    "              0.00000313   0.013611     0.01546\n"
    "              0.00000000   0.013615     0.01544 /\n"
    "\n"
    "    110.00    0.00000798   0.011072     0.01609\n"
    "              0.00000399   0.011076     0.01607\n"
    "              0.00000000   0.011081     0.01605 /\n"
    "\n"
    "    130.00    0.00001041   0.009340     0.01677\n"
    "              0.00000520   0.009346     0.01674\n"
    "              0.00000000   0.009352     0.01671 /\n"
    "\n"
    "    150.00    0.00001365   0.008092     0.01752\n"
    "              0.00000683   0.008099     0.01748\n"
    "              0.00000000   0.008106     0.01743 /\n"
    "\n"
    "    170.00    0.00001786   0.007156     0.01834\n"
    "              0.00000893   0.007164     0.01827\n"
    "              0.00000000   0.007172     0.01819 /\n"
    "\n"
    "    190.00    0.00002316   0.006433     0.01923\n"
    "              0.00001158   0.006442     0.01912\n"
    "              0.00000000   0.006451     0.01900 /\n"
    "\n"
    "    210.00    0.00002972   0.005861     0.02019\n"
    "              0.00001486   0.005871     0.02001\n"
    "              0.00000000   0.005881     0.01984 /\n"
    "\n"
    "    230.00    0.00003767   0.005402     0.02121\n"
    "              0.00001883   0.005412     0.02095\n"
    "              0.00000000   0.005422     0.02071 /\n"
    "\n"
    "    250.80    0.00004756   0.005013     0.02234\n"
    "              0.00002378   0.005022     0.02197\n"
    "              0.00000000   0.005032     0.02162 /\n"
    "\n"
    "    268.42    0.00005757   0.004737     0.02335\n"
    "              0.00002878   0.004746     0.02287\n"
    "              0.00000000   0.004754     0.02240 /\n"
    "\n"
    "    285.33    0.00006853   0.004511     0.02438\n"
    "              0.00003427   0.004518     0.02375\n"
    "              0.00000000   0.004525     0.02315 /\n"
    "\n"
    "    301.59    0.00008041   0.004323     0.02542\n"
    "              0.00004020   0.004327     0.02463\n"
    "              0.00000000   0.004332     0.02387 /\n"
    "\n"
    "    317.23    0.00009313   0.004165     0.02648\n"
    "              0.00004657   0.004166     0.02549\n"
    "              0.00000000   0.004169     0.02456 /\n"
    "\n"
    "    332.29    0.00010668   0.004031     0.02755\n"
    "              0.00005334   0.004029     0.02634\n"
    "              0.00000000   0.004028     0.02522 /\n"
    "\n"
    "    346.80    0.00012100   0.003917     0.02863\n"
    "              0.00006050   0.003911     0.02719\n"
    "              0.00000000   0.003906     0.02585 /\n"
    "\n"
    "    360.80    0.00013607   0.003819     0.02974\n"
    "              0.00006803   0.003808     0.02803\n"
    "              0.00000000   0.003799     0.02645 /\n"
    "\n"
    "    374.31    0.00015188   0.003735     0.03087\n"
    "              0.00007594   0.003718     0.02887\n"
    "              0.00000000   0.003705     0.02703 /\n"
    "\n"
    "    387.36    0.00016843   0.003662     0.03202\n"
    "              0.00008421   0.003639     0.02970\n"
    "              0.00000000   0.003621     0.02758 /\n"
    "\n"
    "    399.99    0.00018571   0.003598     0.03320\n"
    "              0.00009286   0.003570     0.03053\n"
    "              0.00000000   0.003545     0.02810 /\n"
    "\n"
    "    412.21    0.00020375   0.003543     0.03442\n"
    "              0.00010188   0.003508     0.03137\n"
    "              0.00000000   0.003477     0.02861 /\n"
    "\n"
    "    424.05    0.00022256   0.003496     0.03566\n"
    "              0.00011128   0.003453     0.03220\n"
    "              0.00000000   0.003416     0.02909 /\n"
    "\n"
    "    435.53    0.00024218   0.003454     0.03695\n"
    "              0.00012109   0.003404     0.03305\n"
    "              0.00000000   0.003360     0.02956 /\n"
    "\n"
    "    446.68    0.00026266   0.003419     0.03828\n"
    "              0.00013133   0.003360     0.03390\n"
    "              0.00000000   0.003309     0.03000 /\n"
    "\n"
    "    457.51    0.00028404   0.003388     0.03967\n"
    "              0.00014202   0.003320     0.03477\n"
    "              0.00000000   0.003262     0.03043 /\n"
    "\n"
    "    468.04    0.00030639   0.003362     0.04110\n"
    "              0.00015319   0.003285     0.03566\n"
    "              0.00000000   0.003218     0.03085 /\n"
    "\n"
    "    478.30    0.00032980   0.003341     0.04261\n"
    "              0.00016490   0.003253     0.03656\n"
    "              0.00000000   0.003178     0.03125 /\n"
    "\n"
    "    488.30    0.00035436   0.003323     0.04418\n"
    "              0.00017718   0.003225     0.03749\n"
    "              0.00000000   0.003140     0.03164 /\n"
    "\n"
    "    498.06    0.00038020   0.003310     0.04583\n"
    "              0.00019010   0.003200     0.03845\n"
    "              0.00000000   0.003105     0.03202 /\n"
    "\n"
    "    507.59    0.00040745   0.003300     0.04758\n"
    "              0.00020373   0.003178     0.03944\n"
    "              0.00000000   0.003073     0.03238 /\n"
    "\n"
    "    516.92    0.00043630   0.003293     0.04943\n"
    "              0.00021815   0.003158     0.04048\n"
    "              0.00000000   0.003042     0.03273 /\n"
    "\n"
    "    526.06    0.00046694   0.003290     0.05141\n"
    "              0.00023347   0.003141     0.04156\n"
    "              0.00000000   0.003013     0.03308 /\n"
    "\n"
    "    535.02    0.00049963   0.003291     0.05353\n"
    "              0.00024981   0.003126     0.04271\n"
    "              0.00000000   0.002986     0.03342 /\n"
    "\n"
    "    543.83    0.00053469   0.003295     0.05582\n"
    "              0.00026734   0.003114     0.04393\n"
    "              0.00000000   0.002960     0.03374 /\n"
    "\n"
    "    552.49    0.00057251   0.003303     0.05830\n"
    "              0.00028625   0.003105     0.04523\n"
    "              0.00000000   0.002935     0.03407 /\n"
    "\n"
    "    561.04    0.00061359   0.003315     0.06103\n"
    "              0.00030679   0.003098     0.04664\n"
    "              0.00000000   0.002912     0.03438 /\n"
    "\n"
    "    569.48    0.00065855   0.003332     0.06405\n"
    "              0.00032928   0.003093     0.04818\n"
    "              0.00000000   0.002890     0.03469 /\n"
    "\n"
    "    577.82    0.00070820   0.003354     0.06744\n"
    "              0.00035410   0.003092     0.04988\n"
    "              0.00000000   0.002868     0.03500 /\n"
    "\n"
    "    586.09    0.00076355   0.003382     0.07127\n"
    "              0.00038178   0.003094     0.05178\n"
    "              0.00000000   0.002847     0.03530 /\n"
    "\n"
    "    594.29    0.00082592   0.003418     0.07567\n"
    "              0.00041296   0.003099     0.05394\n"
    "              0.00000000   0.002828     0.03560 /\n"
    "/\n"
    "\n"
    "-- PVT region 2 --\n"
    "\n"
    "     80.00    0.00000485   0.015151     0.01506\n"
    "              0.00000242   0.015154     0.01505\n"
    "              0.00000000   0.015157     0.01504 /\n"
    "\n"
    "\n"
    "    100.00    0.00000621   0.012032     0.01566\n"
    "              0.00000310   0.012036     0.01564\n"
    "              0.00000000   0.012040     0.01563 /\n"
    "\n"
    "\n"
    "    120.00    0.00000821   0.009980     0.01632\n"
    "              0.00000411   0.009985     0.01630\n"
    "              0.00000000   0.009990     0.01628 /\n"
    "\n"
    "\n"
    "    140.00    0.00001096   0.008537     0.01706\n"
    "              0.00000548   0.008544     0.01702\n"
    "              0.00000000   0.008550     0.01698 /\n"
    "\n"
    "\n"
    "    160.00    0.00001457   0.007476     0.01786\n"
    "              0.00000728   0.007484     0.01780\n"
    "              0.00000000   0.007492     0.01774 /\n"
    "\n"
    "\n"
    "    180.00    0.00001918   0.006669     0.01873\n"
    "              0.00000959   0.006678     0.01864\n"
    "              0.00000000   0.006687     0.01854 /\n"
    "\n"
    "\n"
    "    200.00    0.00002493   0.006038     0.01967\n"
    "              0.00001247   0.006048     0.01952\n"
    "              0.00000000   0.006058     0.01939 /\n"
    "\n"
    "\n"
    "    216.50    0.00003061   0.005616     0.02049\n"
    "              0.00001530   0.005626     0.02029\n"
    "              0.00000000   0.005636     0.02010 /\n"
    "/\n"
    "\n"
    "PVTO\n"
    "\n"
    "--  RSO    PRESSURE    B-OIL     VISCOSITY\n"
    "--          (BAR)                  (CP)\n"
    "\n"
    "   20.59     50.00    1.10615     1.180\n"
    "             75.00    1.10164     1.247\n"
    "            100.00    1.09744     1.315\n"
    "            125.00    1.09351     1.384\n"
    "            150.00    1.08984     1.453  /\n"
    "\n"
    "   28.19     70.00    1.12522     1.066\n"
    "             95.00    1.12047     1.124\n"
    "            120.00    1.11604     1.182\n"
    "            145.00    1.11191     1.241\n"
    "            170.00    1.10804     1.300  /\n"
    "\n"
    "   36.01     90.00    1.14458     0.964\n"
    "            115.00    1.13959     1.014\n"
    "            140.00    1.13494     1.064\n"
    "            165.00    1.13060     1.115\n"
    "            190.00    1.12653     1.166  /\n"
    "\n"
    "   44.09    110.00    1.16437     0.880\n"
    "            135.00    1.15915     0.924\n"
    "            160.00    1.15428     0.968\n"
    "            185.00    1.14973     1.012\n"
    "            210.00    1.14547     1.056  /\n"
    "\n"
    "   52.46    130.00    1.18467     0.805\n"
    "            155.00    1.17921     0.843\n"
    "            180.00    1.17413     0.882\n"
    "            205.00    1.16937     0.920\n"
    "            230.00    1.16491     0.959  /\n"
    "\n"
    "   61.13    150.00    1.20555     0.746\n"
    "            175.00    1.19985     0.780\n"
    "            200.00    1.19454     0.814\n"
    "            225.00    1.18958     0.849\n"
    "            250.00    1.18492     0.883  /\n"
    "\n"
    "   70.14    170.00    1.22704     0.698\n"
    "            195.00    1.22111     0.729\n"
    "            220.00    1.21558     0.759\n"
    "            245.00    1.21040     0.790\n"
    "            270.00    1.20555     0.821  /\n"
    "\n"
    "   79.50    190.00    1.24922     0.658\n"
    "            215.00    1.24305     0.686\n"
    "            240.00    1.23729     0.714\n"
    "            265.00    1.23190     0.742\n"
    "            290.00    1.22685     0.770  /\n"
    "\n"
    "   89.24    210.00    1.27214     0.637\n"
    "            235.00    1.26573     0.664\n"
    "            260.00    1.25974     0.693\n"
    "            285.00    1.25414     0.725\n"
    "            310.00    1.24888     0.760  /\n"
    "\n"
    "   99.39    230.00    1.29586     0.622\n"
    "            255.00    1.28921     0.641\n"
    "            280.00    1.28300     0.661\n"
    "            305.00    1.27718     0.680\n"
    "            330.00    1.27171     0.699  /\n"
    "\n"
    "  110.41    250.80    1.32148     0.610\n"
    "            275.80    1.31457     0.628\n"
    "            300.80    1.30812     0.647\n"
    "            325.80    1.30207     0.665\n"
    "            350.80    1.29638     0.682  /\n"
    "\n"
    "  120.32    268.42    1.34449     0.576\n"
    "            293.42    1.33735     0.593\n"
    "            318.42    1.33068     0.609\n"
    "            343.42    1.32442     0.626\n"
    "            368.42    1.31853     0.642  /\n"
    "\n"
    "  130.23    285.33    1.36737     0.5335\n"
    "            310.33    1.36001     0.5487\n"
    "            335.33    1.35313     0.5638\n"
    "            360.33    1.34667     0.5787\n"
    "            385.33    1.34059     0.5934  /\n"
    "\n"
    "  140.12    301.59    1.39015     0.4956\n"
    "            326.59    1.38257     0.5094\n"
    "            351.59    1.37548     0.5230\n"
    "            376.59    1.36882     0.5365\n"
    "            401.59    1.36255     0.5498  /\n"
    "\n"
    "  150.01    317.23    1.41282     0.4614\n"
    "            342.23    1.40503     0.4739\n"
    "            367.23    1.39773     0.4863\n"
    "            392.23    1.39088     0.4986\n"
    "            417.23    1.38443     0.5107  /\n"
    "\n"
    "  159.89    332.29    1.43539     0.43042\n"
    "            357.29    1.42739     0.44183\n"
    "            382.29    1.41990     0.45312\n"
    "            407.29    1.41286     0.46430\n"
    "            432.29    1.40622     0.47537 /\n"
    "\n"
    "  169.76    346.80    1.45788     0.41191\n"
    "            371.80    1.44967     0.42260\n"
    "            396.80    1.44198     0.43318\n"
    "            421.80    1.43475     0.44365\n"
    "            446.80    1.42794     0.45402 /\n"
    "\n"
    "  179.63    360.80    1.48028     0.39503\n"
    "            385.80    1.47187     0.40508\n"
    "            410.80    1.46398     0.41502\n"
    "            435.80    1.45657     0.42487\n"
    "            460.80    1.44958     0.43461 /\n"
    "\n"
    "  189.48    374.31    1.50260     0.37959\n"
    "            399.31    1.49399     0.38907\n"
    "            424.31    1.48591     0.39845\n"
    "            449.31    1.47832     0.40773\n"
    "            474.31    1.47116     0.41692 /\n"
    "\n"
    "  199.34    387.36    1.52484     0.36543\n"
    "            412.36    1.51603     0.37439\n"
    "            437.36    1.50777     0.38326\n"
    "            462.36    1.50000     0.39203\n"
    "            487.36    1.49267     0.40072 /\n"
    "\n"
    "  209.18    399.99    1.54700     0.35239\n"
    "            424.99    1.53800     0.36089\n"
    "            449.99    1.52956     0.36929\n"
    "            474.99    1.52161     0.37762\n"
    "            499.99    1.51411     0.38585 /\n"
    "\n"
    "  219.02    412.21    1.56910     0.34035\n"
    "            437.21    1.55991     0.34843\n"
    "            462.21    1.55128     0.35642\n"
    "            487.21    1.54316     0.36433\n"
    "            512.21    1.53549     0.37216 /\n"
    "\n"
    "  228.85    424.05    1.59112     0.32921\n"
    "            449.05    1.58174     0.33691\n"
    "            474.05    1.57294     0.34453\n"
    "            499.05    1.56464     0.35206\n"
    "            524.05    1.55681     0.35952 /\n"
    "\n"
    "  238.67    435.53    1.61307     0.31888\n"
    "            460.53    1.60351     0.32623\n"
    "            485.53    1.59453     0.33350\n"
    "            510.53    1.58606     0.34070\n"
    "            535.53    1.57807     0.34782 /\n"
    "\n"
    "  248.48    446.68    1.63496     0.30927\n"
    "            471.68    1.62522     0.31630\n"
    "            496.68    1.61606     0.32326\n"
    "            521.68    1.60743     0.33014\n"
    "            546.68    1.59927     0.33695 /\n"
    "\n"
    "  258.29    457.51    1.65678     0.30032\n"
    "            482.51    1.64686     0.30706\n"
    "            507.51    1.63753     0.31373\n"
    "            532.51    1.62873     0.32032\n"
    "            557.51    1.62042     0.32685 /\n"
    "\n"
    "  268.09    468.04    1.67853     0.29196\n"
    "            493.04    1.66843     0.29843\n"
    "            518.04    1.65893     0.30483\n"
    "            543.04    1.64997     0.31117\n"
    "            568.04    1.64150     0.31743 /\n"
    "\n"
    "  277.89    478.30    1.70022     0.28414\n"
    "            503.30    1.68994     0.29037\n"
    "            528.30    1.68028     0.29652\n"
    "            553.30    1.67116     0.30261\n"
    "            578.30    1.66253     0.30864 /\n"
    "\n"
    "  287.68    488.30    1.72184     0.27681\n"
    "            513.30    1.71139     0.28281\n"
    "            538.30    1.70156     0.28874\n"
    "            563.30    1.69228     0.29460\n"
    "            588.30    1.68350     0.30040 /\n"
    "\n"
    "  297.46    498.06    1.74339     0.26994\n"
    "            523.06    1.73277     0.27572\n"
    "            548.06    1.72278     0.28144\n"
    "            573.06    1.71334     0.28709\n"
    "            598.06    1.70442     0.29269 /\n"
    "\n"
    "  307.23    507.59    1.76487     0.26347\n"
    "            532.59    1.75409     0.26906\n"
    "            557.59    1.74393     0.27458\n"
    "            582.59    1.73434     0.28004\n"
    "            607.59    1.72527     0.28544 /\n"
    "\n"
    "  317.00    516.92    1.78628     0.25738\n"
    "            541.92    1.77533     0.26279\n"
    "            566.92    1.76502     0.26812\n"
    "            591.92    1.75528     0.27340\n"
    "            616.92    1.74606     0.27863 /\n"
    "\n"
    "  326.76    526.06    1.80761     0.25165\n"
    "            551.06    1.79651     0.25688\n"
    "            576.06    1.78604     0.26204\n"
    "            601.06    1.77615     0.26716\n"
    "            626.06    1.76679     0.27221 /\n"
    "\n"
    "  336.51    535.02    1.82887     0.24623\n"
    "            560.02    1.81761     0.25130\n"
    "            585.02    1.80699     0.25631\n"
    "            610.02    1.79696     0.26126\n"
    "            635.02    1.78746     0.26616 /\n"
    "\n"
    "  346.26    543.83    1.85005     0.24112\n"
    "            568.83    1.83864     0.24603\n"
    "            593.83    1.82787     0.25089\n"
    "            618.83    1.81770     0.25570\n"
    "            643.83    1.80806     0.26045 /\n"
    "\n"
    "  356.00    552.49    1.87115     0.23628\n"
    "            577.49    1.85959     0.24105\n"
    "            602.49    1.84868     0.24577\n"
    "            627.49    1.83836     0.25043\n"
    "            652.49    1.82858     0.25505 /\n"
    "\n"
    "  365.73    561.04    1.89217     0.23170\n"
    "            586.04    1.88046     0.23634\n"
    "            611.04    1.86940     0.24092\n"
    "            636.04    1.85895     0.24546\n"
    "            661.04    1.84904     0.24994 /\n"
    "\n"
    "  375.46    569.48    1.91309     0.22736\n"
    "            594.48    1.90124     0.23187\n"
    "            619.48    1.89004     0.23633\n"
    "            644.48    1.87946     0.24074\n"
    "            669.48    1.86942     0.24510 /\n"
    "\n"
    "  385.18    577.82    1.93391     0.22325\n"
    "            602.82    1.92192     0.22764\n"
    "            627.82    1.91060     0.23198\n"
    "            652.82    1.89988     0.23627\n"
    "            677.82    1.88971     0.24052 /\n"
    "\n"
    "  394.89    586.09    1.95464     0.21934\n"
    "            611.09    1.94252     0.22362\n"
    "            636.09    1.93106     0.22785\n"
    "            661.09    1.92021     0.23204\n"
    "            686.09    1.90993     0.23617 /\n"
    "\n"
    "  404.60    594.29    1.97527     0.21564\n"
    "            619.29    1.96301     0.21981\n"
    "            644.29    1.95143     0.22393\n"
    "            669.29    1.94046     0.22801\n"
    "            694.29    1.93005     0.23204 /\n"
    "/\n"
    "\n"
    "-- PVT region 2\n"
    "\n"
    "   32.91     80.00    1.13304     1.04537\n"
    "            114.00    1.12837     1.10009\n"
    "            148.00    1.12401     1.15521\n"
    "            182.00    1.11994     1.21063 /\n"
    "\n"
    "\n"
    "   40.99    100.00    1.15276     0.97219\n"
    "            134.00    1.14786     1.02086\n"
    "            168.00    1.14328     1.06983\n"
    "            202.00    1.13900     1.11901 /\n"
    "\n"
    "\n"
    "   49.36    120.00    1.17297     0.91124\n"
    "            154.00    1.16783     0.95496\n"
    "            188.00    1.16303     0.99891\n"
    "            222.00    1.15854     1.04301 /\n"
    "\n"
    "\n"
    "   58.04    140.00    1.19374     0.85942\n"
    "            174.00    1.18836     0.89902\n"
    "            208.00    1.18334     0.93878\n"
    "            242.00    1.17864     0.97864 /\n"
    "\n"
    "\n"
    "   67.04    160.00    1.21512     0.81456\n"
    "            194.00    1.20951     0.85065\n"
    "            228.00    1.20426     0.88686\n"
    "            262.00    1.19935     0.92313 /\n"
    "\n"
    "\n"
    "   76.39    180.00    1.23718     0.77508\n"
    "            214.00    1.23132     0.80815\n"
    "            248.00    1.22585     0.84130\n"
    "            282.00    1.22073     0.87448 /\n"
    "\n"
    "\n"
    "   86.11    200.00    1.25996     0.73928\n"
    "            234.00    1.25386     0.76989\n"
    "            268.00    1.24816     0.80050\n"
    "            302.00    1.24283     0.83108 /\n"
    "\n"
    "\n"
    "   94.44    216.50    1.27934     0.67686\n"
    "            250.50    1.27304     0.70995\n"
    "            284.50    1.26716     0.74427\n"
    "            318.50    1.26164     0.77857 /\n"
    "/\n"
    "\n"
    "PVTW\n"
    "    277.0      1.038      4.67E-5    0.318       0.0 /\n"
    "    277.0      1.038      4.67E-5    0.318       0.0 /\n"
    "\n"
    "DENSITY\n"
    "      859.5  1033.0    0.854  /\n"
    "      860.04 1033.0    0.853  /\n"
    "\n";

using Types = boost::mpl::list<double,Opm::DenseAd::Evaluation<double,2>>;

BOOST_AUTO_TEST_CASE_TEMPLATE(BlackOil, Evaluation, Types)
{
    // test the black-oil specific methods of BlackOilFluidSystem. The generic methods
    // for fluid systems are already tested by the generic test for all fluidsystems.

    using Scalar = typename Opm::MathToolbox<Evaluation>::Scalar;
    using FluidSystem = Opm::BlackOilFluidSystem<double>;

    static constexpr int numPhases = FluidSystem::numPhases;

    static constexpr int gasPhaseIdx = FluidSystem::gasPhaseIdx;
    static constexpr int oilPhaseIdx = FluidSystem::oilPhaseIdx;
    static constexpr int waterPhaseIdx = FluidSystem::waterPhaseIdx;

    static constexpr int gasCompIdx = FluidSystem::gasCompIdx;
    static constexpr int oilCompIdx = FluidSystem::oilCompIdx;
    static constexpr int waterCompIdx = FluidSystem::waterCompIdx;

    Opm::Parser parser;

    auto deck = parser.parseString(deckString1);
    auto python = std::make_shared<Opm::Python>();
    Opm::EclipseState eclState(deck);
    Opm::Schedule schedule(deck, eclState, python);

    FluidSystem::initFromState(eclState, schedule);

    // create a parameter cache
    using ParamCache = typename FluidSystem::template ParameterCache<Scalar>;
    ParamCache paramCache(/*maxOilSat=*/0.5, /*regionIdx=*/1);
    BOOST_CHECK_EQUAL(paramCache.regionIndex(), 1);

    BOOST_CHECK_SMALL(Opm::abs(paramCache.maxOilSat() - 0.5), 1e-10);
    BOOST_CHECK_SMALL(Opm::abs(FluidSystem::reservoirTemperature() - (273.15 + 15.555)), 1e-10);
    BOOST_CHECK(FluidSystem::enableDissolvedGas());
    BOOST_CHECK(FluidSystem::enableVaporizedOil());
    BOOST_CHECK_EQUAL(FluidSystem::numRegions(), 2);
    BOOST_CHECK_EQUAL(FluidSystem::numActivePhases(), 3);

    for (unsigned phaseIdx = 0; phaseIdx < numPhases; ++ phaseIdx)
        BOOST_CHECK(FluidSystem::phaseIsActive(phaseIdx));

    BOOST_CHECK_EQUAL(FluidSystem::solventComponentIndex(oilPhaseIdx), oilCompIdx);
    BOOST_CHECK_EQUAL(FluidSystem::solventComponentIndex(gasPhaseIdx), gasCompIdx);
    BOOST_CHECK_EQUAL(FluidSystem::solventComponentIndex(waterPhaseIdx), waterCompIdx);

    BOOST_CHECK_EQUAL(FluidSystem::soluteComponentIndex(oilPhaseIdx), gasCompIdx);
    BOOST_CHECK_EQUAL(FluidSystem::soluteComponentIndex(gasPhaseIdx), oilCompIdx);

    BOOST_CHECK_SMALL(std::abs(FluidSystem::referenceDensity(oilPhaseIdx, /*regionIdx=*/1) - 860.04), 1e-10);
    BOOST_CHECK_SMALL(std::abs(FluidSystem::referenceDensity(gasPhaseIdx, /*regionIdx=*/1) - 0.853), 1e-10);
    BOOST_CHECK_SMALL(std::abs(FluidSystem::referenceDensity(waterPhaseIdx, /*regionIdx=*/1) - 1033), 1e-10);

    Opm::BlackOilFluidState<Scalar, FluidSystem> fluidState;
    Opm::Valgrind::SetUndefined(fluidState);

    static const Scalar eps = std::sqrt(std::numeric_limits<Scalar>::epsilon());
    unsigned regionIdx = paramCache.regionIndex();
    for (unsigned i = 0; i < 1000; ++i) {
        Scalar p = Scalar(i)/1000*350e5 + 100e5;

        for (unsigned phaseIdx = 0; phaseIdx < numPhases; ++phaseIdx) {
            fluidState.setPressure(phaseIdx, p);
            fluidState.setSaturation(phaseIdx, 1e-3);
        }

        Scalar RsSat = FluidSystem::saturatedDissolutionFactor(fluidState, oilPhaseIdx, regionIdx);
        fluidState.setRs(RsSat);

        Scalar RvSat = FluidSystem::saturatedDissolutionFactor(fluidState, gasPhaseIdx, regionIdx);
        fluidState.setRv(RvSat);

        paramCache.updateAll(fluidState);

        for (unsigned phaseIdx = 0; phaseIdx < numPhases; ++phaseIdx) {
            // ensure that the black-oil specific variants of the methods to compute
            // thermdynamic properties return the same value as the generic ones and that
            // the generic methods return the same value as the ones for the saturated
            // quantities (we specify the fluid state to be on the saturation line)
            BOOST_CHECK_SMALL(Opm::abs(FluidSystem::density(fluidState, paramCache, phaseIdx) -
                                       FluidSystem::density(fluidState, phaseIdx, regionIdx)), eps);

            BOOST_CHECK_SMALL(Opm::abs(FluidSystem::density(fluidState, paramCache, phaseIdx) -
                                       FluidSystem::saturatedDensity(fluidState, phaseIdx, regionIdx)), eps);

            Scalar b = FluidSystem::inverseFormationVolumeFactor(fluidState, phaseIdx, regionIdx);
            Scalar bSat = FluidSystem::saturatedInverseFormationVolumeFactor(fluidState, phaseIdx, regionIdx);
            BOOST_CHECK_SMALL(Opm::abs(b - bSat), eps);

            BOOST_CHECK_SMALL(Opm::abs(FluidSystem::viscosity(fluidState, paramCache, phaseIdx) -
                                       FluidSystem::viscosity(fluidState, phaseIdx, regionIdx)), 1e-10);

            Scalar R = FluidSystem::saturatedDissolutionFactor(fluidState, phaseIdx, regionIdx);
            Scalar R2 = FluidSystem::saturatedDissolutionFactor(fluidState, phaseIdx, regionIdx);
            BOOST_CHECK_SMALL(Opm::abs(R - R2), eps);

            // water is immiscible and thus there is no saturation pressure
            if (phaseIdx != waterPhaseIdx) {
                BOOST_CHECK_SMALL(Opm::abs(FluidSystem::saturationPressure(fluidState, phaseIdx, regionIdx) - p), eps*p);
            }
        }

        BOOST_CHECK_SMALL(Opm::abs(FluidSystem::bubblePointPressure(fluidState, regionIdx) - p), eps*p);
        BOOST_CHECK_SMALL(Opm::abs(FluidSystem::dewPointPressure(fluidState, regionIdx) - p), eps*p);
    }

    // make sure that the {oil,gas,water}Pvt() methods are available
    [[maybe_unused]] const auto& gPvt = FluidSystem::gasPvt();
    [[maybe_unused]] const auto& oPvt = FluidSystem::oilPvt();
    [[maybe_unused]] const auto& wPvt = FluidSystem::waterPvt();
}