ilupackmacros.h

数学算法的实现库。可以实现常见的线性计算。

#define QQSORTS      Cqqsorts
#define QQSORTS2     Cqqsorts2
#define CLEAR        Cclear
#define QSPLIT       Cqsplit
#define QSPLIT2      Cqsplit2
#define QSEPARATE    Cqseparate
#define INDSET       Cindset
#define PQPERM       CPQpermF
#define PPERM        CPpermF
#define PPPERM       CPPpermF
#define WDIAG        CwDiagRF
#define WPDIAG       CwPDiagRF
#define SPDWDIAG     CHPDwDiagRF
#define QSORTR2I     sqsortr2i
#define SWAPM        sswapm
#define BISINIT      Cbisinit
#define STOPBIS      Cstopbis
#define TIDYCG       Ctidycg
#define MGSRO        Cmgsro
#define GIVENS       Cgivens
#define BRKDN        Cbrkdn

#define CSRCSC       Ccsrcsc
#define CSRCSC2      Ccsrcsc2
#define READMTC      Creadmtc
#define WRITEMTC     Cwritemtc
#define READVECTORS  Creadvectors
#define WRITEVECTORS Cwritevectors
#define PCG          Cpcg
#define BCG          Cbcg
#define SBCGS        CSYMbcg
#define SBCG         CHERbcg
#define SQMR         CHERqmr
#define SQMRS        CSYMqmr
#define GMRES        Cgmres
#define FGMRES       Cfgmres

#define DISTDOT      cdotc
#define DISTDOTU     cdotu
#define DISTDOT2     cdotc2
#define COPY         ccopy
#define AXPY         caxpy
#define ROTG         crotg
#define ROT          crot
#define NRM          scnrm2
#define I_AMAX       icamax
#define ASUM         scasum
#define SCAL         cscal

#define MC64I        mc64ic
#define MC64A        mc64ac
#define MC64IR       mc64is
#define MC64AR       mc64as

#define FABS(A)      sqrt((double)((A).r*(A).r+(A).i*(A).i))
#define ABS(A)       cabs(A)
#define FLOATABS(A)  abs(A)
#define FLOATABSNAME abs
#define SQRT(A)      sqrt(A)
#define LOG(A)       log(A)
#define LOG10(A)     log10(A)
#define IREAL(A)     real(A)
#define ABSFNAME     cabs
#define SQRTFNAME    sqrt
#define LOGFNAME     log
#define LOGTENFNAME  log10
#define CTOD(A)      real(A)
#define CONJ(A)      conjg(A)
#define FNULL        cmplx(0.0e0,0.0e0)
#define FONE         cmplx(1.0e0,0.0e0)
#define RONE         1.0e0
#define RTWO         2.0e0
#define RFOUR        4.0e0
#define RZERO        0.0e0
#define REAL_MAX     1.0e38
#define SIGNUM       signum
#define SIGNUM2      signum2
#define TWICE        2



#else
#define GETEPS       dgeteps
#define _REAL_MAX_   _D_REAL_MAX_
#define REALS        doubleprecision
#define FLOAT        doublecomplex
#define CSRMAT       Zmat
#define AMGLEVELMAT  ZAMGlevelmat
#define ILUPACKPARAM ZILUPACKparam
#define GNLAGGRESSIVEDROPPING ZGNLaggressivedropping 
#define SPDAGGRESSIVEDROPPING ZHPDaggressivedropping 

#define SPTRF        zsptrf
#define SPTRS        zsptrs

#define ROWSCALE     Zrowscale
#define COLSCALE     Zcolscale
#define SPDSCALE     ZHPDscale
#define SYMSCALE     ZSYMscale
#define CNRMS        Zcnrms
#define RNRMS        Zrnrms
#define DIAMUA       ZDiagMultA
#define AMUDIA       ZAMultDiag

