Skip to content

Commit

Permalink
apply formatter
Browse files Browse the repository at this point in the history
  • Loading branch information
@patrickersing
patrickersing committed Apr 25, 2024
1 parent 592a91a commit eb55c00
Showing 1 changed file with 93 additions and 73 deletions.
166 changes: 93 additions & 73 deletions src/equations/shallow_water_multilayer_2d.jl
View file
Original file line number Diff line number Diff line change
Expand Up @@ -168,85 +168,105 @@ in non-periodic domains).
# this manufactured solution velocity is taken to be constant
ω = 2 * pi * sqrt(2.0)

du1 = (-0.1 * cos(t + x[2] * ω) - 0.1 * sin(t + x[1] * ω) - 0.1 * sin(t + x[2] * ω) -
0.1 * cos(t + x[1] * ω) - 0.1 * cos(t + x[2] * ω) * ω +
0.8 * (-0.1 * sin(t + x[1] * ω) * ω - 0.1 * cos(t + x[1] * ω) * ω) -
0.1 * sin(t + x[2] * ω) * ω)
du1 = (-0.1 * cos(t + x[2] * ω) - 0.1 * sin(t + x[1] * ω) -
0.1 * sin(t + x[2] * ω) -
0.1 * cos(t + x[1] * ω) - 0.1 * cos(t + x[2] * ω) * ω +
0.8 * (-0.1 * sin(t + x[1] * ω) * ω - 0.1 * cos(t + x[1] * ω) * ω) -
0.1 * sin(t + x[2] * ω) * ω)

du2 = (0.1 * cos(t + x[2] * ω) + 0.1 * sin(t + x[1] * ω) + 0.1 * sin(t + x[2] * ω) +
0.1 * cos(t + x[1] * ω) + 0.1 * cos(t + x[2] * ω) * ω +
0.8 * (0.1 * sin(t + x[1] * ω) * ω + 0.1 * cos(t + x[1] * ω) * ω) +
0.1 * sin(t + x[2] * ω) * ω)

du3 = (-0.1 * sin(t + x[1] * ω) - 0.1 * sin(t + x[2] * ω) + 0.1 * sin(x[2] * ω) * ω +
0.8 * (0.1 * sin(x[1] * ω) * ω - 0.1 * sin(t + x[1] * ω) * ω) -
0.1 * sin(t + x[2] * ω) * ω)

du4 = (0.8 * (-0.1 * cos(t + x[2] * ω) - 0.1 * sin(t + x[1] * ω) - 0.1 * sin(t + x[2] * ω) -
0.1 * cos(t + x[1] * ω)) +
0.8 * (-0.1 * cos(t + x[2] * ω) * ω - 0.1 * sin(t + x[2] * ω) * ω) +
0.8^2 * (-0.1 * sin(t + x[1] * ω) * ω - 0.1 * cos(t + x[1] * ω) * ω) +
10.0 * (2.0 + 0.1 * cos(t + x[2] * ω) - 0.1 * sin(t + x[1] * ω) -
0.1 * sin(t + x[2] * ω) + 0.1 * cos(t + x[1] * ω)) *
(-0.1 * sin(t + x[1] * ω) * ω - 0.1 * cos(t + x[1] * ω) * ω) +
(2.0 + 0.1 * cos(t + x[2] * ω) - 0.1 * sin(t + x[1] * ω) - 0.1 * sin(t + x[2] * ω) +
0.1 * cos(t + x[1] * ω)) * cos(t + x[1] * ω) * ω)

du5 = (0.8 * (0.1 * cos(t + x[2] * ω) + 0.1 * sin(t + x[1] * ω) + 0.1 * sin(t + x[2] * ω) +
0.1 * cos(t + x[1] * ω)) +
0.8 * (0.1 * cos(t + x[2] * ω) * ω + 0.1 * sin(t + x[2] * ω) * ω) +
0.8^2 * (0.1 * sin(t + x[1] * ω) * ω + 0.1 * cos(t + x[1] * ω) * ω) +
10.0 * (0.5 - 0.1 * cos(t + x[2] * ω) + 0.1 * sin(t + x[1] * ω) +
0.1 * sin(t + x[2] * ω) - 0.1 * cos(t + x[1] * ω)) *
(0.1 * sin(t + x[1] * ω) * ω + 0.1 * cos(t + x[1] * ω) * ω) +
10.0 * (0.5 - 0.1 * cos(t + x[2] * ω) + 0.1 * sin(t + x[1] * ω) +
0.1 * sin(t + x[2] * ω) - 0.1 * cos(t + x[1] * ω)) *
(-0.1 * sin(t + x[1] * ω) * ω +
0.9 * (-0.1 * sin(t + x[1] * ω) * ω - 0.1 * cos(t + x[1] * ω) * ω)))
0.1 * cos(t + x[1] * ω) + 0.1 * cos(t + x[2] * ω) * ω +
0.8 * (0.1 * sin(t + x[1] * ω) * ω + 0.1 * cos(t + x[1] * ω) * ω) +
0.1 * sin(t + x[2] * ω) * ω)

