1e5c89e4eSSatish Balay #define PETSC_DLL 2e5c89e4eSSatish Balay /* 3e5c89e4eSSatish Balay Provides utility routines for manulating any type of PETSc object. 4e5c89e4eSSatish Balay */ 5e5c89e4eSSatish Balay #include "petsc.h" /*I "petsc.h" I*/ 6e5c89e4eSSatish Balay 7e5c89e4eSSatish Balay PetscCookie PETSC_OBJECT_COOKIE = 0; 8e5c89e4eSSatish Balay 9e5c89e4eSSatish Balay struct _p_Object { 10e5c89e4eSSatish Balay PETSCHEADER(int); 11e5c89e4eSSatish Balay }; 12e5c89e4eSSatish Balay 13e5c89e4eSSatish Balay PetscErrorCode PetscObjectDestroy_Private(PetscObject obj) 14e5c89e4eSSatish Balay { 15e5c89e4eSSatish Balay return 0; 16e5c89e4eSSatish Balay } 17e5c89e4eSSatish Balay 18e5c89e4eSSatish Balay #undef __FUNCT__ 19e5c89e4eSSatish Balay #define __FUNCT__ "PetscObjectCreate" 20e5c89e4eSSatish Balay /*@C 21e5c89e4eSSatish Balay PetscObjectCreate - Creates a PetscObject 22e5c89e4eSSatish Balay 23e5c89e4eSSatish Balay Collective on PetscObject 24e5c89e4eSSatish Balay 25e5c89e4eSSatish Balay Input Parameter: 26e5c89e4eSSatish Balay . comm - An MPI communicator 27e5c89e4eSSatish Balay 28e5c89e4eSSatish Balay Output Parameter: 29e5c89e4eSSatish Balay . obj - The object 30e5c89e4eSSatish Balay 31e5c89e4eSSatish Balay Level: beginner 32e5c89e4eSSatish Balay 33e5c89e4eSSatish Balay Concepts: destroying object 34e5c89e4eSSatish Balay Concepts: freeing object 35e5c89e4eSSatish Balay Concepts: deleting object 36e5c89e4eSSatish Balay 37e5c89e4eSSatish Balay @*/ 38e5c89e4eSSatish Balay PetscErrorCode PETSC_DLLEXPORT PetscObjectCreate(MPI_Comm comm, PetscObject *obj) 39e5c89e4eSSatish Balay { 40e5c89e4eSSatish Balay PetscObject o; 41e5c89e4eSSatish Balay PetscErrorCode ierr; 42e5c89e4eSSatish Balay 43e5c89e4eSSatish Balay PetscFunctionBegin; 44e5c89e4eSSatish Balay PetscValidPointer(obj,2); 45e5c89e4eSSatish Balay 46e5c89e4eSSatish Balay #if !defined(PETSC_USE_DYNAMIC_LIBRARIES) 47e5c89e4eSSatish Balay ierr = PetscInitializePackage(PETSC_NULL);CHKERRQ(ierr); 48e5c89e4eSSatish Balay #endif 49e5c89e4eSSatish Balay ierr = PetscHeaderCreate(o,_p_PetscObject,-1,PETSC_OBJECT_COOKIE,0,"PetscObject",comm,PetscObjectDestroy_Private,0);CHKERRQ(ierr); 50e5c89e4eSSatish Balay /* records not yet defined in PetscObject 51e5c89e4eSSatish Balay o->data = 0; 52e5c89e4eSSatish Balay o->setupcalled = 0; 53e5c89e4eSSatish Balay */ 54e5c89e4eSSatish Balay *obj = o; 55e5c89e4eSSatish Balay PetscFunctionReturn(0); 56e5c89e4eSSatish Balay } 57e5c89e4eSSatish Balay 58e5c89e4eSSatish Balay #undef __FUNCT__ 59e5c89e4eSSatish Balay #define __FUNCT__ "PetscObjectDestroy" 60*e30d2299SSatish Balay /*@ 61e5c89e4eSSatish Balay PetscObjectDestroy - Destroys any PetscObject, regardless of the type. 62e5c89e4eSSatish Balay 63e5c89e4eSSatish Balay Collective on PetscObject 64e5c89e4eSSatish Balay 65e5c89e4eSSatish Balay Input Parameter: 66e5c89e4eSSatish Balay . obj - any PETSc object, for example a Vec, Mat or KSP. 67e5c89e4eSSatish Balay This must be cast with a (PetscObject), for example, 68e5c89e4eSSatish Balay PetscObjectDestroy((PetscObject)mat); 69e5c89e4eSSatish Balay 70e5c89e4eSSatish Balay Level: beginner 71e5c89e4eSSatish Balay 72e5c89e4eSSatish Balay Concepts: destroying object 73e5c89e4eSSatish Balay Concepts: freeing object 74e5c89e4eSSatish Balay Concepts: deleting object 75e5c89e4eSSatish Balay 76e5c89e4eSSatish Balay @*/ 77e5c89e4eSSatish Balay PetscErrorCode PETSC_DLLEXPORT PetscObjectDestroy(PetscObject obj) 78e5c89e4eSSatish Balay { 79e5c89e4eSSatish Balay PetscErrorCode ierr; 80e5c89e4eSSatish Balay 81e5c89e4eSSatish Balay PetscFunctionBegin; 82e5c89e4eSSatish Balay PetscValidHeader(obj,1); 83e5c89e4eSSatish Balay 84e5c89e4eSSatish Balay if (obj->bops->destroy) { 85e5c89e4eSSatish Balay ierr = (*obj->bops->destroy)(obj);CHKERRQ(ierr); 86e5c89e4eSSatish Balay } else { 87e5c89e4eSSatish Balay SETERRQ1(PETSC_ERR_PLIB,"This PETSc object of class %s does not have a generic destroy routine",obj->class_name); 88e5c89e4eSSatish Balay } 89e5c89e4eSSatish Balay PetscFunctionReturn(0); 90e5c89e4eSSatish Balay } 91e5c89e4eSSatish Balay 92e5c89e4eSSatish Balay #undef __FUNCT__ 93e5c89e4eSSatish Balay #define __FUNCT__ "PetscObjectView" 94e5c89e4eSSatish Balay /*@C 95e5c89e4eSSatish Balay PetscObjectView - Views any PetscObject, regardless of the type. 96e5c89e4eSSatish Balay 97e5c89e4eSSatish Balay Collective on PetscObject 98e5c89e4eSSatish Balay 99e5c89e4eSSatish Balay Input Parameters: 100e5c89e4eSSatish Balay + obj - any PETSc object, for example a Vec, Mat or KSP. 101e5c89e4eSSatish Balay This must be cast with a (PetscObject), for example, 102e5c89e4eSSatish Balay PetscObjectView((PetscObject)mat,viewer); 103e5c89e4eSSatish Balay - viewer - any PETSc viewer 104e5c89e4eSSatish Balay 105e5c89e4eSSatish Balay Level: intermediate 106e5c89e4eSSatish Balay 107e5c89e4eSSatish Balay @*/ 108e5c89e4eSSatish Balay PetscErrorCode PETSC_DLLEXPORT PetscObjectView(PetscObject obj,PetscViewer viewer) 109e5c89e4eSSatish Balay { 110e5c89e4eSSatish Balay PetscErrorCode ierr; 111e5c89e4eSSatish Balay 112e5c89e4eSSatish Balay PetscFunctionBegin; 113e5c89e4eSSatish Balay PetscValidHeader(obj,1); 114e5c89e4eSSatish Balay if (!viewer) viewer = PETSC_VIEWER_STDOUT_(obj->comm); 115e5c89e4eSSatish Balay PetscValidHeaderSpecific(viewer,PETSC_VIEWER_COOKIE,2); 116e5c89e4eSSatish Balay 117e5c89e4eSSatish Balay if (obj->bops->view) { 118e5c89e4eSSatish Balay ierr = (*obj->bops->view)(obj,viewer);CHKERRQ(ierr); 119e5c89e4eSSatish Balay } else { 120e5c89e4eSSatish Balay SETERRQ(PETSC_ERR_SUP,"This PETSc object does not have a generic viewer routine"); 121e5c89e4eSSatish Balay } 122e5c89e4eSSatish Balay PetscFunctionReturn(0); 123e5c89e4eSSatish Balay } 124e5c89e4eSSatish Balay 125e5c89e4eSSatish Balay #undef __FUNCT__ 126e5c89e4eSSatish Balay #define __FUNCT__ "PetscTypeCompare" 127e5c89e4eSSatish Balay /*@C 128e5c89e4eSSatish Balay PetscTypeCompare - Determines whether a PETSc object is of a particular type. 129e5c89e4eSSatish Balay 130e5c89e4eSSatish Balay Not Collective 131e5c89e4eSSatish Balay 132e5c89e4eSSatish Balay Input Parameters: 133e5c89e4eSSatish Balay + obj - any PETSc object, for example a Vec, Mat or KSP. 134e5c89e4eSSatish Balay This must be cast with a (PetscObject), for example, 135e5c89e4eSSatish Balay PetscObjectDestroy((PetscObject)mat); 136e5c89e4eSSatish Balay - type_name - string containing a type name 137e5c89e4eSSatish Balay 138e5c89e4eSSatish Balay Output Parameter: 139e5c89e4eSSatish Balay . same - PETSC_TRUE if they are the same, else PETSC_FALSE 140e5c89e4eSSatish Balay 141e5c89e4eSSatish Balay Level: intermediate 142e5c89e4eSSatish Balay 143e5c89e4eSSatish Balay .seealso: VecGetType(), KSPGetType(), PCGetType(), SNESGetType() 144e5c89e4eSSatish Balay 145e5c89e4eSSatish Balay Concepts: comparing^object types 146e5c89e4eSSatish Balay Concepts: types^comparing 147e5c89e4eSSatish Balay Concepts: object type^comparing 148e5c89e4eSSatish Balay 149e5c89e4eSSatish Balay @*/ 150e5c89e4eSSatish Balay PetscErrorCode PETSC_DLLEXPORT PetscTypeCompare(PetscObject obj,const char type_name[],PetscTruth *same) 151e5c89e4eSSatish Balay { 152e5c89e4eSSatish Balay PetscErrorCode ierr; 153e5c89e4eSSatish Balay 154e5c89e4eSSatish Balay PetscFunctionBegin; 155e5c89e4eSSatish Balay if (!obj) { 156e5c89e4eSSatish Balay *same = PETSC_FALSE; 157e5c89e4eSSatish Balay } else if (!type_name && !obj->type_name) { 158e5c89e4eSSatish Balay *same = PETSC_TRUE; 159e5c89e4eSSatish Balay } else if (!type_name || !obj->type_name) { 160e5c89e4eSSatish Balay *same = PETSC_FALSE; 161e5c89e4eSSatish Balay } else { 162e5c89e4eSSatish Balay PetscValidHeader(obj,1); 163e5c89e4eSSatish Balay PetscValidCharPointer(type_name,2); 164e5c89e4eSSatish Balay PetscValidPointer(same,3); 165e5c89e4eSSatish Balay ierr = PetscStrcmp((char*)(obj->type_name),type_name,same);CHKERRQ(ierr); 166e5c89e4eSSatish Balay } 167e5c89e4eSSatish Balay PetscFunctionReturn(0); 168e5c89e4eSSatish Balay } 169e5c89e4eSSatish Balay 170e5c89e4eSSatish Balay static int PetscObjectRegisterDestroy_Count = 0; 171e5c89e4eSSatish Balay static PetscObject PetscObjectRegisterDestroy_Objects[256]; 172e5c89e4eSSatish Balay 173e5c89e4eSSatish Balay #undef __FUNCT__ 174e5c89e4eSSatish Balay #define __FUNCT__ "PetscObjectRegisterDestroy" 175e5c89e4eSSatish Balay /*@C 176e5c89e4eSSatish Balay PetscObjectRegisterDestroy - Registers a PETSc object to be destroyed when 177e5c89e4eSSatish Balay PetscFinalize() is called. 178e5c89e4eSSatish Balay 179e5c89e4eSSatish Balay Collective on PetscObject 180e5c89e4eSSatish Balay 181e5c89e4eSSatish Balay Input Parameter: 182e5c89e4eSSatish Balay . obj - any PETSc object, for example a Vec, Mat or KSP. 183e5c89e4eSSatish Balay This must be cast with a (PetscObject), for example, 184e5c89e4eSSatish Balay PetscObjectRegisterDestroy((PetscObject)mat); 185e5c89e4eSSatish Balay 186e5c89e4eSSatish Balay Level: developer 187e5c89e4eSSatish Balay 188e5c89e4eSSatish Balay Notes: 189e5c89e4eSSatish Balay This is used by, for example, PETSC_VIEWER_XXX_() routines to free the viewer 190e5c89e4eSSatish Balay when PETSc ends. 191e5c89e4eSSatish Balay 192e5c89e4eSSatish Balay .seealso: PetscObjectRegisterDestroyAll() 193e5c89e4eSSatish Balay @*/ 194e5c89e4eSSatish Balay PetscErrorCode PETSC_DLLEXPORT PetscObjectRegisterDestroy(PetscObject obj) 195e5c89e4eSSatish Balay { 196e5c89e4eSSatish Balay PetscFunctionBegin; 197e5c89e4eSSatish Balay PetscValidHeader(obj,1); 198e5c89e4eSSatish Balay PetscObjectRegisterDestroy_Objects[PetscObjectRegisterDestroy_Count++] = obj; 199e5c89e4eSSatish Balay PetscFunctionReturn(0); 200e5c89e4eSSatish Balay } 201e5c89e4eSSatish Balay 202e5c89e4eSSatish Balay #undef __FUNCT__ 203e5c89e4eSSatish Balay #define __FUNCT__ "PetscObjectRegisterDestroyAll" 204e5c89e4eSSatish Balay /*@C 205e5c89e4eSSatish Balay PetscObjectRegisterDestroyAll - Frees all the PETSc objects that have been registered 206e5c89e4eSSatish Balay with PetscObjectRegisterDestroy(). Called by PetscFinalize() 207e5c89e4eSSatish Balay PetscFinalize() is called. 208e5c89e4eSSatish Balay 209e5c89e4eSSatish Balay Collective on individual PetscObjects 210e5c89e4eSSatish Balay 211e5c89e4eSSatish Balay Level: developer 212e5c89e4eSSatish Balay 213e5c89e4eSSatish Balay .seealso: PetscObjectRegisterDestroy() 214e5c89e4eSSatish Balay @*/ 215e5c89e4eSSatish Balay PetscErrorCode PETSC_DLLEXPORT PetscObjectRegisterDestroyAll(void) 216e5c89e4eSSatish Balay { 217e5c89e4eSSatish Balay PetscErrorCode ierr; 218e5c89e4eSSatish Balay int i; 219e5c89e4eSSatish Balay 220e5c89e4eSSatish Balay PetscFunctionBegin; 221e5c89e4eSSatish Balay for (i=0; i<PetscObjectRegisterDestroy_Count; i++) { 222e5c89e4eSSatish Balay ierr = PetscObjectDestroy(PetscObjectRegisterDestroy_Objects[i]);CHKERRQ(ierr); 223e5c89e4eSSatish Balay } 224e5c89e4eSSatish Balay PetscFunctionReturn(0); 225e5c89e4eSSatish Balay } 226e5c89e4eSSatish Balay 227e5c89e4eSSatish Balay 228