1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160 161 162 163
|
/*
ARPACK++ v1.2 2/20/2000
c++ interface to ARPACK code.
MODULE ARBSNSym.h.
Arpack++ class ARluNonSymStdEig definition
(band matrix version).
ARPACK Authors
Richard Lehoucq
Danny Sorensen
Chao Yang
Dept. of Computational & Applied Mathematics
Rice University
Houston, Texas
*/
#ifndef ARBSNSYM_H
#define ARBSNSYM_H
#include <cstddef>
#include "arch.h"
#include "arsnsym.h"
#include "arbnsmat.h"
template<class ARFLOAT>
class ARluNonSymStdEig:
public virtual ARNonSymStdEig<ARFLOAT, ARbdNonSymMatrix<ARFLOAT, ARFLOAT> > {
public:
// a) Public functions:
// a.1) Functions that allow changes in problem parameters.
virtual void ChangeShift(ARFLOAT sigmaRp);
virtual void SetRegularMode();
virtual void SetShiftInvertMode(ARFLOAT sigmap);
// a.2) Constructors and destructor.
ARluNonSymStdEig() { }
// Short constructor.
ARluNonSymStdEig(int nevp, ARbdNonSymMatrix<ARFLOAT, ARFLOAT>& A,
char* whichp = "LM", int ncvp = 0,
ARFLOAT tolp = 0.0, int maxitp = 0,
ARFLOAT* residp = NULL, bool ishiftp = true);
// Long constructor (regular mode).
ARluNonSymStdEig(int nevp, ARbdNonSymMatrix<ARFLOAT, ARFLOAT>& A,
ARFLOAT sigma, char* whichp = "LM", int ncvp = 0,
ARFLOAT tolp = 0.0, int maxitp = 0,
ARFLOAT* residp = NULL, bool ishiftp = true);
// Long constructor (shift and invert mode).
ARluNonSymStdEig(const ARluNonSymStdEig& other) { Copy(other); }
// Copy constructor.
virtual ~ARluNonSymStdEig() { }
// Destructor.
// b) Operators.
ARluNonSymStdEig& operator=(const ARluNonSymStdEig& other);
// Assignment operator.
}; // class ARluNonSymStdEig.
// ------------------------------------------------------------------------ //
// ARluNonSymStdEig member functions definition. //
// ------------------------------------------------------------------------ //
template<class ARFLOAT>
inline void ARluNonSymStdEig<ARFLOAT>::
ChangeShift(ARFLOAT sigmaRp)
{
this->sigmaR = sigmaRp;
this->sigmaI = 0.0;
this->mode = 3;
this->iparam[7] = this->mode;
this->objOP->FactorAsI(this->sigmaR);
this->Restart();
} // ChangeShift.
template<class ARFLOAT>
inline void ARluNonSymStdEig<ARFLOAT>::SetRegularMode()
{
ARStdEig<ARFLOAT, ARFLOAT, ARbdNonSymMatrix<ARFLOAT, ARFLOAT> >::
SetRegularMode(this->objOP, &ARbdNonSymMatrix<ARFLOAT, ARFLOAT>::MultMv);
} // SetRegularMode.
template<class ARFLOAT>
inline void ARluNonSymStdEig<ARFLOAT>::SetShiftInvertMode(ARFLOAT sigmap)
{
ARStdEig<ARFLOAT, ARFLOAT, ARbdNonSymMatrix<ARFLOAT, ARFLOAT> >::
SetShiftInvertMode(sigmap, this->objOP,
&ARbdNonSymMatrix<ARFLOAT, ARFLOAT>::MultInvv);
} // SetShiftInvertMode.
template<class ARFLOAT>
inline ARluNonSymStdEig<ARFLOAT>::
ARluNonSymStdEig(int nevp, ARbdNonSymMatrix<ARFLOAT, ARFLOAT>& A,
char* whichp, int ncvp, ARFLOAT tolp,
int maxitp, ARFLOAT* residp, bool ishiftp)
{
this->NoShift();
this->DefineParameters(A.ncols(), nevp, &A,
&(ARbdNonSymMatrix<ARFLOAT, ARFLOAT>::MultMv),
whichp, ncvp, tolp, maxitp, residp, ishiftp);
} // Long constructor (regular mode).
template<class ARFLOAT>
inline ARluNonSymStdEig<ARFLOAT>::
ARluNonSymStdEig(int nevp, ARbdNonSymMatrix<ARFLOAT, ARFLOAT>& A,
ARFLOAT sigmap, char* whichp, int ncvp, ARFLOAT tolp,
int maxitp, ARFLOAT* residp, bool ishiftp)
{
this->DefineParameters(A.ncols(), nevp, &A,
&ARbdNonSymMatrix<ARFLOAT, ARFLOAT>::MultInvv,
whichp, ncvp, tolp, maxitp, residp, ishiftp);
this->ChangeShift(sigmap);
} // Long constructor (shift and invert mode).
template<class ARFLOAT>
ARluNonSymStdEig<ARFLOAT>& ARluNonSymStdEig<ARFLOAT>::
operator=(const ARluNonSymStdEig<ARFLOAT>& other)
{
if (this != &other) { // Stroustrup suggestion.
this->ClearMem();
Copy(other);
}
return *this;
} // operator=.
#endif // ARBSNSYM_H
|