Class xlifepp::OperatorOnKernel#

class OperatorOnKernel#

Collaboration diagram for xlifepp::OperatorOnKernel:

digraph { graph [bgcolor="#00000000"] node [shape=rectangle style=filled fillcolor="#FFFFFF" font=Helvetica padding=2] edge [color="#1414CE"] "35" [label="xlifepp::Collection< string_t >" tooltip="xlifepp::Collection< string_t >"] "29" [label="xlifepp::Matrix< real_t >" tooltip="xlifepp::Matrix< real_t >"] "32" [label="xlifepp::Vector< real_t >" tooltip="xlifepp::Vector< real_t >"] "42" [label="xlifepp::Vector< xlifepp::Vector< real_t > >" tooltip="xlifepp::Vector< xlifepp::Vector< real_t > >"] "44" [label="std::list< std::pair< number_t, real_t > >" tooltip="std::list< std::pair< number_t, real_t > >"] "25" [label="std::map< number_t, std::vector< number_t > >" tooltip="std::map< number_t, std::vector< number_t > >"] "24" [label="std::map< number_t, xlifepp::Geometry * >" tooltip="std::map< number_t, xlifepp::Geometry * >"] "7" [label="std::map< string_t, number_t >" tooltip="std::map< string_t, number_t >"] "11" [label="std::map< string_t, std::pair< ValueType, StrucType > >" tooltip="std::map< string_t, std::pair< ValueType, StrucType > >"] "51" [label="std::map< xlifepp::GeomElement *, std::set< number_t > >" tooltip="std::map< xlifepp::GeomElement *, std::set< number_t > >"] "45" [label="std::list< T >" tooltip="std::list< T >"] "8" [label="std::map< K, T >" tooltip="std::map< K, T >"] "6" [label="std::vector< T >" tooltip="std::vector< T >"] "38" [label="std::vector< T >" tooltip="std::vector< T >"] "31" [label="std::vector< K >" tooltip="std::vector< K >"] "21" [label="std::vector< RealPair >" tooltip="std::vector< RealPair >"] "50" [label="std::vector< const xlifepp::GeomDomain * >" tooltip="std::vector< const xlifepp::GeomDomain * >"] "34" [label="std::vector< const xlifepp::Transformation * >" tooltip="std::vector< const xlifepp::Transformation * >"] "39" [label="std::vector< number_t >" tooltip="std::vector< number_t >"] "10" [label="std::vector< real_t >" tooltip="std::vector< real_t >"] "49" [label="std::vector< std::vector< GeoNumPair > >" tooltip="std::vector< std::vector< GeoNumPair > >"] "26" [label="std::vector< std::vector< int_t > >" tooltip="std::vector< std::vector< int_t > >"] "36" [label="std::vector< string_t >" tooltip="std::vector< string_t >"] "14" [label="std::vector< xlifepp::DifferentialOperator * >" tooltip="std::vector< xlifepp::DifferentialOperator * >"] "48" [label="std::vector< xlifepp::GeomDomain * >" tooltip="std::vector< xlifepp::GeomDomain * >"] "47" [label="std::vector< xlifepp::GeomElement * >" tooltip="std::vector< xlifepp::GeomElement * >"] "5" [label="std::vector< xlifepp::Parameter * >" tooltip="std::vector< xlifepp::Parameter * >"] "23" [label="std::vector< xlifepp::Point >" tooltip="std::vector< xlifepp::Point >"] "40" [label="std::vector< xlifepp::Transformation * >" tooltip="std::vector< xlifepp::Transformation * >"] "43" [label="std::vector< xlifepp::Vector< real_t > >" tooltip="std::vector< xlifepp::Vector< real_t > >"] "20" [label="xlifepp::BoundingBox" tooltip="xlifepp::BoundingBox"] "37" [label="xlifepp::Collection< T >" tooltip="xlifepp::Collection< T >"] "13" [label="xlifepp::DifferentialOperator" tooltip="xlifepp::DifferentialOperator"] "17" [label="xlifepp::DomainInfo" tooltip="xlifepp::DomainInfo"] "15" [label="xlifepp::Extension" tooltip="xlifepp::Extension"] "27" [label="xlifepp::ExtrusionData" tooltip="xlifepp::ExtrusionData"] "3" [label="xlifepp::Function" tooltip="xlifepp::Function"] "46" [label="xlifepp::GeoNode" tooltip="xlifepp::GeoNode"] "16" [label="xlifepp::GeomDomain" tooltip="xlifepp::GeomDomain"] "19" [label="xlifepp::Geometry" tooltip="xlifepp::Geometry"] "2" [label="xlifepp::Kernel" tooltip="xlifepp::Kernel"] "12" [label="xlifepp::KernelExpansion" tooltip="xlifepp::KernelExpansion"] "30" [label="xlifepp::Matrix< K >" tooltip="xlifepp::Matrix< K >"] "18" [label="xlifepp::Mesh" tooltip="xlifepp::Mesh"] "22" [label="xlifepp::MinimalBox" tooltip="xlifepp::MinimalBox"] "1" [label="xlifepp::OperatorOnKernel" tooltip="xlifepp::OperatorOnKernel" fillcolor="#BFBFBF"] "4" [label="xlifepp::Parameters" tooltip="xlifepp::Parameters"] "41" [label="xlifepp::Parametrization" tooltip="xlifepp::Parametrization"] "9" [label="xlifepp::Point" tooltip="xlifepp::Point"] "28" [label="xlifepp::Transformation" tooltip="xlifepp::Transformation"] "33" [label="xlifepp::Vector< K >" tooltip="xlifepp::Vector< K >"] "35" -> "36" [dir=forward tooltip="public-inheritance"] "35" -> "37" [dir=forward tooltip="template-instance"] "29" -> "10" [dir=forward tooltip="public-inheritance"] "29" -> "30" [dir=forward tooltip="template-instance"] "32" -> "10" [dir=forward tooltip="public-inheritance"] "32" -> "33" [dir=forward tooltip="template-instance"] "42" -> "43" [dir=forward tooltip="public-inheritance"] "42" -> "33" [dir=forward tooltip="template-instance"] "44" -> "45" [dir=forward tooltip="template-instance"] "25" -> "8" [dir=forward tooltip="template-instance"] "24" -> "8" [dir=forward tooltip="template-instance"] "7" -> "8" [dir=forward tooltip="template-instance"] "11" -> "8" [dir=forward tooltip="template-instance"] "51" -> "8" [dir=forward tooltip="template-instance"] "31" -> "6" [dir=forward tooltip="template-instance"] "21" -> "6" [dir=forward tooltip="template-instance"] "50" -> "6" [dir=forward tooltip="template-instance"] "34" -> "6" [dir=forward tooltip="template-instance"] "39" -> "6" [dir=forward tooltip="template-instance"] "10" -> "6" [dir=forward tooltip="template-instance"] "49" -> "6" [dir=forward tooltip="template-instance"] "26" -> "6" [dir=forward tooltip="template-instance"] "36" -> "6" [dir=forward tooltip="template-instance"] "14" -> "6" [dir=forward tooltip="template-instance"] "48" -> "6" [dir=forward tooltip="template-instance"] "47" -> "6" [dir=forward tooltip="template-instance"] "5" -> "6" [dir=forward tooltip="template-instance"] "23" -> "6" [dir=forward tooltip="template-instance"] "40" -> "6" [dir=forward tooltip="template-instance"] "43" -> "6" [dir=forward tooltip="template-instance"] "20" -> "21" [dir=forward tooltip="usage"] "37" -> "38" [dir=forward tooltip="public-inheritance"] "13" -> "14" [dir=forward tooltip="usage"] "17" -> "18" [dir=forward tooltip="usage"] "15" -> "16" [dir=forward tooltip="usage"] "15" -> "51" [dir=forward tooltip="usage"] "27" -> "28" [dir=forward tooltip="usage"] "27" -> "35" [dir=forward tooltip="usage"] "27" -> "39" [dir=forward tooltip="usage"] "27" -> "10" [dir=forward tooltip="usage"] "27" -> "32" [dir=forward tooltip="usage"] "27" -> "40" [dir=forward tooltip="usage"] "27" -> "9" [dir=forward tooltip="usage"] "3" -> "4" [dir=forward tooltip="usage"] "3" -> "3" [dir=forward tooltip="usage"] "3" -> "9" [dir=forward tooltip="usage"] "3" -> "11" [dir=forward tooltip="usage"] "46" -> "46" [dir=forward tooltip="usage"] "46" -> "19" [dir=forward tooltip="usage"] "16" -> "17" [dir=forward tooltip="usage"] "16" -> "16" [dir=forward tooltip="usage"] "16" -> "19" [dir=forward tooltip="usage"] "16" -> "50" [dir=forward tooltip="usage"] "19" -> "20" [dir=forward tooltip="usage"] "19" -> "22" [dir=forward tooltip="usage"] "19" -> "24" [dir=forward tooltip="usage"] "19" -> "25" [dir=forward tooltip="usage"] "19" -> "26" [dir=forward tooltip="usage"] "19" -> "27" [dir=forward tooltip="usage"] "19" -> "41" [dir=forward tooltip="usage"] "19" -> "19" [dir=forward tooltip="usage"] "19" -> "46" [dir=forward tooltip="usage"] "2" -> "3" [dir=forward tooltip="usage"] "2" -> "2" [dir=forward tooltip="usage"] "2" -> "12" [dir=forward tooltip="usage"] "2" -> "4" [dir=forward tooltip="usage"] "2" -> "9" [dir=forward tooltip="usage"] "12" -> "3" [dir=forward tooltip="usage"] "30" -> "31" [dir=forward tooltip="public-inheritance"] "18" -> "19" [dir=forward tooltip="usage"] "18" -> "23" [dir=forward tooltip="usage"] "18" -> "47" [dir=forward tooltip="usage"] "18" -> "48" [dir=forward tooltip="usage"] "18" -> "39" [dir=forward tooltip="usage"] "18" -> "49" [dir=forward tooltip="usage"] "18" -> "18" [dir=forward tooltip="usage"] "22" -> "23" [dir=forward tooltip="usage"] "1" -> "2" [dir=forward tooltip="usage"] "1" -> "13" [dir=forward tooltip="usage"] "1" -> "15" [dir=forward tooltip="usage"] "4" -> "5" [dir=forward tooltip="usage"] "4" -> "7" [dir=forward tooltip="usage"] "41" -> "19" [dir=forward tooltip="usage"] "41" -> "32" [dir=forward tooltip="usage"] "41" -> "42" [dir=forward tooltip="usage"] "41" -> "18" [dir=forward tooltip="usage"] "41" -> "4" [dir=forward tooltip="usage"] "41" -> "44" [dir=forward tooltip="usage"] "9" -> "10" [dir=forward tooltip="public-inheritance"] "28" -> "29" [dir=forward tooltip="usage"] "28" -> "32" [dir=forward tooltip="usage"] "28" -> "34" [dir=forward tooltip="usage"] "33" -> "31" [dir=forward tooltip="public-inheritance"] }

