123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210211212213214215216217218219220221222223224225226227228229230231232233234235236237238239240241242243244245246247248249250251252253254255256257258259260261262263264265266267268269270271272273274275276277278279280281282283284285286287288289290291292293294295296297298299300301302303304305306307308309310311312313314315316317318319320321322323324325326327328329330331332333334335336337338339340341342343344345346347348349350351352353354355 |
-
-
-
-
- #ifndef _HANKELTRANSFORMGAUSSIANQUADRATURE_h_INC
- #define _HANKELTRANSFORMGAUSSIANQUADRATURE_h_INC
-
- #include "HankelTransform.h"
- #include "KernelEM1DBase.h"
-
- #include "boost/math/special_functions.hpp"
-
-
- namespace Lemma {
-
-
-
-
-
-
-
-
- class GQChave : public HankelTransform {
-
- friend std::ostream &operator<<(std::ostream &stream, const GQChave &ob);
-
- struct ctor_key{};
-
- public:
-
-
-
-
- GQChave ( const ctor_key& );
-
-
- GQChave ( const YAML::Node& node, const ctor_key& );
-
-
- ~GQChave ();
-
-
-
- static std::shared_ptr<GQChave> NewSP();
-
-
-
- YAML::Node Serialize() const;
-
-
-
- static std::shared_ptr< GQChave > DeSerialize(const YAML::Node& node);
-
-
-
-
-
-
-
-
-
-
-
-
-
- Complex Zgauss(const int &ikk, const EMMODE &imode,
- const int &itype, const Real &rho,
- const Real &wavef, KernelEM1DBase* Kernel);
-
-
-
-
-
-
- virtual inline std::string GetName() const {
- return CName;
- }
-
-
-
- protected:
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
- void Besautn(Real &besr, Real &besi, const int &iorder,
- const int &nl, const int &nu, const Real &rho,
- const Real &rerr, const Real &aerr,
- const int &npcs, int &inew, const Real &aorb,
- KernelEM1DBase* Kernel);
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
- void Bestrn( Real &BESR, Real &BESI, const int &iorder,
- const int &NG, const Real &R,
- const Real &RERR, const Real &AERR, const int &npcs,
- VectorXi &XSUM, int &NSUM, int &NEW,
- int &IERR, int &NCNTRL, const Real &AORB,
- KernelEM1DBase* Kernel);
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
- void Besqud(const Real &A, const Real &B, Real &BESR, Real &BESI,
- const int &NG, const int &NEW, const int &iorder,
- const Real &R, KernelEM1DBase* Kernel);
-
-
-
-
-
-
-
-
-
-
-
-
- void Padecf(Real &SUMR, Real &SUMI, const int &N);
-
-
-
-
-
-
- void CF( Real& RESR, Real &RESI,
- Eigen::Matrix<Real, 100, 1> &CFCOR,
- Eigen::Matrix<Real, 100, 1> &CFCOI,
- const int &N);
-
-
-
-
-
-
- Real ZeroJ(const int &ZERO, const int &IORDER);
-
-
-
-
-
-
- Real Jbess(const Real &X, const int &IORDER);
-
-
-
-
-
-
- Real _dot(const int&N,
- const Eigen::Matrix<Real, Eigen::Dynamic, Eigen::Dynamic> &X1,
- const int &INC1,
- const Eigen::Matrix<Real, Eigen::Dynamic, Eigen::Dynamic> &X2,
- const int &INC2);
-
-
-
- static const Real PI2;
- static const Real X01P;
- static const Real XMAX;
- static const Real XSMALL;
- static const Real J0_X01;
- static const Real J0_X02;
- static const Real J0_X11;
- static const Real J0_X12;
- static const Real FUDGE;
- static const Real FUDGEX;
- static const Real TWOPI1;
- static const Real TWOPI2;
- static const Real RTPI2;
- static const Real XMIN;
- static const Real J1_X01;
- static const Real J1_X02;
- static const Real J1_X11;
- static const Real J1_X12;
-
-
- int HighestGaussOrder;
-
-
- int NumberPartialIntegrals;
-
-
- int NumberFunctionEvals;
-
- int np;
- int nps;
-
-
-
-
-
- static const VectorXr WT;
- static const VectorXr WA;
-
- Eigen::Matrix<int, 100, 1> Nk;
-
-
- Eigen::Matrix<Real, Eigen::Dynamic, Eigen::Dynamic> karg;
- Eigen::Matrix<Real, Eigen::Dynamic, Eigen::Dynamic> kern;
-
-
- Eigen::Matrix<Real, 100, 1> Xr;
- Eigen::Matrix<Real, 100, 1> Xi;
- Eigen::Matrix<Real, 100, 1> Dr;
- Eigen::Matrix<Real, 100, 1> Di;
- Eigen::Matrix<Real, 100, 1> Sr;
- Eigen::Matrix<Real, 100, 1> Si;
- Eigen::Matrix<Real, 100, 1> Cfcor;
- Eigen::Matrix<Real, 100, 1> Cfcoi;
-
- private:
-
-
- static constexpr auto CName = "FHTKey51";
-
- };
-
-
-
-
-
-
- class LowerGaussLimitGreaterThanUpperGaussLimit :
- public std::runtime_error {
-
-
- public: LowerGaussLimitGreaterThanUpperGaussLimit();
- };
-
- }
-
- #endif
|