du3 = (-0.1 * sin(t + x[1] * ω) - 0.1 * sin(t + x[2] * ω) +
0.1 * sin(x[2] * ω) * ω +
0.8 * (0.1 * sin(x[1] * ω) * ω - 0.1 * sin(t + x[1] * ω) * ω) -
0.1 * sin(t + x[2] * ω) * ω)

du4 = (0.8 * (-0.1 * cos(t + x[2] * ω) - 0.1 * sin(t + x[1] * ω) -
0.1 * sin(t + x[2] * ω) -
0.1 * cos(t + x[1] * ω)) +
0.8 * (-0.1 * cos(t + x[2] * ω) * ω - 0.1 * sin(t + x[2] * ω) * ω) +
0.8^2 * (-0.1 * sin(t + x[1] * ω) * ω - 0.1 * cos(t + x[1] * ω) * ω) +
10.0 *
(2.0 + 0.1 * cos(t + x[2] * ω) - 0.1 * sin(t + x[1] * ω) -
0.1 * sin(t + x[2] * ω) + 0.1 * cos(t + x[1] * ω)) *
(-0.1 * sin(t + x[1] * ω) * ω - 0.1 * cos(t + x[1] * ω) * ω) +
(2.0 + 0.1 * cos(t + x[2] * ω) - 0.1 * sin(t + x[1] * ω) -
0.1 * sin(t + x[2] * ω) +
0.1 * cos(t + x[1] * ω)) * cos(t + x[1] * ω) * ω)

du5 = (0.8 * (0.1 * cos(t + x[2] * ω) + 0.1 * sin(t + x[1] * ω) +
0.1 * sin(t + x[2] * ω) +
0.1 * cos(t + x[1] * ω)) +
0.8 * (0.1 * cos(t + x[2] * ω) * ω + 0.1 * sin(t + x[2] * ω) * ω) +
0.8^2 * (0.1 * sin(t + x[1] * ω) * ω + 0.1 * cos(t + x[1] * ω) * ω) +
10.0 *
(0.5 - 0.1 * cos(t + x[2] * ω) + 0.1 * sin(t + x[1] * ω) +
0.1 * sin(t + x[2] * ω) - 0.1 * cos(t + x[1] * ω)) *
(0.1 * sin(t + x[1] * ω) * ω + 0.1 * cos(t + x[1] * ω) * ω) +
10.0 *
(0.5 - 0.1 * cos(t + x[2] * ω) + 0.1 * sin(t + x[1] * ω) +
0.1 * sin(t + x[2] * ω) - 0.1 * cos(t + x[1] * ω)) *
(-0.1 * sin(t + x[1] * ω) * ω +
0.9 * (-0.1 * sin(t + x[1] * ω) * ω - 0.1 * cos(t + x[1] * ω) * ω)))

du6 = (0.8 * (-0.1 * sin(t + x[1] * ω) - 0.1 * sin(t + x[2] * ω)) +
0.8 * (0.1 * sin(x[2] * ω) * ω - 0.1 * sin(t + x[2] * ω) * ω) +
0.8^2 * (0.1 * sin(x[1] * ω) * ω - 0.1 * sin(t + x[1] * ω) * ω) +
10.0 * (0.5 - 0.1 * cos(x[1] * ω) + 0.1 * cos(t + x[2] * ω) - 0.1 * cos(x[2] * ω) +
0.1 * cos(t + x[1] * ω)) *
(0.1 * sin(x[1] * ω) * ω - 0.1 * sin(t + x[1] * ω) * ω) +
10.0 * (0.5 - 0.1 * cos(x[1] * ω) + 0.1 * cos(t + x[2] * ω) - 0.1 * cos(x[2] * ω) +
0.1 * cos(t + x[1] * ω)) * (-0.1 * sin(x[1] * ω) * ω +
0.9 / 1.1 * (-0.1 * sin(t + x[1] * ω) * ω - 0.1 * cos(t + x[1] * ω) * ω) +
1.0 / 1.1 * (0.1 * sin(t + x[1] * ω) * ω + 0.1 * cos(t + x[1] * ω) * ω)))