describes the operations to process on kernel.

ker_p=0 means K(x,y)=1 !

Public Functions

inline OperatorOnKernel()#

default constructor

OperatorOnKernel(const Kernel&, DiffOpType = _id, DiffOpType = _id)#

basic constructor

OperatorOnKernel(const Kernel*, DiffOpType, DiffOpType, DiffOpType, ValueType vt = _real, StrucType st = _scalar, dimPair = dimPair(1, 1))#

basic constructor

OperatorOnKernel(const Kernel*, DiffOpType = _id, DiffOpType = _id, ValueType vt = _real, StrucType st = _scalar, dimPair = dimPair(1, 1))#

basic constructor

OperatorOnKernel(const OperatorOnKernel&)#

copy constructor

~OperatorOnKernel()#

destructor

string_t asString() const#

return as symbolic string

inline void changeKernel(Kernel *ker)#

change Kernel in operator

inline dimPair dims() const#

return the dimension of returned values

inline bool elementRequired() const#

true if geom element is required by kernel

template<typename T>
Matrix<T> &eval(const Point&, const Point&, Matrix<T>&, const Vector<real_t>* = nullptr, const Vector<real_t>* = nullptr) const#

evaluate matrix kernel operator at (x,y)

template<typename T>
T &eval(const Point&, const Point&, T&, const Vector<real_t>* = nullptr, const Vector<real_t>* = nullptr) const#

