/* Provides utility routines for manulating any type of PETSc object. */ #include /*I "petscsys.h" I*/ #include PetscErrorCode PetscComposedQuantitiesDestroy(PetscObject obj) { PetscErrorCode ierr; PetscInt i; PetscFunctionBegin; if (obj->intstar_idmax>0) { for (i=0; iintstar_idmax; i++) { ierr = PetscFree(obj->intstarcomposeddata[i]);CHKERRQ(ierr); } ierr = PetscFree2(obj->intstarcomposeddata,obj->intstarcomposedstate);CHKERRQ(ierr); } if (obj->realstar_idmax>0) { for (i=0; irealstar_idmax; i++) { ierr = PetscFree(obj->realstarcomposeddata[i]);CHKERRQ(ierr); } ierr = PetscFree2(obj->realstarcomposeddata,obj->realstarcomposedstate);CHKERRQ(ierr); } if (obj->scalarstar_idmax>0) { for (i=0; iscalarstar_idmax; i++) { ierr = PetscFree(obj->scalarstarcomposeddata[i]);CHKERRQ(ierr); } ierr = PetscFree2(obj->scalarstarcomposeddata,obj->scalarstarcomposedstate);CHKERRQ(ierr); } ierr = PetscFree2(obj->intcomposeddata,obj->intcomposedstate);CHKERRQ(ierr); ierr = PetscFree2(obj->realcomposeddata,obj->realcomposedstate);CHKERRQ(ierr); ierr = PetscFree2(obj->scalarcomposeddata,obj->scalarcomposedstate);CHKERRQ(ierr); PetscFunctionReturn(0); } /*@ PetscObjectDestroy - Destroys any PetscObject, regardless of the type. Collective on PetscObject Input Parameter: . obj - any PETSc object, for example a Vec, Mat or KSP. This must be cast with a (PetscObject*), for example, PetscObjectDestroy((PetscObject*)&mat); Level: beginner @*/ PetscErrorCode PetscObjectDestroy(PetscObject *obj) { PetscErrorCode ierr; PetscFunctionBegin; if (!*obj) PetscFunctionReturn(0); PetscValidHeader(*obj,1); if (*obj && (*obj)->bops->destroy) { ierr = (*(*obj)->bops->destroy)(obj);CHKERRQ(ierr); } else if (*obj) SETERRQ1(PETSC_COMM_SELF,PETSC_ERR_PLIB,"This PETSc object of class %s does not have a generic destroy routine",(*obj)->class_name); PetscFunctionReturn(0); } /*@C PetscObjectView - Views any PetscObject, regardless of the type. Collective on PetscObject Input Parameters: + obj - any PETSc object, for example a Vec, Mat or KSP. This must be cast with a (PetscObject), for example, PetscObjectView((PetscObject)mat,viewer); - viewer - any PETSc viewer Level: intermediate @*/ PetscErrorCode PetscObjectView(PetscObject obj,PetscViewer viewer) { PetscErrorCode ierr; PetscFunctionBegin; PetscValidHeader(obj,1); if (!viewer) { ierr = PetscViewerASCIIGetStdout(obj->comm,&viewer);CHKERRQ(ierr); } PetscValidHeaderSpecific(viewer,PETSC_VIEWER_CLASSID,2); if (obj->bops->view) { ierr = (*obj->bops->view)(obj,viewer);CHKERRQ(ierr); } else SETERRQ(PETSC_COMM_SELF,PETSC_ERR_SUP,"This PETSc object does not have a generic viewer routine"); PetscFunctionReturn(0); } #define CHKERRQI(incall,ierr) if (ierr) {incall = PETSC_FALSE; CHKERRQ(ierr);} /*@C PetscObjectViewFromOptions - Processes command line options to determine if/how a PetscObject is to be viewed. Collective on PetscObject Input Parameters: + obj - the object . bobj - optional other object that provides prefix (if NULL then the prefix in obj is used) - optionname - option to activate viewing Level: intermediate @*/ PetscErrorCode PetscObjectViewFromOptions(PetscObject obj,PetscObject bobj,const char optionname[]) { PetscErrorCode ierr; PetscViewer viewer; PetscBool flg; static PetscBool incall = PETSC_FALSE; PetscViewerFormat format; const char *prefix; PetscFunctionBegin; if (incall) PetscFunctionReturn(0); incall = PETSC_TRUE; prefix = bobj ? bobj->prefix : obj->prefix; ierr = PetscOptionsGetViewer(PetscObjectComm((PetscObject)obj),obj->options,prefix,optionname,&viewer,&format,&flg);CHKERRQI(incall,ierr); if (flg) { ierr = PetscViewerPushFormat(viewer,format);CHKERRQI(incall,ierr); ierr = PetscObjectView(obj,viewer);CHKERRQI(incall,ierr); ierr = PetscViewerPopFormat(viewer);CHKERRQI(incall,ierr); ierr = PetscViewerDestroy(&viewer);CHKERRQI(incall,ierr); } incall = PETSC_FALSE; PetscFunctionReturn(0); } /*@C PetscObjectTypeCompare - Determines whether a PETSc object is of a particular type. Not Collective Input Parameters: + obj - any PETSc object, for example a Vec, Mat or KSP. This must be cast with a (PetscObject), for example, PetscObjectTypeCompare((PetscObject)mat); - type_name - string containing a type name Output Parameter: . same - PETSC_TRUE if they are the same, else PETSC_FALSE Level: intermediate .seealso: VecGetType(), KSPGetType(), PCGetType(), SNESGetType(), PetscObjectBaseTypeCompare(), PetscObjectTypeCompareAny() @*/ PetscErrorCode PetscObjectTypeCompare(PetscObject obj,const char type_name[],PetscBool *same) { PetscErrorCode ierr; PetscFunctionBegin; PetscValidPointer(same,3); if (!obj) *same = PETSC_FALSE; else if (!type_name && !obj->type_name) *same = PETSC_TRUE; else if (!type_name || !obj->type_name) *same = PETSC_FALSE; else { PetscValidHeader(obj,1); PetscValidCharPointer(type_name,2); ierr = PetscStrcmp((char*)(obj->type_name),type_name,same);CHKERRQ(ierr); } PetscFunctionReturn(0); } /*@C PetscObjectBaseTypeCompare - Determines whether a PetscObject is of a given base type. For example the base type of MATSEQAIJPERM is MATSEQAIJ Not Collective Input Parameters: + mat - the matrix - type_name - string containing a type name Output Parameter: . same - PETSC_TRUE if it is of the same base type Level: intermediate .seealso: PetscObjectTypeCompare(), PetscObjectTypeCompareAny() @*/ PetscErrorCode PetscObjectBaseTypeCompare(PetscObject obj,const char type_name[],PetscBool *same) { PetscErrorCode ierr; PetscFunctionBegin; PetscValidPointer(same,3); if (!obj) *same = PETSC_FALSE; else if (!type_name && !obj->type_name) *same = PETSC_TRUE; else if (!type_name || !obj->type_name) *same = PETSC_FALSE; else { PetscValidHeader(obj,1); PetscValidCharPointer(type_name,2); ierr = PetscStrbeginswith((char*)(obj->type_name),type_name,same);CHKERRQ(ierr); } PetscFunctionReturn(0); } /*@C PetscObjectTypeCompareAny - Determines whether a PETSc object is of any of a list of types. Not Collective Input Parameters: + obj - any PETSc object, for example a Vec, Mat or KSP. This must be cast with a (PetscObject), for example, PetscObjectTypeCompareAny((PetscObject)mat,...); - type_name - string containing a type name, pass the empty string "" to terminate the list Output Parameter: . match - PETSC_TRUE if the type of obj matches any in the list, else PETSC_FALSE Level: intermediate .seealso: VecGetType(), KSPGetType(), PCGetType(), SNESGetType(), PetscObjectTypeCompare(), PetscObjectBaseTypeCompare() @*/ PetscErrorCode PetscObjectTypeCompareAny(PetscObject obj,PetscBool *match,const char type_name[],...) { PetscErrorCode ierr; va_list Argp; PetscFunctionBegin; PetscValidPointer(match,3); *match = PETSC_FALSE; va_start(Argp,type_name); while (type_name && type_name[0]) { PetscBool found; ierr = PetscObjectTypeCompare(obj,type_name,&found);CHKERRQ(ierr); if (found) { *match = PETSC_TRUE; break; } type_name = va_arg(Argp,const char*); } va_end(Argp); PetscFunctionReturn(0); } #define MAXREGDESOBJS 256 static int PetscObjectRegisterDestroy_Count = 0; static PetscObject PetscObjectRegisterDestroy_Objects[MAXREGDESOBJS]; /*@C PetscObjectRegisterDestroy - Registers a PETSc object to be destroyed when PetscFinalize() is called. Logically Collective on PetscObject Input Parameter: . obj - any PETSc object, for example a Vec, Mat or KSP. This must be cast with a (PetscObject), for example, PetscObjectRegisterDestroy((PetscObject)mat); Level: developer Notes: This is used by, for example, PETSC_VIEWER_XXX_() routines to free the viewer when PETSc ends. .seealso: PetscObjectRegisterDestroyAll() @*/ PetscErrorCode PetscObjectRegisterDestroy(PetscObject obj) { PetscFunctionBegin; PetscValidHeader(obj,1); if (PetscObjectRegisterDestroy_Count < MAXREGDESOBJS) PetscObjectRegisterDestroy_Objects[PetscObjectRegisterDestroy_Count++] = obj; else SETERRQ1(PETSC_COMM_SELF,PETSC_ERR_PLIB,"No more room in array, limit %d \n recompile src/sys/objects/destroy.c with larger value for MAXREGDESOBJS\n",MAXREGDESOBJS); PetscFunctionReturn(0); } /*@C PetscObjectRegisterDestroyAll - Frees all the PETSc objects that have been registered with PetscObjectRegisterDestroy(). Called by PetscFinalize() Logically Collective on individual PetscObjects Level: developer .seealso: PetscObjectRegisterDestroy() @*/ PetscErrorCode PetscObjectRegisterDestroyAll(void) { PetscErrorCode ierr; PetscInt i; PetscFunctionBegin; for (i=0; i