Setting fill in and treshold for ILUT

Hi all.
I'm having some problems to solve accurately the pressure poisson equation 
in a CFD solver, and I think maybe I'm using the wrong preconditioner. My 
matrix type is defined bellow.

typedef mtl::matrix< double, mtl::rectangle<>, mtl::array< mtl::compressed<> 
>, mtl::row_major >::type MTLmatrix;

I'm using ILU<MTLmatrix> as a preconditioner, and I want to switch to ILUT. 
How can I set precisely the values for fill in and treshold? Is it a good 
guess to set fill in based on the largest value in the pressure matrix? In 
that case, how can I find the largest value in the matrix with the max() 
algorithm? Is there a function that returns the i-th row of the matrix?

Thanks in Advance,
Márcio Ricardo
-- 
Márcio Ricardo Pivello
Mechanical Engineer, MSc.
LTCM and CFD Lab
Mechanical Engineering School
Federal University of Uberlândia 
Uberlândia - MG - Brazil

_______________________________________________
This list is archived at http://www.osl.iu.edu/MailArchives/mtl-devel/