Subdivision.h

Go to the documentation of this file.

// SPDX-FileCopyrightText: Copyright (c) OpenGeoSys Community (opengeosys.org)
// SPDX-License-Identifier: BSD-3-Clause
 
#pragma once
 
#include <cstddef>
#include <vector>
 
namespace BaseLib
{
class ISubdivision
{
public:
    virtual std::vector<double> operator()() const = 0;
 
    virtual ~ISubdivision() = default;
};

class UniformSubdivision : public ISubdivision
{
public:
    UniformSubdivision(double length, std::size_t n_subdivision)
        : length_(length), n_subdivision_(n_subdivision)
    {
    }

    std::vector<double> operator()() const override
    {
        std::vector<double> x;
        x.reserve(n_subdivision_ + 1);
        const double dL = length_ / static_cast<double>(n_subdivision_);
        for (std::size_t i = 0; i < n_subdivision_ + 1; i++)
        {
            x.push_back(i * dL);
        }
        return x;
    }
 
private:
    const double length_;
    const std::size_t n_subdivision_;
};

class GradualSubdivision : public ISubdivision
{
public:
    GradualSubdivision(const double L,
                       const double dL0,
                       const double max_dL,
                       const double multiplier);

    std::vector<double> operator()() const override;
 
private:
    const double length_;
    const double dL0_;
    const double max_dL_;
    const double multiplier_;
};

class GradualSubdivisionFixedNum : public ISubdivision
{
public:
    GradualSubdivisionFixedNum(const double L,
                               const std::size_t num_subdivisions,
                               const double multiplier);

    std::vector<double> operator()() const override;
 
private:
    const double length_;
    const std::size_t num_subdivisions_;
    const double multiplier_;
};
 
}  // namespace BaseLib