Whn using two square matrices, A and B, each of size NxN, with A being the Cholesky factorization of a positive definite matrix with a zeroed lower triangle, and B being an identity matrix, the program segfaults. I checked the obvious, and it still fails. It looks like anything beyond tile size either segfaults or gives nondeterministic results. The code is an adapted version from the available examples in the official documentation.
Do you have any code snippets?