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
|
#include <iostream>
#include <random>
#include <Eigen/Core>
#include <Eigen/SparseCore>
#include <Spectra/JDSymEigsBase.h>
#include <Spectra/MatOp/DenseGenMatProd.h>
#include "catch.hpp"
using namespace Spectra;
using Matrix = Eigen::MatrixXd;
using Vector = Eigen::VectorXd;
using ComplexMatrix = Eigen::MatrixXcd;
using ComplexVector = Eigen::VectorXcd;
using SpMatrix = Eigen::SparseMatrix<double>;
using Index = Eigen::Index;
template <typename OpType>
class JDMock : public JDSymEigsBase<JDMock<OpType>, OpType>
{
public:
JDMock(OpType& op, Index nev) :
JDSymEigsBase<JDMock<OpType>, OpType>(op, nev) {}
Matrix setup_initial_search_space(SortRule) const
{
return Matrix::Zero(0, 0);
}
Matrix calculate_correction_vector() const
{
return Matrix::Zero(0, 0);
}
};
TEST_CASE("Constructing JDSymObject", "[eigs_gen]")
{
const Matrix A = Eigen::MatrixXd::Random(10, 10);
DenseGenMatProd<double> op(A);
JDMock<DenseGenMatProd<double>> eigs(op, 5);
REQUIRE(eigs.num_iterations() == 0);
REQUIRE(eigs.info() == Spectra::CompInfo::NotComputed);
}
|
