I think declare(h,scalar) in mnewton causes compiled mnewton to not run. Of course, compiling mnewton doesn't make it faster (likely most users are running uncompiled lu_factor). In general, undeclared arrays makes code more breakable, I think. --Barton