removing_nan.hpp

Go to the documentation of this file.

// Copyright (c) 2024 Andrea Enrico Franzoni (andreaenrico.franzoni@gmail.com)
//
// This file is part of PPCKO
//
// Permission is hereby granted, free of charge, to any person obtaining a copy
// of PPCKO and associated documentation files (the PPCKO software), to deal
// PPCKO without restriction, including without limitation the rights
// to use, copy, modify, merge, publish, distribute, sublicense, and/or sell
// copies of PPCKO, and to permit persons to whom PPCKO is
// furnished to do so, subject to the following conditions:
//
// PPCKO 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 PPCKO OR THE USE OR OTHER DEALINGS IN
// PPCKO.
 
#ifndef KO_REMOVE_NAN_HPP
#define KO_REMOVE_NAN_HPP
 
#include <Eigen/Dense>
#include <iostream>
#include <algorithm>
#include <iterator>
#include <numeric>
#include <functional>
#include <vector>
#include <set>
#include <type_traits>
#include <string>
#include <stdexcept>
 
#include "parameters_wrapper.hpp"
#include "traits_ko.hpp"
 
 

 
 

template <REM_NAN MA_t>
using MAT = std::integral_constant<REM_NAN, MA_t>;
 

template<typename T,REM_NAN MA_t>
class removing_nan
{
  
private:
  KO_Traits::StoringMatrix m_data;
  std::size_t m_m;
  std::size_t m_n;

  void row_removal(Eigen::Block<Eigen::Matrix<T,-1,-1>,1>& row, MAT<REM_NAN::MR>);
  void row_removal(Eigen::Block<Eigen::Matrix<T,-1,-1>,1>& row, MAT<REM_NAN::ZR>);
  
public:
  
  template<typename STOR_OBJ>
  removing_nan(STOR_OBJ&& data)
    :
    m_data{std::forward<STOR_OBJ>(data)}
    { 
      m_m = m_data.rows();
      m_n = m_data.cols();
    }
  
  inline KO_Traits::StoringMatrix data() const {return m_data;};
  
  void row_removal(Eigen::Block<Eigen::Matrix<T,-1,-1>,1>& row) { return row_removal(row, MAT<MA_t>{});};
  
  inline void remove_nan(){   for(auto row : m_data.rowwise()){   row_removal(row);} };
};
 
 
#include "removing_nan_imp.hpp"
#include "removing_nan_cleaner_imp.hpp"
 
#endif /*KO_REMOVE_NAN_HPP*/