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__ 65f86e92e0SMatthew Knepley #define __FUNCT__ "PetscObjectCreateGeneric" 66f86e92e0SMatthew Knepley /*@C 67f86e92e0SMatthew Knepley PetscObjectCreateGeneric - Creates a PetscObject 68f86e92e0SMatthew Knepley 69f86e92e0SMatthew Knepley Collective on PetscObject 70f86e92e0SMatthew Knepley 71f86e92e0SMatthew Knepley Input Parameter: 72f86e92e0SMatthew Knepley + comm - An MPI communicator 73f86e92e0SMatthew Knepley . cookie - The class cookie 74f86e92e0SMatthew Knepley - name - The class name 75f86e92e0SMatthew Knepley 76f86e92e0SMatthew Knepley Output Parameter: 77f86e92e0SMatthew Knepley . obj - The object 78f86e92e0SMatthew Knepley 79f86e92e0SMatthew Knepley Level: developer 80f86e92e0SMatthew Knepley 81f86e92e0SMatthew Knepley Notes: This is a template intended as a starting point to cut and paste with PetscObjectDestroy_PetscObject() 82f86e92e0SMatthew Knepley to make new object classes. 83f86e92e0SMatthew Knepley 84f86e92e0SMatthew Knepley Concepts: destroying object 85f86e92e0SMatthew Knepley Concepts: freeing object 86f86e92e0SMatthew Knepley Concepts: deleting object 87f86e92e0SMatthew Knepley 88f86e92e0SMatthew Knepley @*/ 89f86e92e0SMatthew Knepley PetscErrorCode PETSC_DLLEXPORT PetscObjectCreateGeneric(MPI_Comm comm, PetscCookie cookie, const char name[], PetscObject *obj) 90f86e92e0SMatthew Knepley { 91f86e92e0SMatthew Knepley PetscObject o; 92f86e92e0SMatthew Knepley PetscErrorCode ierr; 93f86e92e0SMatthew Knepley 94f86e92e0SMatthew Knepley PetscFunctionBegin; 95f86e92e0SMatthew Knepley PetscValidPointer(obj,2); 96f86e92e0SMatthew Knepley 97f86e92e0SMatthew Knepley #if !defined(PETSC_USE_DYNAMIC_LIBRARIES) 98f86e92e0SMatthew Knepley ierr = PetscInitializePackage(PETSC_NULL);CHKERRQ(ierr); 99f86e92e0SMatthew Knepley #endif 100f86e92e0SMatthew Knepley ierr = PetscHeaderCreate(o,_p_PetscObject,-1,cookie,0,name,comm,PetscObjectDestroy_PetscObject,0);CHKERRQ(ierr); 101f86e92e0SMatthew Knepley /* records not yet defined in PetscObject 102f86e92e0SMatthew Knepley o->data = 0; 103f86e92e0SMatthew Knepley o->setupcalled = 0; 104f86e92e0SMatthew Knepley */ 105f86e92e0SMatthew Knepley *obj = o; 106f86e92e0SMatthew Knepley PetscFunctionReturn(0); 107f86e92e0SMatthew Knepley } 108f86e92e0SMatthew Knepley 109f86e92e0SMatthew 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); 165*3050cee2SBarry Smith if (!viewer) { 166*3050cee2SBarry Smith ierr = PetscViewerASCIIGetStdout(obj->comm,&viewer);CHKERRQ(ierr); 167*3050cee2SBarry Smith } 168e5c89e4eSSatish Balay PetscValidHeaderSpecific(viewer,PETSC_VIEWER_COOKIE,2); 169e5c89e4eSSatish Balay 170e5c89e4eSSatish Balay if (obj->bops->view) { 171e5c89e4eSSatish Balay ierr = (*obj->bops->view)(obj,viewer);CHKERRQ(ierr); 172e5c89e4eSSatish Balay } else { 173e5c89e4eSSatish Balay SETERRQ(PETSC_ERR_SUP,"This PETSc object does not have a generic viewer routine"); 174e5c89e4eSSatish Balay } 175e5c89e4eSSatish Balay PetscFunctionReturn(0); 176e5c89e4eSSatish Balay } 177e5c89e4eSSatish Balay 178e5c89e4eSSatish Balay #undef __FUNCT__ 179e5c89e4eSSatish Balay #define __FUNCT__ "PetscTypeCompare" 180e5c89e4eSSatish Balay /*@C 181e5c89e4eSSatish Balay PetscTypeCompare - Determines whether a PETSc object is of a particular type. 182e5c89e4eSSatish Balay 183e5c89e4eSSatish Balay Not Collective 184e5c89e4eSSatish Balay 185e5c89e4eSSatish Balay Input Parameters: 186e5c89e4eSSatish Balay + obj - any PETSc object, for example a Vec, Mat or KSP. 187e5c89e4eSSatish Balay This must be cast with a (PetscObject), for example, 188e5c89e4eSSatish Balay PetscObjectDestroy((PetscObject)mat); 189e5c89e4eSSatish Balay - type_name - string containing a type name 190e5c89e4eSSatish Balay 191e5c89e4eSSatish Balay Output Parameter: 192e5c89e4eSSatish Balay . same - PETSC_TRUE if they are the same, else PETSC_FALSE 193e5c89e4eSSatish Balay 194e5c89e4eSSatish Balay Level: intermediate 195e5c89e4eSSatish Balay 196e5c89e4eSSatish Balay .seealso: VecGetType(), KSPGetType(), PCGetType(), SNESGetType() 197e5c89e4eSSatish Balay 198e5c89e4eSSatish Balay Concepts: comparing^object types 199e5c89e4eSSatish Balay Concepts: types^comparing 200e5c89e4eSSatish Balay Concepts: object type^comparing 201e5c89e4eSSatish Balay 202e5c89e4eSSatish Balay @*/ 203e5c89e4eSSatish Balay PetscErrorCode PETSC_DLLEXPORT PetscTypeCompare(PetscObject obj,const char type_name[],PetscTruth *same) 204e5c89e4eSSatish Balay { 205e5c89e4eSSatish Balay PetscErrorCode ierr; 206e5c89e4eSSatish Balay 207e5c89e4eSSatish Balay PetscFunctionBegin; 208e5c89e4eSSatish Balay if (!obj) { 209e5c89e4eSSatish Balay *same = PETSC_FALSE; 210e5c89e4eSSatish Balay } else if (!type_name && !obj->type_name) { 211e5c89e4eSSatish Balay *same = PETSC_TRUE; 212e5c89e4eSSatish Balay } else if (!type_name || !obj->type_name) { 213e5c89e4eSSatish Balay *same = PETSC_FALSE; 214e5c89e4eSSatish Balay } else { 215e5c89e4eSSatish Balay PetscValidHeader(obj,1); 216e5c89e4eSSatish Balay PetscValidCharPointer(type_name,2); 217e5c89e4eSSatish Balay PetscValidPointer(same,3); 218e5c89e4eSSatish Balay ierr = PetscStrcmp((char*)(obj->type_name),type_name,same);CHKERRQ(ierr); 219e5c89e4eSSatish Balay } 220e5c89e4eSSatish Balay PetscFunctionReturn(0); 221e5c89e4eSSatish Balay } 222e5c89e4eSSatish Balay 223e5c89e4eSSatish Balay static int PetscObjectRegisterDestroy_Count = 0; 224e5c89e4eSSatish Balay static PetscObject PetscObjectRegisterDestroy_Objects[256]; 225e5c89e4eSSatish Balay 226e5c89e4eSSatish Balay #undef __FUNCT__ 227e5c89e4eSSatish Balay #define __FUNCT__ "PetscObjectRegisterDestroy" 228e5c89e4eSSatish Balay /*@C 229e5c89e4eSSatish Balay PetscObjectRegisterDestroy - Registers a PETSc object to be destroyed when 230e5c89e4eSSatish Balay PetscFinalize() is called. 231e5c89e4eSSatish Balay 232e5c89e4eSSatish Balay Collective on PetscObject 233e5c89e4eSSatish Balay 234e5c89e4eSSatish Balay Input Parameter: 235e5c89e4eSSatish Balay . obj - any PETSc object, for example a Vec, Mat or KSP. 236e5c89e4eSSatish Balay This must be cast with a (PetscObject), for example, 237e5c89e4eSSatish Balay PetscObjectRegisterDestroy((PetscObject)mat); 238e5c89e4eSSatish Balay 239e5c89e4eSSatish Balay Level: developer 240e5c89e4eSSatish Balay 241e5c89e4eSSatish Balay Notes: 242e5c89e4eSSatish Balay This is used by, for example, PETSC_VIEWER_XXX_() routines to free the viewer 243e5c89e4eSSatish Balay when PETSc ends. 244e5c89e4eSSatish Balay 245e5c89e4eSSatish Balay .seealso: PetscObjectRegisterDestroyAll() 246e5c89e4eSSatish Balay @*/ 247e5c89e4eSSatish Balay PetscErrorCode PETSC_DLLEXPORT PetscObjectRegisterDestroy(PetscObject obj) 248e5c89e4eSSatish Balay { 249e5c89e4eSSatish Balay PetscFunctionBegin; 250e5c89e4eSSatish Balay PetscValidHeader(obj,1); 251e5c89e4eSSatish Balay PetscObjectRegisterDestroy_Objects[PetscObjectRegisterDestroy_Count++] = obj; 252e5c89e4eSSatish Balay PetscFunctionReturn(0); 253e5c89e4eSSatish Balay } 254e5c89e4eSSatish Balay 255e5c89e4eSSatish Balay #undef __FUNCT__ 256e5c89e4eSSatish Balay #define __FUNCT__ "PetscObjectRegisterDestroyAll" 257e5c89e4eSSatish Balay /*@C 258e5c89e4eSSatish Balay PetscObjectRegisterDestroyAll - Frees all the PETSc objects that have been registered 259e5c89e4eSSatish Balay with PetscObjectRegisterDestroy(). Called by PetscFinalize() 260e5c89e4eSSatish Balay PetscFinalize() is called. 261e5c89e4eSSatish Balay 262e5c89e4eSSatish Balay Collective on individual PetscObjects 263e5c89e4eSSatish Balay 264e5c89e4eSSatish Balay Level: developer 265e5c89e4eSSatish Balay 266e5c89e4eSSatish Balay .seealso: PetscObjectRegisterDestroy() 267e5c89e4eSSatish Balay @*/ 268e5c89e4eSSatish Balay PetscErrorCode PETSC_DLLEXPORT PetscObjectRegisterDestroyAll(void) 269e5c89e4eSSatish Balay { 270e5c89e4eSSatish Balay PetscErrorCode ierr; 271e5c89e4eSSatish Balay int i; 272e5c89e4eSSatish Balay 273e5c89e4eSSatish Balay PetscFunctionBegin; 274e5c89e4eSSatish Balay for (i=0; i<PetscObjectRegisterDestroy_Count; i++) { 275e5c89e4eSSatish Balay ierr = PetscObjectDestroy(PetscObjectRegisterDestroy_Objects[i]);CHKERRQ(ierr); 276e5c89e4eSSatish Balay } 277e5c89e4eSSatish Balay PetscFunctionReturn(0); 278e5c89e4eSSatish Balay } 279e5c89e4eSSatish Balay 280e5c89e4eSSatish Balay 281