1 2 #include <petsc/private/matimpl.h> /*I "petscmat.h" I*/ 3 #include <petscpkg_version.h> 4 5 PETSC_EXTERN PetscErrorCode MatCreate_MFFD(Mat); 6 PETSC_EXTERN PetscErrorCode MatCreate_MAIJ(Mat); 7 PETSC_EXTERN PetscErrorCode MatCreate_KAIJ(Mat); 8 PETSC_EXTERN PetscErrorCode MatCreate_IS(Mat); 9 10 PETSC_EXTERN PetscErrorCode MatCreate_SeqAIJ(Mat); 11 PETSC_EXTERN PetscErrorCode MatCreate_MPIAIJ(Mat); 12 13 PETSC_EXTERN PetscErrorCode MatCreate_SeqBAIJ(Mat); 14 PETSC_EXTERN PetscErrorCode MatCreate_MPIBAIJ(Mat); 15 16 PETSC_EXTERN PetscErrorCode MatCreate_SeqSBAIJ(Mat); 17 PETSC_EXTERN PetscErrorCode MatCreate_MPISBAIJ(Mat); 18 19 PETSC_EXTERN PetscErrorCode MatCreate_SeqDense(Mat); 20 PETSC_EXTERN PetscErrorCode MatCreate_MPIDense(Mat); 21 #if defined(PETSC_HAVE_CUDA) 22 PETSC_EXTERN PetscErrorCode MatCreate_SeqDenseCUDA(Mat); 23 PETSC_EXTERN PetscErrorCode MatCreate_MPIDenseCUDA(Mat); 24 #endif 25 26 PETSC_EXTERN PetscErrorCode MatCreate_MPIAdj(Mat); 27 PETSC_EXTERN PetscErrorCode MatCreate_Shell(Mat); 28 PETSC_EXTERN PetscErrorCode MatCreate_Composite(Mat); 29 30 PETSC_EXTERN PetscErrorCode MatCreate_SeqAIJPERM(Mat); 31 PETSC_EXTERN PetscErrorCode MatCreate_MPIAIJPERM(Mat); 32 33 PETSC_EXTERN PetscErrorCode MatCreate_SeqAIJSELL(Mat); 34 PETSC_EXTERN PetscErrorCode MatCreate_MPIAIJSELL(Mat); 35 36 #if defined(PETSC_HAVE_MKL_SPARSE) 37 PETSC_EXTERN PetscErrorCode MatCreate_SeqAIJMKL(Mat); 38 PETSC_EXTERN PetscErrorCode MatCreate_MPIAIJMKL(Mat); 39 #endif 40 41 #if defined(PETSC_HAVE_MKL_SPARSE_OPTIMIZE) 42 PETSC_EXTERN PetscErrorCode MatCreate_SeqBAIJMKL(Mat); 43 PETSC_EXTERN PetscErrorCode MatCreate_MPIBAIJMKL(Mat); 44 #endif 45 46 PETSC_EXTERN PetscErrorCode MatCreate_SeqAIJCRL(Mat); 47 PETSC_EXTERN PetscErrorCode MatCreate_MPIAIJCRL(Mat); 48 49 PETSC_EXTERN PetscErrorCode MatCreate_Scatter(Mat); 50 PETSC_EXTERN PetscErrorCode MatCreate_BlockMat(Mat); 51 PETSC_EXTERN PetscErrorCode MatCreate_Nest(Mat); 52 53 PETSC_EXTERN PetscErrorCode MatCreate_SeqSELL(Mat); 54 PETSC_EXTERN PetscErrorCode MatCreate_MPISELL(Mat); 55 56 #if defined(PETSC_HAVE_CUDA) 57 PETSC_EXTERN PetscErrorCode MatCreate_SeqAIJCUSPARSE(Mat); 58 PETSC_EXTERN PetscErrorCode MatCreate_MPIAIJCUSPARSE(Mat); 59 #endif 60 61 #if defined(PETSC_HAVE_VIENNACL) 62 PETSC_EXTERN PetscErrorCode MatCreate_SeqAIJViennaCL(Mat); 63 PETSC_EXTERN PetscErrorCode MatCreate_MPIAIJViennaCL(Mat); 64 #endif 65 66 #if defined(PETSC_HAVE_FFTW) 67 PETSC_EXTERN PetscErrorCode MatCreate_FFTW(Mat); 68 #endif 69 #if defined(PETSC_HAVE_ELEMENTAL) 70 PETSC_EXTERN PetscErrorCode MatCreate_Elemental(Mat); 71 #endif 72 #if defined(PETSC_HAVE_SCALAPACK) 73 PETSC_EXTERN PetscErrorCode MatCreate_ScaLAPACK(Mat); 74 #endif 75 76 PETSC_EXTERN PetscErrorCode MatCreate_Preallocator(Mat); 77 PETSC_EXTERN PetscErrorCode MatCreate_Dummy(Mat); 78 79 #if defined(PETSC_HAVE_HYPRE) 80 PETSC_EXTERN PetscErrorCode MatCreate_HYPRE(Mat); 81 #endif 82 83 PETSC_EXTERN PetscErrorCode MatCreate_ConstantDiagonal(Mat); 84 85 #if defined(PETSC_HAVE_HARA) 86 PETSC_EXTERN PetscErrorCode MatCreate_HARA(Mat); 87 #endif 88 89 /*@C 90 MatRegisterAll - Registers all of the matrix types in PETSc 91 92 Not Collective 93 94 Level: advanced 95 96 .seealso: MatRegister() 97 @*/ 98 PetscErrorCode MatRegisterAll(void) 99 { 100 PetscErrorCode ierr; 101 102 PetscFunctionBegin; 103 if (MatRegisterAllCalled) PetscFunctionReturn(0); 104 MatRegisterAllCalled = PETSC_TRUE; 105 106 ierr = MatRegister(MATMFFD, MatCreate_MFFD);CHKERRQ(ierr); 107 108 ierr = MatRegister(MATMPIMAIJ, MatCreate_MAIJ);CHKERRQ(ierr); 109 ierr = MatRegister(MATSEQMAIJ, MatCreate_MAIJ);CHKERRQ(ierr); 110 ierr = MatRegister(MATMAIJ, MatCreate_MAIJ);CHKERRQ(ierr); 111 112 ierr = MatRegister(MATMPIKAIJ, MatCreate_KAIJ);CHKERRQ(ierr); 113 ierr = MatRegister(MATSEQKAIJ, MatCreate_KAIJ);CHKERRQ(ierr); 114 ierr = MatRegister(MATKAIJ, MatCreate_KAIJ);CHKERRQ(ierr); 115 116 ierr = MatRegister(MATIS, MatCreate_IS);CHKERRQ(ierr); 117 ierr = MatRegister(MATSHELL, MatCreate_Shell);CHKERRQ(ierr); 118 ierr = MatRegister(MATCOMPOSITE, MatCreate_Composite);CHKERRQ(ierr); 119 120 ierr = MatRegisterRootName(MATAIJ,MATSEQAIJ,MATMPIAIJ);CHKERRQ(ierr); 121 ierr = MatRegister(MATMPIAIJ, MatCreate_MPIAIJ);CHKERRQ(ierr); 122 ierr = MatRegister(MATSEQAIJ, MatCreate_SeqAIJ);CHKERRQ(ierr); 123 124 ierr = MatRegisterRootName(MATAIJPERM,MATSEQAIJPERM,MATMPIAIJPERM);CHKERRQ(ierr); 125 ierr = MatRegister(MATMPIAIJPERM, MatCreate_MPIAIJPERM);CHKERRQ(ierr); 126 ierr = MatRegister(MATSEQAIJPERM, MatCreate_SeqAIJPERM);CHKERRQ(ierr); 127 128 ierr = MatRegisterRootName(MATAIJSELL,MATSEQAIJSELL,MATMPIAIJSELL);CHKERRQ(ierr); 129 ierr = MatRegister(MATMPIAIJSELL, MatCreate_MPIAIJSELL);CHKERRQ(ierr); 130 ierr = MatRegister(MATSEQAIJSELL, MatCreate_SeqAIJSELL);CHKERRQ(ierr); 131 132 #if defined(PETSC_HAVE_MKL_SPARSE) 133 ierr = MatRegisterRootName(MATAIJMKL, MATSEQAIJMKL,MATMPIAIJMKL);CHKERRQ(ierr); 134 ierr = MatRegister(MATMPIAIJMKL, MatCreate_MPIAIJMKL);CHKERRQ(ierr); 135 ierr = MatRegister(MATSEQAIJMKL, MatCreate_SeqAIJMKL);CHKERRQ(ierr); 136 #endif 137 138 #if defined(PETSC_HAVE_MKL_SPARSE_OPTIMIZE) 139 ierr = MatRegisterRootName(MATBAIJMKL,MATSEQBAIJMKL,MATMPIBAIJMKL);CHKERRQ(ierr); 140 ierr = MatRegister(MATMPIBAIJMKL, MatCreate_MPIBAIJMKL);CHKERRQ(ierr); 141 ierr = MatRegister(MATSEQBAIJMKL, MatCreate_SeqBAIJMKL);CHKERRQ(ierr); 142 #endif 143 144 ierr = MatRegisterRootName(MATAIJCRL,MATSEQAIJCRL,MATMPIAIJCRL);CHKERRQ(ierr); 145 ierr = MatRegister(MATSEQAIJCRL, MatCreate_SeqAIJCRL);CHKERRQ(ierr); 146 ierr = MatRegister(MATMPIAIJCRL, MatCreate_MPIAIJCRL);CHKERRQ(ierr); 147 148 ierr = MatRegisterRootName(MATBAIJ,MATSEQBAIJ,MATMPIBAIJ);CHKERRQ(ierr); 149 ierr = MatRegister(MATMPIBAIJ, MatCreate_MPIBAIJ);CHKERRQ(ierr); 150 ierr = MatRegister(MATSEQBAIJ, MatCreate_SeqBAIJ);CHKERRQ(ierr); 151 152 ierr = MatRegisterRootName(MATSBAIJ,MATSEQSBAIJ,MATMPISBAIJ);CHKERRQ(ierr); 153 ierr = MatRegister(MATMPISBAIJ, MatCreate_MPISBAIJ);CHKERRQ(ierr); 154 ierr = MatRegister(MATSEQSBAIJ, MatCreate_SeqSBAIJ);CHKERRQ(ierr); 155 156 ierr = MatRegisterRootName(MATDENSE,MATSEQDENSE,MATMPIDENSE);CHKERRQ(ierr); 157 ierr = MatRegister(MATMPIDENSE, MatCreate_MPIDense);CHKERRQ(ierr); 158 ierr = MatRegister(MATSEQDENSE, MatCreate_SeqDense);CHKERRQ(ierr); 159 #if defined(PETSC_HAVE_CUDA) 160 ierr = MatRegisterRootName(MATDENSECUDA,MATSEQDENSECUDA,MATMPIDENSECUDA);CHKERRQ(ierr); 161 ierr = MatRegister(MATSEQDENSECUDA, MatCreate_SeqDenseCUDA);CHKERRQ(ierr); 162 ierr = MatRegister(MATMPIDENSECUDA, MatCreate_MPIDenseCUDA);CHKERRQ(ierr); 163 #endif 164 165 ierr = MatRegister(MATMPIADJ, MatCreate_MPIAdj);CHKERRQ(ierr); 166 ierr = MatRegister(MATSCATTER, MatCreate_Scatter);CHKERRQ(ierr); 167 ierr = MatRegister(MATBLOCKMAT, MatCreate_BlockMat);CHKERRQ(ierr); 168 ierr = MatRegister(MATNEST, MatCreate_Nest);CHKERRQ(ierr); 169 170 ierr = MatRegisterRootName(MATSELL,MATSEQSELL,MATMPISELL);CHKERRQ(ierr); 171 ierr = MatRegister(MATMPISELL, MatCreate_MPISELL);CHKERRQ(ierr); 172 ierr = MatRegister(MATSEQSELL, MatCreate_SeqSELL);CHKERRQ(ierr); 173 174 #if defined(PETSC_HAVE_CUDA) 175 #if PETSC_PKG_CUDA_VERSION_LT(11,0,0) 176 ierr = MatRegisterRootName(MATAIJCUSPARSE,MATSEQAIJCUSPARSE,MATMPIAIJCUSPARSE);CHKERRQ(ierr); 177 ierr = MatRegister(MATSEQAIJCUSPARSE, MatCreate_SeqAIJCUSPARSE);CHKERRQ(ierr); 178 ierr = MatRegister(MATMPIAIJCUSPARSE, MatCreate_MPIAIJCUSPARSE);CHKERRQ(ierr); 179 #endif 180 #endif 181 182 #if defined(PETSC_HAVE_VIENNACL) 183 ierr = MatRegisterRootName(MATAIJVIENNACL,MATSEQAIJVIENNACL,MATMPIAIJVIENNACL);CHKERRQ(ierr); 184 ierr = MatRegister(MATSEQAIJVIENNACL, MatCreate_SeqAIJViennaCL);CHKERRQ(ierr); 185 ierr = MatRegister(MATMPIAIJVIENNACL, MatCreate_MPIAIJViennaCL);CHKERRQ(ierr); 186 #endif 187 188 #if defined(PETSC_HAVE_FFTW) 189 ierr = MatRegister(MATFFTW, MatCreate_FFTW);CHKERRQ(ierr); 190 #endif 191 #if defined(PETSC_HAVE_ELEMENTAL) 192 ierr = MatRegister(MATELEMENTAL, MatCreate_Elemental);CHKERRQ(ierr); 193 #endif 194 #if defined(PETSC_HAVE_SCALAPACK) 195 ierr = MatRegister(MATSCALAPACK, MatCreate_ScaLAPACK);CHKERRQ(ierr); 196 #endif 197 198 ierr = MatRegister(MATPREALLOCATOR, MatCreate_Preallocator);CHKERRQ(ierr); 199 ierr = MatRegister(MATDUMMY, MatCreate_Dummy);CHKERRQ(ierr); 200 201 ierr = MatRegister(MATCONSTANTDIAGONAL,MatCreate_ConstantDiagonal);CHKERRQ(ierr); 202 203 #if defined(PETSC_HAVE_HYPRE) 204 ierr = MatRegister(MATHYPRE, MatCreate_HYPRE);CHKERRQ(ierr); 205 #endif 206 207 #if defined(PETSC_HAVE_HARA) 208 ierr = MatRegister(MATHARA, MatCreate_HARA);CHKERRQ(ierr); 209 #endif 210 PetscFunctionReturn(0); 211 } 212