Skip to content

Commit 046b666

Browse files
Harsh SinghHarsh Singh
authored andcommitted
fix: runic formatting for H matrix literals in Rosenbrock23/32 tableaus
1 parent 7552c49 commit 046b666

File tree

1 file changed

+8
-4
lines changed

1 file changed

+8
-4
lines changed

lib/OrdinaryDiffEqRosenbrock/src/rosenbrock_tableaus.jl

Lines changed: 8 additions & 4 deletions
Original file line numberDiff line numberDiff line change
@@ -45,8 +45,10 @@ function Rosenbrock23RodasTableau(T, T2)
4545
# Algebraically equals the MATLAB ODE Suite c1·k̃₁ + c2·k̃₂ dense output
4646
# with c1=Θ(1-Θ)/(1-2γ), c2=Θ(Θ-2γ)/(1-2γ); K[2]=0 because the Θ² term
4747
# from c2 is already absorbed into K[1] via the y₁ endpoint.
48-
H = T[zero(T) convert(T, -1/(gamma*(1 - 2*gamma))) zero(T)
49-
zero(T) zero(T) zero(T)]
48+
H = T[
49+
zero(T) convert(T, -1 / (gamma * (1 - 2 * gamma))) zero(T)
50+
zero(T) zero(T) zero(T)
51+
]
5052

5153
return RodasTableau(A, C, gamma, c, d, H, b, btilde)
5254
end
@@ -84,8 +86,10 @@ function Rosenbrock32RodasTableau(T, T2)
8486
]
8587

8688
# Same Shampine interpolation as Rosenbrock23 (identical γ, same stencil)
87-
H = T[zero(T) convert(T, -1/(gamma*(1 - 2*gamma))) zero(T)
88-
zero(T) zero(T) zero(T)]
89+
H = T[
90+
zero(T) convert(T, -1 / (gamma * (1 - 2 * gamma))) zero(T)
91+
zero(T) zero(T) zero(T)
92+
]
8993

9094
return RodasTableau(A, C, gamma, c, d, H, b, btilde)
9195
end

0 commit comments

Comments
 (0)