30 #include <boost/numeric/ublas/vector.hpp> 32 #include <o2scl/constants.h> 33 #include <o2scl/multi_funct.h> 34 #include <o2scl/mmin.h> 35 #include <o2scl/mmin_simp2.h> 36 #include <o2scl/nucmass.h> 37 #include <o2scl/nucmass_ame.h> 38 #include <o2scl/nucdist.h> 40 #ifndef DOXYGEN_NO_O2NS 124 O2SCL_ERR2(
"Tried to give zero uncertainties in nucmass_fit::",
138 ubvector &qual,
int max_iso=20);
143 ubvector &qual,
int min_iso=-8,
int max_iso=60);
147 virtual double min_fun(
size_t nv,
const ubvector &x);
151 #ifndef DOXYGEN_NO_O2NS 169 #ifndef DOXYGEN_NO_O2NS void set_uncerts(vec_t &u)
Set the fit uncertainties (in MeV)
static const int chi_squared_be
Chi-squared for binding energy using specified uncertainties.
mmin_base * mm
The pointer to the minimizer.
virtual void eval(nucmass &n, double &res)
Evaluate quality without fitting.
Nuclear mass formula base [abstract base].
static const int chi_squared_me
Chi-squared for mass excess using specified uncertainties.
Fit a nuclear mass formula.
int minZ
Minimum proton number to fit (default 8)
void set_uncerts(size_t nv, vec_t &u)
Set the fit uncertainties (in MeV) from the first nv elements of u.
ubvector uncs
Uncertainties.
bool even_even
If true, then only fit doubly-even nuclei (default false)
void set_mmin(mmin_base<> &umm)
Change the minimizer for use in the fit.
virtual void fit(nucmass_fit_base &n, double &res)
Fit the nuclear mass formula.
void vector_copy(const vec_t &src, vec2_t &dest)
nucmass_fit_base * nmf
The nuclear mass formula to fit to.
static const int rms_binding_energy
RMS deviation in binding_energy.
#define O2SCL_ERR2(d, d2, n)
mmin_simp2 def_mmin
The default minimizer.
void eval_isospin_beta(nucmass &n, ubvector_int &n_qual, ubvector &qual, int max_iso=20)
Evaluate isospin dependence of fit quality.
static const int rms_mass_excess
RMS deviation in mass excess.
virtual double min_fun(size_t nv, const ubvector &x)
The function to minimize.
void eval_isospin(nucmass &n, ubvector_int &n_qual, ubvector &qual, int min_iso=-8, int max_iso=60)
Evaluate isospin dependence of fit quality.
int minN
Minimum neutron number to fit (default 8)
Fittable mass formula [abstract base].
int fit_method
Current fitting method.
std::vector< nucleus > dist
Select the experimental nuclei to fit.