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 52e5c89e4eSSatish Balay #if !defined(PETSC_USE_DYNAMIC_LIBRARIES) 53e5c89e4eSSatish Balay ierr = PetscInitializePackage(PETSC_NULL);CHKERRQ(ierr); 54e5c89e4eSSatish Balay #endif 55d63b26baSBarry Smith ierr = PetscHeaderCreate(o,_p_PetscObject,-1,PETSC_OBJECT_COOKIE,0,"PetscObject",comm,PetscObjectDestroy_PetscObject,0);CHKERRQ(ierr); 56e5c89e4eSSatish Balay /* records not yet defined in PetscObject 57e5c89e4eSSatish Balay o->data = 0; 58e5c89e4eSSatish Balay o->setupcalled = 0; 59e5c89e4eSSatish Balay */ 60e5c89e4eSSatish Balay *obj = o; 61e5c89e4eSSatish Balay PetscFunctionReturn(0); 62e5c89e4eSSatish Balay } 63e5c89e4eSSatish Balay 64e5c89e4eSSatish Balay #undef __FUNCT__ 65*f86e92e0SMatthew Knepley #define __FUNCT__ "PetscObjectCreateGeneric" 66*f86e92e0SMatthew Knepley /*@C 67*f86e92e0SMatthew Knepley PetscObjectCreateGeneric - Creates a PetscObject 68*f86e92e0SMatthew Knepley 69*f86e92e0SMatthew Knepley Collective on PetscObject 70*f86e92e0SMatthew Knepley 71*f86e92e0SMatthew Knepley Input Parameter: 72*f86e92e0SMatthew Knepley + comm - An MPI communicator 73*f86e92e0SMatthew Knepley . cookie - The class cookie 74*f86e92e0SMatthew Knepley - name - The class name 75*f86e92e0SMatthew Knepley 76*f86e92e0SMatthew Knepley Output Parameter: 77*f86e92e0SMatthew Knepley . obj - The object 78*f86e92e0SMatthew Knepley 79*f86e92e0SMatthew Knepley Level: developer 80*f86e92e0SMatthew Knepley 81*f86e92e0SMatthew Knepley Notes: This is a template intended as a starting point to cut and paste with PetscObjectDestroy_PetscObject() 82*f86e92e0SMatthew Knepley to make new object classes. 83*f86e92e0SMatthew Knepley 84*f86e92e0SMatthew Knepley Concepts: destroying object 85*f86e92e0SMatthew Knepley Concepts: freeing object 86*f86e92e0SMatthew Knepley Concepts: deleting object 87*f86e92e0SMatthew Knepley 88*f86e92e0SMatthew Knepley @*/ 89*f86e92e0SMatthew Knepley PetscErrorCode PETSC_DLLEXPORT PetscObjectCreateGeneric(MPI_Comm comm, PetscCookie cookie, const char name[], PetscObject *obj) 90*f86e92e0SMatthew Knepley { 91*f86e92e0SMatthew Knepley PetscObject o; 92*f86e92e0SMatthew Knepley PetscErrorCode ierr; 93*f86e92e0SMatthew Knepley 94*f86e92e0SMatthew Knepley PetscFunctionBegin; 95*f86e92e0SMatthew Knepley PetscValidPointer(obj,2); 96*f86e92e0SMatthew Knepley 97*f86e92e0SMatthew Knepley #if !defined(PETSC_USE_DYNAMIC_LIBRARIES) 98*f86e92e0SMatthew Knepley ierr = PetscInitializePackage(PETSC_NULL);CHKERRQ(ierr); 99*f86e92e0SMatthew Knepley #endif 100*f86e92e0SMatthew Knepley ierr = PetscHeaderCreate(o,_p_PetscObject,-1,cookie,0,name,comm,PetscObjectDestroy_PetscObject,0);CHKERRQ(ierr); 101*f86e92e0SMatthew Knepley /* records not yet defined in PetscObject 102*f86e92e0SMatthew Knepley o->data = 0; 103*f86e92e0SMatthew Knepley o->setupcalled = 0; 104*f86e92e0SMatthew Knepley */ 105*f86e92e0SMatthew Knepley *obj = o; 106*f86e92e0SMatthew Knepley PetscFunctionReturn(0); 107*f86e92e0SMatthew Knepley } 108*f86e92e0SMatthew Knepley 109*f86e92e0SMatthew Knepley #undef __FUNCT__ 110e5c89e4eSSatish Balay #define __FUNCT__ "PetscObjectDestroy" 111e30d2299SSatish Balay /*@ 112e5c89e4eSSatish Balay PetscObjectDestroy - Destroys any PetscObject, regardless of the type. 113e5c89e4eSSatish Balay 114e5c89e4eSSatish Balay Collective on PetscObject 115e5c89e4eSSatish Balay 116e5c89e4eSSatish Balay Input Parameter: 117e5c89e4eSSatish Balay . obj - any PETSc object, for example a Vec, Mat or KSP. 118e5c89e4eSSatish Balay This must be cast with a (PetscObject), for example, 119e5c89e4eSSatish Balay PetscObjectDestroy((PetscObject)mat); 120e5c89e4eSSatish Balay 121e5c89e4eSSatish Balay Level: beginner 122e5c89e4eSSatish Balay 123e5c89e4eSSatish Balay Concepts: destroying object 124e5c89e4eSSatish Balay Concepts: freeing object 125e5c89e4eSSatish Balay Concepts: deleting object 126e5c89e4eSSatish Balay 127e5c89e4eSSatish Balay @*/ 128e5c89e4eSSatish Balay PetscErrorCode PETSC_DLLEXPORT PetscObjectDestroy(PetscObject obj) 129e5c89e4eSSatish Balay { 130e5c89e4eSSatish Balay PetscErrorCode ierr; 131e5c89e4eSSatish Balay 132e5c89e4eSSatish Balay PetscFunctionBegin; 133e5c89e4eSSatish Balay PetscValidHeader(obj,1); 134e5c89e4eSSatish Balay 135e5c89e4eSSatish Balay if (obj->bops->destroy) { 136e5c89e4eSSatish Balay ierr = (*obj->bops->destroy)(obj);CHKERRQ(ierr); 137e5c89e4eSSatish Balay } else { 138e5c89e4eSSatish Balay SETERRQ1(PETSC_ERR_PLIB,"This PETSc object of class %s does not have a generic destroy routine",obj->class_name); 139e5c89e4eSSatish Balay } 140e5c89e4eSSatish Balay PetscFunctionReturn(0); 141e5c89e4eSSatish Balay } 142e5c89e4eSSatish Balay 143e5c89e4eSSatish Balay #undef __FUNCT__ 144e5c89e4eSSatish Balay #define __FUNCT__ "PetscObjectView" 145e5c89e4eSSatish Balay /*@C 146e5c89e4eSSatish Balay PetscObjectView - Views any PetscObject, regardless of the type. 147e5c89e4eSSatish Balay 148e5c89e4eSSatish Balay Collective on PetscObject 149e5c89e4eSSatish Balay 150e5c89e4eSSatish Balay Input Parameters: 151e5c89e4eSSatish Balay + obj - any PETSc object, for example a Vec, Mat or KSP. 152e5c89e4eSSatish Balay This must be cast with a (PetscObject), for example, 153e5c89e4eSSatish Balay PetscObjectView((PetscObject)mat,viewer); 154e5c89e4eSSatish Balay - viewer - any PETSc viewer 155e5c89e4eSSatish Balay 156e5c89e4eSSatish Balay Level: intermediate 157e5c89e4eSSatish Balay 158e5c89e4eSSatish Balay @*/ 159e5c89e4eSSatish Balay PetscErrorCode PETSC_DLLEXPORT PetscObjectView(PetscObject obj,PetscViewer viewer) 160e5c89e4eSSatish Balay { 161e5c89e4eSSatish Balay PetscErrorCode ierr; 162e5c89e4eSSatish Balay 163e5c89e4eSSatish Balay PetscFunctionBegin; 164e5c89e4eSSatish Balay PetscValidHeader(obj,1); 165e5c89e4eSSatish Balay if (!viewer) viewer = PETSC_VIEWER_STDOUT_(obj->comm); 166e5c89e4eSSatish Balay PetscValidHeaderSpecific(viewer,PETSC_VIEWER_COOKIE,2); 167e5c89e4eSSatish Balay 168e5c89e4eSSatish Balay if (obj->bops->view) { 169e5c89e4eSSatish Balay ierr = (*obj->bops->view)(obj,viewer);CHKERRQ(ierr); 170e5c89e4eSSatish Balay } else { 171e5c89e4eSSatish Balay SETERRQ(PETSC_ERR_SUP,"This PETSc object does not have a generic viewer routine"); 172e5c89e4eSSatish Balay } 173e5c89e4eSSatish Balay PetscFunctionReturn(0); 174e5c89e4eSSatish Balay } 175e5c89e4eSSatish Balay 176e5c89e4eSSatish Balay #undef __FUNCT__ 177e5c89e4eSSatish Balay #define __FUNCT__ "PetscTypeCompare" 178e5c89e4eSSatish Balay /*@C 179e5c89e4eSSatish Balay PetscTypeCompare - Determines whether a PETSc object is of a particular type. 180e5c89e4eSSatish Balay 181e5c89e4eSSatish Balay Not Collective 182e5c89e4eSSatish Balay 183e5c89e4eSSatish Balay Input Parameters: 184e5c89e4eSSatish Balay + obj - any PETSc object, for example a Vec, Mat or KSP. 185e5c89e4eSSatish Balay This must be cast with a (PetscObject), for example, 186e5c89e4eSSatish Balay PetscObjectDestroy((PetscObject)mat); 187e5c89e4eSSatish Balay - type_name - string containing a type name 188e5c89e4eSSatish Balay 189e5c89e4eSSatish Balay Output Parameter: 190e5c89e4eSSatish Balay . same - PETSC_TRUE if they are the same, else PETSC_FALSE 191e5c89e4eSSatish Balay 192e5c89e4eSSatish Balay Level: intermediate 193e5c89e4eSSatish Balay 194e5c89e4eSSatish Balay .seealso: VecGetType(), KSPGetType(), PCGetType(), SNESGetType() 195e5c89e4eSSatish Balay 196e5c89e4eSSatish Balay Concepts: comparing^object types 197e5c89e4eSSatish Balay Concepts: types^comparing 198e5c89e4eSSatish Balay Concepts: object type^comparing 199e5c89e4eSSatish Balay 200e5c89e4eSSatish Balay @*/ 201e5c89e4eSSatish Balay PetscErrorCode PETSC_DLLEXPORT PetscTypeCompare(PetscObject obj,const char type_name[],PetscTruth *same) 202e5c89e4eSSatish Balay { 203e5c89e4eSSatish Balay PetscErrorCode ierr; 204e5c89e4eSSatish Balay 205e5c89e4eSSatish Balay PetscFunctionBegin; 206e5c89e4eSSatish Balay if (!obj) { 207e5c89e4eSSatish Balay *same = PETSC_FALSE; 208e5c89e4eSSatish Balay } else if (!type_name && !obj->type_name) { 209e5c89e4eSSatish Balay *same = PETSC_TRUE; 210e5c89e4eSSatish Balay } else if (!type_name || !obj->type_name) { 211e5c89e4eSSatish Balay *same = PETSC_FALSE; 212e5c89e4eSSatish Balay } else { 213e5c89e4eSSatish Balay PetscValidHeader(obj,1); 214e5c89e4eSSatish Balay PetscValidCharPointer(type_name,2); 215e5c89e4eSSatish Balay PetscValidPointer(same,3); 216e5c89e4eSSatish Balay ierr = PetscStrcmp((char*)(obj->type_name),type_name,same);CHKERRQ(ierr); 217e5c89e4eSSatish Balay } 218e5c89e4eSSatish Balay PetscFunctionReturn(0); 219e5c89e4eSSatish Balay } 220e5c89e4eSSatish Balay 221e5c89e4eSSatish Balay static int PetscObjectRegisterDestroy_Count = 0; 222e5c89e4eSSatish Balay static PetscObject PetscObjectRegisterDestroy_Objects[256]; 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); 249e5c89e4eSSatish Balay PetscObjectRegisterDestroy_Objects[PetscObjectRegisterDestroy_Count++] = obj; 250e5c89e4eSSatish Balay PetscFunctionReturn(0); 251e5c89e4eSSatish Balay } 252e5c89e4eSSatish Balay 253e5c89e4eSSatish Balay #undef __FUNCT__ 254e5c89e4eSSatish Balay #define __FUNCT__ "PetscObjectRegisterDestroyAll" 255e5c89e4eSSatish Balay /*@C 256e5c89e4eSSatish Balay PetscObjectRegisterDestroyAll - Frees all the PETSc objects that have been registered 257e5c89e4eSSatish Balay with PetscObjectRegisterDestroy(). Called by PetscFinalize() 258e5c89e4eSSatish Balay PetscFinalize() is called. 259e5c89e4eSSatish Balay 260e5c89e4eSSatish Balay Collective on individual PetscObjects 261e5c89e4eSSatish Balay 262e5c89e4eSSatish Balay Level: developer 263e5c89e4eSSatish Balay 264e5c89e4eSSatish Balay .seealso: PetscObjectRegisterDestroy() 265e5c89e4eSSatish Balay @*/ 266e5c89e4eSSatish Balay PetscErrorCode PETSC_DLLEXPORT PetscObjectRegisterDestroyAll(void) 267e5c89e4eSSatish Balay { 268e5c89e4eSSatish Balay PetscErrorCode ierr; 269e5c89e4eSSatish Balay int i; 270e5c89e4eSSatish Balay 271e5c89e4eSSatish Balay PetscFunctionBegin; 272e5c89e4eSSatish Balay for (i=0; i<PetscObjectRegisterDestroy_Count; i++) { 273e5c89e4eSSatish Balay ierr = PetscObjectDestroy(PetscObjectRegisterDestroy_Objects[i]);CHKERRQ(ierr); 274e5c89e4eSSatish Balay } 275e5c89e4eSSatish Balay PetscFunctionReturn(0); 276e5c89e4eSSatish Balay } 277e5c89e4eSSatish Balay 278e5c89e4eSSatish Balay 279