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 struct _p_Object { 8e5c89e4eSSatish Balay PETSCHEADER(int); 9e5c89e4eSSatish Balay }; 10e5c89e4eSSatish Balay 11d63b26baSBarry Smith PetscErrorCode PetscObjectDestroy_PetscObject(PetscObject obj) 12e5c89e4eSSatish Balay { 13d63b26baSBarry Smith PetscErrorCode ierr; 14d63b26baSBarry Smith PetscFunctionBegin; 15d63b26baSBarry Smith PetscValidHeader(obj,1); 16d63b26baSBarry Smith if (--obj->refct > 0) PetscFunctionReturn(0); 17d63b26baSBarry Smith ierr = PetscHeaderDestroy(obj);CHKERRQ(ierr); 18d63b26baSBarry Smith PetscFunctionReturn(0); 19e5c89e4eSSatish Balay } 20e5c89e4eSSatish Balay 21e5c89e4eSSatish Balay #undef __FUNCT__ 22e5c89e4eSSatish Balay #define __FUNCT__ "PetscObjectCreate" 23e5c89e4eSSatish Balay /*@C 24e5c89e4eSSatish Balay PetscObjectCreate - Creates a PetscObject 25e5c89e4eSSatish Balay 26e5c89e4eSSatish Balay Collective on PetscObject 27e5c89e4eSSatish Balay 28e5c89e4eSSatish Balay Input Parameter: 29e5c89e4eSSatish Balay . comm - An MPI communicator 30e5c89e4eSSatish Balay 31e5c89e4eSSatish Balay Output Parameter: 32e5c89e4eSSatish Balay . obj - The object 33e5c89e4eSSatish Balay 34d63b26baSBarry Smith Level: developer 35d63b26baSBarry Smith 36d63b26baSBarry Smith Notes: This is a template intended as a starting point to cut and paste with PetscObjectDestroy_PetscObject() 37d63b26baSBarry Smith to make new object classes. 38e5c89e4eSSatish Balay 39e5c89e4eSSatish Balay Concepts: destroying object 40e5c89e4eSSatish Balay Concepts: freeing object 41e5c89e4eSSatish Balay Concepts: deleting object 42e5c89e4eSSatish Balay 43e5c89e4eSSatish Balay @*/ 44e5c89e4eSSatish Balay PetscErrorCode PETSC_DLLEXPORT PetscObjectCreate(MPI_Comm comm, PetscObject *obj) 45e5c89e4eSSatish Balay { 46e5c89e4eSSatish Balay PetscObject o; 47e5c89e4eSSatish Balay PetscErrorCode ierr; 48e5c89e4eSSatish Balay 49e5c89e4eSSatish Balay PetscFunctionBegin; 50e5c89e4eSSatish Balay PetscValidPointer(obj,2); 51e5c89e4eSSatish Balay #if !defined(PETSC_USE_DYNAMIC_LIBRARIES) 52e5c89e4eSSatish Balay ierr = PetscInitializePackage(PETSC_NULL);CHKERRQ(ierr); 53e5c89e4eSSatish Balay #endif 54d63b26baSBarry Smith ierr = PetscHeaderCreate(o,_p_PetscObject,-1,PETSC_OBJECT_COOKIE,0,"PetscObject",comm,PetscObjectDestroy_PetscObject,0);CHKERRQ(ierr); 55e5c89e4eSSatish Balay /* records not yet defined in PetscObject 56e5c89e4eSSatish Balay o->data = 0; 57e5c89e4eSSatish Balay o->setupcalled = 0; 58e5c89e4eSSatish Balay */ 59e5c89e4eSSatish Balay *obj = o; 60e5c89e4eSSatish Balay PetscFunctionReturn(0); 61e5c89e4eSSatish Balay } 62e5c89e4eSSatish Balay 63e5c89e4eSSatish Balay #undef __FUNCT__ 64f86e92e0SMatthew Knepley #define __FUNCT__ "PetscObjectCreateGeneric" 65f86e92e0SMatthew Knepley /*@C 66f86e92e0SMatthew Knepley PetscObjectCreateGeneric - Creates a PetscObject 67f86e92e0SMatthew Knepley 68f86e92e0SMatthew Knepley Collective on PetscObject 69f86e92e0SMatthew Knepley 70f86e92e0SMatthew Knepley Input Parameter: 71f86e92e0SMatthew Knepley + comm - An MPI communicator 72f86e92e0SMatthew Knepley . cookie - The class cookie 73f86e92e0SMatthew Knepley - name - The class name 74f86e92e0SMatthew Knepley 75f86e92e0SMatthew Knepley Output Parameter: 76f86e92e0SMatthew Knepley . obj - The object 77f86e92e0SMatthew Knepley 78f86e92e0SMatthew Knepley Level: developer 79f86e92e0SMatthew Knepley 80f86e92e0SMatthew Knepley Notes: This is a template intended as a starting point to cut and paste with PetscObjectDestroy_PetscObject() 81f86e92e0SMatthew Knepley to make new object classes. 82f86e92e0SMatthew Knepley 83f86e92e0SMatthew Knepley Concepts: destroying object 84f86e92e0SMatthew Knepley Concepts: freeing object 85f86e92e0SMatthew Knepley Concepts: deleting object 86f86e92e0SMatthew Knepley 87f86e92e0SMatthew Knepley @*/ 88f86e92e0SMatthew Knepley PetscErrorCode PETSC_DLLEXPORT PetscObjectCreateGeneric(MPI_Comm comm, PetscCookie cookie, const char name[], PetscObject *obj) 89f86e92e0SMatthew Knepley { 90f86e92e0SMatthew Knepley PetscObject o; 91f86e92e0SMatthew Knepley PetscErrorCode ierr; 92f86e92e0SMatthew Knepley 93f86e92e0SMatthew Knepley PetscFunctionBegin; 94f86e92e0SMatthew Knepley PetscValidPointer(obj,2); 95f86e92e0SMatthew Knepley #if !defined(PETSC_USE_DYNAMIC_LIBRARIES) 96f86e92e0SMatthew Knepley ierr = PetscInitializePackage(PETSC_NULL);CHKERRQ(ierr); 97f86e92e0SMatthew Knepley #endif 98f86e92e0SMatthew Knepley ierr = PetscHeaderCreate(o,_p_PetscObject,-1,cookie,0,name,comm,PetscObjectDestroy_PetscObject,0);CHKERRQ(ierr); 99f86e92e0SMatthew Knepley /* records not yet defined in PetscObject 100f86e92e0SMatthew Knepley o->data = 0; 101f86e92e0SMatthew Knepley o->setupcalled = 0; 102f86e92e0SMatthew Knepley */ 103f86e92e0SMatthew Knepley *obj = o; 104f86e92e0SMatthew Knepley PetscFunctionReturn(0); 105f86e92e0SMatthew Knepley } 106f86e92e0SMatthew Knepley 107f86e92e0SMatthew Knepley #undef __FUNCT__ 108e5c89e4eSSatish Balay #define __FUNCT__ "PetscObjectDestroy" 109e30d2299SSatish Balay /*@ 110e5c89e4eSSatish Balay PetscObjectDestroy - Destroys any PetscObject, regardless of the type. 111e5c89e4eSSatish Balay 112e5c89e4eSSatish Balay Collective on PetscObject 113e5c89e4eSSatish Balay 114e5c89e4eSSatish Balay Input Parameter: 115e5c89e4eSSatish Balay . obj - any PETSc object, for example a Vec, Mat or KSP. 116e5c89e4eSSatish Balay This must be cast with a (PetscObject), for example, 117e5c89e4eSSatish Balay PetscObjectDestroy((PetscObject)mat); 118e5c89e4eSSatish Balay 119e5c89e4eSSatish Balay Level: beginner 120e5c89e4eSSatish Balay 121e5c89e4eSSatish Balay Concepts: destroying object 122e5c89e4eSSatish Balay Concepts: freeing object 123e5c89e4eSSatish Balay Concepts: deleting object 124e5c89e4eSSatish Balay 125e5c89e4eSSatish Balay @*/ 126e5c89e4eSSatish Balay PetscErrorCode PETSC_DLLEXPORT PetscObjectDestroy(PetscObject obj) 127e5c89e4eSSatish Balay { 128e5c89e4eSSatish Balay PetscErrorCode ierr; 129e5c89e4eSSatish Balay 130e5c89e4eSSatish Balay PetscFunctionBegin; 131e5c89e4eSSatish Balay PetscValidHeader(obj,1); 132e5c89e4eSSatish Balay if (obj->bops->destroy) { 133e5c89e4eSSatish Balay ierr = (*obj->bops->destroy)(obj);CHKERRQ(ierr); 134e5c89e4eSSatish Balay } else { 135e5c89e4eSSatish Balay SETERRQ1(PETSC_ERR_PLIB,"This PETSc object of class %s does not have a generic destroy routine",obj->class_name); 136e5c89e4eSSatish Balay } 137e5c89e4eSSatish Balay PetscFunctionReturn(0); 138e5c89e4eSSatish Balay } 139e5c89e4eSSatish Balay 140e5c89e4eSSatish Balay #undef __FUNCT__ 141e5c89e4eSSatish Balay #define __FUNCT__ "PetscObjectView" 142e5c89e4eSSatish Balay /*@C 143e5c89e4eSSatish Balay PetscObjectView - Views any PetscObject, regardless of the type. 144e5c89e4eSSatish Balay 145e5c89e4eSSatish Balay Collective on PetscObject 146e5c89e4eSSatish Balay 147e5c89e4eSSatish Balay Input Parameters: 148e5c89e4eSSatish Balay + obj - any PETSc object, for example a Vec, Mat or KSP. 149e5c89e4eSSatish Balay This must be cast with a (PetscObject), for example, 150e5c89e4eSSatish Balay PetscObjectView((PetscObject)mat,viewer); 151e5c89e4eSSatish Balay - viewer - any PETSc viewer 152e5c89e4eSSatish Balay 153e5c89e4eSSatish Balay Level: intermediate 154e5c89e4eSSatish Balay 155e5c89e4eSSatish Balay @*/ 156e5c89e4eSSatish Balay PetscErrorCode PETSC_DLLEXPORT PetscObjectView(PetscObject obj,PetscViewer viewer) 157e5c89e4eSSatish Balay { 158e5c89e4eSSatish Balay PetscErrorCode ierr; 159e5c89e4eSSatish Balay 160e5c89e4eSSatish Balay PetscFunctionBegin; 161e5c89e4eSSatish Balay PetscValidHeader(obj,1); 1623050cee2SBarry Smith if (!viewer) { 1633050cee2SBarry Smith ierr = PetscViewerASCIIGetStdout(obj->comm,&viewer);CHKERRQ(ierr); 1643050cee2SBarry Smith } 165e5c89e4eSSatish Balay PetscValidHeaderSpecific(viewer,PETSC_VIEWER_COOKIE,2); 166e5c89e4eSSatish Balay 167e5c89e4eSSatish Balay if (obj->bops->view) { 168e5c89e4eSSatish Balay ierr = (*obj->bops->view)(obj,viewer);CHKERRQ(ierr); 169e5c89e4eSSatish Balay } else { 170e5c89e4eSSatish Balay SETERRQ(PETSC_ERR_SUP,"This PETSc object does not have a generic viewer routine"); 171e5c89e4eSSatish Balay } 172e5c89e4eSSatish Balay PetscFunctionReturn(0); 173e5c89e4eSSatish Balay } 174e5c89e4eSSatish Balay 175e5c89e4eSSatish Balay #undef __FUNCT__ 176e5c89e4eSSatish Balay #define __FUNCT__ "PetscTypeCompare" 177e5c89e4eSSatish Balay /*@C 178e5c89e4eSSatish Balay PetscTypeCompare - Determines whether a PETSc object is of a particular type. 179e5c89e4eSSatish Balay 180e5c89e4eSSatish Balay Not Collective 181e5c89e4eSSatish Balay 182e5c89e4eSSatish Balay Input Parameters: 183e5c89e4eSSatish Balay + obj - any PETSc object, for example a Vec, Mat or KSP. 184e5c89e4eSSatish Balay This must be cast with a (PetscObject), for example, 185e5c89e4eSSatish Balay PetscObjectDestroy((PetscObject)mat); 186e5c89e4eSSatish Balay - type_name - string containing a type name 187e5c89e4eSSatish Balay 188e5c89e4eSSatish Balay Output Parameter: 189e5c89e4eSSatish Balay . same - PETSC_TRUE if they are the same, else PETSC_FALSE 190e5c89e4eSSatish Balay 191e5c89e4eSSatish Balay Level: intermediate 192e5c89e4eSSatish Balay 193e5c89e4eSSatish Balay .seealso: VecGetType(), KSPGetType(), PCGetType(), SNESGetType() 194e5c89e4eSSatish Balay 195e5c89e4eSSatish Balay Concepts: comparing^object types 196e5c89e4eSSatish Balay Concepts: types^comparing 197e5c89e4eSSatish Balay Concepts: object type^comparing 198e5c89e4eSSatish Balay 199e5c89e4eSSatish Balay @*/ 200e5c89e4eSSatish Balay PetscErrorCode PETSC_DLLEXPORT PetscTypeCompare(PetscObject obj,const char type_name[],PetscTruth *same) 201e5c89e4eSSatish Balay { 202e5c89e4eSSatish Balay PetscErrorCode ierr; 203e5c89e4eSSatish Balay 204e5c89e4eSSatish Balay PetscFunctionBegin; 205e5c89e4eSSatish Balay if (!obj) { 206e5c89e4eSSatish Balay *same = PETSC_FALSE; 207e5c89e4eSSatish Balay } else if (!type_name && !obj->type_name) { 208e5c89e4eSSatish Balay *same = PETSC_TRUE; 209e5c89e4eSSatish Balay } else if (!type_name || !obj->type_name) { 210e5c89e4eSSatish Balay *same = PETSC_FALSE; 211e5c89e4eSSatish Balay } else { 212e5c89e4eSSatish Balay PetscValidHeader(obj,1); 213e5c89e4eSSatish Balay PetscValidCharPointer(type_name,2); 214e5c89e4eSSatish Balay PetscValidPointer(same,3); 215e5c89e4eSSatish Balay ierr = PetscStrcmp((char*)(obj->type_name),type_name,same);CHKERRQ(ierr); 216e5c89e4eSSatish Balay } 217e5c89e4eSSatish Balay PetscFunctionReturn(0); 218e5c89e4eSSatish Balay } 219e5c89e4eSSatish Balay 220*3cfa8680SLisandro Dalcin #define MAXREGDESOBJS 256 221e5c89e4eSSatish Balay static int PetscObjectRegisterDestroy_Count = 0; 222*3cfa8680SLisandro Dalcin static PetscObject PetscObjectRegisterDestroy_Objects[MAXREGDESOBJS]; 223e5c89e4eSSatish Balay 224e5c89e4eSSatish Balay #undef __FUNCT__ 225e5c89e4eSSatish Balay #define __FUNCT__ "PetscObjectRegisterDestroy" 226e5c89e4eSSatish Balay /*@C 227e5c89e4eSSatish Balay PetscObjectRegisterDestroy - Registers a PETSc object to be destroyed when 228e5c89e4eSSatish Balay PetscFinalize() is called. 229e5c89e4eSSatish Balay 230e5c89e4eSSatish Balay Collective on PetscObject 231e5c89e4eSSatish Balay 232e5c89e4eSSatish Balay Input Parameter: 233e5c89e4eSSatish Balay . obj - any PETSc object, for example a Vec, Mat or KSP. 234e5c89e4eSSatish Balay This must be cast with a (PetscObject), for example, 235e5c89e4eSSatish Balay PetscObjectRegisterDestroy((PetscObject)mat); 236e5c89e4eSSatish Balay 237e5c89e4eSSatish Balay Level: developer 238e5c89e4eSSatish Balay 239e5c89e4eSSatish Balay Notes: 240e5c89e4eSSatish Balay This is used by, for example, PETSC_VIEWER_XXX_() routines to free the viewer 241e5c89e4eSSatish Balay when PETSc ends. 242e5c89e4eSSatish Balay 243e5c89e4eSSatish Balay .seealso: PetscObjectRegisterDestroyAll() 244e5c89e4eSSatish Balay @*/ 245e5c89e4eSSatish Balay PetscErrorCode PETSC_DLLEXPORT PetscObjectRegisterDestroy(PetscObject obj) 246e5c89e4eSSatish Balay { 247e5c89e4eSSatish Balay PetscFunctionBegin; 248e5c89e4eSSatish Balay PetscValidHeader(obj,1); 249*3cfa8680SLisandro Dalcin if (PetscObjectRegisterDestroy_Count < MAXREGDESOBJS) { 250e5c89e4eSSatish Balay PetscObjectRegisterDestroy_Objects[PetscObjectRegisterDestroy_Count++] = obj; 251*3cfa8680SLisandro Dalcin } else { 252*3cfa8680SLisandro Dalcin SETERRQ1(PETSC_ERR_PLIB,"No more room in array, limit %d \n recompile src/sys/objects/destroy.c with larger value for MAXREGDESOBJS\n",MAXREGDESOBJS); 253*3cfa8680SLisandro Dalcin 254*3cfa8680SLisandro Dalcin } 255e5c89e4eSSatish Balay PetscFunctionReturn(0); 256e5c89e4eSSatish Balay } 257e5c89e4eSSatish Balay 258e5c89e4eSSatish Balay #undef __FUNCT__ 259e5c89e4eSSatish Balay #define __FUNCT__ "PetscObjectRegisterDestroyAll" 260e5c89e4eSSatish Balay /*@C 261e5c89e4eSSatish Balay PetscObjectRegisterDestroyAll - Frees all the PETSc objects that have been registered 262e5c89e4eSSatish Balay with PetscObjectRegisterDestroy(). Called by PetscFinalize() 263e5c89e4eSSatish Balay PetscFinalize() is called. 264e5c89e4eSSatish Balay 265e5c89e4eSSatish Balay Collective on individual PetscObjects 266e5c89e4eSSatish Balay 267e5c89e4eSSatish Balay Level: developer 268e5c89e4eSSatish Balay 269e5c89e4eSSatish Balay .seealso: PetscObjectRegisterDestroy() 270e5c89e4eSSatish Balay @*/ 271e5c89e4eSSatish Balay PetscErrorCode PETSC_DLLEXPORT PetscObjectRegisterDestroyAll(void) 272e5c89e4eSSatish Balay { 273e5c89e4eSSatish Balay PetscErrorCode ierr; 274e5c89e4eSSatish Balay int i; 275e5c89e4eSSatish Balay 276e5c89e4eSSatish Balay PetscFunctionBegin; 277e5c89e4eSSatish Balay for (i=0; i<PetscObjectRegisterDestroy_Count; i++) { 278e5c89e4eSSatish Balay ierr = PetscObjectDestroy(PetscObjectRegisterDestroy_Objects[i]);CHKERRQ(ierr); 279e5c89e4eSSatish Balay } 280*3cfa8680SLisandro Dalcin PetscObjectRegisterDestroy_Count = 0; 281e5c89e4eSSatish Balay PetscFunctionReturn(0); 282e5c89e4eSSatish Balay } 283e5c89e4eSSatish Balay 284e5c89e4eSSatish Balay 285