petscdmlibmeshimpl.C File Reference

Go to the source code of this file.

Classes
struct	DM_libMesh
struct	DMVec_libMesh
struct	token

Functions
PetscErrorCode	DMlibMeshGetVec_Private (DM, const char *, Vec *)
PetscErrorCode	DMlibMeshSetUpName_Private (DM dm)
PetscErrorCode	DMlibMeshSetSystem_libMesh (DM dm, NonlinearImplicitSystem &sys)
PetscErrorCode	DMlibMeshGetSystem_libMesh (DM dm, NonlinearImplicitSystem *&sys)
PetscErrorCode	DMlibMeshGetBlocks (DM dm, PetscInt *n, char ***blocknames)
PetscErrorCode	DMlibMeshGetVariables (DM dm, PetscInt *n, char ***varnames)
static PetscErrorCode	DMCreateFieldDecomposition_libMesh (DM dm, PetscInt *len, char ***namelist, IS **islist, DM **dmlist)
static PetscErrorCode	DMCreateDomainDecomposition_libMesh (DM dm, PetscInt *len, char ***namelist, IS **innerislist, IS **outerislist, DM **dmlist)
PetscErrorCode	DMlibMeshCreateFieldDecompositionDM (DM dm, PetscInt dnumber, PetscInt *dsizes, char ***dvarlists, DM *ddm)
PetscErrorCode	DMlibMeshCreateDomainDecompositionDM (DM dm, PetscInt dnumber, PetscInt *dsizes, char ***dblocklists, DM *ddm)
static PetscErrorCode	DMlibMeshParseDecompositionDescriptor_Private (DM dm, const char *ddesc, PetscInt *dtype, PetscInt *dcount, PetscInt **dsizes, char ****dlists)
static PetscErrorCode	DMCreateFieldDecompositionDM_libMesh (DM dm, const char *ddesc, DM *ddm)
static PetscErrorCode	DMCreateDomainDecompositionDM_libMesh (DM dm, const char *ddesc, DM *ddm)
static PetscErrorCode	DMlibMeshFunction (DM dm, Vec x, Vec r)
static PetscErrorCode	SNESFunction_DMlibMesh (SNES, Vec x, Vec r, void *ctx)
static PetscErrorCode	DMlibMeshJacobian (#if PETSC_RELEASE_LESS_THAN(3, 5, 0) DM dm, Vec x, Mat jac, Mat pc, MatStructure *msflag #else DM dm, Vec x, Mat jac, Mat pc #endif)
static PetscErrorCode	SNESJacobian_DMlibMesh (#if PETSC_RELEASE_LESS_THAN(3, 5, 0) SNES, Vec x, Mat *jac, Mat *pc, MatStructure *flag, void *ctx #else SNES, Vec x, Mat jac, Mat pc, void *ctx #endif)
static PetscErrorCode	DMVariableBounds_libMesh (DM dm, Vec xl, Vec xu)
static PetscErrorCode	DMCreateGlobalVector_libMesh (DM dm, Vec *x)
static PetscErrorCode	DMCreateMatrix_libMesh (DM dm, const MatType, Mat *A) static PetscErrorCode DMCreateMatrix_libMesh(DM dm
	CHKERRQ (ierr)
	if (!eq) SETERRQ2(((PetscObject) dm) -> comm, PETSC_ERR_ARG_WRONG, "DM of type %s, not of type %s",((PetscObject) dm) ->type, DMLIBMESH)
	if (!dlm->sys) SETERRQ(PETSC_COMM_WORLD
	PetscFunctionReturn (0)
static PetscErrorCode	DMView_libMesh (DM dm, PetscViewer viewer)
static PetscErrorCode	DMSetUp_libMesh (DM dm)
static PetscErrorCode	DMDestroy_libMesh (DM dm)
PetscErrorCode	DMCreate_libMesh (DM dm)

Variables
static PetscErrorCode Mat *	A
PetscErrorCode	ierr = PetscObjectTypeCompare((PetscObject)dm, DMLIBMESH, &eq)
DM_libMesh *	dlm = (DM_libMesh *)(dm->data)
PetscBool	eq
	PETSC_ERR_ARG_WRONGSTATE
No libMesh system set for	DM_libMesh

Function Documentation

CHKERRQ()

CHKERRQ

(

ierr

)

Referenced by DMlibMeshGetSystem(), and DMlibMeshSetSystem().

DMCreate_libMesh()

PetscErrorCode DMCreate_libMesh

(

DM

dm

)

Definition at line 1170 of file petscdmlibmeshimpl.C.

{
  PetscErrorCode ierr;
  DM_libMesh     * dlm;
 
  PetscFunctionBegin;
  PetscValidHeaderSpecific(dm,DM_CLASSID,1);
#if PETSC_RELEASE_LESS_THAN(3,5,0)
  ierr = PetscNewLog(dm,DM_libMesh,&dlm);CHKERRQ(ierr);
#else
  ierr = PetscNewLog(dm,&dlm);CHKERRQ(ierr);
#endif
  dm->data = dlm;
 
  /* DMlibMesh impl */
  dlm->varids     = new(std::map<std::string, unsigned int>);
  dlm->blockids   = new(std::map<std::string, unsigned int>);
  dlm->varnames   = new(std::map<unsigned int, std::string>);
  dlm->blocknames = new(std::map<unsigned int, std::string>);
  dlm->decomposition   = PETSC_NULL;
  dlm->decomposition_type  = DMLIBMESH_NO_DECOMPOSITION;
 
  /* DM API */
  dm->ops->createglobalvector = DMCreateGlobalVector_libMesh;
  dm->ops->createlocalvector  = 0; // DMCreateLocalVector_libMesh;
  dm->ops->getcoloring        = 0; // DMGetColoring_libMesh;
  dm->ops->creatematrix       = DMCreateMatrix_libMesh;
  dm->ops->createinterpolation= 0; // DMCreateInterpolation_libMesh;
 
  dm->ops->refine             = 0; // DMRefine_libMesh;
  dm->ops->coarsen            = 0; // DMCoarsen_libMesh;
 
  // * dm->ops->getinjection was renamed to dm->ops->createinjection in PETSc 5a84ad338 (5 Jul 2019)
  // * dm->ops-getaggregates was removed in PETSc 97779f9a (5 Jul 2019)
  // * Both changes were merged into PETSc master in 94aad3ce (7 Jul 2019).
#if PETSC_RELEASE_LESS_THAN(3,12,0)
  dm->ops->getinjection       = 0; // DMGetInjection_libMesh;
  dm->ops->getaggregates      = 0; // DMGetAggregates_libMesh;
#else
  dm->ops->createinjection = 0;
#endif
 
 
#if PETSC_RELEASE_LESS_THAN(3,3,1)
  dm->ops->createfielddecompositiondm  = DMCreateFieldDecompositionDM_libMesh;
  dm->ops->createdomaindecompositiondm = DMCreateDomainDecompositionDM_libMesh;
#endif
  dm->ops->createfielddecomposition    = DMCreateFieldDecomposition_libMesh;
  dm->ops->createdomaindecomposition   = DMCreateDomainDecomposition_libMesh;
 
  dm->ops->destroy            = DMDestroy_libMesh;
  dm->ops->view               = DMView_libMesh;
  dm->ops->setfromoptions     = 0; // DMSetFromOptions_libMesh;
  dm->ops->setup              = DMSetUp_libMesh;
 
  /* DMlibMesh API */
#if PETSC_RELEASE_LESS_THAN(3,4,0)
  ierr = PetscObjectComposeFunction((PetscObject)dm,"DMlibMeshSetSystem_C",PETSC_NULL,(PetscVoidFunction)DMlibMeshSetSystem_libMesh);CHKERRQ(ierr);
  ierr = PetscObjectComposeFunction((PetscObject)dm,"DMlibMeshGetSystem_C",PETSC_NULL,(PetscVoidFunction)DMlibMeshGetSystem_libMesh);CHKERRQ(ierr);
#else
  ierr = PetscObjectComposeFunction((PetscObject)dm,"DMlibMeshSetSystem_C",DMlibMeshSetSystem_libMesh);CHKERRQ(ierr);
  ierr = PetscObjectComposeFunction((PetscObject)dm,"DMlibMeshGetSystem_C",DMlibMeshGetSystem_libMesh);CHKERRQ(ierr);
#endif
 
  PetscFunctionReturn(0);
}

References DM_libMesh::blockids, DM_libMesh::blocknames, libMesh::CHKERRQ(), DM_libMesh::decomposition, DM_libMesh::decomposition_type, dlm, DMCreateDomainDecomposition_libMesh(), DMCreateDomainDecompositionDM_libMesh(), DMCreateFieldDecomposition_libMesh(), DMCreateFieldDecompositionDM_libMesh(), DMCreateGlobalVector_libMesh(), DMCreateMatrix_libMesh(), DMDestroy_libMesh(), DMlibMeshGetSystem_libMesh(), DMlibMeshSetSystem_libMesh(), DMSetUp_libMesh(), DMView_libMesh(), libMesh::ierr, PetscFunctionReturn(), DM_libMesh::varids, and DM_libMesh::varnames.

Referenced by libMesh::LibMeshInit::LibMeshInit().

DMCreateDomainDecomposition_libMesh()

static PetscErrorCode DMCreateDomainDecomposition_libMesh ( DM  dm, PetscInt *  len, char ***  namelist, IS **  innerislist, IS **  outerislist, DM **  dmlist  )

static

Definition at line 378 of file petscdmlibmeshimpl.C.

{
  PetscFunctionBegin;
  PetscErrorCode ierr;
  DM_libMesh     * dlm = (DM_libMesh *)(dm->data);
  NonlinearImplicitSystem * sys = dlm->sys;
  IS emb;
  if (dlm->decomposition_type != DMLIBMESH_DOMAIN_DECOMPOSITION) PetscFunctionReturn(0);
  *len = cast_int<unsigned int>(dlm->decomposition->size());
  if (namelist)      {ierr = PetscMalloc(*len*sizeof(char *), namelist);  CHKERRQ(ierr);}
  if (innerislist)   {ierr = PetscMalloc(*len*sizeof(IS),    innerislist);    CHKERRQ(ierr);}
  if (outerislist)   *outerislist = PETSC_NULL; /* FIX: allow mesh-based overlap. */
  if (dmlist)        {ierr = PetscMalloc(*len*sizeof(DM),    dmlist);    CHKERRQ(ierr);}
  for (auto d : index_range(*dlm->decomposition)) {
    std::set<numeric_index_type>               dindices;
    std::string                          dname;
    std::map<std::string, unsigned int>  dblockids;
    std::map<unsigned int,std::string>   dblocknames;
    unsigned int                         dbcount = 0;
    for (const auto & b : (*dlm->decomposition)[d]) {
      std::string bname = (*dlm->blocknames)[b];
      dblockids.insert(std::pair<std::string, unsigned int>(bname,b));
      dblocknames.insert(std::pair<unsigned int,std::string>(b,bname));
      if (!dbcount) dname = bname;
      else   dname += "_" + bname;
      ++dbcount;
      if (!innerislist) continue;
      const subdomain_id_type b_sbd_id = cast_int<subdomain_id_type>(b);
      MeshBase::const_element_iterator       el     = sys->get_mesh().active_local_subdomain_elements_begin(b_sbd_id);
      const MeshBase::const_element_iterator end_el = sys->get_mesh().active_local_subdomain_elements_end(b_sbd_id);
      for ( ; el != end_el; ++el) {
        const Elem * elem = *el;
        std::vector<numeric_index_type> evindices;
        // Iterate only over this DM's variables.
        for (const auto & pr : *(dlm->varids)) {
          // Get the degree of freedom indices for the given variable off the current element.
          sys->get_dof_map().dof_indices(elem, evindices, pr.second);
          for (const auto & dof : evindices) {
            if (dof >= sys->get_dof_map().first_dof() && dof < sys->get_dof_map().end_dof()) // might want to use variable_first/last_local_dof instead
              dindices.insert(dof);
          }
        }
      }
    }
    if (namelist) {
      ierr = PetscStrallocpy(dname.c_str(),(*namelist)+d);            CHKERRQ(ierr);
    }
    if (innerislist) {
      PetscInt * darray;
      IS dis;
      ierr = PetscMalloc(sizeof(PetscInt)*dindices.size(), &darray); CHKERRQ(ierr);
      numeric_index_type i = 0;
      for (const auto & id : dindices) {
        darray[i] = id;
        ++i;
      }
      ierr = ISCreateGeneral(((PetscObject)dm)->comm,
                             cast_int<PetscInt>(dindices.size()),
                             darray, PETSC_OWN_POINTER, &dis);
      CHKERRQ(ierr);
      if (dlm->embedding) {
        /* Create a relative embedding into the parent's index space. */
#if PETSC_RELEASE_LESS_THAN(3,3,1)
        ierr = ISMapFactorRight(dis,dlm->embedding, PETSC_TRUE, &emb); CHKERRQ(ierr);
#else
        ierr = ISEmbed(dis,dlm->embedding, PETSC_TRUE, &emb); CHKERRQ(ierr);
#endif
        PetscInt elen, dlen;
        ierr = ISGetLocalSize(emb, &elen); CHKERRQ(ierr);
        ierr = ISGetLocalSize(dis, &dlen);  CHKERRQ(ierr);
        if (elen != dlen) SETERRQ1(((PetscObject)dm)->comm, PETSC_ERR_PLIB, "Failed to embed field %D", d);
        ierr = ISDestroy(&dis); CHKERRQ(ierr);
        dis = emb;
      }
      else {
        emb = dis;
      }
      if (innerislist) {
        ierr = PetscObjectReference((PetscObject)dis); CHKERRQ(ierr);
        (*innerislist)[d] = dis;
      }
      ierr = ISDestroy(&dis); CHKERRQ(ierr);
    }
    if (dmlist) {
      DM ddm;
      ierr = DMCreate(((PetscObject)dm)->comm, &ddm); CHKERRQ(ierr);
      ierr = DMSetType(ddm, DMLIBMESH);               CHKERRQ(ierr);
      DM_libMesh * ddlm = (DM_libMesh *)(ddm->data);
      ddlm->sys = dlm->sys;
      /* copy over the varids and varnames */
      *ddlm->varids    = *dlm->varids;
      *ddlm->varnames  = *dlm->varnames;
      /* set the blocks from the d-th part of the decomposition. */
      *ddlm->blockids    = dblockids;
      *ddlm->blocknames  = dblocknames;
      ierr = PetscObjectReference((PetscObject)emb); CHKERRQ(ierr);
      ddlm->embedding = emb;
      ddlm->embedding_type = DMLIBMESH_DOMAIN_EMBEDDING;
 
      ierr = DMlibMeshSetUpName_Private(ddm); CHKERRQ(ierr);
      ierr = DMSetFromOptions(ddm);           CHKERRQ(ierr);
      (*dmlist)[d] = ddm;
    }
  }
  PetscFunctionReturn(0);
}

References libMesh::MeshBase::active_local_subdomain_elements_begin(), libMesh::MeshBase::active_local_subdomain_elements_end(), DM_libMesh::blockids, DM_libMesh::blocknames, libMesh::CHKERRQ(), DM_libMesh::decomposition, DM_libMesh::decomposition_type, dlm, DMlibMeshSetUpName_Private(), libMesh::DofMap::dof_indices(), DM_libMesh::embedding, DM_libMesh::embedding_type, libMesh::DofMap::end_dof(), libMesh::DofMap::first_dof(), libMesh::System::get_dof_map(), libMesh::System::get_mesh(), libMesh::ierr, libMesh::index_range(), PETSC_OWN_POINTER, PetscFunctionReturn(), DM_libMesh::sys, DM_libMesh::varids, and DM_libMesh::varnames.

Referenced by DMCreate_libMesh().

DMCreateDomainDecompositionDM_libMesh()

static PetscErrorCode DMCreateDomainDecompositionDM_libMesh ( DM  dm, const char *  ddesc, DM *  ddm  )

static

Definition at line 764 of file petscdmlibmeshimpl.C.

{
  PetscFunctionBegin;
  PetscErrorCode ierr;
  PetscInt dtype, dcount;
  PetscInt * dsizes;
  char *** dlists;
  PetscFunctionBegin;
  *ddm = PETSC_NULL;
  ierr = DMlibMeshParseDecompositionDescriptor_Private(dm,ddesc,&dtype,&dcount,&dsizes,&dlists); CHKERRQ(ierr);
  if (dtype == DMLIBMESH_DOMAIN_DECOMPOSITION) {
    ierr = DMlibMeshCreateDomainDecompositionDM(dm,dcount,dsizes,dlists,ddm); CHKERRQ(ierr);
  }
  else SETERRQ1(((PetscObject)dm)->comm, PETSC_ERR_PLIB, "Unexpected unknown decomposition type for domain decomposition descriptor %s", ddesc);
  PetscFunctionReturn(0);
}

References libMesh::CHKERRQ(), DMlibMeshCreateDomainDecompositionDM(), DMlibMeshParseDecompositionDescriptor_Private(), libMesh::ierr, and PetscFunctionReturn().

Referenced by DMCreate_libMesh().

DMCreateFieldDecomposition_libMesh()

static PetscErrorCode DMCreateFieldDecomposition_libMesh ( DM  dm, PetscInt *  len, char ***  namelist, IS **  islist, DM **  dmlist  )

static

Definition at line 272 of file petscdmlibmeshimpl.C.

{
  PetscFunctionBegin;
  PetscErrorCode ierr;
  DM_libMesh     * dlm = (DM_libMesh *)(dm->data);
  NonlinearImplicitSystem * sys = dlm->sys;
  IS emb;
  if (dlm->decomposition_type != DMLIBMESH_FIELD_DECOMPOSITION) PetscFunctionReturn(0);
 
  *len = cast_int<unsigned int>(dlm->decomposition->size());
  if (namelist) {ierr = PetscMalloc(*len*sizeof(char *), namelist);  CHKERRQ(ierr);}
  if (islist)   {ierr = PetscMalloc(*len*sizeof(IS),    islist);    CHKERRQ(ierr);}
  if (dmlist)   {ierr = PetscMalloc(*len*sizeof(DM),    dmlist);    CHKERRQ(ierr);}
  DofMap & dofmap = dlm->sys->get_dof_map();
  for (auto d : index_range(*dlm->decomposition)) {
    std::set<numeric_index_type>         dindices;
    std::string                          dname;
    std::map<std::string, unsigned int>  dvarids;
    std::map<unsigned int, std::string>  dvarnames;
    unsigned int                         dvcount = 0;
    for (const auto & v : (*dlm->decomposition)[d]) {
      std::string vname = (*dlm->varnames)[v];
      dvarids.insert(std::pair<std::string, unsigned int>(vname,v));
      dvarnames.insert(std::pair<unsigned int,std::string>(v,vname));
      if (!dvcount) dname = vname;
      else   dname += "_" + vname;
      ++dvcount;
      if (!islist) continue;
      // Iterate only over this DM's blocks.
      for (const auto & pr : *(dlm->blockids)) {
        const subdomain_id_type sbd_id = cast_int<subdomain_id_type>(pr.second);
        for (const auto & elem :
               as_range(sys->get_mesh().active_local_subdomain_elements_begin(sbd_id),
                        sys->get_mesh().active_local_subdomain_elements_end(sbd_id))) {
          //unsigned int e_subdomain = elem->subdomain_id();
          std::vector<numeric_index_type> evindices;
          // Get the degree of freedom indices for the given variable off the current element.
          dofmap.dof_indices(elem, evindices, v);
          for (numeric_index_type dof : evindices) {
            if (dof >= dofmap.first_dof() && dof < dofmap.end_dof()) // might want to use variable_first/last_local_dof instead
              dindices.insert(dof);
          }
        }
      }
    }
    if (namelist) {
      ierr = PetscStrallocpy(dname.c_str(),(*namelist)+d);            CHKERRQ(ierr);
    }
    if (islist) {
      IS dis;
      PetscInt * darray;
      ierr = PetscMalloc(sizeof(PetscInt)*dindices.size(), &darray); CHKERRQ(ierr);
      numeric_index_type i = 0;
      for (const auto & id : dindices) {
        darray[i] = id;
        ++i;
      }
      ierr = ISCreateGeneral(((PetscObject)dm)->comm,
                             cast_int<PetscInt>(dindices.size()),
                             darray, PETSC_OWN_POINTER, &dis);
      CHKERRQ(ierr);
      if (dlm->embedding) {
        /* Create a relative embedding into the parent's index space. */
#if PETSC_RELEASE_LESS_THAN(3,3,1)
        ierr = ISMapFactorRight(dis,dlm->embedding, PETSC_TRUE, &emb); CHKERRQ(ierr);
#else
        ierr = ISEmbed(dis,dlm->embedding, PETSC_TRUE, &emb); CHKERRQ(ierr);
#endif
        PetscInt elen, dlen;
        ierr = ISGetLocalSize(emb, &elen); CHKERRQ(ierr);
        ierr = ISGetLocalSize(dis, &dlen); CHKERRQ(ierr);
        if (elen != dlen) SETERRQ1(((PetscObject)dm)->comm, PETSC_ERR_PLIB, "Failed to embed subdomain %D", d);
        ierr = ISDestroy(&dis); CHKERRQ(ierr);
        dis = emb;
      }
      else {
        emb = dis;
      }
      (*islist)[d] = dis;
    }
    if (dmlist) {
      DM ddm;
      ierr = DMCreate(((PetscObject)dm)->comm, &ddm); CHKERRQ(ierr);
      ierr = DMSetType(ddm, DMLIBMESH);               CHKERRQ(ierr);
      DM_libMesh * ddlm = (DM_libMesh *)(ddm->data);
      ddlm->sys = dlm->sys;
      /* copy over the block ids and names */
      *ddlm->blockids = *dlm->blockids;
      *ddlm->blocknames = *dlm->blocknames;
      /* set the vars from the d-th part of the decomposition. */
      *ddlm->varids     = dvarids;
      *ddlm->varnames   = dvarnames;
      ierr = PetscObjectReference((PetscObject)emb); CHKERRQ(ierr);
      ddlm->embedding = emb;
      ddlm->embedding_type = DMLIBMESH_FIELD_EMBEDDING;
 
      ierr = DMlibMeshSetUpName_Private(ddm); CHKERRQ(ierr);
      ierr = DMSetFromOptions(ddm);           CHKERRQ(ierr);
      (*dmlist)[d] = ddm;
    }
  }
  PetscFunctionReturn(0);
}

References libMesh::MeshBase::active_local_subdomain_elements_begin(), libMesh::MeshBase::active_local_subdomain_elements_end(), libMesh::as_range(), DM_libMesh::blockids, DM_libMesh::blocknames, libMesh::CHKERRQ(), DM_libMesh::decomposition, DM_libMesh::decomposition_type, dlm, DMlibMeshSetUpName_Private(), libMesh::DofMap::dof_indices(), DM_libMesh::embedding, DM_libMesh::embedding_type, libMesh::DofMap::end_dof(), libMesh::DofMap::first_dof(), libMesh::System::get_dof_map(), libMesh::System::get_mesh(), libMesh::ierr, libMesh::index_range(), PETSC_OWN_POINTER, PetscFunctionReturn(), DM_libMesh::sys, DM_libMesh::varids, and DM_libMesh::varnames.

Referenced by DMCreate_libMesh().

DMCreateFieldDecompositionDM_libMesh()

static PetscErrorCode DMCreateFieldDecompositionDM_libMesh ( DM  dm, const char *  ddesc, DM *  ddm  )

static

Definition at line 745 of file petscdmlibmeshimpl.C.

{
  PetscFunctionBegin;
  PetscErrorCode ierr;
  PetscInt dtype, dcount;
  PetscInt * dsizes;
  char *** dlists;
  PetscFunctionBegin;
  *ddm = PETSC_NULL;
  ierr = DMlibMeshParseDecompositionDescriptor_Private(dm,ddesc,&dtype,&dcount,&dsizes,&dlists); CHKERRQ(ierr);
  if (dtype == DMLIBMESH_FIELD_DECOMPOSITION){
    ierr = DMlibMeshCreateFieldDecompositionDM(dm,dcount,dsizes,dlists,ddm); CHKERRQ(ierr);
  }
  else SETERRQ1(((PetscObject)dm)->comm, PETSC_ERR_PLIB, "Unexpected unknown decomposition type for field decomposition descriptor %s", ddesc);
  PetscFunctionReturn(0);
}

References libMesh::CHKERRQ(), DMlibMeshCreateFieldDecompositionDM(), DMlibMeshParseDecompositionDescriptor_Private(), libMesh::ierr, and PetscFunctionReturn().

Referenced by DMCreate_libMesh().

DMCreateGlobalVector_libMesh()

static PetscErrorCode DMCreateGlobalVector_libMesh ( DM  dm, Vec *  x  )

static

Definition at line 980 of file petscdmlibmeshimpl.C.

{
  PetscFunctionBegin;
  PetscErrorCode ierr;
  DM_libMesh     * dlm = (DM_libMesh *)(dm->data);
  PetscBool eq;
 
  ierr = PetscObjectTypeCompare((PetscObject)dm, DMLIBMESH, &eq); CHKERRQ(ierr);
 
  if (!eq)
    SETERRQ2(((PetscObject)dm)->comm, PETSC_ERR_ARG_WRONG, "DM of type %s, not of type %s", ((PetscObject)dm)->type, DMLIBMESH);
 
  if (!dlm->sys)
    SETERRQ(PETSC_COMM_WORLD, PETSC_ERR_ARG_WRONGSTATE, "No libMesh system set for DM_libMesh");
 
  NumericVector<Number> * nv = (dlm->sys->solution).get();
  PetscVector<Number> *   pv = dynamic_cast<PetscVector<Number> *>(nv);
  Vec                    v  = pv->vec();
  /* Unfortunately, currently this does not produce a ghosted vector, so nonlinear subproblem solves aren't going to be easily available.
     Should work fine for getting vectors out for linear subproblem solvers. */
  if (dlm->embedding) {
    PetscInt n;
    ierr = VecCreate(((PetscObject)v)->comm, x);       CHKERRQ(ierr);
    ierr = ISGetLocalSize(dlm->embedding, &n);         CHKERRQ(ierr);
    ierr = VecSetSizes(*x,n,PETSC_DETERMINE);           CHKERRQ(ierr);
    ierr = VecSetType(*x,((PetscObject)v)->type_name); CHKERRQ(ierr);
    ierr = VecSetFromOptions(*x);                      CHKERRQ(ierr);
    ierr = VecSetUp(*x);                               CHKERRQ(ierr);
  }
  else {
    ierr = VecDuplicate(v,x); CHKERRQ(ierr);
  }
  ierr = PetscObjectCompose((PetscObject)*x,"DM",(PetscObject)dm); CHKERRQ(ierr);
  PetscFunctionReturn(0);
}

References libMesh::CHKERRQ(), dlm, DM_libMesh::embedding, eq, libMesh::ReferenceElem::get(), libMesh::ierr, PETSC_ERR_ARG_WRONGSTATE, PetscFunctionReturn(), libMesh::System::solution, DM_libMesh::sys, and libMesh::PetscVector< T >::vec().

Referenced by DMCreate_libMesh().

DMCreateMatrix_libMesh()

static PetscErrorCode DMCreateMatrix_libMesh ( DM  dm, const  MatType, Mat *  A  )

static

Referenced by DMCreate_libMesh().

DMDestroy_libMesh()

static PetscErrorCode DMDestroy_libMesh ( DM  dm )

static

Definition at line 1151 of file petscdmlibmeshimpl.C.

{
  DM_libMesh * dlm = (DM_libMesh *)(dm->data);
  PetscErrorCode ierr;
  PetscFunctionBegin;
  delete dlm->varids;
  delete dlm->varnames;
  delete dlm->blockids;
  delete dlm->blocknames;
  delete dlm->decomposition;
  ierr = ISDestroy(&dlm->embedding); CHKERRQ(ierr);
  ierr = PetscFree(dm->data); CHKERRQ(ierr);
 
  PetscFunctionReturn(0);
}

References DM_libMesh::blockids, DM_libMesh::blocknames, libMesh::CHKERRQ(), DM_libMesh::decomposition, dlm, DM_libMesh::embedding, libMesh::ierr, PetscFunctionReturn(), DM_libMesh::varids, and DM_libMesh::varnames.

Referenced by DMCreate_libMesh().

DMlibMeshCreateDomainDecompositionDM()

PetscErrorCode DMlibMeshCreateDomainDecompositionDM	(	DM	dm,
		PetscInt	dnumber,
		PetscInt *	dsizes,
		char ***	dblocklists,
		DM *	ddm
	)

Definition at line 540 of file petscdmlibmeshimpl.C.

{
  PetscErrorCode ierr;
  PetscBool islibmesh;
  PetscFunctionBegin;
  PetscValidHeaderSpecific(dm,DM_CLASSID,1);
  ierr = PetscObjectTypeCompare((PetscObject)dm, DMLIBMESH,&islibmesh);
  if (!islibmesh) SETERRQ2(((PetscObject)dm)->comm, PETSC_ERR_ARG_WRONG, "Got DM of type %s, not of type %s", ((PetscObject)dm)->type_name, DMLIBMESH);
  if (dnumber < 0) SETERRQ1(((PetscObject)dm)->comm, PETSC_ERR_ARG_WRONG, "Negative number %D of decomposition parts", dnumber);
  PetscValidPointer(ddm,5);
  DM_libMesh * dlm = (DM_libMesh *)(dm->data);
  ierr = DMCreate(((PetscObject)dm)->comm, ddm); CHKERRQ(ierr);
  ierr = DMSetType(*ddm, DMLIBMESH);             CHKERRQ(ierr);
  DM_libMesh * ddlm = (DM_libMesh *)((*ddm)->data);
  ddlm->sys = dlm->sys;
  ddlm->varids   = dlm->varids;
  ddlm->varnames = dlm->varnames;
  ddlm->blockids   = dlm->blockids;
  ddlm->blocknames = dlm->blocknames;
  ddlm->decomposition = new(std::vector<std::set<unsigned int>>);
  ddlm->decomposition_type = DMLIBMESH_DOMAIN_DECOMPOSITION;
  if (dnumber) {
    for (PetscInt d = 0; d < dnumber; ++d) {
      if (dsizes[d] < 0) SETERRQ2(((PetscObject)dm)->comm, PETSC_ERR_ARG_WRONG, "Negative size %D of decomposition part %D", dsizes[d],d);
      ddlm->decomposition->push_back(std::set<unsigned int>());
      for (PetscInt b = 0; b < dsizes[d]; ++b) {
        std::string bname(dblocklists[d][b]);
        std::map<std::string, unsigned int>::const_iterator bit = dlm->blockids->find(bname);
        if (bit == dlm->blockids->end())
          SETERRQ3(((PetscObject)dm)->comm, PETSC_ERR_ARG_WRONG, "Block %D on the %D-th list with name %s is not owned by this DM", b, d, dblocklists[d][b]);
        unsigned int bid = bit->second;
        (*ddlm->decomposition)[d].insert(bid);
      }
    }
  }
  else { // Empty splits indicate default: split all blocks with one per split.
    PetscInt d = 0;
    for (const auto & pr : (*ddlm->blockids)) {
      ddlm->decomposition->push_back(std::set<unsigned int>());
      unsigned int bid = pr.second;
      (*ddlm->decomposition)[d].insert(bid);
      ++d;
    }
  }
  ierr = DMlibMeshSetUpName_Private(*ddm); CHKERRQ(ierr);
  ierr = DMSetFromOptions(*ddm);           CHKERRQ(ierr);
  ierr = DMSetUp(*ddm);                    CHKERRQ(ierr);
  PetscFunctionReturn(0);
}

References DM_libMesh::blockids, DM_libMesh::blocknames, libMesh::CHKERRQ(), DM_libMesh::decomposition, DM_libMesh::decomposition_type, dlm, DMlibMeshSetUpName_Private(), libMesh::ierr, PetscFunctionReturn(), DM_libMesh::sys, DM_libMesh::varids, and DM_libMesh::varnames.

Referenced by DMCreateDomainDecompositionDM_libMesh().

DMlibMeshCreateFieldDecompositionDM()

PetscErrorCode DMlibMeshCreateFieldDecompositionDM	(	DM	dm,
		PetscInt	dnumber,
		PetscInt *	dsizes,
		char ***	dvarlists,
		DM *	ddm
	)

Definition at line 488 of file petscdmlibmeshimpl.C.

{
  PetscErrorCode ierr;
  PetscBool islibmesh;
  PetscFunctionBegin;
  PetscValidHeaderSpecific(dm,DM_CLASSID,1);
  ierr = PetscObjectTypeCompare((PetscObject)dm, DMLIBMESH,&islibmesh);
  if (!islibmesh) SETERRQ2(((PetscObject)dm)->comm, PETSC_ERR_ARG_WRONG, "Got DM of type %s, not of type %s", ((PetscObject)dm)->type_name, DMLIBMESH);
  if (dnumber < 0) SETERRQ1(((PetscObject)dm)->comm, PETSC_ERR_ARG_WRONG, "Negative number %D of decomposition parts", dnumber);
  PetscValidPointer(ddm,5);
  DM_libMesh * dlm = (DM_libMesh *)(dm->data);
  ierr = DMCreate(((PetscObject)dm)->comm, ddm); CHKERRQ(ierr);
  ierr = DMSetType(*ddm, DMLIBMESH);             CHKERRQ(ierr);
  DM_libMesh * ddlm = (DM_libMesh *)((*ddm)->data);
  ddlm->sys = dlm->sys;
  ddlm->varids = dlm->varids;
  ddlm->varnames = dlm->varnames;
  ddlm->blockids = dlm->blockids;
  ddlm->blocknames = dlm->blocknames;
  ddlm->decomposition = new(std::vector<std::set<unsigned int>>);
  ddlm->decomposition_type = DMLIBMESH_FIELD_DECOMPOSITION;
  if (dnumber) {
    for (PetscInt d = 0; d < dnumber; ++d) {
      if (dsizes[d] < 0) SETERRQ2(((PetscObject)dm)->comm, PETSC_ERR_ARG_WRONG, "Negative size %D of decomposition part %D", dsizes[d],d);
      ddlm->decomposition->push_back(std::set<unsigned int>());
      for (PetscInt v = 0; v < dsizes[d]; ++v) {
        std::string vname(dvarlists[d][v]);
        std::map<std::string, unsigned int>::const_iterator vit = dlm->varids->find(vname);
        if (vit == dlm->varids->end())
          SETERRQ3(((PetscObject)dm)->comm, PETSC_ERR_ARG_WRONG, "Variable %D on the %D-th list with name %s is not owned by this DM", v, d, dvarlists[d][v]);
        unsigned int vid = vit->second;
        (*ddlm->decomposition)[d].insert(vid);
      }
    }
  }
  else { // Empty splits indicate default: split all variables with one per split.
    PetscInt d = 0;
    for (const auto & pr : (*ddlm->varids)) {
      ddlm->decomposition->push_back(std::set<unsigned int>());
      unsigned int vid = pr.second;
      (*ddlm->decomposition)[d].insert(vid);
      ++d;
    }
  }
  ierr = DMlibMeshSetUpName_Private(*ddm); CHKERRQ(ierr);
  ierr = DMSetFromOptions(*ddm);           CHKERRQ(ierr);
  ierr = DMSetUp(*ddm);                    CHKERRQ(ierr);
  PetscFunctionReturn(0);
}

References DM_libMesh::blockids, DM_libMesh::blocknames, libMesh::CHKERRQ(), DM_libMesh::decomposition, DM_libMesh::decomposition_type, dlm, DMlibMeshSetUpName_Private(), libMesh::ierr, PetscFunctionReturn(), DM_libMesh::sys, DM_libMesh::varids, and DM_libMesh::varnames.

Referenced by DMCreateFieldDecompositionDM_libMesh().

DMlibMeshFunction()

static PetscErrorCode DMlibMeshFunction ( DM  dm, Vec  x, Vec  r  )

static

Definition at line 786 of file petscdmlibmeshimpl.C.

{
  PetscErrorCode ierr;
  PetscFunctionBegin;
  libmesh_assert(x);
  libmesh_assert(r);
 
  NonlinearImplicitSystem * _sys;
  ierr = DMlibMeshGetSystem(dm, _sys);CHKERRQ(ierr);
  NonlinearImplicitSystem & sys = *_sys;
  PetscVector<Number> & X_sys = *libmesh_cast_ptr<PetscVector<Number> *>(sys.solution.get());
  PetscVector<Number> & R_sys = *libmesh_cast_ptr<PetscVector<Number> *>(sys.rhs);
  PetscVector<Number> X_global(x, _sys->comm()), R(r, _sys->comm());
 
  // Use the systems update() to get a good local version of the parallel solution
  X_global.swap(X_sys);
  R.swap(R_sys);
 
  _sys->get_dof_map().enforce_constraints_exactly(*_sys);
  _sys->update();
 
  // Swap back
  X_global.swap(X_sys);
  R.swap(R_sys);
  R.zero();
 
  // if the user has provided both function pointers and objects only the pointer
  // will be used, so catch that as an error
  if (_sys->nonlinear_solver->residual && _sys->nonlinear_solver->residual_object)
    libmesh_error_msg("ERROR: cannot specify both a function and object to compute the Residual!");
 
  if (_sys->nonlinear_solver->matvec && _sys->nonlinear_solver->residual_and_jacobian_object)
    libmesh_error_msg("ERROR: cannot specify both a function and object to compute the combined Residual & Jacobian!");
 
  if (_sys->nonlinear_solver->residual != nullptr)
    _sys->nonlinear_solver->residual(*(_sys->current_local_solution.get()), R, *_sys);
 
  else if (_sys->nonlinear_solver->residual_object != nullptr)
    _sys->nonlinear_solver->residual_object->residual(*(_sys->current_local_solution.get()), R, *_sys);
 
  else if (_sys->nonlinear_solver->matvec   != nullptr)
    _sys->nonlinear_solver->matvec(*(_sys->current_local_solution.get()), &R, nullptr, *_sys);
 
  else if (_sys->nonlinear_solver->residual_and_jacobian_object != nullptr)
    _sys->nonlinear_solver->residual_and_jacobian_object->residual_and_jacobian(*(_sys->current_local_solution.get()), &R, nullptr, *_sys);
 
  else
    libmesh_error_msg("Error! Unable to compute residual and/or Jacobian!");
 
  R.close();
  X_global.close();
  PetscFunctionReturn(0);
}

References libMesh::CHKERRQ(), libMesh::ParallelObject::comm(), libMesh::System::current_local_solution, DMlibMeshGetSystem(), libMesh::DofMap::enforce_constraints_exactly(), libMesh::System::get_dof_map(), libMesh::ierr, libMesh::libmesh_assert(), libMesh::NonlinearImplicitSystem::nonlinear_solver, PetscFunctionReturn(), libMesh::ExplicitSystem::rhs, libMesh::System::solution, libMesh::PetscVector< T >::swap(), and libMesh::System::update().

Referenced by DMSetUp_libMesh(), and SNESFunction_DMlibMesh().

DMlibMeshGetBlocks()

PetscErrorCode DMlibMeshGetBlocks	(	DM	dm,
		PetscInt *	n,
		char ***	blocknames
	)

Definition at line 168 of file petscdmlibmeshimpl.C.

{
  PetscErrorCode ierr;
  PetscInt i;
  PetscFunctionBegin;
  PetscValidHeaderSpecific(dm,DM_CLASSID,1);
  PetscBool islibmesh;
  ierr = PetscObjectTypeCompare((PetscObject)dm, DMLIBMESH,&islibmesh);
  if (!islibmesh) SETERRQ2(((PetscObject)dm)->comm, PETSC_ERR_ARG_WRONG, "Got DM of type %s, not of type %s", ((PetscObject)dm)->type_name, DMLIBMESH);
  DM_libMesh * dlm = (DM_libMesh *)(dm->data);
  PetscValidPointer(n,2);
  *n = cast_int<unsigned int>(dlm->blockids->size());
  if (!blocknames) PetscFunctionReturn(0);
  ierr = PetscMalloc(*n*sizeof(char *), blocknames); CHKERRQ(ierr);
  i = 0;
  for (const auto & pr : *(dlm->blockids))
    {
      ierr = PetscStrallocpy(pr.first.c_str(), *blocknames+i); CHKERRQ(ierr);
      ++i;
    }
  PetscFunctionReturn(0);
}

References DM_libMesh::blockids, libMesh::CHKERRQ(), dlm, libMesh::ierr, and PetscFunctionReturn().

DMlibMeshGetSystem_libMesh()

PetscErrorCode DMlibMeshGetSystem_libMesh	(	DM	dm,
		NonlinearImplicitSystem *&	sys
	)

Definition at line 151 of file petscdmlibmeshimpl.C.

{
  PetscErrorCode ierr;
 
  PetscFunctionBegin;
  PetscValidHeaderSpecific(dm,DM_CLASSID,1);
  PetscBool islibmesh;
  ierr = PetscObjectTypeCompare((PetscObject)dm, DMLIBMESH,&islibmesh);CHKERRQ(ierr);
  if (!islibmesh) SETERRQ2(((PetscObject)dm)->comm, PETSC_ERR_ARG_WRONG, "Got DM of type %s, not of type %s", ((PetscObject)dm)->type_name, DMLIBMESH);
  DM_libMesh * dlm = (DM_libMesh *)(dm->data);
  sys = dlm->sys;
  PetscFunctionReturn(0);
}

References libMesh::CHKERRQ(), dlm, libMesh::ierr, PetscFunctionReturn(), and DM_libMesh::sys.

Referenced by DMCreate_libMesh().

DMlibMeshGetVariables()

PetscErrorCode DMlibMeshGetVariables	(	DM	dm,
		PetscInt *	n,
		char ***	varnames
	)

Definition at line 193 of file petscdmlibmeshimpl.C.

{
  PetscErrorCode ierr;
  PetscFunctionBegin;
  PetscValidHeaderSpecific(dm,DM_CLASSID,1);
  PetscBool islibmesh;
  PetscInt i;
  ierr = PetscObjectTypeCompare((PetscObject)dm, DMLIBMESH,&islibmesh);
  if (!islibmesh) SETERRQ2(((PetscObject)dm)->comm, PETSC_ERR_ARG_WRONG, "Got DM of type %s, not of type %s", ((PetscObject)dm)->type_name, DMLIBMESH);
  DM_libMesh * dlm = (DM_libMesh *)(dm->data);
  PetscValidPointer(n,2);
  *n = cast_int<unsigned int>(dlm->varids->size());
  if (!varnames) PetscFunctionReturn(0);
  ierr = PetscMalloc(*n*sizeof(char *), varnames); CHKERRQ(ierr);
  i = 0;
  for (const auto & pr : *(dlm->varids))
    {
      ierr = PetscStrallocpy(pr.first.c_str(), *varnames+i); CHKERRQ(ierr);
      ++i;
    }
  PetscFunctionReturn(0);
}

References libMesh::CHKERRQ(), dlm, libMesh::ierr, PetscFunctionReturn(), and DM_libMesh::varids.

DMlibMeshGetVec_Private()

PetscErrorCode DMlibMeshGetVec_Private	(	DM	,
		const char *	,
		Vec *
	)

Definition at line 77 of file petscdmlibmeshimpl.C.

{
  PetscFunctionBegin;
 
  PetscFunctionReturn(0);
}

References PetscFunctionReturn().

DMlibMeshJacobian()

static PetscErrorCode DMlibMeshJacobian ( #if PETSC_RELEASE_LESS_THAN(3, 5, 0) DM  dm, Vec  x, Mat  jac, Mat  pc, MatStructure *msflag #else DM  dm, Vec  x, Mat  jac, Mat pc #  endif  )

static

Definition at line 856 of file petscdmlibmeshimpl.C.

{
  PetscErrorCode ierr;
  PetscFunctionBegin;
  NonlinearImplicitSystem * _sys;
  ierr = DMlibMeshGetSystem(dm, _sys); CHKERRQ(ierr);
  NonlinearImplicitSystem & sys = *_sys;
 
  PetscMatrix<Number> the_pc(pc,sys.comm());
  PetscMatrix<Number> Jac(jac,sys.comm());
  PetscVector<Number> & X_sys = *libmesh_cast_ptr<PetscVector<Number> *>(sys.solution.get());
  PetscMatrix<Number> & Jac_sys = *libmesh_cast_ptr<PetscMatrix<Number> *>(sys.matrix);
  PetscVector<Number> X_global(x, sys.comm());
 
  // Set the dof maps
  the_pc.attach_dof_map(sys.get_dof_map());
  Jac.attach_dof_map(sys.get_dof_map());
 
  // Use the systems update() to get a good local version of the parallel solution
  X_global.swap(X_sys);
  Jac.swap(Jac_sys);
 
  sys.get_dof_map().enforce_constraints_exactly(sys);
  sys.update();
 
  X_global.swap(X_sys);
  Jac.swap(Jac_sys);
 
  the_pc.zero();
 
  // if the user has provided both function pointers and objects only the pointer
  // will be used, so catch that as an error
  if (sys.nonlinear_solver->jacobian && sys.nonlinear_solver->jacobian_object)
    libmesh_error_msg("ERROR: cannot specify both a function and object to compute the Jacobian!");
 
  if (sys.nonlinear_solver->matvec && sys.nonlinear_solver->residual_and_jacobian_object)
    libmesh_error_msg("ERROR: cannot specify both a function and object to compute the combined Residual & Jacobian!");
 
  if (sys.nonlinear_solver->jacobian != nullptr)
    sys.nonlinear_solver->jacobian(*(sys.current_local_solution.get()), the_pc, sys);
 
  else if (sys.nonlinear_solver->jacobian_object != nullptr)
    sys.nonlinear_solver->jacobian_object->jacobian(*(sys.current_local_solution.get()), the_pc, sys);
 
  else if (sys.nonlinear_solver->matvec != nullptr)
    sys.nonlinear_solver->matvec(*(sys.current_local_solution.get()), nullptr, &the_pc, sys);
 
  else if (sys.nonlinear_solver->residual_and_jacobian_object != nullptr)
    sys.nonlinear_solver->residual_and_jacobian_object->residual_and_jacobian(*(sys.current_local_solution.get()), nullptr, &the_pc, sys);
 
  else
    libmesh_error_msg("Error! Unable to compute residual and/or Jacobian!");
 
  the_pc.close();
  Jac.close();
  X_global.close();
#if PETSC_RELEASE_LESS_THAN(3,5,0)
  *msflag = SAME_NONZERO_PATTERN;
#endif
  PetscFunctionReturn(0);
}

References libMesh::CHKERRQ(), libMesh::ParallelObject::comm(), libMesh::System::current_local_solution, DMlibMeshGetSystem(), libMesh::DofMap::enforce_constraints_exactly(), libMesh::System::get_dof_map(), libMesh::ierr, libMesh::ImplicitSystem::matrix, libMesh::NonlinearImplicitSystem::nonlinear_solver, PetscFunctionReturn(), libMesh::System::solution, libMesh::PetscVector< T >::swap(), and libMesh::System::update().

Referenced by DMSetUp_libMesh(), and SNESJacobian_DMlibMesh().

DMlibMeshParseDecompositionDescriptor_Private()

static PetscErrorCode DMlibMeshParseDecompositionDescriptor_Private ( DM  dm, const char *  ddesc, PetscInt *  dtype, PetscInt *  dcount, PetscInt **  dsizes, char ****  dlists  )

static

Definition at line 601 of file petscdmlibmeshimpl.C.

{
  PetscFunctionBegin;
  PetscErrorCode ierr;
  PetscBool eq;
  char * s0;
  char * s;
  char * ss;
  struct token * llfirst = PETSC_NULL;
  struct token * lllast = PETSC_NULL;
  struct token * tok;
  PetscInt stcount = 0, brcount = 0, d, i;
  size_t len0, count;
 
  /*
    Parse the decomposition descriptor.
    Decomposition names could be of one of two forms:
    var:v1,v2;v3,v4;v4,v5;
    block:b1,b2;b3,b4;b4,b5;
    resulting in an overlapping decomposition that groups
    variables (v1,v2), (v3,v4), (v4,v5) or
    blocks    (b1,b2), (b3,b4), (b4,b5).
  */
  /* Copy the descriptor so that we can manipulate it in place. */
  ierr = PetscStrallocpy(ddesc,&s0);   CHKERRQ(ierr);
  ierr = PetscStrlen(s0, &len0)  ;     CHKERRQ(ierr);
  ierr = PetscStrstr(s0,":",&ss);      CHKERRQ(ierr);
  if (!ss) {
    ss = s0+len0;
  }
  else {
    *ss = 0;
  }
  ierr = PetscStrcmp(s0,"var",&eq);    CHKERRQ(ierr);
  if (eq) {
    *dtype=DMLIBMESH_FIELD_DECOMPOSITION;
  }
  else {
    ierr = PetscStrcmp(s0,"block",&eq);CHKERRQ(ierr);
    if (!eq)
      SETERRQ1(((PetscObject)dm)->comm, PETSC_ERR_ARG_WRONG, "Could not determine decomposition type from descriptor: %s\n", ddesc); CHKERRQ(ierr);
    *dtype=DMLIBMESH_DOMAIN_DECOMPOSITION;
  }
  ierr = PetscStrlen(s0,&count);       CHKERRQ(ierr);
  while (count < len0) {
    struct token * st;
    struct token * br;
    ++ss; ++count;
    s=ss;
    while (*ss && *ss != ',' && *ss != ';') {
      ++ss; ++count;
    }
    st = PETSC_NULL; br = PETSC_NULL;
    if (*ss) {
      /*
        Found a separator, or a break.
        Add an appropriate token to the list.
        A token separator ',' produces no token.
      */
      if (*ss == ';') {
        /* Create a break token: a token with a null string. */
#if PETSC_RELEASE_LESS_THAN(3,5,0)
        ierr = PetscNew(struct token,&br);CHKERRQ(ierr);
#else
        ierr = PetscNew(&br);CHKERRQ(ierr);
#endif
      }
      *ss = 0;
      if (s != ss) {
        /* A nonempty string. */
#if PETSC_RELEASE_LESS_THAN(3,5,0)
        ierr = PetscNew(struct token, &st);CHKERRQ(ierr);
#else
        ierr = PetscNew(&st);CHKERRQ(ierr);
#endif
        st->s = s; /* The string will be properly copied below. */
      }
      /* Add the new tokens to the list. */
      if (st) {
        if (!lllast) {
          llfirst = lllast = st;
        }
        else {
          lllast->next = st; lllast = st;
        }
      }
      if (br) {
        if (!lllast) {
          llfirst = lllast = br;
        }
        else {
          lllast->next = br; lllast = br;
        }
      }
    }
  }
  /* The result of parsing is in the linked list ll. */
  /* Count up the strings and the breaks. */
  tok = llfirst;
  while (tok) {
    if (tok->s)
      ++stcount;
    else
      ++brcount;
    tok = tok->next;
  }
  /* Allocate the space for the output. */
  *dcount = brcount;
  ierr = PetscMalloc(*dcount*sizeof(PetscInt), dsizes); CHKERRQ(ierr);
  ierr = PetscMalloc(*dcount*sizeof(char **),   dlists); CHKERRQ(ierr);
  for (d = 0; d < *dcount; ++d) (*dsizes)[d] = 0;
  tok = llfirst; d = 0;
  while (tok) {
    if (tok->s)
      ++(*dsizes)[d];
    else
      ++d;
    tok = tok->next;
  }
  for (d = 0; d < *dcount; ++d) {
    ierr = PetscMalloc(sizeof(char **)*(*dsizes)[d], (* dlists)+d); CHKERRQ(ierr);
  }
  /* Now copy strings and destroy tokens. */
  tok = llfirst; d = 0; i = 0;
  while (tok) {
    if (tok->s) {
      ierr = PetscStrallocpy(tok->s, (*dlists)[d]+i); CHKERRQ(ierr);
      ++i;
    }
    else {
      ++d;
      i = 0;
    }
    llfirst = tok;
    tok = tok->next;
    ierr = PetscFree(llfirst); CHKERRQ(ierr);
  }
  /* Deallocate workspace. */
  ierr = PetscFree(s0); CHKERRQ(ierr);
  PetscFunctionReturn(0);
}

References libMesh::CHKERRQ(), eq, libMesh::ierr, token::next, PetscFunctionReturn(), and token::s.

Referenced by DMCreateDomainDecompositionDM_libMesh(), and DMCreateFieldDecompositionDM_libMesh().

DMlibMeshSetSystem_libMesh()

PetscErrorCode DMlibMeshSetSystem_libMesh	(	DM	dm,
		NonlinearImplicitSystem &	sys
	)

Definition at line 90 of file petscdmlibmeshimpl.C.

{
  const Parallel::Communicator & comm = sys.comm();
 
  PetscErrorCode ierr;
  PetscFunctionBegin;
  PetscValidHeaderSpecific(dm,DM_CLASSID,1);
  PetscBool islibmesh;
  ierr = PetscObjectTypeCompare((PetscObject)dm, DMLIBMESH,&islibmesh);
  if (!islibmesh) SETERRQ2(((PetscObject)dm)->comm, PETSC_ERR_ARG_WRONG, "Got DM of type %s, not of type %s", ((PetscObject)dm)->type_name, DMLIBMESH);
 
  if (dm->setupcalled) SETERRQ(((PetscObject)dm)->comm, PETSC_ERR_ARG_WRONGSTATE, "Cannot reset the libMesh system after DM has been set up.");
  DM_libMesh * dlm = (DM_libMesh *)(dm->data);
  dlm->sys =&sys;
  /* Initially populate the sets of active blockids and varids using all of the
     existing blocks/variables (only variables are supported at the moment). */
  DofMap & dofmap = dlm->sys->get_dof_map();
  dlm->varids->clear();
  dlm->varnames->clear();
  for (auto v : IntRange<unsigned int>(0, dofmap.n_variables())) {
    std::string vname = dofmap.variable(v).name();
    dlm->varids->insert(std::pair<std::string,unsigned int>(vname,v));
    dlm->varnames->insert(std::pair<unsigned int,std::string>(v,vname));
  }
  const MeshBase & mesh = dlm->sys->get_mesh();
  dlm->blockids->clear();
  dlm->blocknames->clear();
  std::set<subdomain_id_type> blocks;
  /* The following effectively is a verbatim copy of MeshBase::n_subdomains(). */
  // This requires an inspection on every processor
  libmesh_parallel_only(mesh.comm());
  for (const auto & elem : mesh.active_element_ptr_range())
    blocks.insert(elem->subdomain_id());
  // Some subdomains may only live on other processors
  comm.set_union(blocks);
 
  std::set<subdomain_id_type>::iterator bit = blocks.begin();
  std::set<subdomain_id_type>::iterator bend = blocks.end();
  if (bit == bend) SETERRQ(((PetscObject)dm)->comm, PETSC_ERR_PLIB, "No mesh blocks found.");
 
  for (; bit != bend; ++bit) {
    subdomain_id_type bid = *bit;
    std::string bname = mesh.subdomain_name(bid);
    if (!bname.length()) {
      /* Block names are currently implemented for Exodus II meshes
         only, so we might have to make up our own block names and
         maintain our own mapping of block ids to names.
      */
      std::ostringstream ss;
      ss << "dm" << bid;
      bname = ss.str();
    }
    dlm->blockids->insert(std::pair<std::string,unsigned int>(bname,bid));
    dlm->blocknames->insert(std::pair<unsigned int,std::string>(bid,bname));
  }
  ierr = DMlibMeshSetUpName_Private(dm); CHKERRQ(ierr);
  PetscFunctionReturn(0);
}

References libMesh::MeshBase::active_element_ptr_range(), DM_libMesh::blockids, DM_libMesh::blocknames, libMesh::CHKERRQ(), libMesh::ParallelObject::comm(), dlm, DMlibMeshSetUpName_Private(), libMesh::System::get_dof_map(), libMesh::System::get_mesh(), libMesh::ierr, mesh, libMesh::DofMap::n_variables(), libMesh::Variable::name(), PETSC_ERR_ARG_WRONGSTATE, PetscFunctionReturn(), libMesh::MeshBase::subdomain_name(), DM_libMesh::sys, libMesh::DofMap::variable(), DM_libMesh::varids, and DM_libMesh::varnames.

Referenced by DMCreate_libMesh().

DMlibMeshSetUpName_Private()

PetscErrorCode DMlibMeshSetUpName_Private

(

DM

dm

)

Definition at line 218 of file petscdmlibmeshimpl.C.

{
  DM_libMesh * dlm = (DM_libMesh *)dm->data;
  PetscErrorCode ierr;
  PetscFunctionBegin;
  std::string name = dlm->sys->name();
  std::map<unsigned int, std::string> * dnames = PETSC_NULL, * enames = PETSC_NULL;
  if (dlm->decomposition_type == DMLIBMESH_FIELD_DECOMPOSITION) {
    name += ":dec:var:";
    dnames = dlm->varnames;
  }
  if (dlm->decomposition_type == DMLIBMESH_DOMAIN_DECOMPOSITION) {
    name += ":dec:block:";
    dnames = dlm->blocknames;
  }
  if (dnames) {
    for (auto decomp : *dlm->decomposition) {
      for (std::set<unsigned int>::iterator dit_begin = decomp.begin(),
                                            dit = dit_begin,
                                            dit_end   = decomp.end();
           dit != dit_end; ++dit) {
        unsigned int id = *dit;
        if (dit != dit_begin)
          name += ",";
        name += (*dnames)[id];
      }
      name += ";";
    }
  }
  if (dlm->embedding_type == DMLIBMESH_FIELD_EMBEDDING) {
    name += ":emb:var:";
    enames = dlm->varnames;
  }
  if (dlm->embedding_type == DMLIBMESH_DOMAIN_EMBEDDING) {
    name += ":emb:block:";
    enames = dlm->blocknames;
  }
  if (enames) {
    for (auto eit     = enames->begin(),
              eit_end = enames->end(); eit != eit_end; ++eit) {
      std::string & ename = eit->second;
      if (eit != enames->begin())
        name += ",";
      name += ename;
    }
    name += ";";
  }
  ierr = PetscObjectSetName((PetscObject)dm, name.c_str()); CHKERRQ(ierr);
  PetscFunctionReturn(0);
}

References DM_libMesh::blocknames, libMesh::CHKERRQ(), DM_libMesh::decomposition, DM_libMesh::decomposition_type, dlm, DM_libMesh::embedding_type, libMesh::ierr, libMesh::Quality::name(), libMesh::System::name(), PetscFunctionReturn(), DM_libMesh::sys, and DM_libMesh::varnames.

Referenced by DMCreateDomainDecomposition_libMesh(), DMCreateFieldDecomposition_libMesh(), DMlibMeshCreateDomainDecompositionDM(), DMlibMeshCreateFieldDecompositionDM(), and DMlibMeshSetSystem_libMesh().

DMSetUp_libMesh()

static PetscErrorCode DMSetUp_libMesh ( DM  dm )

static

Definition at line 1107 of file petscdmlibmeshimpl.C.

{
  PetscFunctionBegin;
  PetscErrorCode ierr;
  DM_libMesh     * dlm = (DM_libMesh *)(dm->data);
  PetscBool eq;
 
  ierr = PetscObjectTypeCompare((PetscObject)dm, DMLIBMESH, &eq); CHKERRQ(ierr);
 
  if (!eq)
    SETERRQ2(((PetscObject)dm)->comm, PETSC_ERR_ARG_WRONG, "DM of type %s, not of type %s", ((PetscObject)dm)->type, DMLIBMESH);
 
  if (!dlm->sys)
    SETERRQ(PETSC_COMM_WORLD, PETSC_ERR_ARG_WRONGSTATE, "No libMesh system set for DM_libMesh");
  /*
    Do not evaluate function, Jacobian or bounds for an embedded DM -- the subproblem might not have enough information for that.
  */
  if (!dlm->embedding) {
#if PETSC_RELEASE_LESS_THAN(3,3,1)
    ierr = DMSetFunction(dm, DMlibMeshFunction); CHKERRQ(ierr);
    ierr = DMSetJacobian(dm, DMlibMeshJacobian); CHKERRQ(ierr);
#else
    ierr = DMSNESSetFunction(dm, SNESFunction_DMlibMesh, (void *)dm); CHKERRQ(ierr);
    ierr = DMSNESSetJacobian(dm, SNESJacobian_DMlibMesh, (void *)dm); CHKERRQ(ierr);
#endif
    if (dlm->sys->nonlinear_solver->bounds || dlm->sys->nonlinear_solver->bounds_object)
      ierr = DMSetVariableBounds(dm, DMVariableBounds_libMesh); CHKERRQ(ierr);
  }
  else {
    /*
      Fow now we don't implement even these, although a linear "Dirichlet" subproblem is well-defined.
      Creating the submatrix, however, might require extracting the submatrix preallocation from an unassembled matrix.
    */
    dm->ops->createglobalvector = 0;
    dm->ops->creatematrix = 0;
  }
  PetscFunctionReturn(0);
}

References libMesh::CHKERRQ(), dlm, DMlibMeshFunction(), DMlibMeshJacobian(), DMVariableBounds_libMesh(), DM_libMesh::embedding, eq, libMesh::ierr, libMesh::NonlinearImplicitSystem::nonlinear_solver, PETSC_ERR_ARG_WRONGSTATE, PetscFunctionReturn(), SNESFunction_DMlibMesh(), SNESJacobian_DMlibMesh(), and DM_libMesh::sys.

Referenced by DMCreate_libMesh().

DMVariableBounds_libMesh()

static PetscErrorCode DMVariableBounds_libMesh ( DM  dm, Vec  xl, Vec  xu  )

static

Definition at line 949 of file petscdmlibmeshimpl.C.

{
  PetscErrorCode ierr;
  NonlinearImplicitSystem * _sys;
  ierr = DMlibMeshGetSystem(dm, _sys); CHKERRQ(ierr);
  NonlinearImplicitSystem & sys = *_sys;
  PetscVector<Number> XL(xl, sys.comm());
  PetscVector<Number> XU(xu, sys.comm());
  PetscFunctionBegin;
#if PETSC_VERSION_LESS_THAN(3,5,0) && PETSC_VERSION_RELEASE
  ierr = VecSet(xl, SNES_VI_NINF);CHKERRQ(ierr);
  ierr = VecSet(xu, SNES_VI_INF);CHKERRQ(ierr);
#else
  // Workaround for nonstandard Q suffix warning with quad precision
  const PetscReal petsc_inf = std::numeric_limits<PetscReal>::max() / 4;
  ierr = VecSet(xl, -petsc_inf);CHKERRQ(ierr);
  ierr = VecSet(xu, petsc_inf);CHKERRQ(ierr);
#endif
  if (sys.nonlinear_solver->bounds != nullptr)
    sys.nonlinear_solver->bounds(XL,XU,sys);
  else if (sys.nonlinear_solver->bounds_object != nullptr)
    sys.nonlinear_solver->bounds_object->bounds(XL,XU, sys);
  else
    SETERRQ(((PetscObject)dm)->comm, PETSC_ERR_ARG_WRONG, "No bounds calculation in this libMesh object");
 
  PetscFunctionReturn(0);
}

References libMesh::CHKERRQ(), libMesh::ParallelObject::comm(), DMlibMeshGetSystem(), libMesh::ierr, libMesh::NonlinearImplicitSystem::nonlinear_solver, and PetscFunctionReturn().

Referenced by DMSetUp_libMesh().

DMView_libMesh()

static PetscErrorCode DMView_libMesh ( DM  dm, PetscViewer  viewer  )

static

Definition at line 1048 of file petscdmlibmeshimpl.C.

{
  PetscErrorCode ierr;
  PetscBool isascii;
  const char * name, * prefix;
  DM_libMesh * dlm = (DM_libMesh *)dm->data;
  PetscFunctionBegin;
  ierr = PetscObjectTypeCompare((PetscObject)viewer,PETSCVIEWERASCII,&isascii); CHKERRQ(ierr);
  if (isascii) {
    ierr = PetscObjectGetName((PetscObject)dm, &name);     CHKERRQ(ierr);
    ierr = PetscObjectGetOptionsPrefix((PetscObject)dm, &prefix); CHKERRQ(ierr);
    ierr = PetscViewerASCIIPrintf(viewer, "DM libMesh with name %s and prefix %s\n", name, prefix); CHKERRQ(ierr);
    ierr = PetscViewerASCIIPrintf(viewer, "blocks:", name, prefix); CHKERRQ(ierr);
    std::map<std::string,unsigned int>::iterator bit = dlm->blockids->begin();
    std::map<std::string,unsigned int>::const_iterator bend = dlm->blockids->end();
    for (; bit != bend; ++bit) {
      ierr = PetscViewerASCIIPrintf(viewer, "(%s,%D) ", bit->first.c_str(), bit->second); CHKERRQ(ierr);
    }
    ierr = PetscViewerASCIIPrintf(viewer, "\n"); CHKERRQ(ierr);
    ierr = PetscViewerASCIIPrintf(viewer, "variables:", name, prefix); CHKERRQ(ierr);
    std::map<std::string,unsigned int>::iterator vit = dlm->varids->begin();
    std::map<std::string,unsigned int>::const_iterator vend = dlm->varids->end();
    for (; vit != vend; ++vit) {
      ierr = PetscViewerASCIIPrintf(viewer, "(%s,%D) ", vit->first.c_str(), vit->second); CHKERRQ(ierr);
    }
    ierr = PetscViewerASCIIPrintf(viewer, "\n"); CHKERRQ(ierr);
    if (dlm->decomposition_type == DMLIBMESH_NO_DECOMPOSITION) {
      ierr = PetscViewerASCIIPrintf(viewer, "No decomposition\n"); CHKERRQ(ierr);
    }
    else {
      if (dlm->decomposition_type == DMLIBMESH_FIELD_DECOMPOSITION) {
        ierr = PetscViewerASCIIPrintf(viewer, "Field decomposition by variable: "); CHKERRQ(ierr);
      }
      else if (dlm->decomposition_type == DMLIBMESH_DOMAIN_DECOMPOSITION) {
        ierr = PetscViewerASCIIPrintf(viewer, "Domain decomposition by block: "); CHKERRQ(ierr);
      }
      else SETERRQ1(((PetscObject)dm)->comm, PETSC_ERR_PLIB, "Unexpected decomposition type: %D", dlm->decomposition_type);
      /* FIX: decompositions might have different sizes and components on different ranks. */
      for (auto d : index_range(*dlm->decomposition)) {
        std::set<unsigned int>::iterator dbegin  = (*dlm->decomposition)[d].begin();
        std::set<unsigned int>::iterator dit     = (*dlm->decomposition)[d].begin();
        std::set<unsigned int>::iterator dend    = (*dlm->decomposition)[d].end();
        for (; dit != dend; ++dit) {
          if (dit != dbegin) {
            ierr = PetscViewerASCIIPrintf(viewer, ","); CHKERRQ(ierr);
          }
          ierr = PetscViewerASCIIPrintf(viewer, "%D", *dit); CHKERRQ(ierr);
        }
        ierr = PetscViewerASCIIPrintf(viewer, ";"); CHKERRQ(ierr);
      }
      ierr = PetscViewerASCIIPrintf(viewer, "\n"); CHKERRQ(ierr);
    }
  }
 
  PetscFunctionReturn(0);
}

References DM_libMesh::blockids, libMesh::CHKERRQ(), DM_libMesh::decomposition, DM_libMesh::decomposition_type, dlm, end, libMesh::ierr, libMesh::index_range(), libMesh::Quality::name(), PetscFunctionReturn(), and DM_libMesh::varids.

Referenced by DMCreate_libMesh().

if() [1/2]

if

(

!dlm->

sys

)

if() [2/2]

if

(

!

eq

)

-> comm, PETSC_ERR_ARG_WRONG, "DM of type %s, not of type %s",((PetscObject) dm) ->type, DMLIBMESH)

PetscFunctionReturn()

PetscFunctionReturn

(

0

)

Referenced by DMCreate_libMesh(), DMCreateDomainDecomposition_libMesh(), DMCreateDomainDecompositionDM_libMesh(), DMCreateFieldDecomposition_libMesh(), DMCreateFieldDecompositionDM_libMesh(), DMCreateGlobalVector_libMesh(), DMDestroy_libMesh(), DMlibMeshCreateDomainDecompositionDM(), DMlibMeshCreateFieldDecompositionDM(), DMlibMeshFunction(), DMlibMeshGetBlocks(), DMlibMeshGetSystem(), DMlibMeshGetSystem_libMesh(), DMlibMeshGetVariables(), DMlibMeshGetVec_Private(), DMlibMeshJacobian(), DMlibMeshParseDecompositionDescriptor_Private(), DMlibMeshSetSystem(), DMlibMeshSetSystem_libMesh(), DMlibMeshSetUpName_Private(), DMSetUp_libMesh(), DMVariableBounds_libMesh(), DMView_libMesh(), SNESFunction_DMlibMesh(), and SNESJacobian_DMlibMesh().

SNESFunction_DMlibMesh()

static PetscErrorCode SNESFunction_DMlibMesh ( SNES  , Vec  x, Vec  r, void *  ctx  )

static

Definition at line 843 of file petscdmlibmeshimpl.C.

{
  DM dm = (DM)ctx;
  PetscErrorCode ierr;
  PetscFunctionBegin;
  ierr = DMlibMeshFunction(dm,x,r);CHKERRQ(ierr);
  PetscFunctionReturn(0);
}

References libMesh::CHKERRQ(), libMesh::ctx, DMlibMeshFunction(), libMesh::ierr, and PetscFunctionReturn().

Referenced by DMSetUp_libMesh().

SNESJacobian_DMlibMesh()

static PetscErrorCode SNESJacobian_DMlibMesh ( #if PETSC_RELEASE_LESS_THAN(3, 5, 0)  SNES, Vec  x, Mat *  jac, Mat *  pc, MatStructure *  flag, void *ctx #else  SNES, Vec  x, Mat  jac, Mat  pc, void *ctx #  endif  )

static

Definition at line 927 of file petscdmlibmeshimpl.C.

{
  DM dm = (DM)ctx;
  PetscErrorCode ierr;
  PetscFunctionBegin;
#if PETSC_RELEASE_LESS_THAN(3,5,0)
  ierr = DMlibMeshJacobian(dm,x,*jac,*pc,flag); CHKERRQ(ierr);
#else
  ierr = DMlibMeshJacobian(dm,x,jac,pc); CHKERRQ(ierr);
#endif
  PetscFunctionReturn(0);
}

References libMesh::CHKERRQ(), libMesh::ctx, DMlibMeshJacobian(), libMesh::ierr, and PetscFunctionReturn().

Referenced by DMSetUp_libMesh().

Variable Documentation

A

* A

Initial value:

{
  PetscFunctionBegin

Definition at line 1026 of file petscdmlibmeshimpl.C.

Referenced by libMesh::DenseMatrix< Real >::_cholesky_back_substitute(), libMesh::DenseMatrix< Real >::_cholesky_decompose(), GETPOT_NAMESPACE::GetPot::_DBE_expand(), libMesh::DenseMatrix< Real >::_lu_back_substitute(), libMesh::DenseMatrix< Real >::_lu_decompose(), libMesh::DenseMatrix< Real >::_multiply_blas(), libMesh::EpetraVector< T >::add_vector(), libMesh::PetscVector< libMesh::Number >::add_vector(), libMesh::PetscVector< libMesh::Number >::add_vector_transpose(), libMesh::FEMap::compute_inverse_map_second_derivs(), libMesh::FEMap::compute_single_point_map(), libMesh::QMonomial::init_3D(), libMesh::FESubdivision::init_shape_functions(), libMesh::FESubdivision::init_subdivision_matrix(), libMesh::TypeTensor< T >::inverse(), libMesh::DenseMatrix< Real >::left_multiply_transpose(), libMesh::VariationalMeshSmoother::minJ(), libMesh::RadialBasisInterpolation< KDDim, RBF >::prepare_for_use(), libMesh::Hex::quality(), libMesh::Quad::quality(), libMesh::DenseMatrix< Real >::right_multiply_transpose(), DenseMatrixTest::testComplexSVD(), DenseMatrixTest::testEVD_helper(), DenseMatrixTest::testEVDcomplex(), DenseMatrixTest::testEVDreal(), DenseMatrixTest::testSubMatrix(), DenseMatrixTest::testSVD(), and libMesh::Edge3::volume().

dlm

DM_libMesh* dlm = (DM_libMesh *)(dm->data)

Definition at line 1029 of file petscdmlibmeshimpl.C.

Referenced by DMCreate_libMesh(), DMCreateDomainDecomposition_libMesh(), DMCreateFieldDecomposition_libMesh(), DMCreateGlobalVector_libMesh(), DMDestroy_libMesh(), DMlibMeshCreateDomainDecompositionDM(), DMlibMeshCreateFieldDecompositionDM(), DMlibMeshGetBlocks(), DMlibMeshGetSystem_libMesh(), DMlibMeshGetVariables(), DMlibMeshSetSystem_libMesh(), DMlibMeshSetUpName_Private(), DMSetUp_libMesh(), and DMView_libMesh().

DM_libMesh

No libMesh system set for DM_libMesh

Definition at line 1038 of file petscdmlibmeshimpl.C.

eq

PetscBool eq

Definition at line 1030 of file petscdmlibmeshimpl.C.

Referenced by DMCreateGlobalVector_libMesh(), DMlibMeshParseDecompositionDescriptor_Private(), and DMSetUp_libMesh().

ierr

ierr = PetscObjectTypeCompare((PetscObject)dm, DMLIBMESH, &eq)

Definition at line 1028 of file petscdmlibmeshimpl.C.

Referenced by libMesh::PetscMatrix< libMesh::Number >::_get_submatrix(), libMesh::PetscLinearSolver< Number >::_petsc_shell_matrix_get_diagonal(), libMesh::SlepcEigenSolver< libMesh::Number >::_petsc_shell_matrix_get_diagonal(), libMesh::PetscLinearSolver< Number >::_petsc_shell_matrix_mult(), libMesh::SlepcEigenSolver< libMesh::Number >::_petsc_shell_matrix_mult(), libMesh::PetscLinearSolver< Number >::_petsc_shell_matrix_mult_add(), libMesh::SlepcEigenSolver< libMesh::Number >::_solve_generalized_helper(), libMesh::SlepcEigenSolver< libMesh::Number >::_solve_standard_helper(), libMesh::PetscMatrix< libMesh::Number >::add(), libMesh::PetscMatrix< libMesh::Number >::add_block_matrix(), libMesh::PetscMatrix< libMesh::Number >::add_matrix(), libMesh::PetscLinearSolver< Number >::adjoint_solve(), libMesh::SlepcEigenSolver< libMesh::Number >::attach_deflation_space(), libMesh::PetscNonlinearSolver< Number >::build_mat_null_space(), libMesh::SlepcEigenSolver< libMesh::Number >::clear(), libMesh::PetscNonlinearSolver< Number >::clear(), libMesh::PetscLinearSolver< Number >::clear(), libMesh::PetscVector< libMesh::Number >::clear(), libMesh::PetscMatrix< libMesh::Number >::clear(), libMesh::PetscVector< libMesh::Number >::close(), libMesh::PetscMatrix< libMesh::Number >::close(), libMesh::PetscMatrix< libMesh::Number >::closed(), DMlibMeshGetSystem(), DMlibMeshSetSystem(), libMesh::PetscVector< libMesh::Number >::first_local_index(), libMesh::PetscMatrix< libMesh::Number >::flush(), libMesh::PetscNonlinearSolver< Number >::get_converged_reason(), libMesh::PetscMatrix< libMesh::Number >::get_diagonal(), libMesh::SlepcEigenSolver< libMesh::Number >::get_eigenpair(), libMesh::SlepcEigenSolver< libMesh::Number >::get_eigenvalue(), libMesh::PetscLinearSolver< Number >::get_initial_residual(), libMesh::PetscMatrix< libMesh::Number >::get_local_size(), libMesh::SlepcEigenSolver< libMesh::Number >::get_relative_error(), libMesh::PetscLinearSolver< Number >::get_residual_history(), libMesh::PetscMatrix< libMesh::Number >::get_row(), libMesh::PetscMatrix< libMesh::Number >::get_transpose(), libMesh::SlepcEigenSolver< libMesh::Number >::init(), libMesh::PetscMatrix< libMesh::Number >::init(), libMesh::PetscNonlinearSolver< Number >::init(), libMesh::PetscLinearSolver< Number >::init(), libMesh::PetscVector< libMesh::Number >::init(), libMesh::PetscMatrix< libMesh::Number >::l1_norm(), libMesh::PetscVector< libMesh::Number >::last_local_index(), libMesh::PetscMatrix< libMesh::Number >::linfty_norm(), libMesh::PetscMatrix< libMesh::Number >::local_m(), libMesh::PetscMatrix< libMesh::Number >::local_n(), libMesh::PetscVector< libMesh::Number >::local_size(), libMesh::PetscMatrix< libMesh::Number >::m(), libMesh::PetscVector< libMesh::Number >::map_global_to_local_index(), libMesh::PetscVector< libMesh::Number >::max(), libMesh::PetscVector< libMesh::Number >::min(), libMesh::PetscMatrix< libMesh::Number >::n(), libMesh::PetscMatrix< libMesh::Number >::operator()(), libMesh::PetscVector< libMesh::Number >::PetscVector(), libMesh::PetscMatrix< libMesh::Number >::print_matlab(), libMesh::PetscMatrix< libMesh::Number >::print_personal(), libMesh::PetscMatrix< libMesh::Number >::reset_preallocation(), libMesh::PetscLinearSolver< Number >::restrict_solve_to(), libMesh::PetscMatrix< libMesh::Number >::row_start(), libMesh::PetscMatrix< libMesh::Number >::row_stop(), libMesh::PetscMatrix< libMesh::Number >::set(), libMesh::SlepcEigenSolver< libMesh::Number >::set_initial_space(), libMesh::PetscLinearSolver< Number >::set_petsc_solver_type(), libMesh::SlepcEigenSolver< libMesh::Number >::set_slepc_position_of_spectrum(), libMesh::SlepcEigenSolver< libMesh::Number >::set_slepc_problem_type(), libMesh::SlepcEigenSolver< libMesh::Number >::set_slepc_solver_type(), libMesh::PetscVector< libMesh::Number >::size(), libMesh::PetscNonlinearSolver< Number >::solve(), libMesh::PetscLinearSolver< Number >::solve(), libMesh::SlepcEigenSolver< libMesh::Number >::solve_generalized(), libMesh::SlepcEigenSolver< libMesh::Number >::solve_standard(), libMesh::PetscVector< libMesh::Number >::zero(), libMesh::PetscMatrix< libMesh::Number >::zero(), and libMesh::PetscMatrix< libMesh::Number >::zero_rows().

PETSC_ERR_ARG_WRONGSTATE

PETSC_ERR_ARG_WRONGSTATE

Definition at line 1038 of file petscdmlibmeshimpl.C.

Referenced by DMCreateGlobalVector_libMesh(), DMlibMeshSetSystem_libMesh(), and DMSetUp_libMesh().