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 PETSC_EXTERN PetscErrorCode MatCreate_Elemental(Mat); 69 70 PETSC_EXTERN PetscErrorCode MatCreate_Preallocator(Mat); 71 PETSC_EXTERN PetscErrorCode MatCreate_Dummy(Mat); 72 73 #if defined(PETSC_HAVE_HYPRE) 74 PETSC_EXTERN PetscErrorCode MatCreate_HYPRE(Mat); 75 #endif 76 77 PETSC_EXTERN PetscErrorCode MatCreate_ConstantDiagonal(Mat); 78 79 #if defined(PETSC_HAVE_HARA) 80 PETSC_EXTERN PetscErrorCode MatCreate_HARA(Mat); 81 #endif 82 83 /*@C 84 MatRegisterAll - Registers all of the matrix types in PETSc 85 86 Not Collective 87 88 Level: advanced 89 90 .seealso: MatRegister() 91 @*/ 92 PetscErrorCode MatRegisterAll(void) 93 { 94 PetscErrorCode ierr; 95 96 PetscFunctionBegin; 97 if (MatRegisterAllCalled) PetscFunctionReturn(0); 98 MatRegisterAllCalled = PETSC_TRUE; 99 100 ierr = MatRegister(MATMFFD, MatCreate_MFFD);CHKERRQ(ierr); 101 102 ierr = MatRegister(MATMPIMAIJ, MatCreate_MAIJ);CHKERRQ(ierr); 103 ierr = MatRegister(MATSEQMAIJ, MatCreate_MAIJ);CHKERRQ(ierr); 104 ierr = MatRegister(MATMAIJ, MatCreate_MAIJ);CHKERRQ(ierr); 105 106 ierr = MatRegister(MATMPIKAIJ, MatCreate_KAIJ);CHKERRQ(ierr); 107 ierr = MatRegister(MATSEQKAIJ, MatCreate_KAIJ);CHKERRQ(ierr); 108 ierr = MatRegister(MATKAIJ, MatCreate_KAIJ);CHKERRQ(ierr); 109 110 ierr = MatRegister(MATIS, MatCreate_IS);CHKERRQ(ierr); 111 ierr = MatRegister(MATSHELL, MatCreate_Shell);CHKERRQ(ierr); 112 ierr = MatRegister(MATCOMPOSITE, MatCreate_Composite);CHKERRQ(ierr); 113 114 ierr = MatRegisterRootName(MATAIJ,MATSEQAIJ,MATMPIAIJ);CHKERRQ(ierr); 115 ierr = MatRegister(MATMPIAIJ, MatCreate_MPIAIJ);CHKERRQ(ierr); 116 ierr = MatRegister(MATSEQAIJ, MatCreate_SeqAIJ);CHKERRQ(ierr); 117 118 ierr = MatRegisterRootName(MATAIJPERM,MATSEQAIJPERM,MATMPIAIJPERM);CHKERRQ(ierr); 119 ierr = MatRegister(MATMPIAIJPERM, MatCreate_MPIAIJPERM);CHKERRQ(ierr); 120 ierr = MatRegister(MATSEQAIJPERM, MatCreate_SeqAIJPERM);CHKERRQ(ierr); 121 122 ierr = MatRegisterRootName(MATAIJSELL,MATSEQAIJSELL,MATMPIAIJSELL);CHKERRQ(ierr); 123 ierr = MatRegister(MATMPIAIJSELL, MatCreate_MPIAIJSELL);CHKERRQ(ierr); 124 ierr = MatRegister(MATSEQAIJSELL, MatCreate_SeqAIJSELL);CHKERRQ(ierr); 125 126 #if defined(PETSC_HAVE_MKL_SPARSE) 127 ierr = MatRegisterRootName(MATAIJMKL, MATSEQAIJMKL,MATMPIAIJMKL);CHKERRQ(ierr); 128 ierr = MatRegister(MATMPIAIJMKL, MatCreate_MPIAIJMKL);CHKERRQ(ierr); 129 ierr = MatRegister(MATSEQAIJMKL, MatCreate_SeqAIJMKL);CHKERRQ(ierr); 130 #endif 131 132 #if defined(PETSC_HAVE_MKL_SPARSE_OPTIMIZE) 133 ierr = MatRegisterRootName(MATBAIJMKL,MATSEQBAIJMKL,MATMPIBAIJMKL);CHKERRQ(ierr); 134 ierr = MatRegister(MATMPIBAIJMKL, MatCreate_MPIBAIJMKL);CHKERRQ(ierr); 135 ierr = MatRegister(MATSEQBAIJMKL, MatCreate_SeqBAIJMKL);CHKERRQ(ierr); 136 #endif 137 138 ierr = MatRegisterRootName(MATAIJCRL,MATSEQAIJCRL,MATMPIAIJCRL);CHKERRQ(ierr); 139 ierr = MatRegister(MATSEQAIJCRL, MatCreate_SeqAIJCRL);CHKERRQ(ierr); 140 ierr = MatRegister(MATMPIAIJCRL, MatCreate_MPIAIJCRL);CHKERRQ(ierr); 141 142 ierr = MatRegisterRootName(MATBAIJ,MATSEQBAIJ,MATMPIBAIJ);CHKERRQ(ierr); 143 ierr = MatRegister(MATMPIBAIJ, MatCreate_MPIBAIJ);CHKERRQ(ierr); 144 ierr = MatRegister(MATSEQBAIJ, MatCreate_SeqBAIJ);CHKERRQ(ierr); 145 146 ierr = MatRegisterRootName(MATSBAIJ,MATSEQSBAIJ,MATMPISBAIJ);CHKERRQ(ierr); 147 ierr = MatRegister(MATMPISBAIJ, MatCreate_MPISBAIJ);CHKERRQ(ierr); 148 ierr = MatRegister(MATSEQSBAIJ, MatCreate_SeqSBAIJ);CHKERRQ(ierr); 149 150 ierr = MatRegisterRootName(MATDENSE,MATSEQDENSE,MATMPIDENSE);CHKERRQ(ierr); 151 ierr = MatRegister(MATMPIDENSE, MatCreate_MPIDense);CHKERRQ(ierr); 152 ierr = MatRegister(MATSEQDENSE, MatCreate_SeqDense);CHKERRQ(ierr); 153 #if defined(PETSC_HAVE_CUDA) 154 ierr = MatRegisterRootName(MATDENSECUDA,MATSEQDENSECUDA,MATMPIDENSECUDA);CHKERRQ(ierr); 155 ierr = MatRegister(MATSEQDENSECUDA, MatCreate_SeqDenseCUDA);CHKERRQ(ierr); 156 ierr = MatRegister(MATMPIDENSECUDA, MatCreate_MPIDenseCUDA);CHKERRQ(ierr); 157 #endif 158 159 ierr = MatRegister(MATMPIADJ, MatCreate_MPIAdj);CHKERRQ(ierr); 160 ierr = MatRegister(MATSCATTER, MatCreate_Scatter);CHKERRQ(ierr); 161 ierr = MatRegister(MATBLOCKMAT, MatCreate_BlockMat);CHKERRQ(ierr); 162 ierr = MatRegister(MATNEST, MatCreate_Nest);CHKERRQ(ierr); 163 164 ierr = MatRegisterRootName(MATSELL,MATSEQSELL,MATMPISELL);CHKERRQ(ierr); 165 ierr = MatRegister(MATMPISELL, MatCreate_MPISELL);CHKERRQ(ierr); 166 ierr = MatRegister(MATSEQSELL, MatCreate_SeqSELL);CHKERRQ(ierr); 167 168 #if defined(PETSC_HAVE_CUDA) 169 ierr = MatRegisterRootName(MATAIJCUSPARSE,MATSEQAIJCUSPARSE,MATMPIAIJCUSPARSE);CHKERRQ(ierr); 170 ierr = MatRegister(MATSEQAIJCUSPARSE, MatCreate_SeqAIJCUSPARSE);CHKERRQ(ierr); 171 ierr = MatRegister(MATMPIAIJCUSPARSE, MatCreate_MPIAIJCUSPARSE);CHKERRQ(ierr); 172 #endif 173 174 #if defined(PETSC_HAVE_VIENNACL) 175 ierr = MatRegisterRootName(MATAIJVIENNACL,MATSEQAIJVIENNACL,MATMPIAIJVIENNACL);CHKERRQ(ierr); 176 ierr = MatRegister(MATSEQAIJVIENNACL, MatCreate_SeqAIJViennaCL);CHKERRQ(ierr); 177 ierr = MatRegister(MATMPIAIJVIENNACL, MatCreate_MPIAIJViennaCL);CHKERRQ(ierr); 178 #endif 179 180 #if defined(PETSC_HAVE_FFTW) 181 ierr = MatRegister(MATFFTW, MatCreate_FFTW);CHKERRQ(ierr); 182 #endif 183 #if defined(PETSC_HAVE_ELEMENTAL) 184 ierr = MatRegister(MATELEMENTAL, MatCreate_Elemental);CHKERRQ(ierr); 185 #endif 186 187 ierr = MatRegister(MATPREALLOCATOR, MatCreate_Preallocator);CHKERRQ(ierr); 188 ierr = MatRegister(MATDUMMY, MatCreate_Dummy);CHKERRQ(ierr); 189 190 ierr = MatRegister(MATCONSTANTDIAGONAL,MatCreate_ConstantDiagonal);CHKERRQ(ierr); 191 192 #if defined(PETSC_HAVE_HYPRE) 193 ierr = MatRegister(MATHYPRE, MatCreate_HYPRE);CHKERRQ(ierr); 194 #endif 195 196 #if defined(PETSC_HAVE_HARA) 197 ierr = MatRegister(MATHARA, MatCreate_HARA);CHKERRQ(ierr); 198 #endif 199 PetscFunctionReturn(0); 200 } 201