C++ for Scientists - Technische Universität Dresden

270 CHAPTER 14. NUMERICAL EXERCISES
}
} ;
athens::richardson( poisson scaled(), 0.5∗x, y, 1.e−8, 1000 ) ;
int main() {
typedef glas::dense vector v type ;
v type x( 10 ) ;
glas::random seed seed ;
glas::random( x, seed ) ;
double lambda ;
athens::inverse iteration( solve(), x, 100, lambda ) ;
std::cout ≪ ”lambda = ” ≪ lambda ≪ std::endl ;
std::ofstream xf( ”x.out” ) ;
for ( int i=0; i