40#ifndef ERGO_MAT_ACC_EXTRAPOLATE_HEADER
41#define ERGO_MAT_ACC_EXTRAPOLATE_HEADER
51template<
class Treal,
class Tworker>
56 void Scan(
const Tworker & worker,
61 Treal* errorList_frob_,
62 Treal* errorList_eucl_,
63 Treal* errorList_maxe_,
77template<
class Treal,
class Tworker>
83template<
class Treal,
class Tworker>
110 worker.ComputeMatrix(firstParam, accurateMatrix);
112 Treal currParam = firstParam;
113 for(
int i = 0; i < nSteps; i++)
115 currParam *= stepFactor;
116 time_t startTime, endTime;
118 worker.ComputeMatrix(currParam, otherMatrix);
123 errorMatrix = otherMatrix;
124 errorMatrix += (
ergo_real)(-1) * accurateMatrix;
129 Treal euclAcc = 1e-11;
138template<
class Treal,
class Tworker>
141 Treal* errorList_frob_,
142 Treal* errorList_eucl_,
143 Treal* errorList_maxe_,
Definition mat_acc_extrapolate.h:53
void Scan(const Tworker &worker, Treal firstParam, Treal stepFactor, int nSteps)
Definition mat_acc_extrapolate.h:85
Treal baseThresh
Definition mat_acc_extrapolate.h:68
std::vector< Treal > errorList_frob
Definition mat_acc_extrapolate.h:70
std::vector< Treal > threshList
Definition mat_acc_extrapolate.h:69
std::vector< Treal > errorList_eucl
Definition mat_acc_extrapolate.h:71
mat::SizesAndBlocks matrix_size_block_info
Definition mat_acc_extrapolate.h:66
std::vector< Treal > timeList
Definition mat_acc_extrapolate.h:73
int nScanSteps
Definition mat_acc_extrapolate.h:67
MatAccInvestigator(mat::SizesAndBlocks const &matrix_size_block_info_)
Definition mat_acc_extrapolate.h:78
std::vector< Treal > errorList_maxe
Definition mat_acc_extrapolate.h:72
void GetScanResult(Treal *threshList_, Treal *errorList_frob_, Treal *errorList_eucl_, Treal *errorList_maxe_, Treal *timeList_)
Definition mat_acc_extrapolate.h:140
void resetSizesAndBlocks(SizesAndBlocks const &newRows, SizesAndBlocks const &newCols)
Definition MatrixBase.h:76
Treal eucl(Treal const requestedAccuracy, int maxIter=-1) const
Definition MatrixSymmetric.h:673
Treal frob() const
Definition MatrixSymmetric.h:360
Describes dimensions of matrix and its blocks on all levels.
Definition SizesAndBlocks.h:45
Utilities related to the hierarchical matrix library (HML), including functions for setting up permut...
ergo_real compute_maxabs_sparse(const Tmatrix &M)
Definition matrix_utilities.h:97
double ergo_real
Definition realtype.h:69
MatrixSymmetric< real, matri > symmMatrix
Definition test_LanczosSeveralLargestEig.cc:69