Skip to content

Commit

Permalink
Apply suggestions from code review
Browse files Browse the repository at this point in the history
Co-authored-by: Joshua Lampert <[email protected]>
  • Loading branch information
DanielDoehring and JoshuaLampert authored Dec 19, 2024
1 parent 1001f65 commit e9c988e
Showing 1 changed file with 3 additions and 3 deletions.
6 changes: 3 additions & 3 deletions ext/TrixiConvexECOSExt.jl
Original file line number Diff line number Diff line change
Expand Up @@ -90,12 +90,12 @@ function stability_polynomials_PERK4!(pnoms, num_stage_evals,

# "Fixed" term due to choice of the PERK4 Butcher tableau
# Required to un-do the normalization of the eigenvalues here
pnoms += k1 * dt^5 * normalized_powered_eigvals[:, 5] * factorial(5)
pnoms += k1 * dt^5 * view(normalized_powered_eigvals, :, 5) * factorial(5)

# Contribution from free coefficients
for k in 1:(num_stage_evals - 5)
pnoms += (k2 * dt^(k + 4) * normalized_powered_eigvals[:, k + 4] * gamma[k] +
k1 * dt^(k + 5) * normalized_powered_eigvals[:, k + 5] * gamma[k] *
pnoms += (k2 * dt^(k + 4) * view(normalized_powered_eigvals, :, k + 4) * gamma[k] +
k1 * dt^(k + 5) * view(normalized_powered_eigvals, :, k + 5) * gamma[k] *
(k + 5))
end

Expand Down

0 comments on commit e9c988e

Please sign in to comment.