setuplibceed.c (d310b3d31eeeddd20725517a3a61881a36d919f0) setuplibceed.c (d52d2bab094bb89d57050a07214fd9fdb874afca)
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

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

325 PetscCall(VecGetArrayRead(X_loc, &X_loc_array));
326 CeedVectorSetArray(ceed_data->x_coord, CEED_MEM_HOST, CEED_COPY_VALUES, (PetscScalar *)X_loc_array);
327 PetscCall(VecRestoreArrayRead(X_loc, &X_loc_array));
328
329 // -----------------------------------------------------------------------------
330 // CEED vectors
331 // -----------------------------------------------------------------------------
332 // -- Create CEED vector for geometric data
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

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

325 PetscCall(VecGetArrayRead(X_loc, &X_loc_array));
326 CeedVectorSetArray(ceed_data->x_coord, CEED_MEM_HOST, CEED_COPY_VALUES, (PetscScalar *)X_loc_array);
327 PetscCall(VecRestoreArrayRead(X_loc, &X_loc_array));
328
329 // -----------------------------------------------------------------------------
330 // CEED vectors
331 // -----------------------------------------------------------------------------
332 // -- Create CEED vector for geometric data
333 CeedInt num_qpts_vol;
334 PetscInt loc_num_elem_vol;
335 CeedBasisGetNumQuadraturePoints(ceed_data->basis_q, &num_qpts_vol);
336 CeedElemRestrictionGetNumElements(ceed_data->elem_restr_q, &loc_num_elem_vol);
337 CeedVectorCreate(ceed, q_data_size_vol * loc_num_elem_vol * num_qpts_vol, &ceed_data->q_data);
338
333 CeedElemRestrictionCreateVector(ceed_data->elem_restr_qd_i, &ceed_data->q_data, NULL);
339 CeedElemRestrictionCreateVector(elem_restr_jd_i, &jac_data, NULL);
334 CeedElemRestrictionCreateVector(elem_restr_jd_i, &jac_data, NULL);
335
340 // -----------------------------------------------------------------------------
341 // CEED Operators
342 // -----------------------------------------------------------------------------
343 // -- Create CEED operator for quadrature data
344 CeedOperatorCreate(ceed, ceed_data->qf_setup_vol, NULL, NULL, &ceed_data->op_setup_vol);
345 CeedOperatorSetField(ceed_data->op_setup_vol, "dx", ceed_data->elem_restr_x, ceed_data->basis_x, CEED_VECTOR_ACTIVE);
346 CeedOperatorSetField(ceed_data->op_setup_vol, "weight", CEED_ELEMRESTRICTION_NONE, ceed_data->basis_x, CEED_VECTOR_NONE);
347 CeedOperatorSetField(ceed_data->op_setup_vol, "qdata", ceed_data->elem_restr_qd_i, CEED_BASIS_COLLOCATED, CEED_VECTOR_ACTIVE);

--- 110 unchanged lines hidden --- 		336 // -----------------------------------------------------------------------------
337 // CEED Operators
338 // -----------------------------------------------------------------------------
339 // -- Create CEED operator for quadrature data
340 CeedOperatorCreate(ceed, ceed_data->qf_setup_vol, NULL, NULL, &ceed_data->op_setup_vol);
341 CeedOperatorSetField(ceed_data->op_setup_vol, "dx", ceed_data->elem_restr_x, ceed_data->basis_x, CEED_VECTOR_ACTIVE);
342 CeedOperatorSetField(ceed_data->op_setup_vol, "weight", CEED_ELEMRESTRICTION_NONE, ceed_data->basis_x, CEED_VECTOR_NONE);
343 CeedOperatorSetField(ceed_data->op_setup_vol, "qdata", ceed_data->elem_restr_qd_i, CEED_BASIS_COLLOCATED, CEED_VECTOR_ACTIVE);

--- 110 unchanged lines hidden ---