From 6dd89c30172bb8344a8014a374d867e8d590b309 Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?Luk=C3=A1=C5=A1=20=C4=8Cejka?= <cejkaluk@fjfi.cvut.cz> Date: Thu, 13 Dec 2018 22:47:20 +0100 Subject: [PATCH] Changed Double type casting to roundUpDivision --- src/TNL/Matrices/ChunkedEllpack_impl.h | 3 +-- 1 file changed, 1 insertion(+), 2 deletions(-) diff --git a/src/TNL/Matrices/ChunkedEllpack_impl.h b/src/TNL/Matrices/ChunkedEllpack_impl.h index 29ebfc4157..751fa1c21f 100644 --- a/src/TNL/Matrices/ChunkedEllpack_impl.h +++ b/src/TNL/Matrices/ChunkedEllpack_impl.h @@ -198,8 +198,7 @@ bool ChunkedEllpack< Real, Device, Index >::setSlice( ConstCompressedRowLengthsV // will be zero and make the assertion fail ( https://stackoverflow.com/questions/33273359/in-c-using-the-ceil-a-division-is-not-working ). // To fix this, typecast them to ( float ), instead of ( RealType ) maxChunkInSlice = max( maxChunkInSlice, - ceil( ( double ) rowLengths[ i ] / - ( double ) this->rowToChunkMapping[ i ] ) ); + ceil( roundUpDivision( rowLengths[ i ], this->rowToChunkMapping[ i ] ) ); } TNL_ASSERT( maxChunkInSlice > 0, std::cerr << " maxChunkInSlice = " << maxChunkInSlice << std::endl ); -- GitLab