transport.hpp

#ifndef __GET_CUMULATIVE_PROBABILITY_HPP__
#define __GET_CUMULATIVE_PROBABILITY_HPP__
 
#ifndef NDEBUG
#pragma GCC diagnostic push
#pragma GCC diagnostic ignored "-Wunused-but-set-variable"
#endif 
#include <Eigen/Core>
#ifndef NDEBUG
#pragma GCC diagnostic pop
#endif 
#include <Kokkos_Core.hpp>
#include <cma_read/flowmap.hpp>
#include <cma_read/light_2d_view.hpp>
#include <cma_read/neighbors.hpp>
#include <cma_utils/cache_hydro_state.hpp>
#include <vector>
 
namespace CmaUtils
{
  using CumulativeProbaType = Eigen::Matrix<double, -1, -1, Eigen::RowMajor>;
  CumulativeProbaType
  get_cumulative_probabilities(CmaRead::Neighbors::Neighbors_const_view_t neighbors,
                               const FlowMatrixType& m_transition);

  FlowMatrixType get_transition_matrix(const CmaRead::FlowMap::FlowMap_const_view_t& flows);

  std::vector<double> get_diag_transition(const FlowMatrixType& m_transition);
 
} // namespace CmaUtils
#endif //__GET_CUMULATIVE_PROBABILITY_HPP__