1*98285ab4SZach Atkins // Copyright (c) 2017-2022, Lawrence Livermore National Security, LLC and other CEED contributors. 2*98285ab4SZach Atkins // All Rights Reserved. See the top-level LICENSE and NOTICE files for details. 3*98285ab4SZach Atkins // 4*98285ab4SZach Atkins // SPDX-License-Identifier: BSD-2-Clause 5*98285ab4SZach Atkins // 6*98285ab4SZach Atkins // This file is part of CEED: http://github.com/ceed 7*98285ab4SZach Atkins 8*98285ab4SZach Atkins /// @file 9*98285ab4SZach Atkins /// Utility functions for using PETSc with libCEED 10*98285ab4SZach Atkins 11b7c563b6SJeremy L Thompson #ifndef libceed_petsc_examples_utils_h 12b7c563b6SJeremy L Thompson #define libceed_petsc_examples_utils_h 13e83e87a5Sjeremylt 14e83e87a5Sjeremylt #include <ceed.h> 15e83e87a5Sjeremylt #include <petsc.h> 162b730f8bSJeremy L Thompson 17de1229c5Srezgarshakeri #include "structs.h" 1849a40c8aSKenneth E. Jansen #if PETSC_VERSION_LT(3, 21, 0) 1949a40c8aSKenneth E. Jansen #define DMSetCoordinateDisc(a, b, c) DMProjectCoordinates(a, b) 2049a40c8aSKenneth E. Jansen #endif 21e83e87a5Sjeremylt 22e83e87a5Sjeremylt CeedMemType MemTypeP2C(PetscMemType mtype); 23b6972d74SZach Atkins PetscErrorCode VecP2C(Vec X_petsc, PetscMemType *mem_type, CeedVector x_ceed); 24b6972d74SZach Atkins PetscErrorCode VecC2P(CeedVector x_ceed, PetscMemType mem_type, Vec X_petsc); 25b6972d74SZach Atkins PetscErrorCode VecReadP2C(Vec X_petsc, PetscMemType *mem_type, CeedVector x_ceed); 26b6972d74SZach Atkins PetscErrorCode VecReadC2P(CeedVector x_ceed, PetscMemType mem_type, Vec X_petsc); 279b072555Sjeremylt PetscErrorCode Kershaw(DM dm_orig, PetscScalar eps); 282b730f8bSJeremy L Thompson PetscErrorCode SetupDMByDegree(DM dm, PetscInt p_degree, PetscInt q_extra, PetscInt num_comp_u, PetscInt topo_dim, bool enforce_bc); 292b730f8bSJeremy L Thompson PetscErrorCode CreateRestrictionFromPlex(Ceed ceed, DM dm, CeedInt height, DMLabel domain_label, CeedInt value, CeedElemRestriction *elem_restr); 30de1229c5Srezgarshakeri CeedElemTopology ElemTopologyP2C(DMPolytopeType cell_type); 312b730f8bSJeremy L Thompson PetscErrorCode DMFieldToDSField(DM dm, DMLabel domain_label, PetscInt dm_field, PetscInt *ds_field); 322b730f8bSJeremy L Thompson PetscErrorCode BasisCreateFromTabulation(Ceed ceed, DM dm, DMLabel domain_label, PetscInt label_value, PetscInt height, PetscInt face, PetscFE fe, 332b730f8bSJeremy L Thompson PetscTabulation basis_tabulation, PetscQuadrature quadrature, CeedBasis *basis); 342b730f8bSJeremy L Thompson PetscErrorCode CreateBasisFromPlex(Ceed ceed, DM dm, DMLabel domain_label, CeedInt label_value, CeedInt height, CeedInt dm_field, BPData bp_data, 35f755c37aSrezgarshakeri CeedBasis *basis); 36de1229c5Srezgarshakeri PetscErrorCode CreateDistributedDM(RunParams rp, DM *dm); 37de1229c5Srezgarshakeri 38b7c563b6SJeremy L Thompson #endif // libceed_petsc_examples_utils_h 39