Subversion Repositories gelsvn

Rev

Rev 2 | Rev 45 | Go to most recent revision | Show entire file | Regard whitespace | Details | Blame | Last modification | View Log | RSS feed

Rev 2 Rev 12
Line 23... Line 23...
23 
{
 23 
{
24 
 
 24 
 
25 
	template <class MT>
 25 
  template <class MT>
26 
	int power_eigensolution(const MT& Ap, MT& Q, MT& L, int max_sol)
 26 
  int power_eigensolution(const MT& Ap, MT& Q, MT& L, int max_sol)
27 
	{
 27 
  {
- 
 
 28 
    L = MT(0);
- 
 
 29 
 
28 
		typedef typename MT::VectorType VT;
 30 
    typedef typename MT::VectorType VT;
29 
		MT A = Ap;
 31 
    MT A = Ap;
30 
		int n = s_min(MT::get_v_dim(), max_sol);
 32 
    int n = s_min(MT::get_v_dim(), max_sol);
31 
 
 33 
 
32 
		for(int i=0;i<n;++i)
 34 
    for(int i=0;i<n;++i)