In version 4.4.0 I wold like to use the SolveSteadyNavierStokes() function to solve my system with the Newton method but it always returns a similar output: L2Norm[0] = -nan L2Norm[1] = -nan Is the method working in that version (because I didn't found it in the documentation)? Do I have to initialize in some particular way the m_fields variable? At the moment I initialized each field with zeros but I don't understand if I have to do something else. Thank you, Moreno Pintore