#define CSRMATVEC    ZGNLmatvec
#define CSRMATTVEC   ZGNLmattvec
#define CSRMATHVEC   ZGNLmathvec
#define SYMMATVEC    ZSYMmatvec
#define SSMMATVEC    ZSSMmatvec
#define HERMATVEC    ZHERmatvec
#define SHRMATVEC    ZSHRmatvec
#define SETUPGRAPH   Zsetupgraph
#define SETUPGRAPH_EPSILON   Zsetupgraph_epsilon

#define AMGDELETE       ZGNLAMGdelete
#define SPDAMGDELETE    ZHPDAMGdelete

#define AMGSETUP        ZGNLAMGfactor
#define SPDAMGSETUP     ZHPDAMGfactor
#define AMGFACTOR       ZGNLAMGfactor
#define SPDAMGFACTOR    ZHPDAMGfactor
#define AMGEXTRACT      ZGNLAMGextract
#define SPDAMGEXTRACT   ZHERAMGextract
#define SYMAMGEXTRACT   ZHERAMGextract
#define AMGSOL          ZGNLAMGsol
#define AMGDLSOL        ZGNLAMGdlsol
#define AMGDUSOL        ZGNLAMGdusol
#define AMGLSOL         ZGNLAMGlsol
#define AMGUSOL         ZGNLAMGusol
#define AMGTDLSOL       ZGNLAMGtdlsol
#define AMGTDUSOL       ZGNLAMGtdusol
#define AMGTLSOL        ZGNLAMGtlsol
#define AMGTUSOL        ZGNLAMGtusol
#define AMGTSOL         ZGNLAMGtsol
#define SPDAMGSOL       ZHPDAMGsol
#define AMGDRIVER       ZGNLAMGdriver
#define AMGSOLVER       ZGNLAMGsolver
#define SPDAMGSOLVER    ZHPDAMGsolver
#define SYMAMGSOLVERS    ZSYMAMGsolver
#define SYMAMGSOLVER    ZHERAMGsolver
#define AMGINIT         ZGNLAMGinit
#define AMGGETPARAMS    ZGNLAMGgetparams
#define AMGSETPARAMS    ZGNLAMGsetparams
#define SPDAMGINIT      ZHPDAMGinit
#define SPDAMGGETPARAMS ZHPDAMGgetparams
#define SPDAMGSETPARAMS ZHPDAMGsetparams
#define SYMAMGINITS      ZSYMAMGinit
#define SYMAMGGETPARAMSS ZSYMAMGgetparams
#define SYMAMGSETPARAMSS ZSYMAMGsetparams
#define SYMAMGINIT      ZHERAMGinit
#define SYMAMGGETPARAMS ZHERAMGgetparams
#define SYMAMGSETPARAMS ZHERAMGsetparams

#define LDLP         ZHPDldlp
#define LDLPSOL      ZHPDldlpsol
#define LUPQ         ZGNLlupq
#define LUPQSOL      ZGNLlupqsol
#define LUPQLSOL     ZGNLlupqlsol
#define LUPQUSOL     ZGNLlupqusol
#define LUPQDLSOL    ZGNLlupqdlsol
#define LUPQDUSOL    ZGNLlupqdusol
#define LUPQTSOL     ZGNLlupqtsol
#define LUPQTLSOL    ZGNLlupqtlsol
#define LUPQTUSOL    ZGNLlupqtusol
#define LUPQTDLSOL   ZGNLlupqtdlsol
#define LUPQTDUSOL   ZGNLlupqtdusol


#define ILDLC        ZHERildlc
#define ILDLCSOL     ZHERildlcsol
#define ILUC         ZGNLiluc
#define PILUC        ZGNLpiluc
#define MPILUC       ZGNLmpiluc
#define PILDLC       ZHPDpiluc
#define MPILDLC      ZHPDmpiluc
#define ILUCSOL      ZGNLilucsol
#define ILUCTSOL     ZGNLiluctsol
#define ILUCDLSOL    ZGNLilucdlsol
#define ILUCTDLSOL   ZGNLiluctdlsol
#define ILUCDUSOL    ZGNLilucdusol
#define ILUCTDUSOL   ZGNLiluctdusol
#define ILUCLSOL     ZGNLiluclsol
#define ILUCTLSOL    ZGNLiluctlsol
#define ILUCUSOL     ZGNLilucusol
#define ILUCTUSOL    ZGNLiluctusol

