GSL version 2.7 0 0 1 -nan -nan -nan -nan -nan -nan lu[1,:] = 0 0 1 lu[2,:] = -nan -nan -nan lu[3,:] = -nan -nan -nan p = 0 1 2