lambda_klencki21.h

Go to the documentation of this file.

//
// Created by Giuliano Iorio on 07/12/2021.
// This headed contains the class Lambda Klencki and its derived classes.
// They are used to use the Lambda fit from the paper Klencki+21 (Appendix A,  https://ui.adsabs.harvard.edu/abs/2021A%26A...645A..54K/abstract).
//
 
#ifndef SEVN_LAMBDA_KLENCKI21_H
#define SEVN_LAMBDA_KLENCKI21_H
 
#include <cmath>
#include <iostream>
#include <vector>
#include <sevnlog.h>
#include <set>
#include <lambda_base.h>
 
class Star;
class IO;
 
class Lambda_Klencki  : public Lambda_Base {
 
public:
    Lambda_Klencki(const Star *s);
    Lambda_Klencki(const IO *io);
    virtual ~Lambda_Klencki() = default;
 
    std::vector<double> find_row(double Mzams,double Z);
 
    double operator() (const Star* s) override;
    virtual double operator() (double Mzams, double Z, double R);
 
protected:
 
     sevnstd::SevnLogging svlog;
     double Mzams_cache{0}, Z_cache{0}; 
     std::vector<double> vector_cache; 
     const double Zsun=0.017; 
     static double estimate_lambda(double R, const std::vector<double>& coefficients);
 
     static thread_local std::vector<std::vector<double>> table; 
     static thread_local bool already_loaded; 
};
 
 
 
class Lambda_Klencki_interpolator : public  Lambda_Klencki{
 
public:
    Lambda_Klencki_interpolator(const Star *s);
    Lambda_Klencki_interpolator(const IO *io);
    double operator() (double Mzams, double Z, double R) override;
 
protected:
 
    void fill_interpolators_lists();
 
    void find_interpolators(double Mzams,  double Z);
 
    static thread_local std::vector<double> Mzams_list;  
    static thread_local std::vector<std::vector<double>> Z_list;  
    static thread_local bool already_loaded_lists;  
    std::vector<std::vector<double>> vector_cache; 
private:
 
    double Mzams_interpolators [2] = {0.,0.}; 
    double Z_interpolators [4] = {0.,0.,0.,0.}; 
    double wM [2] = {0., 0.}; 
    double wZ [4] = {0., 0., 0., 0.}; 
};
 
#endif //SEVN_LAMBDA_KLENCKI21_H