/* $Id: bjacobi.h,v 1.24 2001/03/22 20:31:04 bsmith Exp $ */

#if !defined(__BJACOBI_H)
#define __BJACOBI_H
/*
    Private data for block Jacobi and block Gauss-Seidel preconditioner.
*/
#include "petscsles.h"
#include "src/sles/pc/pcimpl.h"

/*
       This data is general for all implementations
*/
typedef struct {
  int        n,n_local;        /* number of blocks (global, local) */
  int        first_local;       /* number of first block on processor */
  PetscTruth use_true_local;    /* use block from true matrix, not preconditioner matrix for local MatMult() */
  SLES       *sles;             /* SLES contexts for blocks */
  void       *data;             /* implementation-specific data */
  PetscTruth same_local_solves; /* flag indicating whether all local solvers are same (used for PCView()) */
  int        *l_lens;           /* lens of each block */
  int        *g_lens;
  Mat        tp_mat,tp_pmat;    /* diagonal block of matrix for this processor */
} PC_BJacobi;

/*
       This data is specific for certain implementations
*/

/*  This is for multiple blocks per processor */

typedef struct {
  Vec              *x,*y;             /* work vectors for solves on each block */
  int              *starts;           /* starting point of each block */
  Mat              *mat,*pmat;        /* submatrices for each block */
  IS               *is;               /* for gathering the submatrices */
} PC_BJacobi_Multiblock;

/*  This is for a single block per processor */
typedef struct {
  Vec  x,y;
} PC_BJacobi_Singleblock;

#endif