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