What happens if you use the /double keyword in INVERT?

Same problem, Does not give a 3x3 identity matrix. 1298510.0 -1.9937495 -0.59229549 -1.9937495 3.7413794e-006 -1.3584707e-007 -0.59229549 -1.3584707e-007 1.8765486e-006

Using the following program: pro ex_invert_test n= [ [ 10.0000, 5.45787e+006, 3.55141e+006], $ [ 5.45787e+006, 2.97883e+012, 1.93831e+012], $ [ 3.55141e+006, 1.93831e+012, 1.26125e+012] ] u=invert(n, /double) print, u#n END outputs the following: IDL> ex_invert_test % Compiled module: EX_INVERT_TEST. 1.0000000 -4.4408921e-016 8.8817842e-016 -0.00035095215 1.0000000 2.3283064e-010 -0.00019073486 -1.1641532e-010 1.0000000 IDL> Do you get the same when you run the program?