du7 = (-0.1 * cos(t + x[2] * ω) - 0.1 * sin(t + x[1] * ω) - 0.1 * sin(t + x[2] * ω) -
0.1 * cos(t + x[1] * ω) - 0.1 * cos(t + x[2] * ω) * ω +
0.8 * (-0.1 * sin(t + x[1] * ω) * ω - 0.1 * cos(t + x[1] * ω) * ω) -
0.1 * sin(t + x[2] * ω) * ω +
10.0 * (2.0 + 0.1 * cos(t + x[2] * ω) - 0.1 * sin(t + x[1] * ω) -
0.1 * sin(t + x[2] * ω) + 0.1 * cos(t + x[1] * ω)) *
(-0.1 * cos(t + x[2] * ω) * ω - 0.1 * sin(t + x[2] * ω) * ω) +
(2.0 + 0.1 * cos(t + x[2] * ω) - 0.1 * sin(t + x[1] * ω) - 0.1 * sin(t + x[2] * ω) +
0.1 * cos(t + x[1] * ω)) * cos(t + x[2] * ω) * ω)
0.8 * (0.1 * sin(x[2] * ω) * ω - 0.1 * sin(t + x[2] * ω) * ω) +
0.8^2 * (0.1 * sin(x[1] * ω) * ω - 0.1 * sin(t + x[1] * ω) * ω) +
10.0 *
(0.5 - 0.1 * cos(x[1] * ω) + 0.1 * cos(t + x[2] * ω) - 0.1 * cos(x[2] * ω) +
0.1 * cos(t + x[1] * ω)) *
(0.1 * sin(x[1] * ω) * ω - 0.1 * sin(t + x[1] * ω) * ω) +
10.0 *
(0.5 - 0.1 * cos(x[1] * ω) + 0.1 * cos(t + x[2] * ω) - 0.1 * cos(x[2] * ω) +
0.1 * cos(t + x[1] * ω)) *
(-0.1 * sin(x[1] * ω) * ω +
0.9 / 1.1 * (-0.1 * sin(t + x[1] * ω) * ω - 0.1 * cos(t + x[1] * ω) * ω) +
1.0 / 1.1 * (0.1 * sin(t + x[1] * ω) * ω + 0.1 * cos(t + x[1] * ω) * ω)))

du7 = (-0.1 * cos(t + x[2] * ω) - 0.1 * sin(t + x[1] * ω) -
0.1 * sin(t + x[2] * ω) -
0.1 * cos(t + x[1] * ω) - 0.1 * cos(t + x[2] * ω) * ω +
0.8 * (-0.1 * sin(t + x[1] * ω) * ω - 0.1 * cos(t + x[1] * ω) * ω) -
0.1 * sin(t + x[2] * ω) * ω +
10.0 *
(2.0 + 0.1 * cos(t + x[2] * ω) - 0.1 * sin(t + x[1] * ω) -
0.1 * sin(t + x[2] * ω) + 0.1 * cos(t + x[1] * ω)) *
(-0.1 * cos(t + x[2] * ω) * ω - 0.1 * sin(t + x[2] * ω) * ω) +
(2.0 + 0.1 * cos(t + x[2] * ω) - 0.1 * sin(t + x[1] * ω) -
0.1 * sin(t + x[2] * ω) +
0.1 * cos(t + x[1] * ω)) * cos(t + x[2] * ω) * ω)

du8 = (0.1 * cos(t + x[2] * ω) + 0.1 * sin(t + x[1] * ω) + 0.1 * sin(t + x[2] * ω) +
0.1 * cos(t + x[1] * ω) + 0.1 * cos(t + x[2] * ω) * ω +
0.8 * (0.1 * sin(t + x[1] * ω) * ω + 0.1 * cos(t + x[1] * ω) * ω) +
0.1 * sin(t + x[2] * ω) * ω +
10.0 * (0.5 - 0.1 * cos(t + x[2] * ω) + 0.1 * sin(t + x[1] * ω) +
0.1 * sin(t + x[2] * ω) - 0.1 * cos(t + x[1] * ω)) *
(0.9 * (-0.1 * cos(t + x[2] * ω) * ω - 0.1 * sin(t + x[2] * ω) * ω) -
0.1 * sin(t + x[2] * ω) * ω) +
10.0 * (0.5 - 0.1 * cos(t + x[2] * ω) + 0.1 * sin(t + x[1] * ω) +
0.1 * sin(t + x[2] * ω) - 0.1 * cos(t + x[1] * ω)) *
(0.1 * cos(t + x[2] * ω) * ω + 0.1 * sin(t + x[2] * ω) * ω))

