axis1d.hpp

Go to the documentation of this file.

/* 
 * This file is part of PLaSK (https://plask.app) by Photonics Group at TUL
 * Copyright (c) 2022 Lodz University of Technology
 * 
 * This program 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, version 3.
 *
 * This program 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.
 */
#ifndef PLASK__AXIS1D_H
#define PLASK__AXIS1D_H
 
#include "mesh.hpp"
#include "interpolation.hpp"
 
#include <memory>
 
namespace plask {
 
class PLASK_API MeshAxis: public MeshD<1> {
 
public:
 
    virtual shared_ptr<MeshAxis> clone() const;
 
    virtual std::size_t findIndex(double to_find) const;
 
    virtual std::size_t findUpIndex(double to_find) const;
 
    virtual std::size_t findNearestIndex(double to_find) const;
 
    virtual shared_ptr<MeshAxis> getMidpointAxis() const;
 
    virtual bool isIncreasing() const = 0;
 
protected:
    void beforeCalcMidpointMesh() const;
};
 
PLASK_API void prepareNearestNeighborInterpolationForAxis(const MeshAxis& axis, const InterpolationFlags& flags, double& wrapped_point_coord, int axis_nr);
 
PLASK_API void prepareInterpolationForAxis(const MeshAxis& axis, const InterpolationFlags& flags, double wrapped_point_coord, int axis_nr, std::size_t& index_lo, std::size_t& index_hi, double& lo, double& hi, bool& invert_lo, bool& invert_hi);
 
 
 
struct OnePointAxis: MeshAxis {
 
    double value;
 
    OnePointAxis(double val): value(val) {}
 
    std::size_t size() const override { return 1; }
 
    double at(std::size_t) const override { return value; }
 
    bool isIncreasing() const override { return true; }
};
 
 
}   // namespace plask
 
#endif // PLASK__AXIS1D_H