diff --git a/src/TNL/Matrices/ChunkedEllpack_impl.h b/src/TNL/Matrices/ChunkedEllpack_impl.h
index 29ebfc4157ffc856f809e10de3b85c65384c02f4..751fa1c21fb4a14956520dc57386c2d73dc72306 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 );