Hi

it seems that there's something wrong with the inverse algorithm that is used to find inverse of a matrix in jacket. if you try following test

>> A=kron(eye(2),rand(N));

>> norm(double(inv(gdouble(A)))-inv(A))

where A is a block diagonal matrix, you'll see that for N<64 norm of difference of the result from gpu and the result from cpu are the same but for N>=65 the results do not match. Actually the result only matches for the first block only!

Does anybody know how to fix this problem please let me know.