TableLookup.h

Go to the documentation of this file.

 
#ifndef RIZER_TABLE_LOOKUP_H
#define RIZER_TABLE_LOOKUP_H
 
#include "numeric/Format.h"
#include "numeric/PropertyTable.h"
#include "numeric/TrialStateError.h"
 
#include <string>
 
namespace rizer
{

    [[noreturn]] inline void throwOutsideTable(const PropertyTable &table, double x,
                                               const char *procedure,
                                               const std::string &what, const char *unit,
                                               const char *remedy)
    {
        throw TrialStateError(
            procedure, what + " is evaluated at " + significant(x) + " " + unit +
                           ", outside its table [" + significant(table.first()) + ", " +
                           significant(table.last()) + "] " + unit + ". " + remedy);
    }

    inline double evalAtState(const PropertyTable &table, double x,
                              const char *procedure, const char *what,
                              const char *unit, const char *remedy)
    {
        if (!table.covers(x))
        {
            throwOutsideTable(table, x, procedure, what, unit, remedy);
        }
        return table.eval(x);
    }
 
} // namespace rizer
 
#endif // RIZER_TABLE_LOOKUP_H