1 2 #include <petsc/private/matimpl.h> 3 4 const char *MatOptions_Shifted[] = {"UNUSED_NONZERO_LOCATION_ERR", 5 "ROW_ORIENTED", 6 "NOT_A_VALID_OPTION", 7 "SYMMETRIC", 8 "STRUCTURALLY_SYMMETRIC", 9 "FORCE_DIAGONAL_ENTRIES", 10 "IGNORE_OFF_PROC_ENTRIES", 11 "USE_HASH_TABLE", 12 "KEEP_NONZERO_PATTERN", 13 "IGNORE_ZERO_ENTRIES", 14 "USE_INODES", 15 "HERMITIAN", 16 "SYMMETRY_ETERNAL", 17 "NEW_NONZERO_LOCATION_ERR", 18 "IGNORE_LOWER_TRIANGULAR", 19 "ERROR_LOWER_TRIANGULAR", 20 "GETROW_UPPERTRIANGULAR", 21 "SPD", 22 "NO_OFF_PROC_ZERO_ROWS", 23 "NO_OFF_PROC_ENTRIES", 24 "NEW_NONZERO_LOCATIONS", 25 "NEW_NONZERO_ALLOCATION_ERR", 26 "SUBSET_OFF_PROC_ENTRIES", 27 "SUBMAT_SINGLEIS", 28 "STRUCTURE_ONLY", 29 "SORTED_FULL", 30 "FORM_EXPLICIT_TRANSPOSE", 31 "MatOption","MAT_",NULL}; 32 const char *const* MatOptions = MatOptions_Shifted+2; 33 const char *const MatFactorShiftTypes[] = {"NONE","NONZERO","POSITIVE_DEFINITE","INBLOCKS","MatFactorShiftType","PC_FACTOR_",NULL}; 34 const char *const MatStructures[] = {"DIFFERENT","SUBSET","SAME","UNKNOWN","MatStructure","MAT_STRUCTURE_",NULL}; 35 const char *const MatFactorShiftTypesDetail[] = {NULL,"diagonal shift to prevent zero pivot","Manteuffel shift","diagonal shift on blocks to prevent zero pivot"}; 36 const char *const MPPTScotchStrategyTypes[] = {"DEFAULT","QUALITY","SPEED","BALANCE","SAFETY","SCALABILITY","MPPTScotchStrategyType","MP_PTSCOTCH_",NULL}; 37 const char *const MPChacoGlobalTypes[] = {"","MULTILEVEL","SPECTRAL","","LINEAR","RANDOM","SCATTERED","MPChacoGlobalType","MP_CHACO_",NULL}; 38 const char *const MPChacoLocalTypes[] = {"","KERNIGHAN","NONE","MPChacoLocalType","MP_CHACO_",NULL}; 39 const char *const MPChacoEigenTypes[] = {"LANCZOS","RQI","MPChacoEigenType","MP_CHACO_",NULL}; 40 41 extern PetscErrorCode MatMFFDInitializePackage(void); 42 extern PetscErrorCode MatSolverTypeDestroy(void); 43 static PetscBool MatPackageInitialized = PETSC_FALSE; 44 /*@C 45 MatFinalizePackage - This function destroys everything in the Petsc interface to the Mat package. It is 46 called from PetscFinalize(). 47 48 Level: developer 49 50 .seealso: PetscFinalize(), MatInitializePackage() 51 @*/ 52 PetscErrorCode MatFinalizePackage(void) 53 { 54 MatRootName nnames,names = MatRootNameList; 55 PetscErrorCode ierr; 56 57 PetscFunctionBegin; 58 ierr = MatSolverTypeDestroy();CHKERRQ(ierr); 59 while (names) { 60 nnames = names->next; 61 ierr = PetscFree(names->rname);CHKERRQ(ierr); 62 ierr = PetscFree(names->sname);CHKERRQ(ierr); 63 ierr = PetscFree(names->mname);CHKERRQ(ierr); 64 ierr = PetscFree(names);CHKERRQ(ierr); 65 names = nnames; 66 } 67 ierr = PetscFunctionListDestroy(&MatList);CHKERRQ(ierr); 68 ierr = PetscFunctionListDestroy(&MatOrderingList);CHKERRQ(ierr); 69 ierr = PetscFunctionListDestroy(&MatColoringList);CHKERRQ(ierr); 70 ierr = PetscFunctionListDestroy(&MatPartitioningList);CHKERRQ(ierr); 71 ierr = PetscFunctionListDestroy(&MatCoarsenList);CHKERRQ(ierr); 72 MatRootNameList = NULL; 73 MatPackageInitialized = PETSC_FALSE; 74 MatRegisterAllCalled = PETSC_FALSE; 75 MatOrderingRegisterAllCalled = PETSC_FALSE; 76 MatColoringRegisterAllCalled = PETSC_FALSE; 77 MatPartitioningRegisterAllCalled = PETSC_FALSE; 78 MatCoarsenRegisterAllCalled = PETSC_FALSE; 79 /* this is not ideal because it exposes SeqAIJ implementation details directly into the base Mat code */ 80 ierr = PetscFunctionListDestroy(&MatSeqAIJList);CHKERRQ(ierr); 81 MatSeqAIJRegisterAllCalled = PETSC_FALSE; 82 PetscFunctionReturn(0); 83 } 84 85 #if defined(PETSC_HAVE_MUMPS) 86 PETSC_EXTERN PetscErrorCode MatSolverTypeRegister_MUMPS(void); 87 #endif 88 #if defined(PETSC_HAVE_CUDA) 89 PETSC_EXTERN PetscErrorCode MatSolverTypeRegister_CUSPARSE(void); 90 #endif 91 #if defined(PETSC_HAVE_KOKKOS_KERNELS) 92 PETSC_EXTERN PetscErrorCode MatSolverTypeRegister_KOKKOS(void); 93 #endif 94 #if defined(PETSC_HAVE_VIENNACL) 95 PETSC_EXTERN PetscErrorCode MatSolverTypeRegister_ViennaCL(void); 96 #endif 97 #if defined(PETSC_HAVE_ELEMENTAL) 98 PETSC_EXTERN PetscErrorCode MatSolverTypeRegister_Elemental(void); 99 #endif 100 #if defined(PETSC_HAVE_SCALAPACK) 101 PETSC_EXTERN PetscErrorCode MatSolverTypeRegister_ScaLAPACK(void); 102 #endif 103 #if defined(PETSC_HAVE_MATLAB_ENGINE) 104 PETSC_EXTERN PetscErrorCode MatSolverTypeRegister_Matlab(void); 105 #endif 106 #if defined(PETSC_HAVE_ESSL) 107 PETSC_EXTERN PetscErrorCode MatSolverTypeRegister_Essl(void); 108 #endif 109 #if defined(PETSC_HAVE_SUPERLU) 110 PETSC_EXTERN PetscErrorCode MatSolverTypeRegister_SuperLU(void); 111 #endif 112 #if defined(PETSC_HAVE_STRUMPACK) 113 PETSC_EXTERN PetscErrorCode MatSolverTypeRegister_STRUMPACK(void); 114 #endif 115 #if defined(PETSC_HAVE_PASTIX) 116 PETSC_EXTERN PetscErrorCode MatSolverTypeRegister_Pastix(void); 117 #endif 118 #if defined(PETSC_HAVE_SUPERLU_DIST) 119 PETSC_EXTERN PetscErrorCode MatSolverTypeRegister_SuperLU_DIST(void); 120 #endif 121 #if defined(PETSC_HAVE_ELEMENTAL) 122 PETSC_EXTERN PetscErrorCode MatSolverTypeRegister_SparseElemental(void); 123 #endif 124 #if defined(PETSC_HAVE_MKL_PARDISO) 125 PETSC_EXTERN PetscErrorCode MatSolverTypeRegister_MKL_Pardiso(void); 126 #endif 127 #if defined(PETSC_HAVE_MKL_CPARDISO) 128 PETSC_EXTERN PetscErrorCode MatSolverTypeRegister_MKL_CPardiso(void); 129 #endif 130 #if defined(PETSC_HAVE_SUITESPARSE) 131 PETSC_EXTERN PetscErrorCode MatSolverTypeRegister_SuiteSparse(void); 132 #endif 133 #if defined(PETSC_HAVE_LUSOL) 134 PETSC_EXTERN PetscErrorCode MatSolverTypeRegister_Lusol(void); 135 #endif 136 137 PETSC_INTERN PetscErrorCode MatGetFactor_seqaij_petsc(Mat,MatFactorType,Mat*); 138 PETSC_INTERN PetscErrorCode MatGetFactor_seqbaij_petsc(Mat,MatFactorType,Mat*); 139 PETSC_INTERN PetscErrorCode MatGetFactor_seqsbaij_petsc(Mat,MatFactorType,Mat*); 140 PETSC_INTERN PetscErrorCode MatGetFactor_seqdense_petsc(Mat,MatFactorType,Mat*); 141 #if defined(PETSC_HAVE_CUDA) 142 PETSC_INTERN PetscErrorCode MatGetFactor_seqdense_cuda(Mat,MatFactorType,Mat*); 143 #endif 144 PETSC_INTERN PetscErrorCode MatGetFactor_constantdiagonal_petsc(Mat,MatFactorType,Mat*); 145 PETSC_INTERN PetscErrorCode MatGetFactor_seqaij_bas(Mat,MatFactorType,Mat*); 146 147 /*@C 148 MatInitializePackage - This function initializes everything in the Mat package. It is called 149 from PetscDLLibraryRegister_petscmat() when using dynamic libraries, and on the first call to MatCreate() 150 when using shared or static libraries. 151 152 Level: developer 153 154 .seealso: PetscInitialize(), MatFinalizePackage() 155 @*/ 156 PetscErrorCode MatInitializePackage(void) 157 { 158 char logList[256]; 159 PetscBool opt,pkg; 160 161 PetscErrorCode ierr; 162 163 PetscFunctionBegin; 164 if (MatPackageInitialized) PetscFunctionReturn(0); 165 MatPackageInitialized = PETSC_TRUE; 166 /* Initialize subpackage */ 167 ierr = MatMFFDInitializePackage();CHKERRQ(ierr); 168 /* Register Classes */ 169 ierr = PetscClassIdRegister("Matrix",&MAT_CLASSID);CHKERRQ(ierr); 170 ierr = PetscClassIdRegister("Matrix FD Coloring",&MAT_FDCOLORING_CLASSID);CHKERRQ(ierr); 171 ierr = PetscClassIdRegister("Matrix Coloring",&MAT_COLORING_CLASSID);CHKERRQ(ierr); 172 ierr = PetscClassIdRegister("Matrix MatTranspose Coloring",&MAT_TRANSPOSECOLORING_CLASSID);CHKERRQ(ierr); 173 ierr = PetscClassIdRegister("Matrix Partitioning",&MAT_PARTITIONING_CLASSID);CHKERRQ(ierr); 174 ierr = PetscClassIdRegister("Matrix Coarsen",&MAT_COARSEN_CLASSID);CHKERRQ(ierr); 175 ierr = PetscClassIdRegister("Matrix Null Space",&MAT_NULLSPACE_CLASSID);CHKERRQ(ierr); 176 /* Register Constructors */ 177 ierr = MatRegisterAll();CHKERRQ(ierr); 178 ierr = MatOrderingRegisterAll();CHKERRQ(ierr); 179 ierr = MatColoringRegisterAll();CHKERRQ(ierr); 180 ierr = MatPartitioningRegisterAll();CHKERRQ(ierr); 181 ierr = MatCoarsenRegisterAll();CHKERRQ(ierr); 182 ierr = MatSeqAIJRegisterAll();CHKERRQ(ierr); 183 /* Register Events */ 184 ierr = PetscLogEventRegister("MatMult", MAT_CLASSID,&MAT_Mult);CHKERRQ(ierr); 185 ierr = PetscLogEventRegister("MatMults", MAT_CLASSID,&MAT_Mults);CHKERRQ(ierr); 186 ierr = PetscLogEventRegister("MatMultAdd", MAT_CLASSID,&MAT_MultAdd);CHKERRQ(ierr); 187 ierr = PetscLogEventRegister("MatMultTranspose", MAT_CLASSID,&MAT_MultTranspose);CHKERRQ(ierr); 188 ierr = PetscLogEventRegister("MatMultTrAdd", MAT_CLASSID,&MAT_MultTransposeAdd);CHKERRQ(ierr); 189 ierr = PetscLogEventRegister("MatSolve", MAT_CLASSID,&MAT_Solve);CHKERRQ(ierr); 190 ierr = PetscLogEventRegister("MatSolves", MAT_CLASSID,&MAT_Solves);CHKERRQ(ierr); 191 ierr = PetscLogEventRegister("MatSolveAdd", MAT_CLASSID,&MAT_SolveAdd);CHKERRQ(ierr); 192 ierr = PetscLogEventRegister("MatSolveTranspos", MAT_CLASSID,&MAT_SolveTranspose);CHKERRQ(ierr); 193 ierr = PetscLogEventRegister("MatSolveTrAdd", MAT_CLASSID,&MAT_SolveTransposeAdd);CHKERRQ(ierr); 194 ierr = PetscLogEventRegister("MatSOR", MAT_CLASSID,&MAT_SOR);CHKERRQ(ierr); 195 ierr = PetscLogEventRegister("MatForwardSolve", MAT_CLASSID,&MAT_ForwardSolve);CHKERRQ(ierr); 196 ierr = PetscLogEventRegister("MatBackwardSolve", MAT_CLASSID,&MAT_BackwardSolve);CHKERRQ(ierr); 197 ierr = PetscLogEventRegister("MatLUFactor", MAT_CLASSID,&MAT_LUFactor);CHKERRQ(ierr); 198 ierr = PetscLogEventRegister("MatLUFactorSym", MAT_CLASSID,&MAT_LUFactorSymbolic);CHKERRQ(ierr); 199 ierr = PetscLogEventRegister("MatLUFactorNum", MAT_CLASSID,&MAT_LUFactorNumeric);CHKERRQ(ierr); 200 ierr = PetscLogEventRegister("MatQRFactor", MAT_CLASSID,&MAT_QRFactor);CHKERRQ(ierr); 201 ierr = PetscLogEventRegister("MatQRFactorSym", MAT_CLASSID,&MAT_QRFactorSymbolic);CHKERRQ(ierr); 202 ierr = PetscLogEventRegister("MatQRFactorNum", MAT_CLASSID,&MAT_QRFactorNumeric);CHKERRQ(ierr); 203 ierr = PetscLogEventRegister("MatCholeskyFctr", MAT_CLASSID,&MAT_CholeskyFactor);CHKERRQ(ierr); 204 ierr = PetscLogEventRegister("MatCholFctrSym", MAT_CLASSID,&MAT_CholeskyFactorSymbolic);CHKERRQ(ierr); 205 ierr = PetscLogEventRegister("MatCholFctrNum", MAT_CLASSID,&MAT_CholeskyFactorNumeric);CHKERRQ(ierr); 206 ierr = PetscLogEventRegister("MatFctrFactSchur", MAT_CLASSID,&MAT_FactorFactS);CHKERRQ(ierr); 207 ierr = PetscLogEventRegister("MatFctrInvSchur", MAT_CLASSID,&MAT_FactorInvS);CHKERRQ(ierr); 208 ierr = PetscLogEventRegister("MatILUFactor", MAT_CLASSID,&MAT_ILUFactor);CHKERRQ(ierr); 209 ierr = PetscLogEventRegister("MatILUFactorSym", MAT_CLASSID,&MAT_ILUFactorSymbolic);CHKERRQ(ierr); 210 ierr = PetscLogEventRegister("MatICCFactorSym", MAT_CLASSID,&MAT_ICCFactorSymbolic);CHKERRQ(ierr); 211 ierr = PetscLogEventRegister("MatCopy", MAT_CLASSID,&MAT_Copy);CHKERRQ(ierr); 212 ierr = PetscLogEventRegister("MatConvert", MAT_CLASSID,&MAT_Convert);CHKERRQ(ierr); 213 ierr = PetscLogEventRegister("MatScale", MAT_CLASSID,&MAT_Scale);CHKERRQ(ierr); 214 ierr = PetscLogEventRegister("MatResidual", MAT_CLASSID,&MAT_Residual);CHKERRQ(ierr); 215 ierr = PetscLogEventRegister("MatAssemblyBegin", MAT_CLASSID,&MAT_AssemblyBegin);CHKERRQ(ierr); 216 ierr = PetscLogEventRegister("MatAssemblyEnd", MAT_CLASSID,&MAT_AssemblyEnd);CHKERRQ(ierr); 217 ierr = PetscLogEventRegister("MatSetValues", MAT_CLASSID,&MAT_SetValues);CHKERRQ(ierr); 218 ierr = PetscLogEventRegister("MatGetValues", MAT_CLASSID,&MAT_GetValues);CHKERRQ(ierr); 219 ierr = PetscLogEventRegister("MatGetRow", MAT_CLASSID,&MAT_GetRow);CHKERRQ(ierr); 220 ierr = PetscLogEventRegister("MatGetRowIJ", MAT_CLASSID,&MAT_GetRowIJ);CHKERRQ(ierr); 221 ierr = PetscLogEventRegister("MatCreateSubMats", MAT_CLASSID,&MAT_CreateSubMats);CHKERRQ(ierr); 222 ierr = PetscLogEventRegister("MatCreateSubMat", MAT_CLASSID,&MAT_CreateSubMat);CHKERRQ(ierr); 223 ierr = PetscLogEventRegister("MatGetOrdering", MAT_CLASSID,&MAT_GetOrdering);CHKERRQ(ierr); 224 ierr = PetscLogEventRegister("MatIncreaseOvrlp", MAT_CLASSID,&MAT_IncreaseOverlap);CHKERRQ(ierr); 225 ierr = PetscLogEventRegister("MatPartitioning", MAT_PARTITIONING_CLASSID,&MAT_Partitioning);CHKERRQ(ierr); 226 ierr = PetscLogEventRegister("MatPartitioningND",MAT_PARTITIONING_CLASSID,&MAT_PartitioningND);CHKERRQ(ierr); 227 ierr = PetscLogEventRegister("MatCoarsen", MAT_COARSEN_CLASSID,&MAT_Coarsen);CHKERRQ(ierr); 228 ierr = PetscLogEventRegister("MatZeroEntries", MAT_CLASSID,&MAT_ZeroEntries);CHKERRQ(ierr); 229 ierr = PetscLogEventRegister("MatLoad", MAT_CLASSID,&MAT_Load);CHKERRQ(ierr); 230 ierr = PetscLogEventRegister("MatView", MAT_CLASSID,&MAT_View);CHKERRQ(ierr); 231 ierr = PetscLogEventRegister("MatAXPY", MAT_CLASSID,&MAT_AXPY);CHKERRQ(ierr); 232 ierr = PetscLogEventRegister("MatFDColorCreate", MAT_FDCOLORING_CLASSID,&MAT_FDColoringCreate);CHKERRQ(ierr); 233 ierr = PetscLogEventRegister("MatFDColorSetUp", MAT_FDCOLORING_CLASSID,&MAT_FDColoringSetUp);CHKERRQ(ierr); 234 ierr = PetscLogEventRegister("MatFDColorApply", MAT_FDCOLORING_CLASSID,&MAT_FDColoringApply);CHKERRQ(ierr); 235 ierr = PetscLogEventRegister("MatFDColorFunc", MAT_FDCOLORING_CLASSID,&MAT_FDColoringFunction);CHKERRQ(ierr); 236 ierr = PetscLogEventRegister("MatTranspose", MAT_CLASSID,&MAT_Transpose);CHKERRQ(ierr); 237 ierr = PetscLogEventRegister("MatMatSolve", MAT_CLASSID,&MAT_MatSolve);CHKERRQ(ierr); 238 ierr = PetscLogEventRegister("MatMatTrSolve", MAT_CLASSID,&MAT_MatTrSolve);CHKERRQ(ierr); 239 ierr = PetscLogEventRegister("MatMatMultSym", MAT_CLASSID,&MAT_MatMultSymbolic);CHKERRQ(ierr); 240 ierr = PetscLogEventRegister("MatMatMultNum", MAT_CLASSID,&MAT_MatMultNumeric);CHKERRQ(ierr); 241 ierr = PetscLogEventRegister("MatMatMatMultSym", MAT_CLASSID,&MAT_MatMatMultSymbolic);CHKERRQ(ierr); 242 ierr = PetscLogEventRegister("MatMatMatMultNum", MAT_CLASSID,&MAT_MatMatMultNumeric);CHKERRQ(ierr); 243 ierr = PetscLogEventRegister("MatPtAPSymbolic", MAT_CLASSID,&MAT_PtAPSymbolic);CHKERRQ(ierr); 244 ierr = PetscLogEventRegister("MatPtAPNumeric", MAT_CLASSID,&MAT_PtAPNumeric);CHKERRQ(ierr); 245 ierr = PetscLogEventRegister("MatRARtSym", MAT_CLASSID,&MAT_RARtSymbolic);CHKERRQ(ierr); 246 ierr = PetscLogEventRegister("MatRARtNum", MAT_CLASSID,&MAT_RARtNumeric);CHKERRQ(ierr); 247 ierr = PetscLogEventRegister("MatMatTrnMultSym", MAT_CLASSID,&MAT_MatTransposeMultSymbolic);CHKERRQ(ierr); 248 ierr = PetscLogEventRegister("MatMatTrnMultNum", MAT_CLASSID,&MAT_MatTransposeMultNumeric);CHKERRQ(ierr); 249 ierr = PetscLogEventRegister("MatTrnMatMultSym", MAT_CLASSID,&MAT_TransposeMatMultSymbolic);CHKERRQ(ierr); 250 ierr = PetscLogEventRegister("MatTrnMatMultNum", MAT_CLASSID,&MAT_TransposeMatMultNumeric);CHKERRQ(ierr); 251 ierr = PetscLogEventRegister("MatTrnColorCreate",MAT_CLASSID,&MAT_TransposeColoringCreate);CHKERRQ(ierr); 252 ierr = PetscLogEventRegister("MatRedundantMat", MAT_CLASSID,&MAT_RedundantMat);CHKERRQ(ierr); 253 ierr = PetscLogEventRegister("MatGetSeqNZStrct", MAT_CLASSID,&MAT_GetSequentialNonzeroStructure);CHKERRQ(ierr); 254 ierr = PetscLogEventRegister("MatGetMultiProcB", MAT_CLASSID,&MAT_GetMultiProcBlock);CHKERRQ(ierr); 255 ierr = PetscLogEventRegister("MatSetRandom", MAT_CLASSID,&MAT_SetRandom);CHKERRQ(ierr); 256 257 /* these may be specific to MPIAIJ matrices */ 258 ierr = PetscLogEventRegister("MatMPISumSeqNumeric",MAT_CLASSID,&MAT_Seqstompinum);CHKERRQ(ierr); 259 ierr = PetscLogEventRegister("MatMPISumSeqSymbolic",MAT_CLASSID,&MAT_Seqstompisym);CHKERRQ(ierr); 260 ierr = PetscLogEventRegister("MatMPISumSeq",MAT_CLASSID,&MAT_Seqstompi);CHKERRQ(ierr); 261 ierr = PetscLogEventRegister("MatMPIConcateSeq",MAT_CLASSID,&MAT_Merge);CHKERRQ(ierr); 262 ierr = PetscLogEventRegister("MatGetLocalMat",MAT_CLASSID,&MAT_Getlocalmat);CHKERRQ(ierr); 263 ierr = PetscLogEventRegister("MatGetLocalMatCondensed",MAT_CLASSID,&MAT_Getlocalmatcondensed);CHKERRQ(ierr); 264 ierr = PetscLogEventRegister("MatGetBrowsOfAcols",MAT_CLASSID,&MAT_GetBrowsOfAcols);CHKERRQ(ierr); 265 ierr = PetscLogEventRegister("MatGetBrAoCol",MAT_CLASSID,&MAT_GetBrowsOfAocols);CHKERRQ(ierr); 266 267 ierr = PetscLogEventRegister("MatApplyPAPt_Symbolic",MAT_CLASSID,&MAT_Applypapt_symbolic);CHKERRQ(ierr); 268 ierr = PetscLogEventRegister("MatApplyPAPt_Numeric",MAT_CLASSID,&MAT_Applypapt_numeric);CHKERRQ(ierr); 269 ierr = PetscLogEventRegister("MatApplyPAPt",MAT_CLASSID,&MAT_Applypapt);CHKERRQ(ierr); 270 271 ierr = PetscLogEventRegister("MatGetSymTrans",MAT_CLASSID,&MAT_Getsymtranspose);CHKERRQ(ierr); 272 ierr = PetscLogEventRegister("MatGetSymTransR",MAT_CLASSID,&MAT_Getsymtransreduced);CHKERRQ(ierr); 273 ierr = PetscLogEventRegister("MatCUSPARSCopyTo",MAT_CLASSID,&MAT_CUSPARSECopyToGPU);CHKERRQ(ierr); 274 ierr = PetscLogEventRegister("MatCUSPARSCopyFr",MAT_CLASSID,&MAT_CUSPARSECopyFromGPU);CHKERRQ(ierr); 275 ierr = PetscLogEventRegister("MatCUSPARSSolAnl",MAT_CLASSID,&MAT_CUSPARSESolveAnalysis);CHKERRQ(ierr); 276 ierr = PetscLogEventRegister("MatCUSPARSGenT",MAT_CLASSID,&MAT_CUSPARSEGenerateTranspose);CHKERRQ(ierr); 277 ierr = PetscLogEventRegister("MatVCLCopyTo", MAT_CLASSID,&MAT_ViennaCLCopyToGPU);CHKERRQ(ierr); 278 ierr = PetscLogEventRegister("MatDenseCopyTo",MAT_CLASSID,&MAT_DenseCopyToGPU);CHKERRQ(ierr); 279 ierr = PetscLogEventRegister("MatDenseCopyFrom",MAT_CLASSID,&MAT_DenseCopyFromGPU);CHKERRQ(ierr); 280 ierr = PetscLogEventRegister("MatSetValBatch",MAT_CLASSID,&MAT_SetValuesBatch);CHKERRQ(ierr); 281 282 ierr = PetscLogEventRegister("MatColoringApply",MAT_COLORING_CLASSID,&MATCOLORING_Apply);CHKERRQ(ierr); 283 ierr = PetscLogEventRegister("MatColoringComm",MAT_COLORING_CLASSID,&MATCOLORING_Comm);CHKERRQ(ierr); 284 ierr = PetscLogEventRegister("MatColoringLocal",MAT_COLORING_CLASSID,&MATCOLORING_Local);CHKERRQ(ierr); 285 ierr = PetscLogEventRegister("MatColoringIS",MAT_COLORING_CLASSID,&MATCOLORING_ISCreate);CHKERRQ(ierr); 286 ierr = PetscLogEventRegister("MatColoringSetUp",MAT_COLORING_CLASSID,&MATCOLORING_SetUp);CHKERRQ(ierr); 287 ierr = PetscLogEventRegister("MatColoringWeights",MAT_COLORING_CLASSID,&MATCOLORING_Weights);CHKERRQ(ierr); 288 289 ierr = PetscLogEventRegister("MatSetPreallCOO",MAT_CLASSID,&MAT_PreallCOO);CHKERRQ(ierr); 290 ierr = PetscLogEventRegister("MatSetValuesCOO",MAT_CLASSID,&MAT_SetVCOO);CHKERRQ(ierr); 291 292 ierr = PetscLogEventRegister("MatH2OpusBuild",MAT_CLASSID,&MAT_H2Opus_Build);CHKERRQ(ierr); 293 ierr = PetscLogEventRegister("MatH2OpusComp", MAT_CLASSID,&MAT_H2Opus_Compress);CHKERRQ(ierr); 294 ierr = PetscLogEventRegister("MatH2OpusOrth", MAT_CLASSID,&MAT_H2Opus_Orthog);CHKERRQ(ierr); 295 ierr = PetscLogEventRegister("MatH2OpusLR", MAT_CLASSID,&MAT_H2Opus_LR);CHKERRQ(ierr); 296 297 /* Mark non-collective events */ 298 ierr = PetscLogEventSetCollective(MAT_SetValues, PETSC_FALSE);CHKERRQ(ierr); 299 ierr = PetscLogEventSetCollective(MAT_SetValuesBatch, PETSC_FALSE);CHKERRQ(ierr); 300 ierr = PetscLogEventSetCollective(MAT_GetRow, PETSC_FALSE);CHKERRQ(ierr); 301 /* Turn off high traffic events by default */ 302 ierr = PetscLogEventSetActiveAll(MAT_SetValues, PETSC_FALSE);CHKERRQ(ierr); 303 ierr = PetscLogEventSetActiveAll(MAT_GetValues, PETSC_FALSE);CHKERRQ(ierr); 304 ierr = PetscLogEventSetActiveAll(MAT_GetRow, PETSC_FALSE);CHKERRQ(ierr); 305 /* Process Info */ 306 { 307 PetscClassId classids[7]; 308 309 classids[0] = MAT_CLASSID; 310 classids[1] = MAT_FDCOLORING_CLASSID; 311 classids[2] = MAT_COLORING_CLASSID; 312 classids[3] = MAT_TRANSPOSECOLORING_CLASSID; 313 classids[4] = MAT_PARTITIONING_CLASSID; 314 classids[5] = MAT_COARSEN_CLASSID; 315 classids[6] = MAT_NULLSPACE_CLASSID; 316 ierr = PetscInfoProcessClass("mat", 7, classids);CHKERRQ(ierr); 317 } 318 319 /* Process summary exclusions */ 320 ierr = PetscOptionsGetString(NULL,NULL,"-log_exclude",logList,sizeof(logList),&opt);CHKERRQ(ierr); 321 if (opt) { 322 ierr = PetscStrInList("mat",logList,',',&pkg);CHKERRQ(ierr); 323 if (pkg) {ierr = PetscLogEventExcludeClass(MAT_CLASSID);CHKERRQ(ierr);} 324 if (pkg) {ierr = PetscLogEventExcludeClass(MAT_FDCOLORING_CLASSID);CHKERRQ(ierr);} 325 if (pkg) {ierr = PetscLogEventExcludeClass(MAT_COLORING_CLASSID);CHKERRQ(ierr);} 326 if (pkg) {ierr = PetscLogEventExcludeClass(MAT_TRANSPOSECOLORING_CLASSID);CHKERRQ(ierr);} 327 if (pkg) {ierr = PetscLogEventExcludeClass(MAT_PARTITIONING_CLASSID);CHKERRQ(ierr);} 328 if (pkg) {ierr = PetscLogEventExcludeClass(MAT_COARSEN_CLASSID);CHKERRQ(ierr);} 329 if (pkg) {ierr = PetscLogEventExcludeClass(MAT_NULLSPACE_CLASSID);CHKERRQ(ierr);} 330 } 331 332 /* Register the PETSc built in factorization based solvers */ 333 ierr = MatSolverTypeRegister(MATSOLVERPETSC, MATSEQAIJ, MAT_FACTOR_LU,MatGetFactor_seqaij_petsc);CHKERRQ(ierr); 334 ierr = MatSolverTypeRegister(MATSOLVERPETSC, MATSEQAIJ, MAT_FACTOR_CHOLESKY,MatGetFactor_seqaij_petsc);CHKERRQ(ierr); 335 ierr = MatSolverTypeRegister(MATSOLVERPETSC, MATSEQAIJ, MAT_FACTOR_ILU,MatGetFactor_seqaij_petsc);CHKERRQ(ierr); 336 ierr = MatSolverTypeRegister(MATSOLVERPETSC, MATSEQAIJ, MAT_FACTOR_ICC,MatGetFactor_seqaij_petsc);CHKERRQ(ierr); 337 338 ierr = MatSolverTypeRegister(MATSOLVERPETSC, MATSEQAIJPERM, MAT_FACTOR_LU,MatGetFactor_seqaij_petsc);CHKERRQ(ierr); 339 ierr = MatSolverTypeRegister(MATSOLVERPETSC, MATSEQAIJPERM, MAT_FACTOR_CHOLESKY,MatGetFactor_seqaij_petsc);CHKERRQ(ierr); 340 ierr = MatSolverTypeRegister(MATSOLVERPETSC, MATSEQAIJPERM, MAT_FACTOR_ILU,MatGetFactor_seqaij_petsc);CHKERRQ(ierr); 341 ierr = MatSolverTypeRegister(MATSOLVERPETSC, MATSEQAIJPERM, MAT_FACTOR_ICC,MatGetFactor_seqaij_petsc);CHKERRQ(ierr); 342 343 ierr = MatSolverTypeRegister(MATSOLVERPETSC, MATCONSTANTDIAGONAL,MAT_FACTOR_LU,MatGetFactor_constantdiagonal_petsc);CHKERRQ(ierr); 344 ierr = MatSolverTypeRegister(MATSOLVERPETSC, MATCONSTANTDIAGONAL,MAT_FACTOR_CHOLESKY,MatGetFactor_constantdiagonal_petsc);CHKERRQ(ierr); 345 ierr = MatSolverTypeRegister(MATSOLVERPETSC, MATCONSTANTDIAGONAL,MAT_FACTOR_ILU,MatGetFactor_constantdiagonal_petsc);CHKERRQ(ierr); 346 ierr = MatSolverTypeRegister(MATSOLVERPETSC, MATCONSTANTDIAGONAL,MAT_FACTOR_ICC,MatGetFactor_constantdiagonal_petsc);CHKERRQ(ierr); 347 348 #if defined(PETSC_HAVE_MKL_SPARSE) 349 ierr = MatSolverTypeRegister(MATSOLVERPETSC, MATSEQAIJMKL, MAT_FACTOR_LU,MatGetFactor_seqaij_petsc);CHKERRQ(ierr); 350 ierr = MatSolverTypeRegister(MATSOLVERPETSC, MATSEQAIJMKL, MAT_FACTOR_CHOLESKY,MatGetFactor_seqaij_petsc);CHKERRQ(ierr); 351 ierr = MatSolverTypeRegister(MATSOLVERPETSC, MATSEQAIJMKL, MAT_FACTOR_ILU,MatGetFactor_seqaij_petsc);CHKERRQ(ierr); 352 ierr = MatSolverTypeRegister(MATSOLVERPETSC, MATSEQAIJMKL, MAT_FACTOR_ICC,MatGetFactor_seqaij_petsc);CHKERRQ(ierr); 353 354 ierr = MatSolverTypeRegister(MATSOLVERPETSC, MATSEQBAIJMKL, MAT_FACTOR_LU,MatGetFactor_seqbaij_petsc);CHKERRQ(ierr); 355 ierr = MatSolverTypeRegister(MATSOLVERPETSC, MATSEQBAIJMKL, MAT_FACTOR_CHOLESKY,MatGetFactor_seqbaij_petsc);CHKERRQ(ierr); 356 ierr = MatSolverTypeRegister(MATSOLVERPETSC, MATSEQBAIJMKL, MAT_FACTOR_ILU,MatGetFactor_seqbaij_petsc);CHKERRQ(ierr); 357 ierr = MatSolverTypeRegister(MATSOLVERPETSC, MATSEQBAIJMKL, MAT_FACTOR_ICC,MatGetFactor_seqbaij_petsc);CHKERRQ(ierr); 358 #endif 359 /* Above, we register the PETSc built-in factorization solvers for MATSEQAIJMKL. In the future, we may want to use 360 * some of the MKL-provided ones instead. */ 361 362 ierr = MatSolverTypeRegister(MATSOLVERPETSC, MATSEQAIJCRL, MAT_FACTOR_LU,MatGetFactor_seqaij_petsc);CHKERRQ(ierr); 363 ierr = MatSolverTypeRegister(MATSOLVERPETSC, MATSEQAIJCRL, MAT_FACTOR_CHOLESKY,MatGetFactor_seqaij_petsc);CHKERRQ(ierr); 364 ierr = MatSolverTypeRegister(MATSOLVERPETSC, MATSEQAIJCRL, MAT_FACTOR_ILU,MatGetFactor_seqaij_petsc);CHKERRQ(ierr); 365 ierr = MatSolverTypeRegister(MATSOLVERPETSC, MATSEQAIJCRL, MAT_FACTOR_ICC,MatGetFactor_seqaij_petsc);CHKERRQ(ierr); 366 367 ierr = MatSolverTypeRegister(MATSOLVERPETSC, MATSEQBAIJ, MAT_FACTOR_LU,MatGetFactor_seqbaij_petsc);CHKERRQ(ierr); 368 ierr = MatSolverTypeRegister(MATSOLVERPETSC, MATSEQBAIJ, MAT_FACTOR_CHOLESKY,MatGetFactor_seqbaij_petsc);CHKERRQ(ierr); 369 ierr = MatSolverTypeRegister(MATSOLVERPETSC, MATSEQBAIJ, MAT_FACTOR_ILU,MatGetFactor_seqbaij_petsc);CHKERRQ(ierr); 370 ierr = MatSolverTypeRegister(MATSOLVERPETSC, MATSEQBAIJ, MAT_FACTOR_ICC,MatGetFactor_seqbaij_petsc);CHKERRQ(ierr); 371 372 ierr = MatSolverTypeRegister(MATSOLVERPETSC, MATSEQSBAIJ, MAT_FACTOR_CHOLESKY,MatGetFactor_seqsbaij_petsc);CHKERRQ(ierr); 373 ierr = MatSolverTypeRegister(MATSOLVERPETSC, MATSEQSBAIJ, MAT_FACTOR_ICC,MatGetFactor_seqsbaij_petsc);CHKERRQ(ierr); 374 375 ierr = MatSolverTypeRegister(MATSOLVERPETSC, MATSEQDENSE, MAT_FACTOR_LU,MatGetFactor_seqdense_petsc);CHKERRQ(ierr); 376 ierr = MatSolverTypeRegister(MATSOLVERPETSC, MATSEQDENSE, MAT_FACTOR_ILU,MatGetFactor_seqdense_petsc);CHKERRQ(ierr); 377 ierr = MatSolverTypeRegister(MATSOLVERPETSC, MATSEQDENSE, MAT_FACTOR_CHOLESKY,MatGetFactor_seqdense_petsc);CHKERRQ(ierr); 378 ierr = MatSolverTypeRegister(MATSOLVERPETSC, MATSEQDENSE, MAT_FACTOR_QR,MatGetFactor_seqdense_petsc);CHKERRQ(ierr); 379 #if defined(PETSC_HAVE_CUDA) 380 ierr = MatSolverTypeRegister(MATSOLVERCUDA, MATSEQDENSE, MAT_FACTOR_LU,MatGetFactor_seqdense_cuda);CHKERRQ(ierr); 381 ierr = MatSolverTypeRegister(MATSOLVERCUDA, MATSEQDENSE, MAT_FACTOR_CHOLESKY,MatGetFactor_seqdense_cuda);CHKERRQ(ierr); 382 ierr = MatSolverTypeRegister(MATSOLVERCUDA, MATSEQDENSE, MAT_FACTOR_QR,MatGetFactor_seqdense_cuda);CHKERRQ(ierr); 383 ierr = MatSolverTypeRegister(MATSOLVERCUDA, MATSEQDENSECUDA, MAT_FACTOR_LU,MatGetFactor_seqdense_cuda);CHKERRQ(ierr); 384 ierr = MatSolverTypeRegister(MATSOLVERCUDA, MATSEQDENSECUDA, MAT_FACTOR_CHOLESKY,MatGetFactor_seqdense_cuda);CHKERRQ(ierr); 385 ierr = MatSolverTypeRegister(MATSOLVERCUDA, MATSEQDENSECUDA, MAT_FACTOR_QR,MatGetFactor_seqdense_cuda);CHKERRQ(ierr); 386 #endif 387 388 ierr = MatSolverTypeRegister(MATSOLVERBAS, MATSEQAIJ, MAT_FACTOR_ICC,MatGetFactor_seqaij_bas);CHKERRQ(ierr); 389 390 /* 391 Register the external package factorization based solvers 392 Eventually we don't want to have these hardwired here at compile time of PETSc 393 */ 394 #if defined(PETSC_HAVE_MUMPS) 395 ierr = MatSolverTypeRegister_MUMPS();CHKERRQ(ierr); 396 #endif 397 #if defined(PETSC_HAVE_CUDA) 398 ierr = MatSolverTypeRegister_CUSPARSE();CHKERRQ(ierr); 399 #endif 400 #if defined(PETSC_HAVE_KOKKOS_KERNELS) 401 ierr = MatSolverTypeRegister_KOKKOS();CHKERRQ(ierr); 402 #endif 403 #if defined(PETSC_HAVE_VIENNACL) 404 ierr = MatSolverTypeRegister_ViennaCL();CHKERRQ(ierr); 405 #endif 406 #if defined(PETSC_HAVE_ELEMENTAL) 407 ierr = MatSolverTypeRegister_Elemental();CHKERRQ(ierr); 408 #endif 409 #if defined(PETSC_HAVE_SCALAPACK) 410 ierr = MatSolverTypeRegister_ScaLAPACK();CHKERRQ(ierr); 411 #endif 412 #if defined(PETSC_HAVE_MATLAB_ENGINE) 413 ierr = MatSolverTypeRegister_Matlab();CHKERRQ(ierr); 414 #endif 415 #if defined(PETSC_HAVE_ESSL) 416 ierr = MatSolverTypeRegister_Essl();CHKERRQ(ierr); 417 #endif 418 #if defined(PETSC_HAVE_SUPERLU) 419 ierr = MatSolverTypeRegister_SuperLU();CHKERRQ(ierr); 420 #endif 421 #if defined(PETSC_HAVE_STRUMPACK) 422 ierr = MatSolverTypeRegister_STRUMPACK();CHKERRQ(ierr); 423 #endif 424 #if defined(PETSC_HAVE_PASTIX) 425 ierr = MatSolverTypeRegister_Pastix();CHKERRQ(ierr); 426 #endif 427 #if defined(PETSC_HAVE_SUPERLU_DIST) 428 ierr = MatSolverTypeRegister_SuperLU_DIST();CHKERRQ(ierr); 429 #endif 430 #if defined(PETSC_HAVE_ELEMENTAL) 431 ierr = MatSolverTypeRegister_SparseElemental();CHKERRQ(ierr); 432 #endif 433 #if defined(PETSC_HAVE_MKL_PARDISO) 434 ierr = MatSolverTypeRegister_MKL_Pardiso();CHKERRQ(ierr); 435 #endif 436 #if defined(PETSC_HAVE_MKL_CPARDISO) 437 ierr = MatSolverTypeRegister_MKL_CPardiso();CHKERRQ(ierr); 438 #endif 439 #if defined(PETSC_HAVE_SUITESPARSE) 440 ierr = MatSolverTypeRegister_SuiteSparse();CHKERRQ(ierr); 441 #endif 442 #if defined(PETSC_HAVE_LUSOL) 443 ierr = MatSolverTypeRegister_Lusol();CHKERRQ(ierr); 444 #endif 445 /* Register package finalizer */ 446 ierr = PetscRegisterFinalize(MatFinalizePackage);CHKERRQ(ierr); 447 PetscFunctionReturn(0); 448 } 449 450 #if defined(PETSC_HAVE_DYNAMIC_LIBRARIES) 451 /* 452 PetscDLLibraryRegister - This function is called when the dynamic library it is in is opened. 453 454 This one registers all the matrix methods that are in the basic PETSc Matrix library. 455 456 */ 457 PETSC_EXTERN PetscErrorCode PetscDLLibraryRegister_petscmat(void) 458 { 459 PetscErrorCode ierr; 460 461 PetscFunctionBegin; 462 ierr = MatInitializePackage();CHKERRQ(ierr); 463 PetscFunctionReturn(0); 464 } 465 466 #endif /* PETSC_HAVE_DYNAMIC_LIBRARIES */ 467