1 /* $Id: mat.h,v 1.176 1999/04/02 15:05:18 bsmith Exp bsmith $ */ 2 /* 3 Include file for the matrix component of PETSc 4 */ 5 #ifndef __MAT_H 6 #define __MAT_H 7 #include "vec.h" 8 9 #define MAT_COOKIE PETSC_COOKIE+5 10 11 typedef struct _p_Mat* Mat; 12 13 #define MAX_MATRIX_TYPES 14 14 /* 15 The default matrix data storage formats and routines to create them. 16 17 MATLASTTYPE is "end-of-list" marker that can be used to check that 18 MAX_MATRIX_TYPES is large enough. The rule is 19 MAX_MATRIX_TYPES >= MATLASTTYPE . 20 21 To do: add a test program that checks the consistency of these values. 22 */ 23 typedef enum { MATSAME=-1, MATSEQDENSE, MATSEQAIJ, MATMPIAIJ, MATSHELL, 24 MATMPIROWBS, MATSEQBDIAG, MATMPIBDIAG, MATMPIDENSE, MATSEQBAIJ, 25 MATMPIBAIJ, MATMPICSN, MATSEQCSN, MATSEQADJ, MATMPIADJ, 26 MATLASTTYPE } MatType; 27 28 extern int MatCreate(MPI_Comm,int,int,int,int,Mat*); 29 extern int MatCreateSeqDense(MPI_Comm,int,int,Scalar*,Mat*); 30 extern int MatCreateMPIDense(MPI_Comm,int,int,int,int,Scalar*,Mat*); 31 extern int MatCreateSeqAIJ(MPI_Comm,int,int,int,int*,Mat*); 32 extern int MatCreateMPIAIJ(MPI_Comm,int,int,int,int,int,int*,int,int*,Mat*); 33 extern int MatCreateMPIRowbs(MPI_Comm,int,int,int,int*,void*,Mat*); 34 extern int MatCreateSeqBDiag(MPI_Comm,int,int,int,int,int*,Scalar**,Mat*); 35 extern int MatCreateMPIBDiag(MPI_Comm,int,int,int,int,int,int*,Scalar**,Mat*); 36 extern int MatCreateSeqBAIJ(MPI_Comm,int,int,int,int,int*,Mat*); 37 extern int MatCreateMPIBAIJ(MPI_Comm,int,int,int,int,int,int,int*,int,int*,Mat*); 38 extern int MatCreateSeqAdj(MPI_Comm,int,int,int*,int*,Mat *); 39 extern int MatCreateMPIAdj(MPI_Comm,int,int,int*,int*,Mat*); 40 41 extern int MatDestroy(Mat); 42 43 extern int MatCreateShell(MPI_Comm,int,int,int,int,void *,Mat*); 44 extern int MatShellGetContext(Mat,void **); 45 46 extern int MatPrintHelp(Mat); 47 extern int MatGetMaps(Mat,Map*,Map*); 48 49 /* ------------------------------------------------------------*/ 50 extern int MatSetValues(Mat,int,int*,int,int*,Scalar*,InsertMode); 51 extern int MatSetValuesBlocked(Mat,int,int*,int,int*,Scalar*,InsertMode); 52 53 typedef enum {MAT_FLUSH_ASSEMBLY=1,MAT_FINAL_ASSEMBLY=0} MatAssemblyType; 54 extern int MatAssemblyBegin(Mat,MatAssemblyType); 55 extern int MatAssemblyEnd(Mat,MatAssemblyType); 56 #define MatSetValue(v,i,j,va,mode) \ 57 {int _ierr,_row = i,_col = j; Scalar _va = va; \ 58 _ierr = MatSetValues(v,1,&_row,1,&_col,&_va,mode);CHKERRQ(_ierr); \ 59 } 60 #define MatGetValue(v,i,j,va) \ 61 {int _ierr,_row = i,_col = j; \ 62 _ierr = MatGetValues(v,1,&_row,1,&_col,&va);CHKERRQ(_ierr); \ 63 } 64 /* 65 Any additions/changes here MUST also be made in include/finclude/mat.h 66 */ 67 typedef enum {MAT_ROW_ORIENTED=1,MAT_COLUMN_ORIENTED=2,MAT_ROWS_SORTED=4, 68 MAT_COLUMNS_SORTED=8,MAT_NO_NEW_NONZERO_LOCATIONS=16, 69 MAT_YES_NEW_NONZERO_LOCATIONS=32,MAT_SYMMETRIC=64, 70 MAT_STRUCTURALLY_SYMMETRIC=65,MAT_NO_NEW_DIAGONALS=66, 71 MAT_YES_NEW_DIAGONALS=67,MAT_INODE_LIMIT_1=68,MAT_INODE_LIMIT_2=69, 72 MAT_INODE_LIMIT_3=70,MAT_INODE_LIMIT_4=71,MAT_INODE_LIMIT_5=72, 73 MAT_IGNORE_OFF_PROC_ENTRIES=73,MAT_ROWS_UNSORTED=74, 74 MAT_COLUMNS_UNSORTED=75,MAT_NEW_NONZERO_LOCATION_ERR=76, 75 MAT_NEW_NONZERO_ALLOCATION_ERR=77,MAT_USE_HASH_TABLE=78} MatOption; 76 extern int MatSetOption(Mat,MatOption); 77 extern int MatGetType(Mat,MatType*,char**); 78 extern int MatGetTypeFromOptions(MPI_Comm,char*,MatType*,PetscTruth*); 79 80 extern int MatGetValues(Mat,int,int*,int,int*,Scalar*); 81 extern int MatGetRow(Mat,int,int *,int **,Scalar**); 82 extern int MatRestoreRow(Mat,int,int *,int **,Scalar**); 83 extern int MatGetColumn(Mat,int,int *,int **,Scalar**); 84 extern int MatRestoreColumn(Mat,int,int *,int **,Scalar**); 85 extern int MatGetColumnVector(Mat,Vec,int); 86 extern int MatGetArray(Mat,Scalar **); 87 extern int MatRestoreArray(Mat,Scalar **); 88 extern int MatGetBlockSize(Mat,int *); 89 90 extern int MatMult(Mat,Vec,Vec); 91 extern int MatMultAdd(Mat,Vec,Vec,Vec); 92 extern int MatMultTrans(Mat,Vec,Vec); 93 extern int MatMultTransAdd(Mat,Vec,Vec,Vec); 94 95 typedef enum {MAT_DO_NOT_COPY_VALUES, MAT_COPY_VALUES} MatDuplicateOption; 96 97 extern int MatConvert(Mat,MatType,Mat*); 98 extern int MatDuplicate(Mat,MatDuplicateOption,Mat*); 99 extern int MatConvertRegister(MatType,MatType,int (*)(Mat,MatType,Mat*)); 100 extern int MatConvertRegisterAll(void); 101 102 typedef enum {SAME_NONZERO_PATTERN,DIFFERENT_NONZERO_PATTERN,SAME_PRECONDITIONER} MatStructure; 103 104 extern int MatCopy(Mat,Mat,MatStructure); 105 extern int MatView(Mat,Viewer); 106 extern int MatLoad(Viewer,MatType,Mat*); 107 extern int MatLoadRegister(MatType,int (*)(Viewer,MatType,Mat*)); 108 extern int MatLoadRegisterAll(void); 109 110 extern int MatGetRowIJ(Mat,int,PetscTruth,int*,int **,int **,PetscTruth *); 111 extern int MatRestoreRowIJ(Mat,int,PetscTruth,int *,int **,int **,PetscTruth *); 112 extern int MatGetColumnIJ(Mat,int,PetscTruth,int*,int **,int **,PetscTruth *); 113 extern int MatRestoreColumnIJ(Mat,int,PetscTruth,int *,int **,int **,PetscTruth *); 114 115 /* 116 Context of matrix information, used with MatGetInfo() 117 Note: If any entries are added to this context, be sure 118 to adjust MAT_INFO_SIZE in finclude/mat.h 119 */ 120 typedef struct { 121 PLogDouble rows_global, columns_global; /* number of global rows and columns */ 122 PLogDouble rows_local, columns_local; /* number of local rows and columns */ 123 PLogDouble block_size; /* block size */ 124 PLogDouble nz_allocated, nz_used, nz_unneeded; /* number of nonzeros */ 125 PLogDouble memory; /* memory allocated */ 126 PLogDouble assemblies; /* number of matrix assemblies */ 127 PLogDouble mallocs; /* number of mallocs during MatSetValues() */ 128 PLogDouble fill_ratio_given, fill_ratio_needed; /* fill ratio for LU/ILU */ 129 PLogDouble factor_mallocs; /* number of mallocs during factorization */ 130 } MatInfo; 131 132 typedef enum {MAT_LOCAL=1,MAT_GLOBAL_MAX=2,MAT_GLOBAL_SUM=3} MatInfoType; 133 extern int MatGetInfo(Mat,MatInfoType,MatInfo*); 134 extern int MatValid(Mat,PetscTruth*); 135 extern int MatGetDiagonal(Mat,Vec); 136 extern int MatTranspose(Mat,Mat*); 137 extern int MatPermute(Mat,IS,IS,Mat *); 138 extern int MatDiagonalScale(Mat,Vec,Vec); 139 extern int MatDiagonalShift(Mat,Vec); 140 extern int MatEqual(Mat,Mat, PetscTruth*); 141 142 extern int MatNorm(Mat,NormType,double *); 143 extern int MatZeroEntries(Mat); 144 extern int MatZeroRows(Mat,IS,Scalar*); 145 extern int MatZeroColumns(Mat,IS,Scalar*); 146 147 extern int MatUseScaledForm(Mat,PetscTruth); 148 extern int MatScaleSystem(Mat,Vec,Vec); 149 extern int MatUnScaleSystem(Mat,Vec,Vec); 150 151 extern int MatGetSize(Mat,int*,int*); 152 extern int MatGetLocalSize(Mat,int*,int*); 153 extern int MatGetOwnershipRange(Mat,int*,int*); 154 155 typedef enum {MAT_INITIAL_MATRIX, MAT_REUSE_MATRIX} MatReuse; 156 extern int MatGetSubMatrices(Mat,int,IS *,IS *,MatReuse,Mat **); 157 extern int MatDestroyMatrices(int, Mat **); 158 extern int MatGetSubMatrix(Mat,IS,IS,int,MatReuse,Mat *); 159 160 extern int MatIncreaseOverlap(Mat,int,IS *,int); 161 162 extern int MatAXPY(Scalar *,Mat,Mat); 163 extern int MatAYPX(Scalar *,Mat,Mat); 164 extern int MatCompress(Mat); 165 166 extern int MatScale(Scalar *,Mat); 167 extern int MatShift(Scalar *,Mat); 168 169 extern int MatSetLocalToGlobalMapping(Mat, ISLocalToGlobalMapping); 170 extern int MatSetLocalToGlobalMappingBlocked(Mat, ISLocalToGlobalMapping); 171 extern int MatZeroRowsLocal(Mat,IS,Scalar*); 172 extern int MatSetValuesLocal(Mat,int,int*,int,int*,Scalar*,InsertMode); 173 extern int MatSetValuesBlockedLocal(Mat,int,int*,int,int*,Scalar*,InsertMode); 174 175 extern int MatSetStashInitialSize(Mat,int, int); 176 177 /* Routines unique to particular data structures */ 178 extern int MatBDiagGetData(Mat,int*,int*,int**,int**,Scalar***); 179 extern int MatSeqAIJSetColumnIndices(Mat,int *); 180 extern int MatSeqBAIJSetColumnIndices(Mat,int *); 181 182 extern int MatStoreValues(Mat); 183 extern int MatRetrieveValues(Mat); 184 185 /* 186 These routines are not usually accessed directly, rather solving is 187 done through the SLES, KSP and PC interfaces. 188 */ 189 190 typedef char* MatOrderingType; 191 #define MATORDERING_NATURAL "natural" 192 #define MATORDERING_ND "nd" 193 #define MATORDERING_1WD "1wd" 194 #define MATORDERING_RCM "rcm" 195 #define MATORDERING_QMD "qmd" 196 #define MATORDERING_ROWLENGTH "rowlength" 197 198 extern int MatGetOrdering(Mat,MatOrderingType,IS*,IS*); 199 extern int MatOrderingRegister_Private(char*,char*,char*,int(*)(Mat,MatOrderingType,IS*,IS*)); 200 #if defined(USE_DYNAMIC_LIBRARIES) 201 #define MatOrderingRegister(a,b,c,d) MatOrderingRegister_Private(a,b,c,0) 202 #else 203 #define MatOrderingRegister(a,b,c,d) MatOrderingRegister_Private(a,b,c,d) 204 #endif 205 extern int MatOrderingRegisterDestroy(void); 206 extern int MatOrderingRegisterAll(char*); 207 extern int MatOrderingRegisterAllCalled; 208 209 extern int MatReorderForNonzeroDiagonal(Mat,double,IS,IS); 210 211 extern int MatCholeskyFactor(Mat,IS,double); 212 extern int MatCholeskyFactorSymbolic(Mat,IS,double,Mat*); 213 extern int MatCholeskyFactorNumeric(Mat,Mat*); 214 215 /* 216 Context of matrix information, used with MatILUFactor() and MatILUFactorSymbolic() 217 Note: If any entries are added to this context, be sure 218 to adjust MAT_ILUINFO_SIZE in finclude/mat.h 219 220 Note: The integer values below are passed in double to allow easy use from 221 Fortran 222 */ 223 typedef struct { 224 double levels; /* ILU(levels) */ 225 double fill; /* expected fill; nonzeros in factored matrix/nonzeros in original matrix*/ 226 double diagonal_fill; /* force diagonal to fill in if initially not filled */ 227 } MatILUInfo; 228 229 extern int MatLUFactor(Mat,IS,IS,double); 230 extern int MatILUFactor(Mat,IS,IS,MatILUInfo*); 231 extern int MatLUFactorSymbolic(Mat,IS,IS,double,Mat*); 232 extern int MatILUFactorSymbolic(Mat,IS,IS,MatILUInfo*,Mat*); 233 extern int MatIncompleteCholeskyFactorSymbolic(Mat,IS,double,int,Mat*); 234 extern int MatLUFactorNumeric(Mat,Mat*); 235 extern int MatILUDTFactor(Mat,double,int,IS,IS,Mat *); 236 237 extern int MatSolve(Mat,Vec,Vec); 238 extern int MatForwardSolve(Mat,Vec,Vec); 239 extern int MatBackwardSolve(Mat,Vec,Vec); 240 extern int MatSolveAdd(Mat,Vec,Vec,Vec); 241 extern int MatSolveTrans(Mat,Vec,Vec); 242 extern int MatSolveTransAdd(Mat,Vec,Vec,Vec); 243 244 extern int MatSetUnfactored(Mat); 245 246 /* MatSORType may be bitwise ORd together, so do not change the numbers */ 247 248 typedef enum {SOR_FORWARD_SWEEP=1,SOR_BACKWARD_SWEEP=2,SOR_SYMMETRIC_SWEEP=3, 249 SOR_LOCAL_FORWARD_SWEEP=4,SOR_LOCAL_BACKWARD_SWEEP=8, 250 SOR_LOCAL_SYMMETRIC_SWEEP=12,SOR_ZERO_INITIAL_GUESS=16, 251 SOR_EISENSTAT=32,SOR_APPLY_UPPER=64,SOR_APPLY_LOWER=128} MatSORType; 252 extern int MatRelax(Mat,Vec,double,MatSORType,double,int,Vec); 253 254 /* 255 These routines are for efficiently computing Jacobians via finite differences. 256 */ 257 258 typedef char* MatColoringType; 259 #define MATCOLORING_NATURAL "natural" 260 #define MATCOLORING_SL "sl" 261 #define MATCOLORING_LF "lf" 262 #define MATCOLORING_ID "id" 263 264 extern int MatGetColoring(Mat,MatColoringType,ISColoring*); 265 extern int MatColoringRegister_Private(char*,char*,char*,int(*)(Mat,MatColoringType,ISColoring *)); 266 #if defined(USE_DYNAMIC_LIBRARIES) 267 #define MatColoringRegister(a,b,c,d) MatColoringRegister_Private(a,b,c,0) 268 #else 269 #define MatColoringRegister(a,b,c,d) MatColoringRegister_Private(a,b,c,d) 270 #endif 271 extern int MatColoringRegisterAll(char *); 272 extern int MatColoringRegisterAllCalled; 273 extern int MatColoringRegisterDestroy(void); 274 extern int MatColoringPatch(Mat,int,int *,ISColoring*); 275 276 /* 277 Data structures used to compute Jacobian vector products 278 efficiently using finite differences. 279 */ 280 #define MAT_FDCOLORING_COOKIE PETSC_COOKIE + 22 281 282 typedef struct _p_MatFDColoring *MatFDColoring; 283 284 extern int MatFDColoringCreate(Mat,ISColoring,MatFDColoring *); 285 extern int MatFDColoringDestroy(MatFDColoring); 286 extern int MatFDColoringView(MatFDColoring,Viewer); 287 extern int MatFDColoringSetFunction(MatFDColoring,int (*)(void),void*); 288 extern int MatFDColoringSetParameters(MatFDColoring,double,double); 289 extern int MatFDColoringSetFrequency(MatFDColoring,int); 290 extern int MatFDColoringGetFrequency(MatFDColoring,int*); 291 extern int MatFDColoringSetFromOptions(MatFDColoring); 292 extern int MatFDColoringPrintHelp(MatFDColoring); 293 extern int MatFDColoringApply(Mat,MatFDColoring,Vec,MatStructure*,void *); 294 extern int MatFDColoringApplyTS(Mat,MatFDColoring,double,Vec,MatStructure*,void *); 295 296 /* 297 These routines are for partitioning matrices: currently used only 298 for adjacency matrix, MatCreateSeqAdj() or MatCreateMPIAdj(). 299 */ 300 #define MATPARTITIONING_COOKIE PETSC_COOKIE + 25 301 302 typedef struct _p_MatPartitioning *MatPartitioning; 303 typedef char* MatPartitioningType; 304 #define MATPARTITIONING_CURRENT "current" 305 #define MATPARTITIONING_PARMETIS "parmetis" 306 307 extern int MatPartitioningCreate(MPI_Comm,MatPartitioning*); 308 extern int MatPartitioningSetType(MatPartitioning,MatPartitioningType); 309 extern int MatPartitioningSetAdjacency(MatPartitioning,Mat); 310 extern int MatPartitioningSetVertexWeights(MatPartitioning,double*); 311 extern int MatPartitioningApply(MatPartitioning,IS*); 312 extern int MatPartitioningDestroy(MatPartitioning); 313 314 extern int MatPartitioningRegister_Private(char*,char*,char*,int(*)(MatPartitioning)); 315 #if defined(USE_DYNAMIC_LIBRARIES) 316 #define MatPartitioningRegister(a,b,c,d) MatPartitioningRegister_Private(a,b,c,0) 317 #else 318 #define MatPartitioningRegister(a,b,c,d) MatPartitioningRegister_Private(a,b,c,d) 319 #endif 320 321 extern int MatPartitioningRegisterAll(char *); 322 extern int MatPartitioningRegisterAllCalled; 323 extern int MatPartitioningRegisterDestroy(void); 324 extern int MatPartitioningView(MatPartitioning,Viewer); 325 extern int MatPartitioningSetFromOptions(MatPartitioning); 326 extern int MatPartitioningPrintHelp(MatPartitioning); 327 extern int MatPartitioningGetType(MatPartitioning,MatPartitioningType*); 328 329 extern int MatPartitioningParmetisSetCoarseSequential(MatPartitioning); 330 331 /* 332 If you add entries here you must also add them to finclude/mat.h 333 */ 334 typedef enum { MATOP_SET_VALUES=0, 335 MATOP_GET_ROW=1, 336 MATOP_RESTORE_ROW=2, 337 MATOP_MULT=3, 338 MATOP_MULT_ADD=4, 339 MATOP_MULT_TRANS=5, 340 MATOP_MULT_TRANS_ADD=6, 341 MATOP_SOLVE=7, 342 MATOP_SOLVE_ADD=8, 343 MATOP_SOLVE_TRANS=9, 344 MATOP_SOLVE_TRANS_ADD=10, 345 MATOP_LUFACTOR=11, 346 MATOP_CHOLESKYFACTOR=12, 347 MATOP_RELAX=13, 348 MATOP_TRANSPOSE=14, 349 MATOP_GETINFO=15, 350 MATOP_EQUAL=16, 351 MATOP_GET_DIAGONAL=17, 352 MATOP_DIAGONAL_SCALE=18, 353 MATOP_NORM=19, 354 MATOP_ASSEMBLY_BEGIN=20, 355 MATOP_ASSEMBLY_END=21, 356 MATOP_COMPRESS=22, 357 MATOP_SET_OPTION=23, 358 MATOP_ZERO_ENTRIES=24, 359 MATOP_ZERO_ROWS=25, 360 MATOP_LUFACTOR_SYMBOLIC=26, 361 MATOP_LUFACTOR_NUMERIC=27, 362 MATOP_CHOLESKY_FACTOR_SYMBOLIC=28, 363 MATOP_CHOLESKY_FACTOR_NUMERIC=29, 364 MATOP_GET_SIZE=30, 365 MATOP_GET_LOCAL_SIZE=31, 366 MATOP_GET_OWNERSHIP_RANGE=32, 367 MATOP_ILUFACTOR_SYMBOLIC=33, 368 MATOP_INCOMPLETECHOLESKYFACTOR_SYMBOLIC=34, 369 MATOP_GET_ARRAY=35, 370 MATOP_RESTORE_ARRAY=36, 371 372 MATOP_CONVERT_SAME_TYPE=37, 373 MATOP_FORWARD_SOLVE=38, 374 MATOP_BACKWARD_SOLVE=39, 375 MATOP_ILUFACTOR=40, 376 MATOP_INCOMPLETECHOLESKYFACTOR=41, 377 MATOP_AXPY=42, 378 MATOP_GET_SUBMATRICES=43, 379 MATOP_INCREASE_OVERLAP=44, 380 MATOP_GET_VALUES=45, 381 MATOP_COPY=46, 382 MATOP_PRINT_HELP=47, 383 MATOP_SCALE=48, 384 MATOP_SHIFT=49, 385 MATOP_DIAGONAL_SHIFT=50, 386 MATOP_ILUDT_FACTOR=51, 387 MATOP_GET_BLOCK_SIZE=52, 388 MATOP_GET_ROW_IJ=53, 389 MATOP_RESTORE_ROW_IJ=54, 390 MATOP_GET_COLUMN_IJ=55, 391 MATOP_RESTORE_COLUMN_IJ=56, 392 MATOP_FDCOLORING_CREATE=57, 393 MATOP_COLORING_PATCH=58, 394 MATOP_SET_UNFACTORED=59, 395 MATOP_PERMUTE=60, 396 MATOP_SET_VALUES_BLOCKED=61, 397 MATOP_DESTROY=250, 398 MATOP_VIEW=251 399 } MatOperation; 400 extern int MatHasOperation(Mat,MatOperation,PetscTruth*); 401 extern int MatShellSetOperation(Mat,MatOperation,void *); 402 extern int MatShellGetOperation(Mat,MatOperation,void **); 403 404 /* 405 Codes for matrices stored on disk. By default they are 406 stored in a universal format. By changing the format with 407 ViewerSetFormat(viewer,VIEWER_FORMAT_BINARY_NATIVE); the matrices will 408 be stored in a way natural for the matrix, for example dense matrices 409 would be stored as dense. Matrices stored this way may only be 410 read into matrices of the same time. 411 */ 412 #define MATRIX_BINARY_FORMAT_DENSE -1 413 414 /* 415 New matrix classes not yet distributed 416 */ 417 /* 418 MatAIJIndices is a data structure for storing the nonzero location information 419 for sparse matrices. Several matrices with identical nonzero structure can share 420 the same MatAIJIndices. 421 */ 422 typedef struct _p_MatAIJIndices* MatAIJIndices; 423 424 extern int MatCreateAIJIndices(int,int,int*,int*,PetscTruth,MatAIJIndices*); 425 extern int MatCreateAIJIndicesEmpty(int,int,int*,PetscTruth,MatAIJIndices*); 426 extern int MatAttachAIJIndices(MatAIJIndices,MatAIJIndices*); 427 extern int MatDestroyAIJIndices(MatAIJIndices); 428 extern int MatCopyAIJIndices(MatAIJIndices,MatAIJIndices*); 429 extern int MatValidateAIJIndices(int,MatAIJIndices); 430 extern int MatShiftAIJIndices(MatAIJIndices); 431 extern int MatShrinkAIJIndices(MatAIJIndices); 432 extern int MatTransposeAIJIndices(MatAIJIndices, MatAIJIndices*); 433 434 extern int MatCreateSeqCSN(MPI_Comm,int,int,int*,int,Mat*); 435 extern int MatCreateSeqCSN_Single(MPI_Comm,int,int,int*,int,Mat*); 436 extern int MatCreateSeqCSNWithPrecision(MPI_Comm,int,int,int*,int,ScalarPrecision,Mat*); 437 438 extern int MatCreateSeqCSNIndices(MPI_Comm,MatAIJIndices,int,Mat *); 439 extern int MatCreateSeqCSNIndices_Single(MPI_Comm,MatAIJIndices,int,Mat *); 440 extern int MatCreateSeqCSNIndicesWithPrecision(MPI_Comm,MatAIJIndices,int,ScalarPrecision,Mat *); 441 442 extern int MatMPIBAIJSetHashTableFactor(Mat,double); 443 extern int MatSeqAIJGetInodeSizes(Mat,int *,int *[],int *); 444 445 446 #endif 447 448 449 450