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
|
/*
ARPACK++ v1.2 2/20/2000
c++ interface to ARPACK code.
MODULE ARSComp.h.
Arpack++ class ARCompStdEig definition.
ARPACK Authors
Richard Lehoucq
Danny Sorensen
Chao Yang
Dept. of Computational & Applied Mathematics
Rice University
Houston, Texas
*/
#ifndef ARSCOMP_H
#define ARSCOMP_H
#include <cstddef>
#include "arch.h"
#include "arseig.h"
#include "arrscomp.h"
template<class ARFLOAT, class ARFOP>
class ARCompStdEig:
virtual public ARStdEig<ARFLOAT, arcomplex<ARFLOAT>, ARFOP>,
virtual public ARrcCompStdEig<ARFLOAT> {
public:
// a) Constructors and destructor.
ARCompStdEig() { }
// Short constructor.
ARCompStdEig(int np, int nevp, ARFOP* objOPp,
void (ARFOP::* MultOPxp)(arcomplex<ARFLOAT>[],arcomplex<ARFLOAT>[]),
char* whichp = "LM", int ncvp = 0,
ARFLOAT tolp = 0.0, int maxitp = 0,
arcomplex<ARFLOAT>* residp = NULL, bool ishiftp = true);
// Long constructor (regular mode).
ARCompStdEig(int np, int nevp, ARFOP* objOPp,
void (ARFOP::* MultOPxp)(arcomplex<ARFLOAT>[],arcomplex<ARFLOAT>[]),
arcomplex<ARFLOAT> sigma, char* whichp = "LM",
int ncvp = 0, ARFLOAT tolp = 0.0, int maxitp = 0,
arcomplex<ARFLOAT>* residp = NULL, bool ishiftp = true);
// Long constructor (shift and invert mode).
ARCompStdEig(const ARCompStdEig& other) { Copy(other); }
// Copy constructor.
virtual ~ARCompStdEig() { }
// Destructor.
// b) Operators.
ARCompStdEig& operator=(const ARCompStdEig& other);
// Assignment operator.
}; // class ARCompStdEig.
// ------------------------------------------------------------------------ //
// ARCompStdEig member functions definition. //
// ------------------------------------------------------------------------ //
template<class ARFLOAT, class ARFOP>
inline ARCompStdEig<ARFLOAT, ARFOP>::
ARCompStdEig(int np, int nevp, ARFOP* objOPp,
void (ARFOP::* MultOPxp)(arcomplex<ARFLOAT>[],arcomplex<ARFLOAT>[]),
char* whichp, int ncvp, ARFLOAT tolp, int maxitp,
arcomplex<ARFLOAT>* residp, bool ishiftp)
{
this->NoShift();
this->DefineParameters(np, nevp, objOPp, MultOPxp, whichp,
ncvp, tolp, maxitp, residp, ishiftp);
} // Long constructor (regular mode).
template<class ARFLOAT, class ARFOP>
inline ARCompStdEig<ARFLOAT, ARFOP>::
ARCompStdEig(int np, int nevp, ARFOP* objOPp,
void (ARFOP::* MultOPxp)(arcomplex<ARFLOAT>[],arcomplex<ARFLOAT>[]),
arcomplex<ARFLOAT> sigmap, char* whichp, int ncvp,
ARFLOAT tolp, int maxitp, arcomplex<ARFLOAT>* residp,
bool ishiftp)
{
this->ChangeShift(sigmap);
this->DefineParameters(np, nevp, objOPp, MultOPxp, whichp,
ncvp, tolp, maxitp, residp, ishiftp);
} // Long constructor (shift and invert mode).
template<class ARFLOAT, class ARFOP>
ARCompStdEig<ARFLOAT, ARFOP>& ARCompStdEig<ARFLOAT, ARFOP>::
operator=(const ARCompStdEig<ARFLOAT, ARFOP>& other)
{
if (this != &other) { // Stroustrup suggestion.
this->ClearMem();
Copy(other);
}
return *this;
} // operator=.
#endif // ARSCOMP_H
|
