28 #ifndef ERGO_MAT_ACC_EXTRAPOLATE_HEADER 29 #define ERGO_MAT_ACC_EXTRAPOLATE_HEADER 39 template<
class Treal,
class Tworker>
44 void Scan(
const Tworker & worker,
49 Treal* errorList_frob_,
50 Treal* errorList_eucl_,
51 Treal* errorList_maxe_,
65 template<
class Treal,
class Tworker>
71 template<
class Treal,
class Tworker>
98 worker.ComputeMatrix(firstParam, accurateMatrix);
100 Treal currParam = firstParam;
101 for(
int i = 0; i < nSteps; i++)
103 currParam *= stepFactor;
104 time_t startTime, endTime;
106 worker.ComputeMatrix(currParam, otherMatrix);
111 errorMatrix = otherMatrix;
112 errorMatrix += (
ergo_real)(-1) * accurateMatrix;
117 Treal euclAcc = 1e-11;
126 template<
class Treal,
class Tworker>
129 Treal* errorList_frob_,
130 Treal* errorList_eucl_,
131 Treal* errorList_maxe_,
double ergo_real
Definition: realtype.h:53
Treal eucl(Treal const requestedAccuracy, int maxIter=-1) const
Definition: MatrixSymmetric.h:661
std::vector< Treal > timeList
Definition: mat_acc_extrapolate.h:61
Describes dimensions of matrix and its blocks on all levels.
Definition: SizesAndBlocks.h:37
Definition: mat_acc_extrapolate.h:40
ergo_real compute_maxabs_sparse(const Tmatrix &M)
Definition: matrix_utilities.h:87
MatAccInvestigator(mat::SizesAndBlocks const &matrix_size_block_info_)
Definition: mat_acc_extrapolate.h:66
void resetSizesAndBlocks(SizesAndBlocks const &newRows, SizesAndBlocks const &newCols)
Definition: MatrixBase.h:74
std::vector< Treal > errorList_eucl
Definition: mat_acc_extrapolate.h:59
std::vector< Treal > threshList
Definition: mat_acc_extrapolate.h:57
std::vector< Treal > errorList_frob
Definition: mat_acc_extrapolate.h:58
void Scan(const Tworker &worker, Treal firstParam, Treal stepFactor, int nSteps)
Definition: mat_acc_extrapolate.h:73
Treal frob() const
Definition: MatrixSymmetric.h:349
int nScanSteps
Definition: mat_acc_extrapolate.h:55
void GetScanResult(Treal *threshList_, Treal *errorList_frob_, Treal *errorList_eucl_, Treal *errorList_maxe_, Treal *timeList_)
Definition: mat_acc_extrapolate.h:128
mat::SizesAndBlocks matrix_size_block_info
Definition: mat_acc_extrapolate.h:54
Treal baseThresh
Definition: mat_acc_extrapolate.h:56
std::vector< Treal > errorList_maxe
Definition: mat_acc_extrapolate.h:60