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
|
/*
ARPACK++ v1.2 2/20/2000
c++ interface to ARPACK code.
MODULE ARGComp.h.
Arpack++ class ARCompGenEig definition.
ARPACK Authors
Richard Lehoucq
Danny Sorensen
Chao Yang
Dept. of Computational & Applied Mathematics
Rice University
Houston, Texas
*/
#ifndef ARGCOMP_H
#define ARGCOMP_H
#include <cstddef>
#include "arch.h"
#include "arscomp.h"
#include "argeig.h"
template<class ARFLOAT, class ARFOP, class ARFB>
class ARCompGenEig:
virtual public ARGenEig<ARFLOAT, arcomplex<ARFLOAT>, ARFOP, ARFB>,
virtual public ARCompStdEig<ARFLOAT, ARFOP> {
public:
// a) Constructors and destructor.
ARCompGenEig() { }
// Short constructor (Does nothing but calling base classes constructors).
ARCompGenEig(int np, int nevp, ARFOP* objOPp,
void (ARFOP::* MultOPxp)(arcomplex<ARFLOAT>[],arcomplex<ARFLOAT>[]),
ARFB* objBp,
void (ARFB::* MultBxp)(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).
ARCompGenEig(int np, int nevp, ARFOP* objOPp,
void (ARFOP::* MultOPxp)(arcomplex<ARFLOAT>[],arcomplex<ARFLOAT>[]),
ARFB* objBp,
void (ARFB::* MultBxp)(arcomplex<ARFLOAT>[],arcomplex<ARFLOAT>[]),
arcomplex<ARFLOAT> sigmap,
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).
ARCompGenEig(const ARCompGenEig& other) { Copy(other); }
// Copy constructor.
virtual ~ARCompGenEig() { }
// Destructor.
// b) Operators.
ARCompGenEig& operator=(const ARCompGenEig& other);
// Assignment operator.
}; // class ARCompGenEig.
// ------------------------------------------------------------------------ //
// ARCompGenEig member functions definition. //
// ------------------------------------------------------------------------ //
template<class ARFLOAT, class ARFOP, class ARFB>
inline ARCompGenEig<ARFLOAT, ARFOP, ARFB>::
ARCompGenEig(int np, int nevp, ARFOP* objOPp,
void (ARFOP::* MultOPxp)(arcomplex<ARFLOAT>[],arcomplex<ARFLOAT>[]),
ARFB* objBp,
void (ARFB::* MultBxp)(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, objBp, MultBxp,
whichp, ncvp, tolp, maxitp, residp, ishiftp);
} // Long constructor (regular mode).
template<class ARFLOAT, class ARFOP, class ARFB>
inline ARCompGenEig<ARFLOAT, ARFOP, ARFB>::
ARCompGenEig(int np, int nevp, ARFOP* objOPp,
void (ARFOP::* MultOPxp)(arcomplex<ARFLOAT>[],arcomplex<ARFLOAT>[]),
ARFB* objBp,
void (ARFB::* MultBxp)(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, objBp, MultBxp,
whichp, ncvp, tolp, maxitp, residp, ishiftp);
} // Long constructor (shift and invert mode).
template<class ARFLOAT, class ARFOP, class ARFB>
ARCompGenEig<ARFLOAT, ARFOP, ARFB>& ARCompGenEig<ARFLOAT, ARFOP, ARFB>::
operator=(const ARCompGenEig<ARFLOAT, ARFOP, ARFB>& other)
{
if (this != &other) { // Stroustrup suggestion.
this->ClearMem();
Copy(other);
}
return *this;
} // operator=.
#endif // ARGCOMP_H
|
