12eac72dbSBarry Smith /* 22eac72dbSBarry Smith Include file for the matrix component of PETSc 32eac72dbSBarry Smith */ 41eb62cbbSBarry Smith #ifndef __MAT_PACKAGE 51eb62cbbSBarry Smith #define __MAT_PACKAGE 68ed539a5SBarry Smith #include "vec.h" 72eac72dbSBarry Smith 8f0479e8cSBarry Smith #define MAT_COOKIE 0x404040 9f0479e8cSBarry Smith 102eac72dbSBarry Smith typedef struct _Mat* Mat; 112eac72dbSBarry Smith typedef struct _MatScatterCtx* MatScatterCtx; 122eac72dbSBarry Smith 132eac72dbSBarry Smith 14*28988994SBarry Smith #define MATDENSESEQ 1 15*28988994SBarry Smith #define MATAIJSEQ 2 16*28988994SBarry Smith #define MATAIJMPI 3 17*28988994SBarry Smith 188ed539a5SBarry Smith extern int MatCreateSequentialDense(int,int,Mat*); 19aabeff55SBarry Smith extern int MatCreateSequentialAIJ(int,int,int,int *,Mat*); 209717bf64SBarry Smith extern int MatCreateMPIAIJ(MPI_Comm,int,int,int,int,int,int*,int,int*,Mat*); 219717bf64SBarry Smith 22aabeff55SBarry Smith extern int MatShellCreate(int,int,void *,Mat*); 23aabeff55SBarry Smith extern int MatShellSetMult(Mat,int (*)(void*,Vec,Vec)); 24aabeff55SBarry Smith extern int MatShellSetMultTrans(Mat,int (*)(void*,Vec,Vec)); 25f0479e8cSBarry Smith extern int MatShellSetMultTransAdd(Mat,int (*)(void*,Vec,Vec,Vec)); 262eac72dbSBarry Smith 278ed539a5SBarry Smith /* ------------------------------------------------------------*/ 288ed539a5SBarry Smith extern int MatValidMatrix(Mat); 292eac72dbSBarry Smith 3020563c6bSBarry Smith extern int MatSetValues(Mat,int,int*,int,int*,Scalar*,InsertMode); 318ed539a5SBarry Smith extern int MatBeginAssembly(Mat); 328ed539a5SBarry Smith extern int MatEndAssembly(Mat); 3320563c6bSBarry Smith extern int MatSetOption(Mat,int); 348ed539a5SBarry Smith #define ROW_ORIENTED 1 358ed539a5SBarry Smith #define COLUMN_ORIENTED 2 3620563c6bSBarry Smith #define ROWS_SORTED 4 3720563c6bSBarry Smith #define COLUMNS_SORTED 8 3820563c6bSBarry Smith #define NO_NEW_NONZERO_LOCATIONS 16 3920563c6bSBarry Smith #define YES_NEW_NONZERO_LOCATIONS 32 402eac72dbSBarry Smith 418ed539a5SBarry Smith extern int MatGetValues(Mat,Scalar*,int,int*,int,int*); 428ed539a5SBarry Smith extern int MatGetRow(Mat,int,int *,int **,Scalar**); 438ed539a5SBarry Smith extern int MatRestoreRow(Mat,int,int *,int **,Scalar**); 448ed539a5SBarry Smith extern int MatGetCol(Mat,int,int *,int **,Scalar**); 458ed539a5SBarry Smith extern int MatRestoreCol(Mat,int,int *,int **,Scalar**); 462eac72dbSBarry Smith 478ed539a5SBarry Smith extern int MatMult(Mat,Vec,Vec); 488ed539a5SBarry Smith extern int MatMultAdd(Mat,Vec,Vec,Vec); 498ed539a5SBarry Smith extern int MatMultTrans(Mat,Vec,Vec); 508ed539a5SBarry Smith extern int MatMultTransAdd(Mat,Vec,Vec,Vec); 512eac72dbSBarry Smith 528ed539a5SBarry Smith #define ORDER_NATURAL 0 538ed539a5SBarry Smith #define ORDER_ND 1 548ed539a5SBarry Smith #define ORDER_1WD 2 558ed539a5SBarry Smith #define ORDER_RCM 3 568ed539a5SBarry Smith #define ORDER_QMD 4 578ed539a5SBarry Smith extern int MatGetReordering(Mat,int,IS*,IS*); 582eac72dbSBarry Smith 598ed539a5SBarry Smith extern int MatLUFactor(Mat,IS,IS); 608ed539a5SBarry Smith extern int MatCholeskyFactor(Mat,IS); 618ed539a5SBarry Smith extern int MatLUFactorSymbolic(Mat,IS,IS,Mat*); 62da3a660dSBarry Smith extern int MatILUFactorSymbolic(Mat,IS,IS,int,Mat*); 638ed539a5SBarry Smith extern int MatCholeskyFactorSymbolic(Mat,IS,Mat*); 64da3a660dSBarry Smith extern int MatIncompleteCholeskyFactorSymbolic(Mat,IS,int,Mat*); 6520563c6bSBarry Smith extern int MatLUFactorNumeric(Mat,Mat*); 6620563c6bSBarry Smith extern int MatCholeskyFactorNumeric(Mat,Mat*); 678ed539a5SBarry Smith 688ed539a5SBarry Smith extern int MatSolve(Mat,Vec,Vec); 698ed539a5SBarry Smith extern int MatSolveAdd(Mat,Vec,Vec,Vec); 70*28988994SBarry Smith extern int MatSolveTrans(Mat,Vec,Vec); 71*28988994SBarry Smith extern int MatSolveTransAdd(Mat,Vec,Vec,Vec); 728ed539a5SBarry Smith 738ed539a5SBarry Smith #define SOR_FORWARD_SWEEP 1 748ed539a5SBarry Smith #define SOR_BACKWARD_SWEEP 2 758ed539a5SBarry Smith #define SOR_SYMMETRIC_SWEEP 3 76d6dfbf8fSBarry Smith #define SOR_LOCAL_FORWARD_SWEEP 4 77d6dfbf8fSBarry Smith #define SOR_LOCAL_BACKWARD_SWEEP 8 78d6dfbf8fSBarry Smith #define SOR_LOCAL_SYMMETRIC_SWEEP 12 79d6dfbf8fSBarry Smith #define SOR_ZERO_INITIAL_GUESS 16 80da3a660dSBarry Smith #define SOR_EISENSTAT 32 81d6dfbf8fSBarry Smith extern int MatRelax(Mat,Vec,double,int,double,int,Vec); 828ed539a5SBarry Smith 838ed539a5SBarry Smith extern int MatCopy(Mat,Mat*); 848ed539a5SBarry Smith extern int MatView(Mat,Viewer); 858ed539a5SBarry Smith #include <stdio.h> 868ed539a5SBarry Smith extern int MatPrintMatlab(Mat,FILE*,char *); 878ed539a5SBarry Smith extern int MatNonZeros(Mat,int*); 888ed539a5SBarry Smith extern int MatMemoryUsed(Mat,int*); 898ed539a5SBarry Smith extern int MatGetDiagonal(Mat,Vec); 908ed539a5SBarry Smith extern int MatTranspose(Mat); 918ed539a5SBarry Smith extern int MatScale(Mat,Vec,Vec); 928ed539a5SBarry Smith extern int MatShrink(Mat,int,int*,int,int*); 938ed539a5SBarry Smith extern int MatEqual(Mat,Mat); 9420563c6bSBarry Smith extern int MatScatterBegin(Mat,IS,IS,Mat,IS,IS,InsertMode,MatScatterCtx*); 9520563c6bSBarry Smith extern int MatScatterEnd(Mat,IS,IS,Mat,IS,IS,InsertMode,MatScatterCtx*); 968ed539a5SBarry Smith extern int MatReOrder(Mat,IS,IS); 972eac72dbSBarry Smith 982eac72dbSBarry Smith #define NORM_1 1 992eac72dbSBarry Smith #define NORM_2 2 1002eac72dbSBarry Smith #define NORM_FROBENIUS 3 1012eac72dbSBarry Smith #define NORM_INFINITY 4 1028ed539a5SBarry Smith extern int MatNorm(Mat,int,double *); 1032eac72dbSBarry Smith 1048ed539a5SBarry Smith extern int MatZeroEntries(Mat); 10520563c6bSBarry Smith extern int MatZeroRows(Mat,IS,Scalar*); 10620563c6bSBarry Smith extern int MatZeroColumns(Mat,IS,Scalar*); 1078ed539a5SBarry Smith 1088ed539a5SBarry Smith extern int MatCompress(Mat); 1098ed539a5SBarry Smith extern int MatDestroy(Mat); 1108ed539a5SBarry Smith 1111eb62cbbSBarry Smith extern int MatGetSize(Mat,int*,int*); 1121eb62cbbSBarry Smith extern int MatGetLocalSize(Mat,int*,int*); 113d6dfbf8fSBarry Smith extern int MatGetOwnershipRange(Mat,int*,int*); 1141eb62cbbSBarry Smith 1158ed539a5SBarry Smith extern int MatCreateInitialMatrix(int,int,Mat*); 1168ed539a5SBarry Smith 1178ed539a5SBarry Smith 1188ed539a5SBarry Smith 1191eb62cbbSBarry Smith 1202eac72dbSBarry Smith #endif 1212eac72dbSBarry Smith 1222eac72dbSBarry Smith 123