Solved by ptrblck in post #2 
              
              
                Linear algebra methods are usually unstable in reduced precision, so you should use float32 for these kind of operations. 
I don’t think float16 support will be added to lu_cuda for the aforementioned stability reasons.