LCOV - code coverage report
Current view: top level - src/utils - PorousFlowFLACrelperm.C (source / functions) Hit Total Coverage
Test: porous_flow Test Coverage Lines: 12 12 100.0 %
Date: 2017-11-21 14:47:27 Functions: 5 5 100.0 %
Legend: Lines: hit not hit

          Line data    Source code
       1             : /****************************************************************/
       2             : /* MOOSE - Multiphysics Object Oriented Simulation Environment  */
       3             : /*                                                              */
       4             : /*          All contents are licensed under LGPL V2.1           */
       5             : /*             See LICENSE for full restrictions                */
       6             : /****************************************************************/
       7             : 
       8             : #include "PorousFlowFLACrelperm.h"
       9             : 
      10             : namespace PorousFlowFLACrelperm
      11             : {
      12             : Real
      13      220214 : relativePermeability(Real seff, Real m)
      14             : {
      15      220214 :   if (seff <= 0.0)
      16             :     return 0.0;
      17      220212 :   else if (seff >= 1.0)
      18             :     return 1.0;
      19        7316 :   return (1.0 + m) * std::pow(seff, m) - m * std::pow(seff, m + 1.0);
      20             : }
      21             : 
      22             : Real
      23      220214 : dRelativePermeability(Real seff, Real m)
      24             : {
      25      220214 :   if (seff <= 0.0 || seff >= 1.0)
      26             :     return 0.0;
      27        7316 :   return (1.0 + m) * m * std::pow(seff, m - 1.0) - m * (m + 1.0) * std::pow(seff, m);
      28             : }
      29             : 
      30             : Real
      31           8 : d2RelativePermeability(Real seff, Real m)
      32             : {
      33           8 :   if (seff <= 0.0 || seff >= 1.0)
      34             :     return 0.0;
      35           4 :   return (1.0 + m) * m * (m - 1.0) * std::pow(seff, m - 2.0) -
      36           4 :          m * (m + 1.0) * m * std::pow(seff, m - 1.0);
      37             : }
      38        2499 : }

Generated by: LCOV version 1.11