Commit d0e3ce8d authored by Jakub Klinkovský's avatar Jakub Klinkovský
Browse files

Fixed bug in the Euler solver 

I don't understand how it could work on the small grids...
parent 986aaeab

src/TNL/Solvers/ODE/Euler_impl.h

+2 −1
Original line number Diff line number Diff line
@@ -193,8 +193,9 @@ void Euler< Problem > :: computeNewTimeLevel( DofVectorPointer& u, 
         const IndexType sharedMemory = cudaBlockSize.x * sizeof( RealType );

         const IndexType gridOffset = gridIdx * threadsPerGrid;

         const IndexType currentSize = min( size - gridOffset, threadsPerGrid );

         const IndexType currentGridSize = Devices::Cuda::getNumberOfBlocks( currentSize, cudaBlockSize.x );



         updateUEuler<<< cudaBlocks, cudaBlockSize, sharedMemory >>>( currentSize,

         updateUEuler<<< currentGridSize, cudaBlockSize, sharedMemory >>>( currentSize,

                                                                      tau,

                                                                      &_k1[ gridOffset ],

                                                                      &_u[ gridOffset ],