Arpack.hpp

/* -----------------------------------------------------------------------------
 *  Diana process modelling, simulation and analysis software
 *  Copyright (c) 2005, Michael Krasnyk
 *  All rights reserved.

 *  This file is a part of Diana process modelling, simulation and analysis software
 *
 *  Diana is free software; you can redistribute it and/or modify it
 *  under the terms of the GNU General Public License as published
 *  by the Free Software Foundation (see accompanying file LICENSE)
 * -----------------------------------------------------------------------------
 *  $Id: Arpack.hpp 9369 2011-05-04 09:33:10Z miha $
 * -----------------------------------------------------------------------------
 *  Description:
 */

#ifndef ARPACK__HPP
#define ARPACK__HPP

#include "Fortran.hpp"

enum ArpackIparam{
  // IPARAM(1) = ISHIFT: method for selecting the implicit shifts.
  ArpackIparam_IShift  = 0,
  // IPARAM(2) = No longer referenced.
  // IPARAM(3) = MXITER
  // On INPUT:  maximum number of Arnoldi update iterations allowed.
  // On OUTPUT: actual number of Arnoldi update iterations taken. 
  ArpackIparam_MaxIter = 2,
  // IPARAM(4) = NB: blocksize to be used in the recurrence.
  // The code currently works only for NB = 1.
  ArpackIparam_NB = 3,
  // IPARAM(5) = NCONV: number of "converged" Ritz values.
  // This represents the number of Ritz values that satisfy the convergence criterion.
  ArpackIparam_NConv = 4,
  // IPARAM(6) = IUPD
  // No longer referenced. Implicit restarting is ALWAYS used.  
  // IPARAM(7) = MODE
  // On INPUT determines what type of eigenproblem is being solved.
  // Must be 1,2,3,4
  ArpackIparam_Mode = 6,
  // IPARAM(8) = NP
  ArpackIparam_NP = 7,
  // IPARAM(9) = NUMOP
  // OUTPUT: NUMOP  = total number of OP*x operations,
  ArpackIparam_NumOP = 8,
  // IPARAM(10) = NUMOPB
  // NUMOPB = total number of B*x operations if BMAT='G',
  ArpackIparam_NumOPB = 9,
  // IPARAM(11) = NUMREO
  // NUMREO = total number of steps of re-orthogonalization.
  ArpackIparam_NumREO = 10,
  ArpackIparam_Size,
};

enum ArpackMode{
  ArpackMode_Regular         = 1,
  ArpackMode_RegularInverse  = 2,
  ArpackMode_ShiftInvert     = 3,
  ArpackMode_Buckling        = 4,
  ArpackMode_Cayley          = 5,
};

// IPNTR   Integer array of length 14.  (OUTPUT)
const int ArpackIpntr_Size=14;

struct arpack_debug_t {
  fortrani logfil, ndigit, mgetv0, msaupd, msaup2, msaitr, mseigt, msapps, 
    msgets, mseupd, mnaupd, mnaup2, mnaitr, mneigh, mnapps, mngets, 
    mneupd, mcaupd, mcaup2, mcaitr, mceigh, mcapps, mcgets, mceupd;
}debug_;

// TODO: configuration of names (non-gcc mangling)
#define DNAUPD (dnaupd_)
#define DNEUPD (dneupd_)
#define DSAUPD (dsaupd_)
#define DSEUPD (dseupd_)

extern "C" {
  void DNAUPD(fortrani& ido, const fortrans bmat, const fortrani& n, const fortrans which, const fortrani& nev,
              const fortrand& tol, fortrand* resid, const fortrani& ncv, fortrand* v, const fortrani& ldv,
              fortrani* iparam, fortrani* ipntr, fortrand* workd, fortrand* workl, const fortrani& lworkl,
              fortrani& info);
  void DNEUPD(const fortranl& rvec, const fortrans howmny, const fortranl* select, fortrand* dr, fortrand* di,
              fortrand* z, const fortrani& ldz, const fortrand& sigmar, const fortrand& sigmai,
              fortrand* workev, const fortrans bmat, const fortrani& n, const fortrans which, const fortrani& nev,
              const fortrand& tol, fortrand* resid, const fortrani& ncv, fortrand* v, const fortrani& ldv,
              fortrani* iparam, fortrani* ipntr, fortrand* workd, fortrand* workl, const fortrani& lworkl,
              fortrani& info);
  void DSAUPD(fortrani& ido, const fortrans bmat, const fortrani& n, const fortrans which, const fortrani& nev,
              const fortrand& tol, fortrand* resid, const fortrani& ncv, fortrand* v, const fortrani& ldv,
              fortrani* iparam, fortrani* ipntr, fortrand* workd, fortrand* workl, const fortrani& lworkl,
              fortrani& info);
  void DSEUPD(const fortranl& rvec, const fortrans howmny, const fortranl* select, fortrand* dr, 
              fortrand* z, const fortrani& ldz, const fortrand& sigmar, 
              const fortrans bmat, const fortrani& n, const fortrans which, const fortrani& nev,
              const fortrand& tol, fortrand* resid, const fortrani& ncv, fortrand* v, const fortrani& ldv,
              fortrani* iparam, fortrani* ipntr, fortrand* workd, fortrand* workl, const fortrani& lworkl,
              fortrani& info);
};

#endif // ARPACK__HPP