* eigenvalue estimates.
*/
bool log_coefficients;
- /**
+
+ /**
* Compute the condition
* number of the projected
* matrix.
* @note Requires LAPACK support.
*/
bool compute_condition_number;
- /**
+
+ /**
* Compute the condition
* number of the projected
* matrix in each step.
// Write eigenvalues or condition number
if (do_eigenvalues)
- {
- TridiagonalMatrix<double> T(diagonal.size(), true);
- for (unsigned int i=0;i<diagonal.size();++i)
- {
- T(i,i) = diagonal[i];
- if (i< diagonal.size()-1)
- T(i,i+1) = offdiagonal[i];
- }
- T.compute_eigenvalues();
- if (additional_data.compute_condition_number
- & ! additional_data.compute_all_condition_numbers)
- deallog << "Condition number estimate: " <<
- T.eigenvalue(T.n()-1)/T.eigenvalue(0) << std::endl;
- if (additional_data.compute_eigenvalues)
- {
- for (unsigned int i=0;i<T.n();++i)
- deallog << ' ' << T.eigenvalue(i);
- deallog << std::endl;
- }
- }
+ {
+ TridiagonalMatrix<double> T(diagonal.size(), true);
+ for (unsigned int i=0;i<diagonal.size();++i)
+ {
+ T(i,i) = diagonal[i];
+ if (i< diagonal.size()-1)
+ T(i,i+1) = offdiagonal[i];
+ }
+ T.compute_eigenvalues();
+ if (additional_data.compute_condition_number
+ & ! additional_data.compute_all_condition_numbers)
+ deallog << "Condition number estimate: " <<
+ T.eigenvalue(T.n()-1)/T.eigenvalue(0) << std::endl;
+ if (additional_data.compute_eigenvalues)
+ {
+ for (unsigned int i=0;i<T.n();++i)
+ deallog << ' ' << T.eigenvalue(i);
+ deallog << std::endl;
+ }
+ }
// Deallocate Memory
cleanup();