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
|
/*
ARPACK++ v1.2 2/20/2000
c++ interface to ARPACK code.
MODULE ARBSSym.h.
Arpack++ class ARluSymStdEig definition
(band matrix version).
ARPACK Authors
Richard Lehoucq
Danny Sorensen
Chao Yang
Dept. of Computational & Applied Mathematics
Rice University
Houston, Texas
*/
#ifndef ARBSSYM_H
#define ARBSSYM_H
#include <cstddef>
#include "arch.h"
#include "arssym.h"
#include "arbsmat.h"
template<class ARFLOAT>
class ARluSymStdEig:
public virtual ARSymStdEig<ARFLOAT, ARbdSymMatrix<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.
ARluSymStdEig() { }
// Short constructor.
ARluSymStdEig(int nevp, ARbdSymMatrix<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).
ARluSymStdEig(int nevp, ARbdSymMatrix<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).
ARluSymStdEig(const ARluSymStdEig& other) { Copy(other); }
// Copy constructor.
virtual ~ARluSymStdEig() { }
// Destructor.
// b) Operators.
ARluSymStdEig& operator=(const ARluSymStdEig& other);
// Assignment operator.
}; // class ARluSymStdEig.
// ------------------------------------------------------------------------ //
// ARluSymStdEig member functions definition. //
// ------------------------------------------------------------------------ //
template<class ARFLOAT>
inline void ARluSymStdEig<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 ARluSymStdEig<ARFLOAT>::SetRegularMode()
{
ARStdEig<ARFLOAT, ARFLOAT, ARbdSymMatrix<ARFLOAT> >::
SetRegularMode(this->objOP, &ARbdSymMatrix<ARFLOAT>::MultMv);
} // SetRegularMode.
template<class ARFLOAT>
inline void ARluSymStdEig<ARFLOAT>::SetShiftInvertMode(ARFLOAT sigmap)
{
ARStdEig<ARFLOAT, ARFLOAT, ARbdSymMatrix<ARFLOAT> >::
SetShiftInvertMode(sigmap, this->objOP, &ARbdSymMatrix<ARFLOAT>::MultInvv);
} // SetShiftInvertMode.
template<class ARFLOAT>
inline ARluSymStdEig<ARFLOAT>::
ARluSymStdEig(int nevp, ARbdSymMatrix<ARFLOAT>& A,
char* whichp, int ncvp, ARFLOAT tolp,
int maxitp, ARFLOAT* residp, bool ishiftp)
{
this->NoShift();
this->DefineParameters(A.ncols(), nevp, &A, &ARbdSymMatrix<ARFLOAT>::MultMv,
whichp, ncvp, tolp, maxitp, residp, ishiftp);
} // Long constructor (regular mode).
template<class ARFLOAT>
inline ARluSymStdEig<ARFLOAT>::
ARluSymStdEig(int nevp, ARbdSymMatrix<ARFLOAT>& A,
ARFLOAT sigmap, char* whichp, int ncvp, ARFLOAT tolp,
int maxitp, ARFLOAT* residp, bool ishiftp)
{
this->DefineParameters(A.ncols(), nevp, &A, &ARbdSymMatrix<ARFLOAT>::MultInvv,
whichp, ncvp, tolp, maxitp, residp, ishiftp);
this->ChangeShift(sigmap);
} // Long constructor (shift and invert mode).
template<class ARFLOAT>
ARluSymStdEig<ARFLOAT>& ARluSymStdEig<ARFLOAT>::
operator=(const ARluSymStdEig<ARFLOAT>& other)
{
if (this != &other) { // Stroustrup suggestion.
this->ClearMem();
Copy(other);
}
return *this;
} // operator=.
#endif // ARBSSYM_H
|