evsl.h

Go to the documentation of this file.

#ifndef EVSL_H
#define EVSL_H

#include <complex.h>
#include <stdio.h>
#include "struct.h"

/*- - - - - - - - - cheblanNr.c */
int ChebLanNr(double *intv, int maxit, double tol, double *vinit,
              polparams *pol, int *nevOut, double **lamo, double **Wo,
              double **reso, FILE *fstats);

/*- - - - - - - - - cheblanTr.c */
int ChebLanTr(int lanm, int nev, double *intv, int maxit, double tol,
              double *vinit, polparams *pol, int *nev2, double **vals,
              double **W, double **resW, FILE *fstats);

/*- - - - - - - - - chebpoly.c */
//
void set_pol_def(polparams *pol);
//
int find_pol(double *intv, polparams *pol);
//
void free_pol(polparams *pol);

/*- - - - - - - - - chebsi.c */
int ChebSI(int nev, double *intv, int maxit, double tol, double *vinit,
           polparams *pol, int *nevo, double **lamo, double **Yo, double **reso,
           FILE *fstats);

/*- - - - - - - - - - dos_utils.c */
void SetupPolRec(int n, int max_deg, double tol, double lmin, double lmax, BSolDataPol *data);
//
void SetupPolSqrt(int n, int max_deg, double tol, double lmin, double lmax, BSolDataPol *data);
//
void FreeBSolPolData(BSolDataPol *data);
//
void BSolPol(double *b, double *x, void *data);
//
void extrDiagCsr(csrMat *B, double *d);
//
void diagScalCsr(csrMat *A, double *d);

/*- - - - - - - - - lanbounds.c */
int LanBounds(int msteps, double *v, double *lmin, double *lmax);

/*- - - - - - - - - - landos.c */
// Computes the density of states (DOS, or spectral density)
int LanDos(const int nvec, int msteps, const int npts, double *xdos,
           double *ydos, double *neig, const double *const intv);

//*- - - - - - - - -  landosG.c - Generalized lanDOS */
int LanDosG(const int nvec, int msteps, const int npts, double *xdos, double *ydos,
            double *neig, const double *const intv);

/*- - - - - - - - - lanTrbounds.c */
int LanTrbounds(int lanm, int maxit, double tol, double *vinit, int bndtype,
                double *lammin, double *lammax, FILE *fstats);

/*- - - -- - - - - - misc_la.c */
int scalEigVec(int n, int nev, double *Y, double* sqrtdiag);


/*- - - - - - - - - ratfilter.c */
//
void set_ratf_def(ratparams *rat);
//
int find_ratf(double *intv, ratparams *rat);
//
int set_ratf_solfunc(ratparams *rat, csrMat *A, csrMat *B, SolFuncC *funcs,
                     void **data);
//
void free_rat(ratparams *rat);

/*- - - - - - - - - ratlanNr.c */
//
int RatLanNr(double *intv, int maxit, double tol, double *vinit, ratparams *rat,
             int *nevOut, double **lamo, double **Wo, double **reso,
             FILE *fstats);

/*- - - - - - - - - ratlanTr.c */
//
int RatLanTr(int lanm, int nev, double *intv, int maxit, double tol,
             double *vinit, ratparams *rat, int *nev2, double **vals,
             double **W, double **resW, FILE *fstats);

/*- - - - - - - - - spmat.c */
void csr_copy(csrMat *A, csrMat *B, int allocB);
// convert a COO matrix to a CSR matrix
int cooMat_to_csrMat(int cooidx, cooMat *coo, csrMat *csr);
// free a CSR
void free_csr(csrMat *csr);
// free a COO
void free_coo(cooMat *coo);
// matvec y = A*x
int matvec(char trans, csrMat *A, double *x, double *y);

/*- - - - - - - - - evsl.c */
/* set an external matvec function */
int SetAMatvec(int n, MVFunc func, void *data);
int SetBMatvec(int n, MVFunc func, void *data);
/* unset an external matvec function */
int UnsetAMatvec();
int UnsetBMatvec();
int SetAMatrix(csrMat *A);
int SetBMatrix(csrMat *B);
int SetBSol(SolFuncR func, void *data);
int SetLTSol(SolFuncR func, void *data);
int SetASigmaBSol(ratparams *rat, SolFuncC *func, SolFuncC allf, void **data);
int SetStdEig();
int SetGenEig();
/* start EVSL */
int EVSLStart();
/* finalize EVSL */
int EVSLFinish();
void SetDiagScal(double *ds);

/*- - - - - - - - - spslicer.c */
//
int spslicer(double *sli, double *mu, int Mdeg, double *intv, int n_int,
             int npts);
//
int kpmdos(int Mdeg, int damping, int nvec, double *ab, double *mu,
           double *ecnt);

/*- - - - - - - - - - spslicer2.c */
void spslicer2(double *xi, double *yi, int n_int, int npts, double *sli);

/*- - - - - - - - - timing.c */
//
double evsl_timer();

/*- - - - - - - - - vect.c */
// generate a random vector
void rand_double(int n, double *v);
// generate a normally distributed random vector
void randn_double(int n, double *v);
// sort a vector
void sort_double(int n, double *v, int *ind);
//
void linspace(double a, double b, int num, double *arr);

void StatsPrint(FILE *fstats);

void StatsReset();

#endif