ilupackmacros.h

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

#ifdef _DOUBLE_REAL_
#define GETEPS       dgeteps
#define _REAL_MAX_   _D_REAL_MAX_
#define REALS        doubleprecision
#define FLOAT        doubleprecision
#define CSRMAT       Dmat
#define AMGLEVELMAT  DAMGlevelmat
#define ILUPACKPARAM DILUPACKparam
#define GNLAGGRESSIVEDROPPING DGNLaggressivedropping 
#define SPDAGGRESSIVEDROPPING DSPDaggressivedropping 
#define SPTRF        dsptrf
#define SPTRS        dsptrs

#define ROWSCALE     Drowscale
#define COLSCALE     Dcolscale
#define SPDSCALE     DSPDscale
#define SYMSCALE     DSYMscale
#define CNRMS        Dcnrms
#define RNRMS        Drnrms
#define DIAMUA       DDiagMultA
#define AMUDIA       DAMultDiag

#define CSRMATVEC    DGNLmatvec
#define CSRMATTVEC   DGNLmattvec
#define CSRMATHVEC   DGNLmathvec
#define SYMMATVEC    DSYMmatvec
#define SSMMATVEC    DSSMmatvec
#define HERMATVEC    DSYMmatvec
#define SETUPGRAPH   Dsetupgraph
#define SETUPGRAPH_EPSILON   Dsetupgraph_epsilon

#define AMGDELETE       DGNLAMGdelete
#define SPDAMGDELETE    DSPDAMGdelete

#define AMGSETUP        DGNLAMGfactor
#define SPDAMGSETUP     DSPDAMGfactor
#define AMGFACTOR       DGNLAMGfactor
#define SPDAMGFACTOR    DSPDAMGfactor
#define AMGEXTRACT      DGNLAMGextract
#define SPDAMGEXTRACT   DSYMAMGextract
#define SYMAMGEXTRACT   DSYMAMGextract
#define AMGSOL          DGNLAMGsol
#define AMGTSOL         DGNLAMGtsol
#define AMGDLSOL        DGNLAMGdlsol
#define AMGDUSOL        DGNLAMGdusol
#define AMGLSOL         DGNLAMGlsol
#define AMGUSOL         DGNLAMGusol
#define AMGTDLSOL       DGNLAMGtdlsol
#define AMGTDUSOL       DGNLAMGtdusol
#define AMGTLSOL        DGNLAMGtlsol
#define AMGTUSOL        DGNLAMGtusol
#define AMGDRIVER       DGNLAMGdriver
#define AMGSOLVER       DGNLAMGsolver
#define SPDAMGSOLVER    DSPDAMGsolver
#define SYMAMGSOLVER    DSYMAMGsolver
#define AMGINIT         DGNLAMGinit
#define AMGGETPARAMS    DGNLAMGgetparams
#define AMGSETPARAMS    DGNLAMGsetparams
#define SPDAMGSOL       DSPDAMGsol
#define SPDAMGINIT      DSPDAMGinit
#define SPDAMGGETPARAMS DSPDAMGgetparams
#define SPDAMGSETPARAMS DSPDAMGsetparams
#define SYMAMGINIT      DSYMAMGinit
#define SYMAMGGETPARAMS DSYMAMGgetparams
#define SYMAMGSETPARAMS DSYMAMGsetparams

#define LDLP         DSPDldlp
#define LDLPSOL      DSPDldlpsol
#define LUPQ         DGNLlupq
#define LUPQSOL      DGNLlupqsol
#define LUPQLSOL     DGNLlupqlsol
#define LUPQUSOL     DGNLlupqusol
#define LUPQDLSOL    DGNLlupqdlsol
#define LUPQDUSOL    DGNLlupqdusol
#define LUPQTSOL     DGNLlupqtsol
#define LUPQTLSOL    DGNLlupqtlsol
#define LUPQTUSOL    DGNLlupqtusol
#define LUPQTDLSOL   DGNLlupqtdlsol
#define LUPQTDUSOL   DGNLlupqtdusol