evaluate scalar kernel operator at (x,y)

evaluate operator on kernel operator at (x,y), k kernel value, nxp, nyp pointer to normal vectors nx and ny (null if not available)

template<typename T>
Vector<T> &eval(const Point&, const Point&, Vector<T>&, const Vector<real_t>* = nullptr, const Vector<real_t>* = nullptr, dimPair* = nullptr, const ExtensionData *extdata = nullptr) const#

evaluate vector kernel operator at (x,y)

general function evaluating operator on kernel, output may be either a scalar, a vector or a matrix but it is always returned as a vector scalar output: set opk[0] = out vector output: set opk[i] = out[i] 0<i<n matrix output: set opk[i*n+j]=out[i,j], 0<i<m; 0<j<n

x, y: points where current operator on kernel is evaluated opk: return value of opk as a vector (may be resized here) dims: dimension of the output (m=n=1 if scalar, n=1 if a vector), computed if pointer is not 0 nxp, nyp: pointers to normal vectors (0 if not allocated) extdata: data related to extension if required, used if not 0

template<typename T>
T &eval(const Point&, T&, const Vector<real_t>* = nullptr, const Vector<real_t>* = nullptr) const#

evaluate kernel operator at (x,P) or (P,x)

evaluate operator on kernel operator at (x,P) or (P,x) regarding xpar, xory=P k kernel value, nxp, nyp pointer to normal vectors nx and ny (null if not available)

