Show
Ignore:
Timestamp:
06/05/08 14:58:48 (6 months ago)
Author:
pernet
Message:

Update the randiter:
* add the non-zero randiter class (following LinBox? structure)
* modify finite fields accordingly
* Update implementation of finite fields to fix bugs, improve efficiency and readability
* minor tweeks

Files:
1 modified

Legend:

Unmodified
Added
Removed
  • include/fflas-ffpack/ffpack_frobenius.inl

    r58 r66  
    4343        // Picking a random noc x N block vector U^T 
    4444        typename Field::RandIter g (F); 
     45        typename Field::NonZeroRandIter nzg (F,g); 
    4546        for (size_t i = 0; i < noc; ++i) 
    4647                for (size_t j = 0; j < N; ++j) 
    4748                        g.random( *(K + i*ldk +j) ); 
    4849        for (size_t i = 0; i < noc; ++i) 
    49                 g.nonzerorandom (*(K + i*ldk +i)); 
     50                nzg.random (*(K + i*ldk +i)); 
    5051         
    5152        // Computing the bloc Krylov matrix [U AU .. A^(c-1) U]^T