#define ILDLC        DSYMildlc
#define ILDLCSOL     DSYMildlcsol
#define ILUC         DGNLiluc
#define PILUC        DGNLpiluc
#define MPILUC       DGNLmpiluc
#define PILDLC       DSPDpiluc
#define MPILDLC      DSPDmpiluc
#define ILUCSOL      DGNLilucsol
#define ILUCTSOL     DGNLiluctsol
#define ILUCDLSOL    DGNLilucdlsol
#define ILUCTDLSOL   DGNLiluctdlsol
#define ILUCDUSOL    DGNLilucdusol
#define ILUCTDUSOL   DGNLiluctdusol
#define ILUCLSOL     DGNLiluclsol
#define ILUCTLSOL    DGNLiluctlsol
#define ILUCUSOL     DGNLilucusol
#define ILUCTUSOL    DGNLiluctusol

#define PILUCDLSOL    DGNLpilucdlsol
#define PILUCTDLSOL   DGNLpiluctdlsol
#define PILUCDUSOL    DGNLpilucdusol
#define PILUCTDUSOL   DGNLpiluctdusol
#define PILUCLSOL     DGNLpiluclsol
#define PILUCTLSOL    DGNLpiluctlsol
#define PILUCUSOL     DGNLpilucusol
#define PILUCTUSOL    DGNLpiluctusol

#define ILUTP        DGNLilutp
#define ILUT         DGNLilut
#define LUSOL        DGNLlusol
#define LUTSOL       DGNLlutsol
#define LULSOL       DGNLlulsol
#define LUTLSOL      DGNLlutlsol
#define LUUSOL       DGNLluusol
#define LUTUSOL      DGNLlutusol
#define LUDLSOL      DGNLludlsol
#define LUTDLSOL     DGNLlutdlsol
#define LUDUSOL      DGNLludusol
#define LUTDUSOL     DGNLlutdusol

#define CPERM        Dcperm
#define RPERM        Drperm

#define PERMNULL     DGNLperm_null
#define PERMND       DGNLperm_nd
#define PERMRCM      DGNLperm_rcm
#define PERMAMF      DGNLperm_amf
#define PERMAMD      DGNLperm_amd
#define PERMMMD      DGNLperm_mmd
#define PERMMETISE   DGNLperm_metis_e
#define PERMMETISN   DGNLperm_metis_n
#define PERMPQ       DGNLperm_pq
#define PERMFC       DGNLperm_fc
#define INDFC        Dindfc
#define SYMPERMFC    DSYMperm_fc
#define SYMINDFC     DSYMindfc
#define SPDPERMFC    DSPDperm_fc
#define SPDINDFC     DSPDindfc
#define PERMP        DGNLperm_p
#define PERMINDSET   DGNLperm_indset

#define SPDPERMNULL     DSPDperm_null
#define SPDPERMND       DSPDperm_nd
#define SPDPERMRCM      DSPDperm_rcm
#define SPDPERMAMF      DSPDperm_amf
#define SPDPERMAMD      DSPDperm_amd
#define SPDPERMMMD      DSPDperm_mmd
#define SPDPERMMETISE   DSPDperm_metis_e
#define SPDPERMMETISN   DSPDperm_metis_n
#define SPDPERMPP       DSPDperm_pp
#define SPDPERMINDSET   DSPDperm_indset

#define SYMPERMNULL     DSYMperm_null
#define SYMPERMND       DSYMperm_nd
#define SYMPERMRCM      DSYMperm_rcm
#define SYMPERMAMF      DSYMperm_amf
#define SYMPERMAMD      DSYMperm_amd
#define SYMPERMMMD      DSYMperm_mmd
#define SYMPERMMETISE   DSYMperm_metis_e
#define SYMPERMMETISN   DSYMperm_metis_n
#define SYMPERMINDSET   DSYMperm_indset

#define PERMMWMRCM    DGNLperm_mwm_rcm
#define PERMMWMMMD    DGNLperm_mwm_mmd
#define PERMMWMAMF    DGNLperm_mwm_amf
#define PERMMWMAMD    DGNLperm_mwm_amd
#define PERMMWMMETISE DGNLperm_mwm_metis_e
#define PERMMWMMETISN DGNLperm_mwm_metis_n

