Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Definition of stieltjes? #39

Open
putianyi889 opened this issue Mar 22, 2024 · 1 comment
Open

Definition of stieltjes? #39

putianyi889 opened this issue Mar 22, 2024 · 1 comment

Comments

@putianyi889
Copy link

The transform is defined on $\mathbb{R}$ and in terms of quasimatrices it's an $\mathbb{R}\times I$ kernel multiplied by an $I\times\mathbb{N}$ basis, so the result should be an $\mathbb{R}\times\mathbb{N}$ basis.

@dlfivefifty
Copy link
Member

dlfivefifty commented Apr 20, 2024

When you use the full broadcast notation it's clear what the domain is:

W = Weighted(ChebyshevT())
x = axes(W,1)
S = inv.(x .- x') * W
@test axes(S) == (Inclusion(-1..1), 1:∞)

So I guess you are referring to the behaviour of stieltjes(W)? I think the convention is that it evaluates on the domain axes(W,1).

This may also be true for complex supported domains but I've just realised if I write

Γ = somecomplexcontour
W = Weighted(chebyshevt(Γ))
z = axes(W,1)
S = inv.(z .- z') * W

The definition follows from the definition of the adjoint/inner product. Which would only sensibly be the $\int_Γ \bar f(z) g(z) ds$ where we would also want a version with $dz$. Probably we would use stieltjes to mean the latter and cauchy (with the $1/(2πi)$ to mean the former.

If you indeed want the whole real line we can support the syntax

W = Weighted(ChebyshevT())
x = Inclusion(ℝ)
t = axes(W,1)
S = inv.(x .- t') * W
@test axes(S) == (x, 1:∞)

which could lower to a call steiltjes(W, ℝ). Would take some thought how to implement this.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants