eulervortex.c (71bbf5a7c845c47c2193d271e8411746270b9cb1) eulervortex.c (a2f6637ed9f7b99e1cdac198632ca3cb6dbe41d4)
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

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

18
19PetscErrorCode NS_EULER_VORTEX(ProblemData *problem, DM dm, void *ctx, SimpleBC bc) {
20 EulerTestType euler_test;
21 User user = *(User *)ctx;
22 StabilizationType stab;
23 MPI_Comm comm = user->comm;
24 Ceed ceed = user->ceed;
25 PetscBool implicit;
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

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

18
19PetscErrorCode NS_EULER_VORTEX(ProblemData *problem, DM dm, void *ctx, SimpleBC bc) {
20 EulerTestType euler_test;
21 User user = *(User *)ctx;
22 StabilizationType stab;
23 MPI_Comm comm = user->comm;
24 Ceed ceed = user->ceed;
25 PetscBool implicit;
26 PetscBool has_curr_time = PETSC_TRUE;
27 PetscBool has_neumann = PETSC_TRUE;
28 EulerContext euler_ctx;
29 CeedQFunctionContext euler_context;
30
31 PetscFunctionBeginUser;
32 PetscCall(PetscCalloc1(1, &euler_ctx));
33
34 // ------------------------------------------------------
35 // SET UP DENSITY_CURRENT

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

118 for (PetscInt i = 0; i < 3; i++) {
119 center[i] *= meter;
120 domain_size[i] *= meter;
121 mean_velocity[i] *= (meter / second);
122 }
123 problem->dm_scale = meter;
124
125 // -- QFunction Context
26 EulerContext euler_ctx;
27 CeedQFunctionContext euler_context;
28
29 PetscFunctionBeginUser;
30 PetscCall(PetscCalloc1(1, &euler_ctx));
31
32 // ------------------------------------------------------
33 // SET UP DENSITY_CURRENT

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

116 for (PetscInt i = 0; i < 3; i++) {
117 center[i] *= meter;
118 domain_size[i] *= meter;
119 mean_velocity[i] *= (meter / second);
120 }
121 problem->dm_scale = meter;
122
123 // -- QFunction Context
126 user->phys->stab = stab;
127 user->phys->euler_test = euler_test;
128 user->phys->implicit = implicit;
124 user->phys->implicit = implicit;
129 user->phys->has_curr_time = has_curr_time;
130 user->phys->has_neumann = has_neumann;
131 euler_ctx->curr_time = 0.;
132 euler_ctx->implicit = implicit;
133 euler_ctx->euler_test = euler_test;
134 euler_ctx->center[0] = center[0];
135 euler_ctx->center[1] = center[1];
136 euler_ctx->center[2] = center[2];
137 euler_ctx->vortex_strength = vortex_strength;
138 euler_ctx->c_tau = c_tau;

--- 37 unchanged lines hidden --- 		125 euler_ctx->curr_time = 0.;
126 euler_ctx->implicit = implicit;
127 euler_ctx->euler_test = euler_test;
128 euler_ctx->center[0] = center[0];
129 euler_ctx->center[1] = center[1];
130 euler_ctx->center[2] = center[2];
131 euler_ctx->vortex_strength = vortex_strength;
132 euler_ctx->c_tau = c_tau;

--- 37 unchanged lines hidden ---