text: /* MIT License
   *
   * Copyright (c) 2020 Aleksandr Zhuravlyov and Zakhar Lanets
   *
   * Permission is hereby granted, free of charge, to any person obtaining a copy
   * of this software and associated documentation files (the "Software"), to deal
   * in the Software without restriction, including without limitation the rights
   * to use, copy, modify, merge, publish, distribute, sublicense, and/or sell
   * copies of the Software, and to permit persons to whom the Software is
   * furnished to do so, subject to the following conditions:
   *
   * The above copyright notice and this permission notice shall be included in all
   * copies or substantial portions of the Software.
   *
   * THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
   * IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
   * FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
   * AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
   * LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
   * OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE
   * SOFTWARE.
   */
  
  #include "Local.h"
  #include "funcs.h"
  #include <algorithm>
  
  Local::Local(std::shared_ptr<Props> props, std::shared_ptr<Sgrid> sgrid) :
          _props(props),
          _sgrid(sgrid),
          _alphas(_sgrid->_cellsN, 0) {}
  
  void Local::calcTimeSteps() {
  
      auto time = std::get<double>(_props->_params["time_period"]);
      auto timeStep = std::get<double>(_props->_params["time_step"]);
      double division = time / timeStep;
      double fullStepsN;
      auto lastStep = std::modf(division, &fullStepsN);
  
      _timeSteps.clear();
      _timeSteps = std::vector<double>(fullStepsN, timeStep);
      if (lastStep > 0)
          _timeSteps.push_back(lastStep * timeStep);
  }
  
  void Local::calcAlphas(Eigen::Ref<Eigen::VectorXd> concs,
                         const double &timeStep) {
  
      auto poroIni = std::get<double>(_props->_params["poro"]);
  
      for (int i = 0; i < _alphas.size(); i++) {
          auto aCoeff = calcAFunc(concs[i], poroIni);
          _alphas[i] = aCoeff * _sgrid->_cellV / timeStep;
      }
  }