#define PILUCDLSOL    ZGNLpilucdlsol
#define PILUCTDLSOL   ZGNLpiluctdlsol
#define PILUCDUSOL    ZGNLpilucdusol
#define PILUCTDUSOL   ZGNLpiluctdusol
#define PILUCLSOL     ZGNLpiluclsol
#define PILUCTLSOL    ZGNLpiluctlsol
#define PILUCUSOL     ZGNLpilucusol
#define PILUCTUSOL    ZGNLpiluctusol

#define ILUTP        ZGNLilutp
#define ILUT         ZGNLilut
#define LUSOL        ZGNLlusol
#define LUTSOL       ZGNLlutsol
#define LUUSOL       ZGNLluusol
#define LUTUSOL      ZGNLlutusol
#define LULSOL       ZGNLlulsol
#define LUTLSOL      ZGNLlutlsol
#define LUDUSOL      ZGNLludusol
#define LUTDUSOL     ZGNLlutdusol
#define LUDLSOL      ZGNLludlsol
#define LUTDLSOL     ZGNLlutdlsol

#define CPERM        Zcperm
#define RPERM        Zrperm

#define PERMNULL     ZGNLperm_null
#define PERMND       ZGNLperm_nd
#define PERMRCM      ZGNLperm_rcm
#define PERMAMF      ZGNLperm_amf
#define PERMAMD      ZGNLperm_amd
#define PERMMMD      ZGNLperm_mmd
#define PERMMETISE   ZGNLperm_metis_e
#define PERMMETISN   ZGNLperm_metis_n
#define PERMPQ       ZGNLperm_pq
#define PERMFC       ZGNLperm_fc
#define INDFC        Zindfc
#define SYMPERMFC    ZSYMperm_fc
#define SYMINDFC     ZSYMindfc
#define SPDPERMFC    ZHPDperm_fc
#define SPDINDFC     ZHPDindfc
#define PERMP        ZGNLperm_p
#define PERMINDSET   ZGNLperm_indset

#define PERMMWMRCM    ZGNLperm_mwm_rcm
#define PERMMWMMMD    ZGNLperm_mwm_mmd
#define PERMMWMAMF    ZGNLperm_mwm_amf
#define PERMMWMAMD    ZGNLperm_mwm_amd
#define PERMMWMMETISE ZGNLperm_mwm_metis_e
#define PERMMWMMETISN ZGNLperm_mwm_metis_n

#define SYMPERMMWMRCM    ZHERperm_mwm_rcm
#define SYMPERMMWMMMD    ZHERperm_mwm_mmd
#define SYMPERMMWMAMF    ZHERperm_mwm_amf
#define SYMPERMMWMAMD    ZHERperm_mwm_amd
#define SYMPERMMWMMETISE ZHERperm_mwm_metis_e
#define SYMPERMMWMMETISN ZHERperm_mwm_metis_n

#define PERMMC64RCM    ZGNLperm_mc64_rcm
#define PERMMC64MMD    ZGNLperm_mc64_mmd
#define PERMMC64AMF    ZGNLperm_mc64_amf
#define PERMMC64AMD    ZGNLperm_mc64_amd
#define PERMMC64METISE ZGNLperm_mc64_metis_e
#define PERMMC64METISN ZGNLperm_mc64_metis_n

#define SPDPERMNULL     ZHPDperm_null
#define SPDPERMND       ZHPDperm_nd
#define SPDPERMRCM      ZHPDperm_rcm
#define SPDPERMAMF      ZHPDperm_amf
#define SPDPERMAMD      ZHPDperm_amd
#define SPDPERMMMD      ZHPDperm_mmd
#define SPDPERMMETISE   ZHPDperm_metis_e
#define SPDPERMMETISN   ZHPDperm_metis_n
#define SPDPERMPP       ZHPDperm_pp
#define SPDPERMINDSET   ZHPDperm_indset

