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