eulervortex.c (ce18bed930e8f3bfebcf709a18844aba97fe4630) eulervortex.c (67490bc6de48fa3615ca135c74fbc2a7c9584dee)
1// Copyright (c) 2017-2022, Lawrence Livermore National Security, LLC and other CEED contributors.
2// All Rights Reserved. See the top-level LICENSE and NOTICE files for details.
3//
4// SPDX-License-Identifier: BSD-2-Clause
5//
6// This file is part of CEED: http://github.com/ceed
7
8/// @file

--- 59 unchanged lines hidden (view full) ---

68 // Create the PETSc context
69 // ------------------------------------------------------
70 PetscScalar meter = 1e-2; // 1 meter in scaled length units
71 PetscScalar second = 1e-2; // 1 second in scaled time units
72
73 // ------------------------------------------------------
74 // Command line Options
75 // ------------------------------------------------------
1// Copyright (c) 2017-2022, Lawrence Livermore National Security, LLC and other CEED contributors.
2// All Rights Reserved. See the top-level LICENSE and NOTICE files for details.
3//
4// SPDX-License-Identifier: BSD-2-Clause
5//
6// This file is part of CEED: http://github.com/ceed
7
8/// @file

--- 59 unchanged lines hidden (view full) ---

68 // Create the PETSc context
69 // ------------------------------------------------------
70 PetscScalar meter = 1e-2; // 1 meter in scaled length units
71 PetscScalar second = 1e-2; // 1 second in scaled time units
72
73 // ------------------------------------------------------
74 // Command line Options
75 // ------------------------------------------------------
76 ierr = PetscOptionsBegin(comm, NULL, "Options for EULER_VORTEX problem",
77 NULL); CHKERRQ(ierr);
76 PetscOptionsBegin(comm, NULL, "Options for EULER_VORTEX problem", NULL);
78 // -- Physics
79 ierr = PetscOptionsScalar("-vortex_strength", "Strength of Vortex",
80 NULL, vortex_strength, &vortex_strength, NULL);
81 CHKERRQ(ierr);
82 PetscInt n = problem->dim;
83 PetscBool user_velocity;
84 ierr = PetscOptionsRealArray("-mean_velocity", "Background velocity vector",
85 NULL, mean_velocity, &n, &user_velocity);

--- 29 unchanged lines hidden (view full) ---

115 }
116 if (user_velocity && (euler_test == EULER_TEST_1
117 || euler_test == EULER_TEST_3)) {
118 ierr = PetscPrintf(comm,
119 "Warning! Background velocity vector for -euler_test t1 and -euler_test t3 is (0,0,0)\n");
120 CHKERRQ(ierr);
121 }
122
77 // -- Physics
78 ierr = PetscOptionsScalar("-vortex_strength", "Strength of Vortex",
79 NULL, vortex_strength, &vortex_strength, NULL);
80 CHKERRQ(ierr);
81 PetscInt n = problem->dim;
82 PetscBool user_velocity;
83 ierr = PetscOptionsRealArray("-mean_velocity", "Background velocity vector",
84 NULL, mean_velocity, &n, &user_velocity);

--- 29 unchanged lines hidden (view full) ---

114 }
115 if (user_velocity && (euler_test == EULER_TEST_1
116 || euler_test == EULER_TEST_3)) {
117 ierr = PetscPrintf(comm,
118 "Warning! Background velocity vector for -euler_test t1 and -euler_test t3 is (0,0,0)\n");
119 CHKERRQ(ierr);
120 }
121
123 ierr = PetscOptionsEnd(); CHKERRQ(ierr);
122 PetscOptionsEnd();
124
125 // ------------------------------------------------------
126 // Set up the PETSc context
127 // ------------------------------------------------------
128 user->units->meter = meter;
129 user->units->second = second;
130
131 // ------------------------------------------------------

--- 85 unchanged lines hidden --- 		123
124 // ------------------------------------------------------
125 // Set up the PETSc context
126 // ------------------------------------------------------
127 user->units->meter = meter;
128 user->units->second = second;
129
130 // ------------------------------------------------------

--- 85 unchanged lines hidden ---