Ignore:
Timestamp:
Jun 20, 2005 5:53:55 PM (14 years ago)
Author:
andreasw
Message:

ran astyle

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/PDSystemImpl/PDFullSpace/IpPardisoSolverInterface.cpp

    r304 r321  
    1919{
    2020  void F77_FUNC(pardisoinit,PARDISOINIT)(void* PT, const ipfint* MTYPE,
    21                                         ipfint* IPARM);
     21                                        ipfint* IPARM);
    2222  void F77_FUNC(pardiso,PARDISO)(void** PT, const ipfint* MAXFCT,
    23                                 const ipfint* MNUM, const ipfint* MTYPE,
    24                                 const ipfint* PHASE, const ipfint* N,
    25                                 const double* A, const ipfint* IA,
    26                                 const ipfint* JA, const ipfint* PERM,
    27                                 const ipfint* NRHS, ipfint* IPARM,
    28                                 const ipfint* MSGLVL, double* B, double* X,
    29                                 ipfint* ERROR);
     23                                const ipfint* MNUM, const ipfint* MTYPE,
     24                                const ipfint* PHASE, const ipfint* N,
     25                                const double* A, const ipfint* IA,
     26                                const ipfint* JA, const ipfint* PERM,
     27                                const ipfint* NRHS, ipfint* IPARM,
     28                                const ipfint* MSGLVL, double* B, double* X,
     29                                ipfint* ERROR);
    3030}
    3131
     
    6767      ipfint ERROR;
    6868      F77_FUNC(pardiso,PARDISO)(PT_, &MAXFCT_, &MNUM_, &MTYPE_, &PHASE, &N,
    69                                 NULL, NULL, NULL, NULL, &NRHS, IPARM_,
    70                                 &MSGLVL_, NULL, NULL, &ERROR);
     69                                NULL, NULL, NULL, NULL, &NRHS, IPARM_,
     70                                &MSGLVL_, NULL, NULL, &ERROR);
    7171      DBG_ASSERT(ERROR==0);
    7272    }
     
    8989      ipfint ERROR;
    9090      F77_FUNC(pardiso,PARDISO)(PT_, &MAXFCT_, &MNUM_, &MTYPE_, &PHASE, &N,
    91                                 NULL, NULL, NULL, NULL, &NRHS, IPARM_,
    92                                 &MSGLVL_, NULL, NULL, &ERROR);
     91                                NULL, NULL, NULL, NULL, &NRHS, IPARM_,
     92                                &MSGLVL_, NULL, NULL, &ERROR);
    9393      DBG_ASSERT(ERROR==0);
    9494    }
     
    102102    // Call Pardiso's initialization routine
    103103    IPARM_[0] = 0;  // Tell it to fill IPARM with default values(?)
    104     F77_FUNC(pardisoinit,PARDISOINIT)(PT_, &MTYPE_, IPARM_);   
     104    F77_FUNC(pardisoinit,PARDISOINIT)(PT_, &MTYPE_, IPARM_);
    105105
    106106    // Set some parameters for Pardiso
     
    119119
    120120  ESymSolverStatus PardisoSolverInterface::MultiSolve(bool new_matrix,
    121                                                       const Index* ia,
    122                                                       const Index* ja,
    123                                                       Index nrhs,
    124                                                       double* rhs_vals,
    125                                                       bool check_NegEVals,
    126                                                       Index numberOfNegEVals)
     121      const Index* ia,
     122      const Index* ja,
     123      Index nrhs,
     124      double* rhs_vals,
     125      bool check_NegEVals,
     126      Index numberOfNegEVals)
    127127  {
    128128    DBG_START_METH("PardisoSolverInterface::MultiSolve",dbg_verbosity);
     
    155155      elements */
    156156  ESymSolverStatus PardisoSolverInterface::InitializeStructure
    157       (Index dim, Index nonzeros,
    158        const Index* ia,
    159        const Index* ja)
     157  (Index dim, Index nonzeros,
     158   const Index* ia,
     159   const Index* ja)
    160160  {
    161161    DBG_START_METH("PardisoSolverInterface::InitializeStructure",dbg_verbosity);
     
    180180  ESymSolverStatus
    181181  PardisoSolverInterface::SymbolicFactorization(const Index* ia,
    182                                                 const Index* ja)
     182      const Index* ja)
    183183  {
    184184    DBG_START_METH("PardisoSolverInterface::SymbolicFactorization",
    185                    dbg_verbosity);
     185                   dbg_verbosity);
    186186
    187187    // Call Pardiso to do the analysis phase
     
    191191    ipfint NRHS = 0;
    192192    double B;  // This should not be accessed by Pardiso in analysis
    193                // phase
     193    // phase
    194194    double X;  // This should not be accessed by Pardiso in analysis
    195                // phase
     195    // phase
    196196    ipfint ERROR;
    197197
    198198    F77_FUNC(pardiso,PARDISO)(PT_, &MAXFCT_, &MNUM_, &MTYPE_,
    199                               &PHASE, &N, a_, ia, ja, &PERM,
    200                               &NRHS, IPARM_, &MSGLVL_, &B, &X,
    201                               &ERROR);
     199                              &PHASE, &N, a_, ia, ja, &PERM,
     200                              &NRHS, IPARM_, &MSGLVL_, &B, &X,
     201                              &ERROR);
    202202    if (ERROR!=0) {
    203203      Jnlst().Printf(J_ERROR, J_LINEAR_ALGEBRA,
    204204                     "Error in Pardiso during analysis phase.  ERROR = %d.\n",
    205                      ERROR);
     205                     ERROR);
    206206      return SYMSOLVER_FATAL_ERROR;
    207207    }
     
    212212  ESymSolverStatus
    213213  PardisoSolverInterface::Factorization(const Index* ia,
    214                                         const Index* ja,
    215                                         bool check_NegEVals,
    216                                         Index numberOfNegEVals)
     214                                        const Index* ja,
     215                                        bool check_NegEVals,
     216                                        Index numberOfNegEVals)
    217217  {
    218218    DBG_START_METH("PardisoSolverInterface::Factorization",dbg_verbosity);
     
    224224    ipfint NRHS = 0;
    225225    double B;  // This should not be accessed by Pardiso in factorization
    226                // phase
     226    // phase
    227227    double X;  // This should not be accessed by Pardiso in factorization
    228                // phase
     228    // phase
    229229    ipfint ERROR;
    230230
    231231    F77_FUNC(pardiso,PARDISO)(PT_, &MAXFCT_, &MNUM_, &MTYPE_,
    232                               &PHASE, &N, a_, ia, ja, &PERM,
    233                               &NRHS, IPARM_, &MSGLVL_, &B, &X,
    234                               &ERROR);
     232                              &PHASE, &N, a_, ia, ja, &PERM,
     233                              &NRHS, IPARM_, &MSGLVL_, &B, &X,
     234                              &ERROR);
    235235    if (ERROR==-4) {
    236236      // I think this means that the matrix is singular
     
    264264
    265265  ESymSolverStatus PardisoSolverInterface::Solve(const Index* ia,
    266                                                 const Index* ja,
    267                                                 Index nrhs,
    268                                                 double *rhs_vals)
     266      const Index* ja,
     267      Index nrhs,
     268      double *rhs_vals)
    269269  {
    270270    DBG_START_METH("PardisoSolverInterface::Solve",dbg_verbosity);
     
    279279
    280280    F77_FUNC(pardiso,PARDISO)(PT_, &MAXFCT_, &MNUM_, &MTYPE_,
    281                               &PHASE, &N, a_, ia, ja, &PERM,
    282                               &NRHS, IPARM_, &MSGLVL_, rhs_vals, X,
    283                               &ERROR);
     281                              &PHASE, &N, a_, ia, ja, &PERM,
     282                              &NRHS, IPARM_, &MSGLVL_, rhs_vals, X,
     283                              &ERROR);
    284284    if (ERROR!=0 ) {
    285285      Jnlst().Printf(J_ERROR, J_LINEAR_ALGEBRA,
Note: See TracChangeset for help on using the changeset viewer.