template<typename T>
std::vector<Matrix<T>> &eval(const std::vector<Point>&, const std::vector<Point>&, std::vector<Matrix<T>>&, number_t = 0, const std::vector<Vector<real_t>>* = nullptr, const std::vector<Vector<real_t>>* = nullptr) const#

evaluate vector kernel operator at (xi,yj)

template<typename T>
std::vector<T> &eval(const std::vector<Point>&, const std::vector<Point>&, std::vector<T>&, number_t = 0, const std::vector<Vector<real_t>>* = nullptr, const std::vector<Vector<real_t>>* = nullptr) const#

evaluate kernel operator at (xi,yj)

evaluate kernel operator at (xi,yj) : op [k(x1,y1), k(x1,y2),…, k(x1,yn), k(x2,y1), k(x2,y2), …, k(xm,yn)] xs, ys: points collection nxs, nys: pointer to normal vectors nx and ny (null if not available) shift: shift index of result (by default 0)

NOTE: be careful, address of ks begins at 1 and not at 0 !!!!

template<typename T>
std::vector<Vector<T>> &eval(const std::vector<Point>&, const std::vector<Point>&, std::vector<Vector<T>>&, number_t = 0, const std::vector<Vector<real_t>>* = nullptr, const std::vector<Vector<real_t>>* = nullptr) const#

evaluate vector kernel operator at (xi,yj)

inline number_t extensionOrder() const#

derivative order in extension

inline const string_t &kernelName() const#

returns kernel shortname

inline const Kernel *kernelp() const#

return the Kernel pointer

inline bool normalRequired() const#

true if x or y normals involved

const OperatorOnFunction &operator()(const Point&, VariableName) const#

interpret operator on kernel as operator on y function

const OperatorOnFunction &operator()(VariableName) const#

interpret operator on kernel as operator on variableName function

const OperatorOnFunction &operator()(VariableName, const Point&) const#

interpret operator on kernel as operator on x function

OperatorOnKernel &operator=(const OperatorOnKernel&)#

assignment operator

void print(std::ostream&) const#

print on stream attributes

inline void setX(const Point &P) const#

set X point (usefull for Kernel)

inline void setY(const Point &P) const#

set Y point (usefull for Kernel)

inline real_t singularOrder() const#

returns singularity order (0 if not singular)

inline SingularityType singularType() const#

returns type of singularity (_notsingular, _r, _logr,_loglogr)

inline StrucType &strucType()#

return the operator returned structure

inline StrucType strucType() const#

return the operator returned structure

inline bool tangentRequired() const#

true if x or y tangents involved

inline ValueType &valueType()#

return the operator returned type

inline ValueType valueType() const#

return the operator returned type

inline bool voidKernel() const#

returns true if kernel pointer is null

inline number_t xdiffOrder() const#

order of differential operator

inline DifferentialOperator &xdifOp() const#

return reference to differential operator

inline DifferentialOperator *&xdifOp_()#

return pointer to differential operator

inline DiffOpType xdifOpType() const#

return type of x differential operator

inline bool xnormalRequired() const#

true if x-normal involved

inline bool xtangentRequired() const#

true if x-tangent involved

inline DifferentialOperator &xydifOp() const#

return reference to differential operator

inline DifferentialOperator *&xydifOp_()#

return pointer to differential operator

inline DiffOpType xydifOpType() const#

return type of y differential operator

inline number_t ydiffOrder() const#

order of differential operator

inline DifferentialOperator &ydifOp() const#

return reference to differential operator

inline DifferentialOperator *&ydifOp_()#

return pointer to differential operator

inline DiffOpType ydifOpType() const#

return type of y differential operator

inline bool ynormalRequired() const#

true if y-normal involved

inline bool ytangentRequired() const#

true if y-tangent involved

Public Members

mutable bool conjugate_#

temporary flag for conjugate operation construction

dimPair dimsRes_#

dimension of result, (0,0) if unset (NOT FULL MANAGED)

Extension *ext_p#

extension operator to apply to

bool in_ext#

if true interpreted as ext(difop(f)) else interpreted as difop(ext(f)); default is true

mutable bool noUpdatedNormal#

true means that the normal is not available in Kernel parameters

StrucType struct_#

structure of returned value (one among _scalar, _vector, _matrix)

mutable bool transpose_#

temporary flag for transpose operation construction

ValueType type_#

type of returned value (one among _real, _complex)

Friends

friend bool operator!=(const OperatorOnKernel&, const OperatorOnKernel&)#

different operator on kernel

friend bool operator==(const OperatorOnKernel&, const OperatorOnKernel&)#

same operator on kernel