DiracHydrogen.hpp

#pragma once
#include "qip/StrongType.hpp"
 
 
namespace DiracHydrogen {
 
//------------------------------------------------------------------------------
// Uses Strong Types:
enum class DiracTypes { DiracQN, AlphaFS, Zeff, PrincipalQN, RaB };
 
using DiracQN = qip::StrongType<DiracTypes::DiracQN, int>;
using AlphaFS = qip::StrongType<DiracTypes::AlphaFS, double>;
using Zeff = qip::StrongType<DiracTypes::Zeff, double>;
// double (allow eff):
using PrincipalQN = qip::StrongType<DiracTypes::PrincipalQN, double>;
using RaB = qip::StrongType<DiracTypes::RaB, double>;
 
//------------------------------------------------------------------------------
 
double gamma(DiracQN k, Zeff z, AlphaFS a);
 
double enk(PrincipalQN n, DiracQN k, Zeff z, AlphaFS a);
 
double Enk(PrincipalQN n, DiracQN k, Zeff z, AlphaFS a);
 
double f(RaB r, PrincipalQN n, DiracQN k, Zeff z, AlphaFS a);
 
double g(RaB r, PrincipalQN n, DiracQN k, Zeff z, AlphaFS a);
 
double gfratio(double r, int k, double z, double a, double e, double m = 1.0);
 
} // namespace DiracHydrogen