fasp_block.h

Go to the documentation of this file.

 
#include "fasp.h"
 
#ifndef __FASPBLOCK_HEADER__ /*-- allow multiple inclusions --*/
#define __FASPBLOCK_HEADER__ 
/*---------------------------*/
/*---   Data structures   ---*/
/*---------------------------*/
 
typedef struct dBSRmat {
 
    INT ROW;
 
    INT COL;
 
    INT NNZ;
 
    INT nb; // NOTE: for the moment, allow nb*nb full block
 
    INT storage_manner; // 0: row-major order, 1: column-major order
 
    REAL* val;
 
    INT* IA;
 
    INT* JA;
 
} dBSRmat; 
typedef struct dBLCmat {
 
    INT brow;
 
    INT bcol;
 
    dCSRmat** blocks;
 
} dBLCmat; 
typedef struct iBLCmat {
 
    INT brow;
 
    INT bcol;
 
    iCSRmat** blocks;
 
} iBLCmat; 
typedef struct block_dvector {
 
    INT brow;
 
    dvector** blocks;
 
} block_dvector; 
typedef struct block_ivector {
 
    INT brow;
 
    ivector** blocks;
 
} block_ivector; 
/*---------------------------*/
/*--- Parameter structures --*/
/*---------------------------*/
 
typedef struct {
 
    INT max_levels;
 
    INT num_levels;
 
    dBSRmat A;
 
    dBSRmat R;
 
    dBSRmat P;
 
    dvector b;
 
    dvector x;
 
    dvector diaginv;
 
    dCSRmat Ac;
 
    void* Numeric;
 
    Pardiso_data pdata;
 
    dCSRmat PP;
 
    AMG_data* mglP;
 
    dCSRmat TT;
 
    AMG_data* mglT;
 
    dBSRmat PT;
 
    REAL* pw;
 
    dBSRmat SS;
 
    REAL* sw;
 
    dvector diaginv_SS;
 
    ILU_data PP_LU;
 
    ivector cfmark;
 
    INT ILU_levels;
 
    ILU_data LU;
 
    INT near_kernel_dim;
 
    REAL** near_kernel_basis;
 
    //-----------------------------------------
    // extra near kernal space for extra solve
 
    dCSRmat* A_nk;
 
    dCSRmat* P_nk;
 
    dCSRmat* R_nk;
    //-----------------------------------------
 
    dvector w;
 
    Mumps_data mumps;
 
} AMG_data_bsr; 
typedef struct {
 
    INT nb;
 
    dvector diag;
 
} precond_diag_bsr; 
typedef struct {
 
    SHORT AMG_type;
 
    SHORT print_level;
 
    INT maxit;
 
    INT max_levels;
 
    REAL tol;
 
    SHORT cycle_type;
 
    SHORT smoother;
 
    SHORT smooth_order;
 
    SHORT presmooth_iter;
 
    SHORT postsmooth_iter;
 
    SHORT coarsening_type;
 
    REAL relaxation;
 
    SHORT coarse_solver;
 
    SHORT coarse_scaling;
 
    SHORT amli_degree;
 
    REAL* amli_coef;
 
    REAL tentative_smooth;
 
    SHORT nl_amli_krylov_type;
 
    AMG_data_bsr* mgl_data;
 
    AMG_data* pres_mgl_data;
 
    ILU_data* LU;
 
    dBSRmat* A;
 
    // extra near kernal space
 
    dCSRmat* A_nk;
 
    dCSRmat* P_nk;
 
    dCSRmat* R_nk;
 
    dvector r;
 
    REAL* w;
 
} precond_data_bsr; 
typedef struct {
 
    /*-------------------------------------*/
    /* Basic data for block preconditioner */
    /*-------------------------------------*/
    dBLCmat* Ablc;   
    dCSRmat* A_diag; 
    dvector r;       
    /*------------------------------*/
    /* Data for the diagonal blocks */
    /*------------------------------*/
 
    /*--- solve by direct solver ---*/
    void** LU_diag; 
    /*--- solve by AMG ---*/
    AMG_data** mgl;      
    AMG_param* amgparam; 
} precond_data_blc;      
typedef struct {
 
    INT NumLayers;        
    dBLCmat* A;           
    dBLCmat* Ai;          
    dCSRmat* local_A;     
    void**   local_LU;    
    ivector* local_index; 
    // temprary work spaces
    dvector r; 
    REAL*   w; 
} precond_data_sweeping; 
#endif                   /* end if for __FASPBLOCK_HEADER__ */
 
/*---------------------------------*/
/*--        End of File          --*/
/*---------------------------------*/