Class xlifepp::DuffyIM#

class DuffyIM : public xlifepp::DoubleIM #

integral over a product of 2D geometric elements with singularity using

for adjacent elements : a method based on Duffy transform
for self influence : an hybrid method based on t^p transform note : use 1D quadrature with quadrature points != 0.5, for instance Gauss-Legendre of order 4k+1 may manage different quadrature rules Elements must be segment! Note : does not address the case of separate elements

Public Functions

inline DuffyIM(const Quadrature &q)#: full constructor from a quadrature object

inline DuffyIM(IntegrationMethodType imt)#: basic constructor

inline DuffyIM(number_t os = 5, number_t oa = 5)#: full constructor from quadrature order

inline DuffyIM(QuadRule qsX, number_t osX, QuadRule qsY, number_t osY, QuadRule qaX, number_t oaX, QuadRule qaY, number_t oaY, number_t so = 5, bool us = true)#: full constructor from quadrules

inline virtual void print(std::ostream &os) const#: print IntegrationMethod on stream

inline virtual std::list<Quadrature*> quadratures() const#: return the list of (single) quadratures in a list

Public Members

number_t ordAdjtY#: order of quadratures on the segment for adjacent elements

number_t ordSelfY#: order of quadratures on the segment for self influence

Quadrature *quadAdjtY#: quadratures on segment [0,1] for adjacent elements

Quadrature *quadSelfY#: quadratures on segment [0,1] for self influence

number_t selfOrd#: order of t^p transform or number of diminishing layers (def=5)

bool useSelfTP#: use t^p improved method for self influence (def=true)