instream_profile_guesser.hpp

Go to the documentation of this file.

/*
 * This file is part of BubbleProfiler.
 *
 * BubbleProfiler 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 3 of the License, or
 * (at your option) any later version.
 *
 * BubbleProfiler 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 BubbleProfiler.  If not, see <http://www.gnu.org/licenses/>.
 */

#ifndef BUBBLEPROFILER_INSTREAM_PROFILE_GUESSER_HPP_INCLUDED
#define BUBBLEPROFILER_INSTREAM_PROFILE_GUESSER_HPP_INCLUDED

#include "basic_logger.hpp"
#include "profile_guesser.hpp"

#include <Eigen/Core>

#include <istream>
#include <string>
#include <vector>

namespace BubbleProfiler {

class Instream_profile_guesser: public Profile_guesser {
public:
   explicit Instream_profile_guesser(std::istream&);

   Instream_profile_guesser(std::istream&, const std::string&);

   Instream_profile_guesser(std::istream&, const std::string&,
                            const std::string&);
   virtual ~Instream_profile_guesser() = default;
   Instream_profile_guesser(const Instream_profile_guesser&) = default;
   Instream_profile_guesser& operator=(const Instream_profile_guesser&) = default;
   Instream_profile_guesser(Instream_profile_guesser&&) = default;
   Instream_profile_guesser& operator=(Instream_profile_guesser&&) = default;

   virtual Field_profiles get_profile_guess(
      const Potential&, const Eigen::VectorXd&, int,
      double, double, double, double) override;

   std::size_t get_number_of_grid_points() const { return spatial_grid.size(); }

   std::size_t get_number_of_fields() const { return field_values.size(); }

   Eigen::VectorXd get_spatial_grid() const;

   Eigen::MatrixXd get_field_values() const;

   void set_spatial_tolerance(double t) { tolerance = t; }

   void set_delimiters(const std::string& d) { delimiters = d; };

   void set_line_comment(const std::string& c) { line_comment_start = c; };

   void read_profiles(std::istream&);

private:
   std::string delimiters{" \t"};
   std::string line_comment_start{"#"};
   bool has_n_fields{false};
   std::size_t n_fields{0};
   double tolerance{1.e-6};
   std::vector<double> spatial_grid{};
   std::vector<std::vector<double> > field_values{};
   logging::Basic_logger logger{};

   void reset();
   void strip_comments(std::string&) const;
   void process_line(const std::string&);
};

} // namespace BubbleProfiler

#endif