diff options
Diffstat (limited to 'src/kernels/level3/invert_diagonal_blocks.opencl')
| -rw-r--r-- | src/kernels/level3/invert_diagonal_blocks.opencl | 4 |
1 files changed, 3 insertions, 1 deletions
diff --git a/src/kernels/level3/invert_diagonal_blocks.opencl b/src/kernels/level3/invert_diagonal_blocks.opencl index e94b4d30..d43b9b7c 100644 --- a/src/kernels/level3/invert_diagonal_blocks.opencl +++ b/src/kernels/level3/invert_diagonal_blocks.opencl @@ -100,7 +100,9 @@ void InvertDiagonalBlock(int n, __global const real* restrict src, const int src if (unit_diagonal == 0) { const real diagonal_value = lm[thread_index][thread_index]; if (!IsZero(diagonal_value)) { // Only for non-singular values and values inside the matrix - DivideReal(inverted_diagonal, inverted_diagonal, diagonal_value); + real constant_one; + SetToOne(constant_one); + DivideReal(inverted_diagonal, constant_one, diagonal_value); } } lm[thread_index][thread_index] = inverted_diagonal; |