/*@C
DMMGDestroy - Destroys a DA based multigrid solver object.
Collective on DMMG
Input Parameter:
. - the context
Level: advanced
.seealso DMMGCreate()
@*/
PetscErrorCode PETSCSNES_DLLEXPORT DMMGDestroy(DMMG *dmmg)
{
PetscErrorCode ierr;
PetscInt i,nlevels = dmmg[0]->nlevels;
PetscFunctionBegin;
if (!dmmg) SETERRQ(PETSC_ERR_ARG_NULL,"Passing null as DMMG");
for (i=1; i<nlevels; i++) {
if (dmmg[i]->R) {ierr = MatDestroy(dmmg[i]->R);CHKERRQ(ierr);}
}
for (i=0; i<nlevels; i++) {
ierr = PetscStrfree(dmmg[i]->prefix);CHKERRQ(ierr);
ierr = PetscStrfree(dmmg[i]->mtype);CHKERRQ(ierr);
if (dmmg[i]->dm) {ierr = DMDestroy(dmmg[i]->dm);CHKERRQ(ierr);}
if (dmmg[i]->x) {ierr = VecDestroy(dmmg[i]->x);CHKERRQ(ierr);}
if (dmmg[i]->b) {ierr = VecDestroy(dmmg[i]->b);CHKERRQ(ierr);}
if (dmmg[i]->r) {ierr = VecDestroy(dmmg[i]->r);CHKERRQ(ierr);}
if (dmmg[i]->work1) {ierr = VecDestroy(dmmg[i]->work1);CHKERRQ(ierr);}
if (dmmg[i]->w) {ierr = VecDestroy(dmmg[i]->w);CHKERRQ(ierr);}
if (dmmg[i]->work2) {ierr = VecDestroy(dmmg[i]->work2);CHKERRQ(ierr);}
if (dmmg[i]->lwork1) {ierr = VecDestroy(dmmg[i]->lwork1);CHKERRQ(ierr);}
if (dmmg[i]->B) {ierr = MatDestroy(dmmg[i]->B);CHKERRQ(ierr);}
if (dmmg[i]->J) {ierr = MatDestroy(dmmg[i]->J);CHKERRQ(ierr);}
if (dmmg[i]->Rscale) {ierr = VecDestroy(dmmg[i]->Rscale);CHKERRQ(ierr);}
if (dmmg[i]->fdcoloring){ierr = MatFDColoringDestroy(dmmg[i]->fdcoloring);CHKERRQ(ierr);}
if (dmmg[i]->ksp && !dmmg[i]->snes) {ierr = KSPDestroy(dmmg[i]->ksp);CHKERRQ(ierr);}
if (dmmg[i]->snes) {ierr = PetscObjectDestroy((PetscObject)dmmg[i]->snes);CHKERRQ(ierr);}
if (dmmg[i]->inject) {ierr = VecScatterDestroy(dmmg[i]->inject);CHKERRQ(ierr);}
ierr = PetscFree(dmmg[i]);CHKERRQ(ierr);
}
ierr = PetscFree(dmmg);CHKERRQ(ierr);
PetscFunctionReturn(0);
}
static PetscErrorCode PCDestroy_LU(PC pc)
{
PC_LU *dir = (PC_LU*)pc->data;
PetscErrorCode ierr;
PetscFunctionBegin;
if (!dir->inplace && ((PC_Factor*)dir)->fact) {ierr = MatDestroy(((PC_Factor*)dir)->fact);CHKERRQ(ierr);}
if (dir->row && dir->col && dir->row != dir->col) {ierr = ISDestroy(dir->row);CHKERRQ(ierr);}
if (dir->col) {ierr = ISDestroy(dir->col);CHKERRQ(ierr);}
ierr = PetscStrfree(((PC_Factor*)dir)->ordering);CHKERRQ(ierr);
ierr = PetscStrfree(((PC_Factor*)dir)->solvertype);CHKERRQ(ierr);
ierr = PetscFree(dir);CHKERRQ(ierr);
PetscFunctionReturn(0);
}
PetscErrorCode PFDestroy_String(void *value)
{
PetscErrorCode ierr;
PetscFunctionBegin;
ierr = PetscStrfree(value);CHKERRQ(ierr);
PetscFunctionReturn(0);
}
/*@C
PetscObjectSetName - Sets a string name associated with a PETSc object.
Not Collective
Input Parameters:
+ obj - the Petsc variable
Thus must be cast with a (PetscObject), for example,
PetscObjectSetName((PetscObject)mat,name);
- name - the name to give obj
Level: advanced
Concepts: object name^setting
.seealso: PetscObjectGetName()
@*/
PetscErrorCode PETSC_DLLEXPORT PetscObjectSetName(PetscObject obj,const char name[])
{
PetscErrorCode ierr;
PetscFunctionBegin;
PetscValidHeader(obj,1);
ierr = PetscStrfree(obj->name);CHKERRQ(ierr);
ierr = PetscStrallocpy(name,&obj->name);CHKERRQ(ierr);
PetscFunctionReturn(0);
}
PetscErrorCode PFDestroy_Matlab(void *value)
{
PetscErrorCode ierr;
PF_Matlab *matlab = (PF_Matlab*)value;
PetscFunctionBegin;
ierr = PetscStrfree(matlab->string);CHKERRQ(ierr);
ierr = PetscMatlabEngineDestroy(matlab->mengine);CHKERRQ(ierr);
ierr = PetscFree(matlab);CHKERRQ(ierr);
PetscFunctionReturn(0);
}
PetscErrorCode PETSC_DLLEXPORT PetscObjectChangeTypeName(PetscObject obj,const char type_name[])
{
PetscErrorCode ierr;
PetscFunctionBegin;
PetscValidHeader(obj,1);
ierr = PetscObjectTakeAccess(obj);CHKERRQ(ierr);
ierr = PetscStrfree(obj->type_name);CHKERRQ(ierr);
ierr = PetscStrallocpy(type_name,&obj->type_name);CHKERRQ(ierr);
ierr = PetscObjectGrantAccess(obj);CHKERRQ(ierr);
PetscFunctionReturn(0);
}
static PetscErrorCode PCDestroy_Shell(PC pc)
{
PC_Shell *shell = (PC_Shell*)pc->data;
PetscErrorCode ierr;
PetscFunctionBegin;
ierr = PetscStrfree(shell->name);CHKERRQ(ierr);
if (shell->destroy) {
ierr = (*shell->destroy)(pc);CHKERRQ(ierr);
}
ierr = PetscFree(shell);CHKERRQ(ierr);
PetscFunctionReturn(0);
}
/*
PFStringCreateFunction - Creates a function from a string
Collective over PF
Input Parameters:
+ pf - the function object
- string - the string that defines the function
Output Parameter:
. f - the function pointer.
.seealso: PFSetFromOptions()
*/
PetscErrorCode PETSCVEC_DLLEXPORT PFStringCreateFunction(PF pf,char *string,void **f)
{
#if defined(PETSC_USE_DYNAMIC_LIBRARIES)
PetscErrorCode ierr;
char task[1024],tmp[256],lib[PETSC_MAX_PATH_LEN],username[64];
FILE *fd;
PetscTruth tmpshared,wdshared,keeptmpfiles = PETSC_FALSE;
MPI_Comm comm;
#endif
PetscFunctionBegin;
#if defined(PETSC_USE_DYNAMIC_LIBRARIES)
ierr = PetscStrfree(pf->data);CHKERRQ(ierr);
ierr = PetscStrallocpy(string,(char**)&pf->data);CHKERRQ(ierr);
/* create the new C function and compile it */
ierr = PetscSharedTmp(((PetscObject)pf)->comm,&tmpshared);CHKERRQ(ierr);
ierr = PetscSharedWorkingDirectory(((PetscObject)pf)->comm,&wdshared);CHKERRQ(ierr);
if (tmpshared) { /* do it in /tmp since everyone has one */
ierr = PetscGetTmp(((PetscObject)pf)->comm,tmp,PETSC_MAX_PATH_LEN);CHKERRQ(ierr);
comm = ((PetscObject)pf)->comm;
} else if (!wdshared) { /* each one does in private /tmp */
ierr = PetscGetTmp(((PetscObject)pf)->comm,tmp,PETSC_MAX_PATH_LEN);CHKERRQ(ierr);
comm = PETSC_COMM_SELF;
} else { /* do it in current directory */
ierr = PetscStrcpy(tmp,".");CHKERRQ(ierr);
comm = ((PetscObject)pf)->comm;
}
ierr = PetscOptionsGetTruth(((PetscObject)pf)->prefix,"-pf_string_keep_files",&keeptmpfiles,PETSC_NULL);CHKERRQ(ierr);
if (keeptmpfiles) {
sprintf(task,"cd %s ; mkdir ${USERNAME} ; cd ${USERNAME} ; \\cp -f ${PETSC_DIR}/src/pf/impls/string/makefile ./makefile ; ke MIN=%d NOUT=%d petscdlib STRINGFUNCTION=\"%s\" ; sync\n",tmp,(int)pf->dimin,(int)pf->dimout,string);
} else {
sprintf(task,"cd %s ; mkdir ${USERNAME} ;cd ${USERNAME} ; \\cp -f ${PETSC_DIR}/src/pf/impls/string/makefile ./makefile ; make MIN=%d NOUT=%d -f makefile petscdlib STRINGFUNCTION=\"%s\" ; \\rm -f makefile petscdlib.c libpetscdlib.a ; sync\n",tmp,(int)pf->dimin,(int)pf->dimout,string);
}
#if defined(PETSC_HAVE_POPEN)
ierr = PetscPOpen(comm,PETSC_NULL,task,"r",&fd);CHKERRQ(ierr);
ierr = PetscPClose(comm,fd);CHKERRQ(ierr);
#else
SETERRQ(PETSC_ERR_SUP_SYS,"Cannot run external programs on this machine");
#endif
ierr = MPI_Barrier(comm);CHKERRQ(ierr);
/* load the apply function from the dynamic library */
ierr = PetscGetUserName(username,64);CHKERRQ(ierr);
sprintf(lib,"%s/%s/libpetscdlib",tmp,username);
ierr = PetscDLLibrarySym(comm,PETSC_NULL,lib,"PFApply_String",f);CHKERRQ(ierr);
#endif
PetscFunctionReturn(0);
}
/*
PetscHeaderDestroy_Private - Destroys a base PETSc object header. Called by
the macro PetscHeaderDestroy().
*/
PetscErrorCode PETSC_DLLEXPORT PetscHeaderDestroy_Private(PetscObject h)
{
PetscErrorCode ierr;
PetscFunctionBegin;
if (PetscMemoryCollectMaximumUsage) {
PetscLogDouble usage;
ierr = PetscMemoryGetCurrentUsage(&usage);CHKERRQ(ierr);
if (usage > PetscMemoryMaximumUsage) PetscMemoryMaximumUsage = usage;
}
/* first destroy things that could execute arbitrary code */
if (h->python_destroy) {
void *python_context = h->python_context;
PetscErrorCode (*python_destroy)(void*) = h->python_destroy;
h->python_context = 0;
h->python_destroy = 0;
ierr = (*python_destroy)(python_context);CHKERRQ(ierr);
}
ierr = PetscOListDestroy(h->olist);CHKERRQ(ierr);
ierr = PetscCommDestroy(&h->comm);CHKERRQ(ierr);
/* next destroy other things */
h->cookie = PETSCFREEDHEADER;
ierr = PetscFree(h->bops);CHKERRQ(ierr);
ierr = PetscFListDestroy(&h->qlist);CHKERRQ(ierr);
ierr = PetscStrfree(h->type_name);CHKERRQ(ierr);
ierr = PetscStrfree(h->name);CHKERRQ(ierr);
ierr = PetscStrfree(h->prefix);CHKERRQ(ierr);
ierr = PetscFree(h->fortran_func_pointers);CHKERRQ(ierr);
ierr = PetscFree(h->intcomposeddata);CHKERRQ(ierr);
ierr = PetscFree(h->intcomposedstate);CHKERRQ(ierr);
ierr = PetscFree(h->realcomposeddata);CHKERRQ(ierr);
ierr = PetscFree(h->realcomposedstate);CHKERRQ(ierr);
ierr = PetscFree(h->scalarcomposeddata);CHKERRQ(ierr);
ierr = PetscFree(h->scalarcomposedstate);CHKERRQ(ierr);
PetscFunctionReturn(0);
}
EXTERN_C_END
EXTERN_C_BEGIN
#undef __FUNCT__
#define __FUNCT__ "PCShellSetName_Shell"
PetscErrorCode PETSCKSP_DLLEXPORT PCShellSetName_Shell(PC pc,const char name[])
{
PC_Shell *shell;
PetscErrorCode ierr;
PetscFunctionBegin;
shell = (PC_Shell*)pc->data;
ierr = PetscStrfree(shell->name);CHKERRQ(ierr);
ierr = PetscStrallocpy(name,&shell->name);CHKERRQ(ierr);
PetscFunctionReturn(0);
}
/*@
DADestroy - Destroys a distributed array.
Collective on DA
Input Parameter:
. da - the distributed array to destroy
Level: beginner
.keywords: distributed array, destroy
.seealso: DACreate1d(), DACreate2d(), DACreate3d()
@*/
PetscErrorCode PETSCDM_DLLEXPORT DADestroy(DA da)
{
PetscErrorCode ierr;
PetscErrorCode i;
PetscTruth done;
PetscFunctionBegin;
PetscValidHeaderSpecific(da,DM_COOKIE,1);
ierr = DMDestroy_Private((DM)da,&done);CHKERRQ(ierr);
if (!done) PetscFunctionReturn(0);
/* destroy the internal part */
if (da->ops->destroy) {
ierr = (*da->ops->destroy)(da);CHKERRQ(ierr);
}
/* destroy the external/common part */
for (i=0; i<DA_MAX_AD_ARRAYS; i++) {
ierr = PetscFree(da->adstartghostedout[i]);CHKERRQ(ierr);
ierr = PetscFree(da->adstartghostedin[i]);CHKERRQ(ierr);
ierr = PetscFree(da->adstartout[i]);CHKERRQ(ierr);
ierr = PetscFree(da->adstartin[i]);CHKERRQ(ierr);
}
for (i=0; i<DA_MAX_AD_ARRAYS; i++) {
ierr = PetscFree(da->admfstartghostedout[i]);CHKERRQ(ierr);
ierr = PetscFree(da->admfstartghostedin[i]);CHKERRQ(ierr);
ierr = PetscFree(da->admfstartout[i]);CHKERRQ(ierr);
ierr = PetscFree(da->admfstartin[i]);CHKERRQ(ierr);
}
for (i=0; i<DA_MAX_WORK_ARRAYS; i++) {
ierr = PetscFree(da->startghostedout[i]);CHKERRQ(ierr);
ierr = PetscFree(da->startghostedin[i]);CHKERRQ(ierr);
ierr = PetscFree(da->startout[i]);CHKERRQ(ierr);
ierr = PetscFree(da->startin[i]);CHKERRQ(ierr);
}
/* if memory was published with AMS then destroy it */
ierr = PetscObjectDepublish(da);CHKERRQ(ierr);
if (da->ltog) {ierr = VecScatterDestroy(da->ltog);CHKERRQ(ierr);}
if (da->gtol) {ierr = VecScatterDestroy(da->gtol);CHKERRQ(ierr);}
if (da->ltol) {ierr = VecScatterDestroy(da->ltol);CHKERRQ(ierr);}
if (da->natural){
ierr = VecDestroy(da->natural);CHKERRQ(ierr);
}
if (da->gton) {
ierr = VecScatterDestroy(da->gton);CHKERRQ(ierr);
}
if (da->ao) {
ierr = AODestroy(da->ao);CHKERRQ(ierr);
}
if (da->ltogmap) {
ierr = ISLocalToGlobalMappingDestroy(da->ltogmap);CHKERRQ(ierr);
}
if (da->ltogmapb) {
ierr = ISLocalToGlobalMappingDestroy(da->ltogmapb);CHKERRQ(ierr);
}
ierr = PetscFree(da->lx);CHKERRQ(ierr);
ierr = PetscFree(da->ly);CHKERRQ(ierr);
ierr = PetscFree(da->lz);CHKERRQ(ierr);
if (da->fieldname) {
for (i=0; i<da->w; i++) {
ierr = PetscStrfree(da->fieldname[i]);CHKERRQ(ierr);
}
ierr = PetscFree(da->fieldname);CHKERRQ(ierr);
}
if (da->localcoloring) {
ierr = ISColoringDestroy(da->localcoloring);CHKERRQ(ierr);
}
if (da->ghostedcoloring) {
ierr = ISColoringDestroy(da->ghostedcoloring);CHKERRQ(ierr);
}
if (da->coordinates) {ierr = VecDestroy(da->coordinates);CHKERRQ(ierr);}
if (da->ghosted_coordinates) {ierr = VecDestroy(da->ghosted_coordinates);CHKERRQ(ierr);}
if (da->da_coordinates && da != da->da_coordinates) {ierr = DADestroy(da->da_coordinates);CHKERRQ(ierr);}
ierr = PetscFree(da->neighbors);CHKERRQ(ierr);
ierr = PetscFree(da->dfill);CHKERRQ(ierr);
ierr = PetscFree(da->ofill);CHKERRQ(ierr);
ierr = PetscFree(da->e);CHKERRQ(ierr);
ierr = PetscHeaderDestroy(da);CHKERRQ(ierr);
PetscFunctionReturn(0);
}