#define SYMPERMMWMRCM    DSYMperm_mwm_rcm
#define SYMPERMMWMMMD    DSYMperm_mwm_mmd
#define SYMPERMMWMAMF    DSYMperm_mwm_amf
#define SYMPERMMWMAMD    DSYMperm_mwm_amd
#define SYMPERMMWMMETISE DSYMperm_mwm_metis_e
#define SYMPERMMWMMETISN DSYMperm_mwm_metis_n

#define PERMMC64RCM    DGNLperm_mc64_rcm
#define PERMMC64MMD    DGNLperm_mc64_mmd
#define PERMMC64AMF    DGNLperm_mc64_amf
#define PERMMC64AMD    DGNLperm_mc64_amd
#define PERMMC64METISE DGNLperm_mc64_metis_e
#define PERMMC64METISN DGNLperm_mc64_metis_n

#define SPARTRAN     Dspartran
#define QSORT        Dqsort
#define QSORT2       Dqsort2
#define QQSORT       Dqqsort
#define QQSORT2      Dqqsort2
#define QQSORTS      Dqqsorts
#define QQSORTS2     Dqqsorts2
#define QSPLIT       Dqsplit
#define QSPLIT2      Dqsplit2
#define QSEPARATE    Dqseparate
#define CLEAR        Dclear
#define INDSET       Dindset
#define PQPERM       DPQpermF
#define PPERM        DPpermF
#define PPPERM       DPPpermF
#define WDIAG        DwDiagRF
#define WPDIAG       DwPDiagRF
#define QSORTR2I     dqsortr2i
#define SWAPM        dswapm
#define SPDWDIAG     DSPDwDiagRF
#define BISINIT      Dbisinit
#define STOPBIS      Dstopbis
#define TIDYCG       Dtidycg
#define MGSRO        Dmgsro
#define GIVENS       Dgivens
#define BRKDN        Dbrkdn

#define CSRCSC       Dcsrcsc
#define CSRCSC2      Dcsrcsc2
#define READMTC      Dreadmtc
#define WRITEMTC     Dwritemtc
#define READVECTORS  Dreadvectors
#define WRITEVECTORS Dwritevectors
#define PCG          Dpcg
#define BCG          Dbcg
#define SBCG         DSYMbcg
#define SQMR         DSYMqmr
#define GMRES        Dgmres
#define FGMRES       Dfgmres

#define DISTDOT      ddot
#define DISTDOTU     ddot
#define DISTDOT2     ddot2
#define COPY         dcopy
#define AXPY         daxpy
#define ROTG         drotg
#define ROT          drot
#define NRM          dnrm2
#define I_AMAX       idamax
#define ASUM         dasum
#define SCAL         dscal

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

#define FABS(A)      (((A)>=0)?(A):(-(A)))
#define ABS(A)       dabs(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     dabs
#define SQRTFNAME    dsqrt
#define LOGFNAME     dlog
#define LOGTENFNAME  dlog10
#define CTOD(A)      (A)
#define CONJ(A)      (A)
#define FNULL        0.0d0
#define FONE         1.0d0
#define RONE         1.0d0
#define RTWO         2.0d0
#define RFOUR        4.0d0
#define RZERO        0.0d0
#define REAL_MAX     1.0d300
#define SIGNUM       1.0d0
#define SIGNUM2      1.0d0
#define TWICE        1



#elif defined _SINGLE_REAL_
#define GETEPS       sgeteps
#define _REAL_MAX_   _S_REAL_MAX_
#define REALS        real
#define FLOAT        real
#define CSRMAT       Smat
#define AMGLEVELMAT  SAMGlevelmat
#define ILUPACKPARAM SILUPACKparam
#define GNLAGGRESSIVEDROPPING SGNLaggressivedropping 
#define SPDAGGRESSIVEDROPPING SSPDaggressivedropping 
#define SPTRF        ssptrf
#define SPTRS        ssptrs

