crystalfields.h

#ifndef CRYSTALFIELDS_H
#define CRYSTALFIELDS_H
#define _USE_MATH_DEFINES

#include <cmath>
#include "array.h"
#include <gsl/gsl_rng.h>
class CrystalFields
{
public:
    CrystalFields(double a, double c, int na, int nc);
    double B_FeSeIntercalated_fm(double m, int cycles);
    double B_FeSeIntercalated_fm_Li(double m, int cycles);
    double* Vzz_FeSe(double q, double alpha, double SeVacancies, int cycles);
    double* Vzz_FeSe_excessFe(double q, double alpha,
                              double SeVacancies, int cycles,
                              double FeExcess, double zFeExcess);
private:
    void add_dipol_B(double x, double y, double z, double phi, double theta, double m, double *B);
    void add_point_charge(double x, double y, double z, double q, double *Vzz);
    double a,c;
    int Na,Nc;
    double mu_0;
    const gsl_rng_type * T;
    gsl_rng * r;

};

#endif // CRYSTALFIELDS_H