1 #define PETSCMAT_DLL 2 3 #include "private/matimpl.h" 4 5 const char *MatOptions[] = {"ROW_ORIENTED","NEW_NONZERO_LOCATIONS", 6 "SYMMETRIC", 7 "STRUCTURALLY_SYMMETRIC", 8 "NEW_DIAGONALS", 9 "IGNORE_OFF_PROC_ENTRIES", 10 "NEW_NONZERO_LOCATION_ERR", 11 "NEW_NONZERO_ALLOCATION_ERR","USE_HASH_TABLE", 12 "KEEP_ZEROED_ROWS","IGNORE_ZERO_ENTRIES","USE_INODES", 13 "HERMITIAN", 14 "SYMMETRY_ETERNAL", 15 "CHECK_COMPRESSED_ROW", 16 "IGNORE_LOWER_TRIANGULAR","ERROR_LOWER_TRIANGULAR","GETROW_UPPERTRIANGULAR","SPD","NO_OFF_PROC_ENTRIES","NO_OFF_PROC_ZERO_ROWS","MatOption","MAT_",0}; 17 const char *MatFactorShiftTypes[] = {"NONE","NONZERO","POSITIVE_DEFINITE","INBLOCKS","MatFactorShiftType","PC_FACTOR_",0}; 18 19 extern PetscErrorCode MatMFFDInitializePackage(const char[]); 20 static PetscBool MatPackageInitialized = PETSC_FALSE; 21 #undef __FUNCT__ 22 #define __FUNCT__ "MatFinalizePackage" 23 /*@C 24 MatFinalizePackage - This function destroys everything in the Petsc interface to the Mat package. It is 25 called from PetscFinalize(). 26 27 Level: developer 28 29 .keywords: Petsc, destroy, package, mathematica 30 .seealso: PetscFinalize() 31 @*/ 32 PetscErrorCode MatFinalizePackage(void) 33 { 34 PetscFunctionBegin; 35 MatPackageInitialized = PETSC_FALSE; 36 MatRegisterAllCalled = PETSC_FALSE; 37 MatList = PETSC_NULL; 38 MatOrderingRegisterAllCalled = PETSC_FALSE; 39 MatOrderingList = PETSC_NULL; 40 MatColoringList = PETSC_NULL; 41 MatColoringRegisterAllCalled = PETSC_FALSE; 42 MatPartitioningList = PETSC_NULL; 43 MatPartitioningRegisterAllCalled = PETSC_FALSE; 44 PetscFunctionReturn(0); 45 } 46 47 #undef __FUNCT__ 48 #define __FUNCT__ "MatInitializePackage" 49 /*@C 50 MatInitializePackage - This function initializes everything in the Mat package. It is called 51 from PetscDLLibraryRegister() when using dynamic libraries, and on the first call to MatCreate() 52 when using static libraries. 53 54 Input Parameter: 55 path - The dynamic library path, or PETSC_NULL 56 57 Level: developer 58 59 .keywords: Mat, initialize, package 60 .seealso: PetscInitialize() 61 @*/ 62 PetscErrorCode MatInitializePackage(const char path[]) 63 { 64 char logList[256]; 65 char *className; 66 PetscBool opt; 67 PetscErrorCode ierr; 68 69 PetscFunctionBegin; 70 if (MatPackageInitialized) PetscFunctionReturn(0); 71 MatPackageInitialized = PETSC_TRUE; 72 /* Inialize subpackage */ 73 ierr = MatMFFDInitializePackage(path);CHKERRQ(ierr); 74 /* Register Classes */ 75 ierr = PetscClassIdRegister("Matrix",&MAT_CLASSID);CHKERRQ(ierr); 76 ierr = PetscClassIdRegister("Matrix FD Coloring",&MAT_FDCOLORING_CLASSID);CHKERRQ(ierr); 77 ierr = PetscClassIdRegister("Matrix Partitioning",&MAT_PARTITIONING_CLASSID);CHKERRQ(ierr); 78 ierr = PetscClassIdRegister("Matrix Null Space",&MAT_NULLSPACE_CLASSID);CHKERRQ(ierr); 79 /* Register Constructors */ 80 ierr = MatRegisterAll(path);CHKERRQ(ierr); 81 ierr = MatOrderingRegisterAll(path);CHKERRQ(ierr); 82 ierr = MatColoringRegisterAll(path);CHKERRQ(ierr); 83 ierr = MatPartitioningRegisterAll(path);CHKERRQ(ierr); 84 /* Register Events */ 85 ierr = PetscLogEventRegister("MatMult", MAT_CLASSID,&MAT_Mult);CHKERRQ(ierr); 86 ierr = PetscLogEventRegister("MatMults", MAT_CLASSID,&MAT_Mults);CHKERRQ(ierr); 87 ierr = PetscLogEventRegister("MatMultConstr", MAT_CLASSID,&MAT_MultConstrained);CHKERRQ(ierr); 88 ierr = PetscLogEventRegister("MatMultAdd", MAT_CLASSID,&MAT_MultAdd);CHKERRQ(ierr); 89 ierr = PetscLogEventRegister("MatMultTranspose", MAT_CLASSID,&MAT_MultTranspose);CHKERRQ(ierr); 90 ierr = PetscLogEventRegister("MatMultTrConstr", MAT_CLASSID,&MAT_MultTransposeConstrained);CHKERRQ(ierr); 91 ierr = PetscLogEventRegister("MatMultTrAdd", MAT_CLASSID,&MAT_MultTransposeAdd);CHKERRQ(ierr); 92 ierr = PetscLogEventRegister("MatSolve", MAT_CLASSID,&MAT_Solve);CHKERRQ(ierr); 93 ierr = PetscLogEventRegister("MatSolves", MAT_CLASSID,&MAT_Solves);CHKERRQ(ierr); 94 ierr = PetscLogEventRegister("MatSolveAdd", MAT_CLASSID,&MAT_SolveAdd);CHKERRQ(ierr); 95 ierr = PetscLogEventRegister("MatSolveTranspos", MAT_CLASSID,&MAT_SolveTranspose);CHKERRQ(ierr); 96 ierr = PetscLogEventRegister("MatSolveTrAdd", MAT_CLASSID,&MAT_SolveTransposeAdd);CHKERRQ(ierr); 97 ierr = PetscLogEventRegister("MatSOR", MAT_CLASSID,&MAT_SOR);CHKERRQ(ierr); 98 ierr = PetscLogEventRegister("MatForwardSolve", MAT_CLASSID,&MAT_ForwardSolve);CHKERRQ(ierr); 99 ierr = PetscLogEventRegister("MatBackwardSolve", MAT_CLASSID,&MAT_BackwardSolve);CHKERRQ(ierr); 100 ierr = PetscLogEventRegister("MatLUFactor", MAT_CLASSID,&MAT_LUFactor);CHKERRQ(ierr); 101 ierr = PetscLogEventRegister("MatLUFactorSym", MAT_CLASSID,&MAT_LUFactorSymbolic);CHKERRQ(ierr); 102 ierr = PetscLogEventRegister("MatLUFactorNum", MAT_CLASSID,&MAT_LUFactorNumeric);CHKERRQ(ierr); 103 ierr = PetscLogEventRegister("MatCholeskyFctr", MAT_CLASSID,&MAT_CholeskyFactor);CHKERRQ(ierr); 104 ierr = PetscLogEventRegister("MatCholFctrSym", MAT_CLASSID,&MAT_CholeskyFactorSymbolic);CHKERRQ(ierr); 105 ierr = PetscLogEventRegister("MatCholFctrNum", MAT_CLASSID,&MAT_CholeskyFactorNumeric);CHKERRQ(ierr); 106 ierr = PetscLogEventRegister("MatILUFactor", MAT_CLASSID,&MAT_ILUFactor);CHKERRQ(ierr); 107 ierr = PetscLogEventRegister("MatILUFactorSym", MAT_CLASSID,&MAT_ILUFactorSymbolic);CHKERRQ(ierr); 108 ierr = PetscLogEventRegister("MatICCFactorSym", MAT_CLASSID,&MAT_ICCFactorSymbolic);CHKERRQ(ierr); 109 ierr = PetscLogEventRegister("MatCopy", MAT_CLASSID,&MAT_Copy);CHKERRQ(ierr); 110 ierr = PetscLogEventRegister("MatConvert", MAT_CLASSID,&MAT_Convert);CHKERRQ(ierr); 111 ierr = PetscLogEventRegister("MatScale", MAT_CLASSID,&MAT_Scale);CHKERRQ(ierr); 112 ierr = PetscLogEventRegister("MatAssemblyBegin", MAT_CLASSID,&MAT_AssemblyBegin);CHKERRQ(ierr); 113 ierr = PetscLogEventRegister("MatAssemblyEnd", MAT_CLASSID,&MAT_AssemblyEnd);CHKERRQ(ierr); 114 ierr = PetscLogEventRegister("MatSetValues", MAT_CLASSID,&MAT_SetValues);CHKERRQ(ierr); 115 ierr = PetscLogEventRegister("MatGetValues", MAT_CLASSID,&MAT_GetValues);CHKERRQ(ierr); 116 ierr = PetscLogEventRegister("MatGetRow", MAT_CLASSID,&MAT_GetRow);CHKERRQ(ierr); 117 ierr = PetscLogEventRegister("MatGetRowIJ", MAT_CLASSID,&MAT_GetRowIJ);CHKERRQ(ierr); 118 ierr = PetscLogEventRegister("MatGetSubMatrice", MAT_CLASSID,&MAT_GetSubMatrices);CHKERRQ(ierr); 119 ierr = PetscLogEventRegister("MatGetColoring", MAT_CLASSID,&MAT_GetColoring);CHKERRQ(ierr); 120 ierr = PetscLogEventRegister("MatGetOrdering", MAT_CLASSID,&MAT_GetOrdering);CHKERRQ(ierr); 121 ierr = PetscLogEventRegister("MatIncreaseOvrlp", MAT_CLASSID,&MAT_IncreaseOverlap);CHKERRQ(ierr); 122 ierr = PetscLogEventRegister("MatPartitioning", MAT_PARTITIONING_CLASSID,&MAT_Partitioning);CHKERRQ(ierr); 123 ierr = PetscLogEventRegister("MatZeroEntries", MAT_CLASSID,&MAT_ZeroEntries);CHKERRQ(ierr); 124 ierr = PetscLogEventRegister("MatLoad", MAT_CLASSID,&MAT_Load);CHKERRQ(ierr); 125 ierr = PetscLogEventRegister("MatView", MAT_CLASSID,&MAT_View);CHKERRQ(ierr); 126 ierr = PetscLogEventRegister("MatAXPY", MAT_CLASSID,&MAT_AXPY);CHKERRQ(ierr); 127 ierr = PetscLogEventRegister("MatFDColorCreate", MAT_FDCOLORING_CLASSID,&MAT_FDColoringCreate);CHKERRQ(ierr); 128 ierr = PetscLogEventRegister("MatFDColorApply", MAT_FDCOLORING_CLASSID,&MAT_FDColoringApply);CHKERRQ(ierr); 129 ierr = PetscLogEventRegister("MatFDColorFunc", MAT_FDCOLORING_CLASSID,&MAT_FDColoringFunction);CHKERRQ(ierr); 130 ierr = PetscLogEventRegister("MatTranspose", MAT_CLASSID,&MAT_Transpose);CHKERRQ(ierr); 131 ierr = PetscLogEventRegister("MatMatMult", MAT_CLASSID,&MAT_MatMult);CHKERRQ(ierr); 132 ierr = PetscLogEventRegister("MatMatSolve", MAT_CLASSID,&MAT_MatSolve);CHKERRQ(ierr); 133 ierr = PetscLogEventRegister("MatMatMultSym", MAT_CLASSID,&MAT_MatMultSymbolic);CHKERRQ(ierr); 134 ierr = PetscLogEventRegister("MatMatMultNum", MAT_CLASSID,&MAT_MatMultNumeric);CHKERRQ(ierr); 135 ierr = PetscLogEventRegister("MatPtAP", MAT_CLASSID,&MAT_PtAP);CHKERRQ(ierr); 136 ierr = PetscLogEventRegister("MatPtAPSymbolic", MAT_CLASSID,&MAT_PtAPSymbolic);CHKERRQ(ierr); 137 ierr = PetscLogEventRegister("MatPtAPNumeric", MAT_CLASSID,&MAT_PtAPNumeric);CHKERRQ(ierr); 138 ierr = PetscLogEventRegister("MatMatMultTrans", MAT_CLASSID,&MAT_MatMultTranspose);CHKERRQ(ierr); 139 ierr = PetscLogEventRegister("MatMatMultTrnSym" ,MAT_CLASSID,&MAT_MatMultTransposeSymbolic);CHKERRQ(ierr); 140 ierr = PetscLogEventRegister("MatMatMultTrnNum", MAT_CLASSID,&MAT_MatMultTransposeNumeric);CHKERRQ(ierr); 141 ierr = PetscLogEventRegister("MatGetRedundant", MAT_CLASSID,&MAT_GetRedundantMatrix);CHKERRQ(ierr); 142 ierr = PetscLogEventRegister("MatGetSeqNZStrct", MAT_CLASSID,&MAT_GetSequentialNonzeroStructure);CHKERRQ(ierr); 143 ierr = PetscLogEventRegister("MatGetMultiProcBlock", MAT_CLASSID,&MAT_GetMultiProcBlock);CHKERRQ(ierr); 144 145 /* these may be specific to MPIAIJ matrices */ 146 ierr = PetscLogEventRegister("MatMerge_SeqsToMPINumeric",MAT_CLASSID,&MAT_Seqstompinum); 147 ierr = PetscLogEventRegister("MatMerge_SeqsToMPISymbolic",MAT_CLASSID,&MAT_Seqstompisym); 148 ierr = PetscLogEventRegister("MatMerge_SeqsToMPI",MAT_CLASSID,&MAT_Seqstompi); 149 ierr = PetscLogEventRegister("MatGetLocalMat",MAT_CLASSID,&MAT_Getlocalmat); 150 ierr = PetscLogEventRegister("MatGetLocalMatCondensed",MAT_CLASSID,&MAT_Getlocalmatcondensed); 151 ierr = PetscLogEventRegister("MatGetBrowsOfAcols",MAT_CLASSID,&MAT_GetBrowsOfAcols); 152 ierr = PetscLogEventRegister("MatGetBrAoCol",MAT_CLASSID,&MAT_GetBrowsOfAocols); 153 154 ierr = PetscLogEventRegister("MatApplyPAPt_Symbolic",MAT_CLASSID,&MAT_Applypapt_symbolic);CHKERRQ(ierr); 155 ierr = PetscLogEventRegister("MatApplyPAPt_Numeric",MAT_CLASSID,&MAT_Applypapt_numeric);CHKERRQ(ierr); 156 ierr = PetscLogEventRegister("MatApplyPAPt",MAT_CLASSID,&MAT_Applypapt);CHKERRQ(ierr); 157 158 ierr = PetscLogEventRegister("MatGetSymTrans",MAT_CLASSID,&MAT_Getsymtranspose);CHKERRQ(ierr); 159 ierr = PetscLogEventRegister("MatGetSymTransR",MAT_CLASSID,&MAT_Getsymtransreduced);CHKERRQ(ierr); 160 ierr = PetscLogEventRegister("MatTranspose_SeqAIJ_FAST",MAT_CLASSID,&MAT_Transpose_SeqAIJ);CHKERRQ(ierr); 161 ierr = PetscLogEventRegister("MatCUDACopyTo",MAT_CLASSID,&MAT_CUDACopyToGPU);CHKERRQ(ierr); 162 163 /* Turn off high traffic events by default */ 164 ierr = PetscLogEventSetActiveAll(MAT_SetValues, PETSC_FALSE);CHKERRQ(ierr); 165 /* Process info exclusions */ 166 ierr = PetscOptionsGetString(PETSC_NULL, "-info_exclude", logList, 256, &opt);CHKERRQ(ierr); 167 if (opt) { 168 ierr = PetscStrstr(logList, "mat", &className);CHKERRQ(ierr); 169 if (className) { 170 ierr = PetscInfoDeactivateClass(MAT_CLASSID);CHKERRQ(ierr); 171 } 172 } 173 /* Process summary exclusions */ 174 ierr = PetscOptionsGetString(PETSC_NULL, "-log_summary_exclude", logList, 256, &opt);CHKERRQ(ierr); 175 if (opt) { 176 ierr = PetscStrstr(logList, "mat", &className);CHKERRQ(ierr); 177 if (className) { 178 ierr = PetscLogEventDeactivateClass(MAT_CLASSID);CHKERRQ(ierr); 179 } 180 } 181 ierr = PetscRegisterFinalize(MatFinalizePackage);CHKERRQ(ierr); 182 PetscFunctionReturn(0); 183 } 184 185 #ifdef PETSC_USE_DYNAMIC_LIBRARIES 186 EXTERN_C_BEGIN 187 #undef __FUNCT__ 188 #define __FUNCT__ "PetscDLLibraryRegister_petscmat" 189 /* 190 PetscDLLibraryRegister - This function is called when the dynamic library it is in is opened. 191 192 This one registers all the matrix methods that are in the basic PETSc Matrix library. 193 194 Input Parameter: 195 path - library path 196 */ 197 PetscErrorCode PetscDLLibraryRegister_petscmat(const char path[]) 198 { 199 PetscErrorCode ierr; 200 201 PetscFunctionBegin; 202 ierr = MatInitializePackage(path);CHKERRQ(ierr); 203 PetscFunctionReturn(0); 204 } 205 EXTERN_C_END 206 207 208 #endif /* PETSC_USE_DYNAMIC_LIBRARIES */ 209