libceedsetup.c (a171b6ef1192a8f5f3031330339aaab2ded4c8fb) libceedsetup.c (4d00b080eb3f95d2e04e55c0ff369c5c847bb288)
1// Copyright (c) 2017-2024, 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#include "../include/libceedsetup.h"

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

45 DM dm_coord;
46 Vec coords;
47 const PetscScalar *coord_array;
48 CeedBasis basis_x, basis_u;
49 CeedElemRestriction elem_restr_x, elem_restr_u, elem_restr_u_i, elem_restr_qd_i;
50 CeedQFunction qf_setup_geo, qf_apply;
51 CeedOperator op_setup_geo, op_apply;
52 CeedVector x_coord, q_data, x_ceed, y_ceed;
1// Copyright (c) 2017-2024, 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#include "../include/libceedsetup.h"

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

45 DM dm_coord;
46 Vec coords;
47 const PetscScalar *coord_array;
48 CeedBasis basis_x, basis_u;
49 CeedElemRestriction elem_restr_x, elem_restr_u, elem_restr_u_i, elem_restr_qd_i;
50 CeedQFunction qf_setup_geo, qf_apply;
51 CeedOperator op_setup_geo, op_apply;
52 CeedVector x_coord, q_data, x_ceed, y_ceed;
53 CeedInt num_qpts, c_start, c_end, num_elem, q_data_size = bp_data.q_data_size;
53 PetscInt c_start, c_end, num_elem;
54 CeedInt num_qpts, q_data_size = bp_data.q_data_size;
54 CeedScalar R = 1; // radius of the sphere
55 CeedScalar l = 1.0 / PetscSqrtReal(3.0); // half edge of the inscribed cube
56
57 PetscFunctionBeginUser;
58 PetscCall(DMGetCoordinateDM(dm, &dm_coord));
59
60 // CEED bases
61 PetscCall(CreateBasisFromPlex(ceed, dm_coord, 0, 0, 0, 0, bp_data, &basis_x));

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

217 DM dm_coord;
218 Vec coords;
219 const PetscScalar *coord_array;
220 CeedBasis basis_x, basis_u;
221 CeedElemRestriction elem_restr_x, elem_restr_u, elem_restr_u_i, elem_restr_qd_i;
222 CeedQFunction qf_setup_geo, qf_setup_rhs, qf_error;
223 CeedOperator op_setup_geo, op_setup_rhs;
224 CeedVector x_coord, q_data, target, rhs;
55 CeedScalar R = 1; // radius of the sphere
56 CeedScalar l = 1.0 / PetscSqrtReal(3.0); // half edge of the inscribed cube
57
58 PetscFunctionBeginUser;
59 PetscCall(DMGetCoordinateDM(dm, &dm_coord));
60
61 // CEED bases
62 PetscCall(CreateBasisFromPlex(ceed, dm_coord, 0, 0, 0, 0, bp_data, &basis_x));

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

218 DM dm_coord;
219 Vec coords;
220 const PetscScalar *coord_array;
221 CeedBasis basis_x, basis_u;
222 CeedElemRestriction elem_restr_x, elem_restr_u, elem_restr_u_i, elem_restr_qd_i;
223 CeedQFunction qf_setup_geo, qf_setup_rhs, qf_error;
224 CeedOperator op_setup_geo, op_setup_rhs;
225 CeedVector x_coord, q_data, target, rhs;
225 CeedInt num_qpts, c_start, c_end, num_elem, q_data_size = bp_data.q_data_size;
226 PetscInt c_start, c_end, num_elem;
227 CeedInt num_qpts, q_data_size = bp_data.q_data_size;
226 CeedScalar R = 1; // radius of the sphere
227 CeedScalar l = 1.0 / PetscSqrtReal(3.0); // half edge of the inscribed cube
228
229 PetscFunctionBeginUser;
230 PetscCall(DMGetCoordinateDM(dm, &dm_coord));
231
232 // CEED bases
233 PetscCall(CreateBasisFromPlex(ceed, dm_coord, 0, 0, 0, 0, bp_data, &basis_x));

--- 103 unchanged lines hidden --- 		228 CeedScalar R = 1; // radius of the sphere
229 CeedScalar l = 1.0 / PetscSqrtReal(3.0); // half edge of the inscribed cube
230
231 PetscFunctionBeginUser;
232 PetscCall(DMGetCoordinateDM(dm, &dm_coord));
233
234 // CEED bases
235 PetscCall(CreateBasisFromPlex(ceed, dm_coord, 0, 0, 0, 0, bp_data, &basis_x));

--- 103 unchanged lines hidden ---