#define SYMPERMNULL     ZHERperm_null
#define SYMPERMND       ZHERperm_nd
#define SYMPERMRCM      ZHERperm_rcm
#define SYMPERMAMF      ZHERperm_amf
#define SYMPERMAMD      ZHERperm_amd
#define SYMPERMMMD      ZHERperm_mmd
#define SYMPERMMETISE   ZHERperm_metis_e
#define SYMPERMMETISN   ZHERperm_metis_n
#define SYMPERMPP       ZHERperm_pp
#define SYMPERMINDSET   ZHERperm_indset

#define SPARTRAN     Zspartran
#define QQSORT       Zqqsort
#define QSORT        Zqsort
#define QQSORT2      Zqqsort2
#define QSORT2       Zqsort2
#define QQSORTS      Zqqsorts
#define QQSORTS2     Zqqsorts2
#define CLEAR        Zclear
#define QSPLIT       Zqsplit
#define QSPLIT2      Zqsplit2
#define QSEPARATE    Zqseparate
#define INDSET       Zindset
#define PQPERM       ZPQpermF
#define PPERM        ZPpermF
#define PPPERM       ZPPpermF
#define WDIAG        ZwDiagRF
#define WPDIAG       ZwPDiagRF
#define SPDWDIAG     ZHPDwDiagRF
#define QSORTR2I     dqsortr2i
#define SWAPM        dswapm
#define BISINIT      Zbisinit
#define STOPBIS      Zstopbis
#define TIDYCG       Ztidycg
#define MGSRO        Zmgsro
#define GIVENS       Zgivens
#define BRKDN        Zbrkdn

#define CSRCSC       Zcsrcsc
#define CSRCSC2      Zcsrcsc2
#define READMTC      Zreadmtc
#define WRITEMTC     Zwritemtc
#define READVECTORS  Zreadvectors
#define WRITEVECTORS Zwritevectors
#define PCG          Zpcg
#define BCG          Zbcg
#define SBCGS        ZSYMbcg
#define SBCG         ZHERbcg
#define SQMR         ZHERqmr
#define SQMRS        ZSYMqmr
#define GMRES        Zgmres
#define FGMRES       Zfgmres

#define DISTDOT      zdotc
#define DISTDOTU     zdotu
#define DISTDOT2     zdotc2
#define COPY         zcopy
#define AXPY         zaxpy
#define ROTG         zrotg
#define ROT          zrot
#define NRM          dznrm2
#define I_AMAX       izamax
#define ASUM         dzasum
#define SCAL         zscal

#define MC64I        mc64iz
#define MC64A        mc64az
#define MC64IR       mc64id
#define MC64AR       mc64ad

#define FABS(A)      sqrt((double)((A).r*(A).r+(A).i*(A).i))
#define ABS(A)       cdabs(A)
#define FLOATABS(A)  dabs(A)
#define FLOATABSNAME dabs
#define SQRT(A)      dsqrt(A)
#define LOG(A)       dlog(A)
#define LOG10(A)     dlog10(A)
#define IREAL(A)     dble(A)
#define ABSFNAME     cdabs
#define SQRTFNAME    dsqrt
#define LOGFNAME     dlog
#define LOGTENFNAME  dlog10
#define CTOD(A)      dreal(A)
#define CONJ(A)      dconjg(A)
#define FNULL        dcmplx(0.0d0,0.0d0)
#define FONE         dcmplx(1.0d0,0.0d0)
#define RONE         1.0d0
#define RTWO         2.0d0
#define RFOUR        4.0d0
#define RZERO        0.0d0
#define REAL_MAX     1.0d300
#define SIGNUM       signum
#define SIGNUM2      signum2
#define TWICE        2
#endif


#define MALLOC        MAlloc
#define REALLOC       ReAlloc
#define MC64THRESHOLD 1.0e-2
#define QQSORTI       qqsorti