We read every piece of feedback, and take your input very seriously.
To see all available qualifiers, see our documentation.
There was an error while loading. Please reload this page.
1 parent cbbfd01 commit 6b3ffa0Copy full SHA for 6b3ffa0
src/tableau_info.jl
@@ -12,15 +12,15 @@ Base.length(tab::ODERKTableau) = tab.stages
12
Calculates the stability function from the tableau at `z`. Stable if <1.
13
14
```math
15
-r(z) = 1 + z bᵀ(I - zA)⁻¹ 𝟙
+r(z) = 1 + z bᵀ(I - zA)⁻¹ e
16
```
17
-where 𝟙 denotes a vector of ones.
+where e denotes a vector of ones.
18
"""
19
function stability_region(z,tab::ODERKTableau; embedded=false)
20
A, c = tab.A, tab.c
21
b = embedded ? tab.αEEst : tab.α
22
- 𝟙 = ones(eltype(A), length(b))
23
- stages = (I - z*A) \ 𝟙
+ e = ones(eltype(A), length(b))
+ stages = (I - z*A) \ e
24
1 + z * (transpose(b) * stages)
25
end
26
0 commit comments