Skip to content

Commit 05d85d2

Browse files
authored
Update jacobi-SH internal doc (#838)
1 parent 77cec2d commit 05d85d2

File tree

1 file changed

+8
-4
lines changed

1 file changed

+8
-4
lines changed

docs/src/internals/multivariate.md

Lines changed: 8 additions & 4 deletions
Original file line numberDiff line numberDiff line change
@@ -14,14 +14,18 @@ In the example above, the resulting space is a [`TensorSpace`](@ref).
1414

1515
## AbstractProductSpace
1616

17-
A `TensorSpace` is a subtype of `AbstractProductSpace`. The purpose of `AbstractProductSpace`s beyond `TensorSpace` is that it is often convenient to use different bases in one-dimension depending on the order. Thus we want to be able to represent functions in the basis `φ_k^j(x)ζ_j(y)`: for example, we could have the basis `(1-x^2)^j P_k^{(j,j)}(x) e^{ijy}` which is related to spherical harmonics, where `P_k^{(a,b)}` are the Jacobi polynomials.
17+
A `TensorSpace` is a subtype of `AbstractProductSpace`. The purpose of `AbstractProductSpace`s beyond `TensorSpace` is that it is often convenient to use different bases in one-dimension depending on the order. Thus we want to be able to represent functions in the basis `φ_k^j(x)ζ_j(y)`: for example, we could have the basis
18+
```math
19+
(1-x^2)^{m/2} \mathrm{P}_k^{(m,m)}(x) e^{imy}
20+
```
21+
(which is related to spherical harmonics), where ``\mathrm{P}_k^{(m,m)}(x)`` are Jacobi polynomials.
1822

19-
To handle this more general setting, an `AbstractProductSpace` implements `ApproxFunBase.columnspace`: In the `(1-x^2)^j P_k^{(j,j)}(x) e^{ijy}` example, we would have
23+
To handle this more general setting, an `AbstractProductSpace` implements `ApproxFunBase.columnspace`: In the `(1-x^2)^(m/2) P_k^{(m,m)}(x) e^{imy}` example, we would have
2024

2125
```julia
2226
columnspace(myproductspace,1) == JacobiWeight(0,0,Legendre()),
23-
columnspace(myproductspace,2) == JacobiWeight(1,1,Jacobi(1,1)),
24-
columnspace(myproductspace,k) == JacobiWeight(k-1,k-1,Jacobi(k-1,k-1))
27+
columnspace(myproductspace,2) == JacobiWeight(0.5,0.5,Jacobi(1,1)),
28+
columnspace(myproductspace,k) == JacobiWeight((k-1)/2,(k-1)/2,Jacobi(k-1,k-1))
2529
```
2630

2731
The basis is then interlaced as in `TensorSpace`.

0 commit comments

Comments
 (0)