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
|
/*
ARPACK++ v1.0 8/1/1997
c++ interface to ARPACK code.
MODULE ARUSComp.h.
Arpack++ class ARluCompStdEig definition
(umfpack version).
ARPACK Authors
Richard Lehoucq
Danny Sorensen
Chao Yang
Dept. of Computational & Applied Mathematics
Rice University
Houston, Texas
*/
#ifndef ARUSCOMP_H
#define ARUSCOMP_H
#include <stddef.h>
#include "arch.h"
#include "arscomp.h"
#include "arunsmat.h"
#include "arrseig.h"
template<class FLOAT>
class ARluCompStdEig:
public virtual ARCompStdEig<FLOAT, ARumNonSymMatrix<arcomplex<FLOAT> > > {
public:
// a) Public functions:
// a.1) Functions that allow changes in problem parameters.
virtual void ChangeShift(arcomplex<FLOAT> sigmaRp);
virtual void SetRegularMode();
virtual void SetShiftInvertMode(arcomplex<FLOAT> sigmap);
// a.2) Constructors and destructor.
ARluCompStdEig() { }
// Short constructor.
ARluCompStdEig(int nevp, ARumNonSymMatrix<arcomplex<FLOAT> >& A,
char* whichp = "LM", int ncvp = 0,
FLOAT tolp = 0.0, int maxitp = 0,
arcomplex<FLOAT>* residp = NULL, bool ishiftp = true);
// Long constructor (regular mode).
ARluCompStdEig(int nevp, ARumNonSymMatrix<arcomplex<FLOAT> >& A,
arcomplex<FLOAT> sigma, char* whichp = "LM",
int ncvp = 0, FLOAT tolp = 0.0, int maxitp = 0,
arcomplex<FLOAT>* residp = NULL, bool ishiftp = true);
// Long constructor (shift and invert mode).
ARluCompStdEig(const ARluCompStdEig& other) { Copy(other); }
// Copy constructor.
virtual ~ARluCompStdEig() { }
// Destructor.
// b) Operators.
ARluCompStdEig& operator=(const ARluCompStdEig& other);
// Assignment operator.
}; // class ARluCompStdEig.
// ------------------------------------------------------------------------ //
// ARluCompStdEig member functions definition. //
// ------------------------------------------------------------------------ //
template<class FLOAT>
inline void ARluCompStdEig<FLOAT>::
ChangeShift(arcomplex<FLOAT> sigmaRp)
{
objOP->FactorAsI(sigmaRp);
ARrcStdEig<FLOAT, arcomplex<FLOAT> >::ChangeShift(sigmaRp);
} // ChangeShift.
template<class FLOAT>
inline void ARluCompStdEig<FLOAT>::SetRegularMode()
{
ARStdEig<FLOAT, arcomplex<FLOAT>, ARumNonSymMatrix<arcomplex<FLOAT> > >::
SetRegularMode(objOP, &ARumNonSymMatrix<arcomplex<FLOAT> >::MultMv);
} // SetRegularMode.
template<class FLOAT>
inline void ARluCompStdEig<FLOAT>::
SetShiftInvertMode(arcomplex<FLOAT> sigmap)
{
ARStdEig<FLOAT, arcomplex<FLOAT>, ARumNonSymMatrix<arcomplex<FLOAT> > >::
SetShiftInvertMode(sigmap, objOP,
&ARumNonSymMatrix<arcomplex<FLOAT> >::MultInvv);
} // SetShiftInvertMode.
template<class FLOAT>
inline ARluCompStdEig<FLOAT>::
ARluCompStdEig(int nevp, ARumNonSymMatrix<arcomplex<FLOAT> >& A,
char* whichp, int ncvp, FLOAT tolp,
int maxitp, arcomplex<FLOAT>* residp, bool ishiftp)
{
NoShift();
DefineParameters(A.ncols(), nevp, &A,
ARumNonSymMatrix<arcomplex<FLOAT> >::MultMv,
whichp, ncvp, tolp, maxitp, residp, ishiftp);
} // Long constructor (regular mode).
template<class FLOAT>
inline ARluCompStdEig<FLOAT>::
ARluCompStdEig(int nevp, ARumNonSymMatrix<arcomplex<FLOAT> >& A,
arcomplex<FLOAT> sigmap, char* whichp, int ncvp,
FLOAT tolp, int maxitp, arcomplex<FLOAT>* residp,
bool ishiftp)
{
DefineParameters(A.ncols(), nevp, &A,
&ARumNonSymMatrix<arcomplex<FLOAT> >::MultInvv,
whichp, ncvp, tolp, maxitp, residp, ishiftp);
ChangeShift(sigmap);
} // Long constructor (shift and invert mode).
template<class FLOAT>
ARluCompStdEig<FLOAT>& ARluCompStdEig<FLOAT>::
operator=(const ARluCompStdEig<FLOAT>& other)
{
if (this != &other) { // Stroustrup suggestion.
ClearMem();
Copy(other);
}
return *this;
} // operator=.
#endif // ARUSCOMP_H
|