EpsilonScalingHandler.h

Go to the documentation of this file.

#ifndef EPSILON_SCALING_HANDLER_H
#define EPSILON_SCALING_HANDLER_H
 
#include <deal.II/base/conditional_ostream.h>
#include <deal.II/base/mpi.h>
#include <vector>
#include <string>
#include <memory>
#include <cmath>
 
using namespace dealii;
 
class EpsilonScalingHandler {
public:
    EpsilonScalingHandler(const MPI_Comm& comm,
                         double initial_epsilon,
                         double scaling_factor,
                         unsigned int num_steps);
 
    std::vector<std::vector<double>> compute_epsilon_distribution(
        unsigned int num_levels);
 
    const std::vector<double>& get_epsilon_values_for_level(unsigned int level_index) const;
 
    void print_epsilon_distribution() const;
 
private:
    MPI_Comm mpi_communicator; 
    const unsigned int this_mpi_process; 
    ConditionalOStream pcout; 
 
    double initial_epsilon; 
    double scaling_factor; 
    unsigned int num_steps; 
 
    std::vector<std::vector<double>> epsilon_distribution; 
 
    std::vector<double> generate_epsilon_sequence() const;
 
    std::vector<std::vector<double>> distribute_epsilon_values(
        const std::vector<double>& epsilon_sequence,
        unsigned int num_levels);
};
 
#endif // EPSILON_SCALING_HANDLER_H