1*ca71c51bSBarry Smith /* $Id: is.h,v 1.38 1997/09/26 02:22:17 bsmith Exp bsmith $ */ 22eac72dbSBarry Smith 32eac72dbSBarry Smith /* 4f8256253SLois Curfman McInnes An index set is a generalization of a subset of integers. Index sets 5f8256253SLois Curfman McInnes are used for defining scatters and gathers. 62eac72dbSBarry Smith */ 72eac72dbSBarry Smith #if !defined(__IS_PACKAGE) 82eac72dbSBarry Smith #define __IS_PACKAGE 92eac72dbSBarry Smith #include "petsc.h" 102eac72dbSBarry Smith 119e25ed09SBarry Smith #define IS_COOKIE PETSC_COOKIE+2 12f0479e8cSBarry Smith 13f09e8eb9SSatish Balay typedef struct _p_IS* IS; 142eac72dbSBarry Smith 15639f9d9dSBarry Smith /* 16639f9d9dSBarry Smith Default index set data structures that PETSc provides. 17639f9d9dSBarry Smith */ 18639f9d9dSBarry Smith typedef enum {IS_GENERAL=0, IS_STRIDE=1, IS_BLOCK = 2} ISType; 19537820f0SBarry Smith extern int ISCreateGeneral(MPI_Comm,int,int *,IS *); 20537820f0SBarry Smith extern int ISCreateBlock(MPI_Comm,int,int,int *,IS *); 21537820f0SBarry Smith extern int ISCreateStride(MPI_Comm,int,int,int,IS *); 224b0e389bSBarry Smith 234b0e389bSBarry Smith extern int ISDestroy(IS); 244b0e389bSBarry Smith 25aabeff55SBarry Smith extern int ISSetPermutation(IS); 2677c4ece6SBarry Smith extern int ISPermutation(IS,PetscTruth*); 2708480c60SBarry Smith extern int ISSetIdentity(IS); 2877c4ece6SBarry Smith extern int ISIdentity(IS,PetscTruth*); 2908480c60SBarry Smith 30aabeff55SBarry Smith extern int ISGetIndices(IS,int **); 31aabeff55SBarry Smith extern int ISRestoreIndices(IS,int **); 32aabeff55SBarry Smith extern int ISGetSize(IS,int *); 33aabeff55SBarry Smith extern int ISInvertPermutation(IS,IS*); 34aabeff55SBarry Smith extern int ISView(IS,Viewer); 3577c4ece6SBarry Smith extern int ISEqual(IS, IS, PetscTruth *); 360f453280SSatish Balay extern int ISSort(IS); 3777c4ece6SBarry Smith extern int ISSorted(IS, PetscTruth *); 38*ca71c51bSBarry Smith extern int ISDifference(IS,IS,IS*); 39612dd529SBarry Smith 40537820f0SBarry Smith extern int ISBlock(IS,PetscTruth*); 41c16cb8f2SBarry Smith extern int ISBlockGetIndices(IS,int **); 42c16cb8f2SBarry Smith extern int ISBlockRestoreIndices(IS,int **); 43537820f0SBarry Smith extern int ISBlockGetSize(IS,int *); 44537820f0SBarry Smith extern int ISBlockGetBlockSize(IS,int *); 45c16cb8f2SBarry Smith 46537820f0SBarry Smith extern int ISStride(IS,PetscTruth*); 47c16cb8f2SBarry Smith extern int ISStrideGetInfo(IS,int *,int*); 48c16cb8f2SBarry Smith 49d4bb536fSBarry Smith extern int ISDuplicate(IS, IS *); 5056cd22aeSBarry Smith /* --------------------------------------------------------------------------*/ 5156cd22aeSBarry Smith 52639f9d9dSBarry Smith /* 5390f02eecSBarry Smith ISLocalToGlobalMappings are mappings from an arbitrary 5490f02eecSBarry Smith local ordering from 0 to n-1 to a global PETSc ordering 55d4bb536fSBarry Smith used by a vector or matrix. 56d4bb536fSBarry Smith 57d4bb536fSBarry Smith Note: mapping from Local to Global is scalable; but Global 58d4bb536fSBarry Smith to local may not be if the range of global values represented locally 59d4bb536fSBarry Smith is very large. 6090f02eecSBarry Smith */ 61d4bb536fSBarry Smith #define IS_LTOGM_COOKIE PETSC_COOKIE+12 62f09e8eb9SSatish Balay typedef struct _p_ISLocalToGlobalMapping* ISLocalToGlobalMapping; 6390f02eecSBarry Smith 64d4bb536fSBarry Smith extern int ISLocalToGlobalMappingCreate(MPI_Comm,int, int*, ISLocalToGlobalMapping*); 6590f02eecSBarry Smith extern int ISLocalToGlobalMappingDestroy(ISLocalToGlobalMapping); 66d4bb536fSBarry Smith extern int ISLocalToGlobalMappingApply(ISLocalToGlobalMapping,int,int*,int *); 6790f02eecSBarry Smith extern int ISLocalToGlobalMappingApplyIS(ISLocalToGlobalMapping,IS,IS*); 68d4bb536fSBarry Smith typedef enum {IS_GTOLM_MASK,IS_GTOLM_DROP} ISGlobalToLocalMappingType; 69d4bb536fSBarry Smith extern int ISGlobalToLocalMappingApply(ISLocalToGlobalMapping,ISGlobalToLocalMappingType, 70d4bb536fSBarry Smith int,int *,int*,int *); 7190f02eecSBarry Smith 7256cd22aeSBarry Smith /* --------------------------------------------------------------------------*/ 7356cd22aeSBarry Smith 7490f02eecSBarry Smith /* 75639f9d9dSBarry Smith ISColorings are sets of IS's that define a coloring 76639f9d9dSBarry Smith of the underlying indices 77639f9d9dSBarry Smith */ 78f09e8eb9SSatish Balay struct _p_ISColoring { 79639f9d9dSBarry Smith int n; 80639f9d9dSBarry Smith IS *is; 81639f9d9dSBarry Smith MPI_Comm comm; 82639f9d9dSBarry Smith }; 83f09e8eb9SSatish Balay typedef struct _p_ISColoring* ISColoring; 84639f9d9dSBarry Smith 85639f9d9dSBarry Smith extern int ISColoringDestroy(ISColoring); 86639f9d9dSBarry Smith extern int ISColoringView(ISColoring,Viewer); 87639f9d9dSBarry Smith extern int ISColoringCreate(MPI_Comm,int,int*,ISColoring*); 88639f9d9dSBarry Smith 89dbef8a1cSBarry Smith /* --------------------------------------------------------------------------*/ 90dbef8a1cSBarry Smith 91dbef8a1cSBarry Smith /* 92dbef8a1cSBarry Smith ISPartitioning are sets of IS's that define a partioning 93dbef8a1cSBarry Smith of the underlying indices. This is the same as a ISColoring. 94dbef8a1cSBarry Smith */ 95dbef8a1cSBarry Smith #define ISPartitioning ISColoring 960752156aSBarry Smith #define ISPartitioningView ISColoringView 970752156aSBarry Smith #define ISPartitioningCreate ISColoringCreate 980752156aSBarry Smith #define ISPartitioningDestroy ISColoringDestroy 990752156aSBarry Smith extern int ISPartitioningToLocalIS(ISPartitioning,IS*); 100dbef8a1cSBarry Smith 101a2ce50c7SBarry Smith #endif 1027588ac45SBarry Smith 1037588ac45SBarry Smith 104639f9d9dSBarry Smith 105639f9d9dSBarry Smith 106