1 2 #include <petsc/private/tsimpl.h> 3 4 static PetscBool TSPackageInitialized = PETSC_FALSE; 5 /*@C 6 TSFinalizePackage - This function destroys everything in the Petsc interface to Mathematica. It is 7 called from PetscFinalize(). 8 9 Level: developer 10 11 .keywords: Petsc, destroy, package, mathematica 12 .seealso: PetscFinalize() 13 @*/ 14 PetscErrorCode TSFinalizePackage(void) 15 { 16 PetscErrorCode ierr; 17 18 PetscFunctionBegin; 19 ierr = PetscFunctionListDestroy(&TSList);CHKERRQ(ierr); 20 ierr = PetscFunctionListDestroy(&TSTrajectoryList);CHKERRQ(ierr); 21 TSPackageInitialized = PETSC_FALSE; 22 TSRegisterAllCalled = PETSC_FALSE; 23 PetscFunctionReturn(0); 24 } 25 26 /*@C 27 TSInitializePackage - This function initializes everything in the TS package. It is called 28 from PetscDLLibraryRegister() when using dynamic libraries, and on the first call to TSCreate() 29 when using static libraries. 30 31 Level: developer 32 33 .keywords: TS, initialize, package 34 .seealso: PetscInitialize() 35 @*/ 36 PetscErrorCode TSInitializePackage(void) 37 { 38 char logList[256]; 39 PetscBool opt,pkg,cls; 40 PetscErrorCode ierr; 41 42 PetscFunctionBegin; 43 if (TSPackageInitialized) PetscFunctionReturn(0); 44 TSPackageInitialized = PETSC_TRUE; 45 /* Inialize subpackages */ 46 ierr = TSAdaptInitializePackage();CHKERRQ(ierr); 47 ierr = TSGLLEInitializePackage();CHKERRQ(ierr); 48 ierr = TSRKInitializePackage();CHKERRQ(ierr); 49 ierr = TSGLEEInitializePackage();CHKERRQ(ierr); 50 ierr = TSARKIMEXInitializePackage();CHKERRQ(ierr); 51 ierr = TSRosWInitializePackage();CHKERRQ(ierr); 52 ierr = TSSSPInitializePackage();CHKERRQ(ierr); 53 ierr = TSGLLEAdaptInitializePackage();CHKERRQ(ierr); 54 /* Register Classes */ 55 ierr = PetscClassIdRegister("TS",&TS_CLASSID);CHKERRQ(ierr); 56 ierr = PetscClassIdRegister("DMTS",&DMTS_CLASSID);CHKERRQ(ierr); 57 ierr = PetscClassIdRegister("TSTrajectory",&TSTRAJECTORY_CLASSID);CHKERRQ(ierr); 58 /* Register Constructors */ 59 ierr = TSRegisterAll();CHKERRQ(ierr); 60 ierr = TSTrajectoryRegisterAll();CHKERRQ(ierr); 61 /* Register Events */ 62 ierr = PetscLogEventRegister("TSStep", TS_CLASSID,&TS_Step);CHKERRQ(ierr); 63 ierr = PetscLogEventRegister("TSFunctionEval", TS_CLASSID,&TS_FunctionEval);CHKERRQ(ierr); 64 ierr = PetscLogEventRegister("TSJacobianEval", TS_CLASSID,&TS_JacobianEval);CHKERRQ(ierr); 65 ierr = PetscLogEventRegister("TSForwardStep", TS_CLASSID,&TS_ForwardStep);CHKERRQ(ierr); 66 ierr = PetscLogEventRegister("TSAdjointStep", TS_CLASSID,&TS_AdjointStep);CHKERRQ(ierr); 67 ierr = PetscLogEventRegister("TSTrajectorySet", TSTRAJECTORY_CLASSID,&TSTrajectory_Set);CHKERRQ(ierr); 68 ierr = PetscLogEventRegister("TSTrajectoryGet", TSTRAJECTORY_CLASSID,&TSTrajectory_Get);CHKERRQ(ierr); 69 ierr = PetscLogEventRegister("TSTrajDiskWrite", TSTRAJECTORY_CLASSID,&TSTrajectory_DiskWrite);CHKERRQ(ierr); 70 ierr = PetscLogEventRegister("TSTrajDiskRead", TSTRAJECTORY_CLASSID,&TSTrajectory_DiskRead);CHKERRQ(ierr); 71 ierr = PetscLogEventRegister("TSPseudoCmptTStp",TS_CLASSID,&TS_PseudoComputeTimeStep);CHKERRQ(ierr); 72 /* Process info exclusions */ 73 ierr = PetscOptionsGetString(NULL,NULL,"-info_exclude",logList,sizeof(logList),&opt);CHKERRQ(ierr); 74 if (opt) { 75 ierr = PetscStrInList("ts",logList,',',&pkg);CHKERRQ(ierr); 76 if (pkg) {ierr = PetscInfoDeactivateClass(TS_CLASSID);CHKERRQ(ierr);} 77 ierr = PetscStrInList("dm",logList,',',&cls);CHKERRQ(ierr); 78 if (pkg || cls) {ierr = PetscInfoDeactivateClass(DMTS_CLASSID);CHKERRQ(ierr);} 79 ierr = PetscStrInList("tsadapt",logList,',',&cls);CHKERRQ(ierr); 80 if (pkg || cls) {ierr = PetscInfoDeactivateClass(TSADAPT_CLASSID);CHKERRQ(ierr);} 81 ierr = PetscStrInList("tstrajectory",logList,',',&cls);CHKERRQ(ierr); 82 if (pkg || cls) {ierr = PetscInfoDeactivateClass(TSTRAJECTORY_CLASSID);CHKERRQ(ierr);} 83 } 84 /* Process summary exclusions */ 85 ierr = PetscOptionsGetString(NULL,NULL,"-log_exclude",logList,sizeof(logList),&opt);CHKERRQ(ierr); 86 if (opt) { 87 ierr = PetscStrInList("ts",logList,',',&pkg);CHKERRQ(ierr); 88 if (pkg) {ierr = PetscLogEventExcludeClass(TS_CLASSID);CHKERRQ(ierr);} 89 ierr = PetscStrInList("dm",logList,',',&cls);CHKERRQ(ierr); 90 if (pkg || cls) {ierr = PetscLogEventExcludeClass(DMTS_CLASSID);CHKERRQ(ierr);} 91 ierr = PetscStrInList("tsadapt",logList,',',&cls);CHKERRQ(ierr); 92 if (pkg || cls) {ierr = PetscLogEventExcludeClass(TSADAPT_CLASSID);CHKERRQ(ierr);} 93 ierr = PetscStrInList("tstrajectory",logList,',',&cls);CHKERRQ(ierr); 94 if (pkg || cls) {ierr = PetscLogEventExcludeClass(TSTRAJECTORY_CLASSID);CHKERRQ(ierr);} 95 } 96 /* Register package finalizer */ 97 ierr = PetscRegisterFinalize(TSFinalizePackage);CHKERRQ(ierr); 98 PetscFunctionReturn(0); 99 } 100 101 #if defined(PETSC_HAVE_DYNAMIC_LIBRARIES) 102 /* 103 PetscDLLibraryRegister - This function is called when the dynamic library it is in is opened. 104 105 This one registers all the TS methods that are in the basic PETSc libpetscts library. 106 107 */ 108 PETSC_EXTERN PetscErrorCode PetscDLLibraryRegister_petscts(void); /*prototype*/ 109 PETSC_EXTERN PetscErrorCode PetscDLLibraryRegister_petscts(void) 110 { 111 PetscErrorCode ierr; 112 113 PetscFunctionBegin; 114 ierr = TSInitializePackage();CHKERRQ(ierr); 115 PetscFunctionReturn(0); 116 } 117 118 119 #endif /* PETSC_HAVE_DYNAMIC_LIBRARIES */ 120