From eb55c008eaff63b0d49026e1157cac842c426a21 Mon Sep 17 00:00:00 2001 From: patrickersing Date: Thu, 25 Apr 2024 13:27:36 +0200 Subject: [PATCH] apply formatter --- src/equations/shallow_water_multilayer_2d.jl | 166 +++++++++++-------- 1 file changed, 93 insertions(+), 73 deletions(-) diff --git a/src/equations/shallow_water_multilayer_2d.jl b/src/equations/shallow_water_multilayer_2d.jl index ab3009c..61dfaad 100644 --- a/src/equations/shallow_water_multilayer_2d.jl +++ b/src/equations/shallow_water_multilayer_2d.jl @@ -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