#define ROWSCALE     Srowscale
#define COLSCALE     Scolscale
#define SPDSCALE     SSPDscale
#define SYMSCALE     SSYMscale
#define CNRMS        Scnrms
#define RNRMS        Srnrms
#define DIAMUA       SDiagMultA
#define AMUDIA       SAMultDiag

#define CSRMATVEC    SGNLmatvec
#define CSRMATTVEC   SGNLmattvec
#define CSRMATHVEC   SGNLmathvec
#define SYMMATVEC    SSYMmatvec
#define SSMMATVEC    SSSMmatvec
#define HERMATVEC    SSYMmatvec
#define SETUPGRAPH   Ssetupgraph
#define SETUPGRAPH_EPSILON   Ssetupgraph_epsilon

#define AMGDELETE       SGNLAMGdelete
#define SPDAMGDELETE    SSPDAMGdelete

#define AMGSETUP        SGNLAMGfactor
#define SPDAMGSETUP     SSPDAMGfactor
#define AMGFACTOR       SGNLAMGfactor
#define SPDAMGFACTOR    SSPDAMGfactor
#define AMGEXTRACT      SGNLAMGextract
#define SPDAMGEXTRACT   SSYMAMGextract
#define SYMAMGEXTRACT   SSYMAMGextract
#define AMGSOL          SGNLAMGsol
#define AMGDLSOL        SGNLAMGdlsol
#define AMGDUSOL        SGNLAMGdusol
#define AMGLSOL         SGNLAMGlsol
#define AMGUSOL         SGNLAMGusol
#define AMGTDLSOL       SGNLAMGtdlsol
#define AMGTDUSOL       SGNLAMGtdusol
#define AMGTLSOL        SGNLAMGtlsol
#define AMGTUSOL        SGNLAMGtusol
#define AMGTSOL         SGNLAMGtsol
#define SPDAMGSOL       SSPDAMGsol
#define AMGDRIVER       SGNLAMGdriver
#define AMGSOLVER       SGNLAMGsolver
#define SPDAMGSOLVER    SSPDAMGsolver
#define SYMAMGSOLVER    SSYMAMGsolver
#define AMGINIT         SGNLAMGinit
#define AMGGETPARAMS    SGNLAMGgetparams
#define AMGSETPARAMS    SGNLAMGsetparams
#define SPDAMGINIT      SSPDAMGinit
#define SPDAMGGETPARAMS SSPDAMGgetparams
#define SPDAMGSETPARAMS SSPDAMGsetparams
#define SYMAMGINIT      SSYMAMGinit
#define SYMAMGGETPARAMS SSYMAMGgetparams
#define SYMAMGSETPARAMS SSYMAMGsetparams

#define LDLP         SSPDldlp
#define LDLPSOL      SSPDldlpsol
#define LUPQ         SGNLlupq
#define LUPQSOL      SGNLlupqsol
#define LUPQLSOL     SGNLlupqlsol
#define LUPQUSOL     SGNLlupqusol
#define LUPQDLSOL    SGNLlupqdlsol
#define LUPQDUSOL    SGNLlupqdusol
#define LUPQTSOL     SGNLlupqtsol
#define LUPQTLSOL    SGNLlupqtlsol
#define LUPQTUSOL    SGNLlupqtusol
#define LUPQTDLSOL   SGNLlupqtdlsol
#define LUPQTDUSOL   SGNLlupqtdusol



#define ILDLC        SSYMildlc
#define ILDLCSOL     SSYMildlcsol
#define ILUC         SGNLiluc
#define PILUC        SGNLpiluc
#define MPILUC       SGNLmpiluc
#define PILDLC       SSPDpiluc
#define MPILDLC      SSPDmpiluc
#define ILUCSOL      SGNLilucsol
#define ILUCTSOL     SGNLiluctsol
#define ILUCDLSOL    SGNLilucdlsol
#define ILUCTDLSOL   SGNLiluctdlsol
#define ILUCDUSOL    SGNLilucdusol
#define ILUCTDUSOL   SGNLiluctdusol
#define ILUCLSOL     SGNLiluclsol
#define ILUCTLSOL    SGNLiluctlsol
#define ILUCUSOL     SGNLilucusol
#define ILUCTUSOL    SGNLiluctusol