123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201 |
-
-
-
-
-
-
- #ifndef QWEKEY_INC
- #define QWEKEY_INC
-
- #include "HankelTransform.h"
-
- #include <Eigen/Eigenvalues>
-
- #ifdef HAVE_BOOST_SPECIAL_FUNCTIONS
- #include "boost/math/special_functions.hpp"
- #include "boost/math/special_functions/bessel.hpp"
- #endif
-
- namespace Lemma {
-
-
-
- enum sZeroType{J0, J1, NPI};
-
-
-
- class QWEKey : public HankelTransform {
-
- friend std::ostream &operator<<(std::ostream &stream, const QWEKey &ob);
-
- public:
-
-
-
-
- QWEKey ( const ctor_key& );
-
-
- QWEKey ( const YAML::Node& node, const ctor_key& );
-
-
- ~QWEKey ();
-
-
-
- static std::shared_ptr<QWEKey> NewSP();
-
-
-
- YAML::Node Serialize() const;
-
-
-
- static std::shared_ptr< QWEKey > DeSerialize(const YAML::Node& node);
-
-
-
- void TestPrivate(const int& N);
-
-
-
- Complex Zgauss(const int &ikk, const EMMODE &imode,
- const int &itype, const Real &rho,
- const Real &wavef, KernelEM1DBase* Kernel);
-
-
-
- void ComputeRelated(const Real& rho, std::shared_ptr<KernelEM1DBase> Kernel);
-
- void ComputeRelated(const Real& rho, std::vector< std::shared_ptr<KernelEM1DBase> > KernelVec);
-
- void ComputeRelated(const Real& rho, std::shared_ptr<KernelEM1DManager> KernelManager);
-
-
-
-
-
-
- virtual std::string GetName() const ;
-
- protected:
-
-
-
-
-
- void GaussQuadWeights(const int& N);
-
-
-
- void BesselWeights( const sZeroType& sType);
-
-
-
- void QWE(const Real& rho);
-
-
-
- void getEyKernel(const int& i, const int& idx, const Real& rho);
-
- private:
-
-
-
-
- Real RelTol;
-
-
- Real AbsTol;
-
-
- int nQuad;
-
-
- int nDelay;
-
-
- int nIntervalsMax;
-
-
- VectorXr GaussWeights;
-
-
- VectorXr GaussAbscissa;
-
-
- VectorXr xInt;
-
-
- VectorXr Bx;
-
-
- VectorXr BJ0;
-
-
- VectorXr BJ1;
-
-
- VectorXr Lambda;
-
-
- VectorXr Intervals;
-
- MatrixXcr TS;
- VectorXi Tn;
- MatrixXcr Textrap;
- MatrixXr TrelErr;
- MatrixXr TabsErr;
-
-
- Eigen::Matrix<Complex, Eigen::Dynamic, Eigen::Dynamic > Zwork;
-
-
- Eigen::Matrix<Complex, Eigen::Dynamic, Eigen::Dynamic > Zans;
-
-
- std::shared_ptr<KernelEM1DManager> KernelManager;
-
-
- static constexpr auto CName = "QWEKey";
-
- };
-
- }
-
- #endif
-
|