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