struct.h

Go to the documentation of this file.

#ifndef STRUCT_H
#define STRUCT_H

#include <complex.h>
#include <stddef.h>
typedef struct _cooMat {
  int nrows,  
      ncols,  
      nnz,    
      *ir,    
      *jc;    
  double *vv; 
} cooMat;


typedef struct _csrMat {
  int owndata, 
      nrows,   
      ncols,   
      *ia,     
      *ja;     
  double *a;   
} csrMat;

typedef struct _polparams {
  int     max_deg;     
  int     min_deg ;    
  int     damping;     
  double  thresh_ext;  
  double  thresh_int;  
  double  tol;         
  double  intvtol;     
  int     type;       
  double *mu;         
  double  cc;         
  double  dd;         
  double  gam;        
  double  bar;        
  int deg ;           
} polparams;

typedef void (*SolFuncC)(int n, double *br, double *bz, double *xr, double *xz, void *data);

typedef void (*SolFuncR)(double *b, double *x, void *data);

typedef void (*MVFunc)(double *x, double *y, void *data);

typedef struct _EVSLASIGMABSol {
  SolFuncC func;       
  void *data;          
} EVSLASIGMABSol;

typedef struct _ratparams {
  /*  */
  int num;            
  int pw;             
  int method;         
  double beta;        
  double aa;          
  double bb;          
  double bar;         
  int *mulp;          
  int pow;            
  complex double *omega; 
  complex double *zk;    
  EVSLASIGMABSol *ASIGBsol; 
} ratparams;


typedef struct _EVSLMatvec {
  MVFunc func;         
  void *data;          
} EVSLMatvec;

typedef struct _EVSLBSol {
  SolFuncR func;       
  void *data;          
} EVSLBSol;

typedef struct _EVSLLTSol {
  SolFuncR func;       
  void *data;          
} EVSLLTSol;


typedef struct _evsldata {
  int n;                    
  int ifGenEv;              
  EVSLMatvec *Amv;          
  EVSLMatvec *Bmv;          
  EVSLBSol *Bsol;           
  EVSLLTSol *LTsol;         
  double *ds;               
} evslData;

/*
 * Define a struct for L-S polynomial approximation to a matrix function
 */
typedef struct _BSolDataPol {
  int max_deg;    
  int deg;        
  double intv[2]; 
  double cc, dd;  
  double tol;     
  double *mu;     
  double *wk;     
} BSolDataPol;


/* global variable: evslData */
extern evslData evsldata;

typedef struct _evslstat {
  /* timing [level-1 funcs] */
  double t_setBsv;
  double t_setASigBsv;
  double t_iter;
  /* timing [level-2 funcs] */
  double t_mvA;
  double t_mvB;
  double t_svB;
  double t_svLT;
  double t_svASigB;
  double t_reorth;
  double t_eig;
  double t_blas;
  double t_ritz;
  double t_polAv;
  double t_ratAv;
  double t_sth;
  size_t n_mvA;
  size_t n_mvB;
  size_t n_svB;
  size_t n_svLT;
  size_t n_svASigB;
  size_t n_polAv;
  size_t n_ratAv;
  /* memory */
  size_t alloced;
  size_t alloced_total;
  size_t alloced_max;
} evslStat;

/* global variable: pevsl_stats */
extern evslStat evslstat;
#endif