du9 = (-0.1 * sin(t + x[1] * ω) - 0.1 * sin(t + x[2] * ω) + 0.1 * sin(x[2] * ω) * ω +
0.8 * (0.1 * sin(x[1] * ω) * ω - 0.1 * sin(t + x[1] * ω) * ω) -
0.1 * sin(t + x[2] * ω) * ω +
10.0 * (0.5 - 0.1 * cos(x[1] * ω) + 0.1 * cos(t + x[2] * ω) - 0.1 * cos(x[2] * ω) +
0.1 * cos(t + x[1] * ω)) *
(0.9 / 1.1 * (-0.1 * cos(t + x[2] * ω) * ω - 0.1 * sin(t + x[2] * ω) * ω) +
1.0 / 1.1 * (0.1 * cos(t + x[2] * ω) * ω + 0.1 * sin(t + x[2] * ω) * ω) -
0.1 * sin(x[2] * ω) * ω) +
10.0 * (0.5 - 0.1 * cos(x[1] * ω) + 0.1 * cos(t + x[2] * ω) - 0.1 * cos(x[2] * ω) +
0.1 * cos(t + x[1] * ω)) * (0.1 * sin(x[2] * ω) * ω - 0.1 * sin(t + x[2] * ω) * ω))
0.1 * cos(t + x[1] * ω) + 0.1 * cos(t + x[2] * ω) * ω +
0.8 * (0.1 * sin(t + x[1] * ω) * ω + 0.1 * cos(t + x[1] * ω) * ω) +
0.1 * sin(t + x[2] * ω) * ω +
10.0 *
(0.5 - 0.1 * cos(t + x[2] * ω) + 0.1 * sin(t + x[1] * ω) +
0.1 * sin(t + x[2] * ω) - 0.1 * cos(t + x[1] * ω)) *
(0.9 * (-0.1 * cos(t + x[2] * ω) * ω - 0.1 * sin(t + x[2] * ω) * ω) -
0.1 * sin(t + x[2] * ω) * ω) +
10.0 *
(0.5 - 0.1 * cos(t + x[2] * ω) + 0.1 * sin(t + x[1] * ω) +
0.1 * sin(t + x[2] * ω) - 0.1 * cos(t + x[1] * ω)) *
(0.1 * cos(t + x[2] * ω) * ω + 0.1 * sin(t + x[2] * ω) * ω))

du9 = (-0.1 * sin(t + x[1] * ω) - 0.1 * sin(t + x[2] * ω) +
0.1 * sin(x[2] * ω) * ω +
0.8 * (0.1 * sin(x[1] * ω) * ω - 0.1 * sin(t + x[1] * ω) * ω) -
0.1 * sin(t + x[2] * ω) * ω +
10.0 *
(0.5 - 0.1 * cos(x[1] * ω) + 0.1 * cos(t + x[2] * ω) - 0.1 * cos(x[2] * ω) +
0.1 * cos(t + x[1] * ω)) *
(0.9 / 1.1 * (-0.1 * cos(t + x[2] * ω) * ω - 0.1 * sin(t + x[2] * ω) * ω) +
1.0 / 1.1 * (0.1 * cos(t + x[2] * ω) * ω + 0.1 * sin(t + x[2] * ω) * ω) -
0.1 * sin(x[2] * ω) * ω) +
10.0 *
(0.5 - 0.1 * cos(x[1] * ω) + 0.1 * cos(t + x[2] * ω) - 0.1 * cos(x[2] * ω) +
0.1 * cos(t + x[1] * ω)) *
(0.1 * sin(x[2] * ω) * ω - 0.1 * sin(t + x[2] * ω) * ω))

return SVector(du1, du2, du3, du4, du5, du6, du7, du8, du9, zero(eltype(u)))
end
Expand Down

0 comments on commit eb55c00

Please sign in to comment.