InhomogenousGreens.hpp

#pragma once
#include <vector>
class DiracSpinor;
 
namespace DiracODE {
 
//==============================================================================
 
 
DiracSpinor solve_inhomog(const int kappa, const double en,
                          const std::vector<double> &v,
                          const std::vector<double> &H_mag, const double alpha,
                          const DiracSpinor &source,
                          const DiracSpinor *const VxFa = nullptr,
                          const DiracSpinor *const Fa0 = nullptr,
                          double zion = 1, double mass = 1.0);
 
 
void solve_inhomog(DiracSpinor &Fa, const double en,
                   const std::vector<double> &v,
                   const std::vector<double> &H_mag, const double alpha,
                   const DiracSpinor &source,
                   const DiracSpinor *const VxFa = nullptr,
                   const DiracSpinor *const Fa0 = nullptr, double zion = 1,
                   double mass = 1.0);
 
 
void solve_inhomog(DiracSpinor &Fa, DiracSpinor &Fzero, DiracSpinor &Finf,
                   const double en, const std::vector<double> &v,
                   const std::vector<double> &H_mag, const double alpha,
                   const DiracSpinor &source,
                   const DiracSpinor *const VxFa = nullptr,
                   const DiracSpinor *const Fa0 = nullptr, double zion = 1,
                   double mass = 1.0);
 
//==============================================================================
 
namespace Internal {
 
// Takes solution regular at infinity (Finf), and that regular at zero (Fzero),
// and the inhomogenous source term, Sr, to find particular solution, Fa.
void GreenSolution(DiracSpinor &Fa, const DiracSpinor &Finf,
                   const DiracSpinor &Fzero, const double alpha,
                   const DiracSpinor &Sr);
 
} // namespace Internal
} // namespace DiracODE