Loading src/Benchmarks/HeatEquation/HeatEquationBenchmarkProblem_impl.h +1 −1 Original line number Diff line number Diff line Loading @@ -646,7 +646,6 @@ HeatEquationBenchmarkProblem< Mesh, BoundaryCondition, RightHandSide, Differenti applyBoundaryConditions( const RealType& time, DofVectorPointer& uDofs ) { const MeshPointer& mesh = this->getMesh(); if( this->cudaKernelType == "templated" ) { this->bindDofs( uDofs ); Loading Loading @@ -698,6 +697,7 @@ applyBoundaryConditions( const RealType& time, #endif userData.boundaryConditions = &this->boundaryConditionPointer.template getData< Devices::Cuda >(); Meshes::Traverser< MeshType, Cell > meshTraverser; const MeshPointer& mesh = this->getMesh(); // */ /*meshTraverser.template processBoundaryEntities< BoundaryEntitiesProcessor > ( mesh, Loading src/Benchmarks/ODESolvers/Merson.hpp +2 −2 Original line number Diff line number Diff line Loading @@ -244,7 +244,7 @@ void Merson< Problem, SolverMonitor >::computeKFunctions( DofVectorPointer& u, RealType* _k2 = k2->getData(); RealType* _k3 = k3->getData(); RealType* _k4 = k4->getData(); RealType* _k5 = k5->getData(); //RealType* _k5 = k5->getData(); RealType* _kAux = kAux->getData(); RealType* _u = u->getData(); Loading Loading @@ -353,7 +353,6 @@ typename Problem :: RealType Merson< Problem, SolverMonitor >::computeError( con const RealType* _k3 = k3->getData(); const RealType* _k4 = k4->getData(); const RealType* _k5 = k5->getData(); RealType* _kAux = kAux->getData(); RealType eps( 0.0 ), maxEps( 0.0 ); if( std::is_same< DeviceType, Devices::Host >::value ) Loading Loading @@ -383,6 +382,7 @@ typename Problem :: RealType Merson< Problem, SolverMonitor >::computeError( con if( std::is_same< DeviceType, Devices::Cuda >::value ) { #ifdef HAVE_CUDA RealType* _kAux = kAux->getData(); dim3 cudaBlockSize( 512 ); const IndexType cudaBlocks = Cuda::getNumberOfBlocks( size, cudaBlockSize.x ); const IndexType cudaGrids = Cuda::getNumberOfGrids( cudaBlocks ); Loading src/Examples/transport-equation/transportEquationProblem_impl.h +1 −1 Original line number Diff line number Diff line Loading @@ -188,7 +188,7 @@ getExplicitUpdate( const RealType& time, */ const MeshPointer& mesh = this->getMesh(); typedef typename MeshType::Cell Cell; int count = ::sqrt(mesh->template getEntitiesCount< Cell >()); //int count = ::sqrt(mesh->template getEntitiesCount< Cell >()); this->bindDofs( _u ); Solvers::PDE::ExplicitUpdater< Mesh, MeshFunctionType, DifferentialOperator, BoundaryCondition, RightHandSide > explicitUpdater; Pointers::SharedPointer< MeshFunctionType > u( mesh, _u ); Loading src/TNL/Containers/Segments/ChunkedEllpackView.hpp +1 −1 Original line number Diff line number Diff line Loading @@ -622,7 +622,7 @@ segmentsReductionKernel( IndexType gridIdx, const IndexType sliceOffset = sliceInfo.pointer; const IndexType chunkSize = sliceInfo.chunkSize; const IndexType chunkIdx = sliceIdx * chunksInSlice + threadIdx.x; //const IndexType chunkIdx = sliceIdx * chunksInSlice + threadIdx.x; bool compute( true ); if( Organization == RowMajorOrder ) Loading src/TNL/Containers/Segments/SlicedEllpackView.hpp +2 −2 Original line number Diff line number Diff line Loading @@ -253,7 +253,7 @@ forSegments( IndexType first, IndexType last, Function& f, Args... args ) const auto l = [=] __cuda_callable__ ( const IndexType segmentIdx, Args... args ) mutable { const IndexType sliceIdx = segmentIdx / SliceSize; const IndexType segmentInSliceIdx = segmentIdx % SliceSize; const IndexType segmentSize = sliceSegmentSizes_view[ sliceIdx ]; //const IndexType segmentSize = sliceSegmentSizes_view[ sliceIdx ]; const IndexType begin = sliceOffsets_view[ sliceIdx ] + segmentInSliceIdx; const IndexType end = sliceOffsets_view[ sliceIdx + 1 ]; IndexType localIdx( 0 ); Loading Loading @@ -312,7 +312,7 @@ segmentsReduction( IndexType first, IndexType last, Fetch& fetch, const Reductio auto l = [=] __cuda_callable__ ( const IndexType segmentIdx, Args... args ) mutable { const IndexType sliceIdx = segmentIdx / SliceSize; const IndexType segmentInSliceIdx = segmentIdx % SliceSize; const IndexType segmentSize = sliceSegmentSizes_view[ sliceIdx ]; //const IndexType segmentSize = sliceSegmentSizes_view[ sliceIdx ]; const IndexType begin = sliceOffsets_view[ sliceIdx ] + segmentInSliceIdx; const IndexType end = sliceOffsets_view[ sliceIdx + 1 ]; RealType aux( zero ); Loading Loading
src/Benchmarks/HeatEquation/HeatEquationBenchmarkProblem_impl.h +1 −1 Original line number Diff line number Diff line Loading @@ -646,7 +646,6 @@ HeatEquationBenchmarkProblem< Mesh, BoundaryCondition, RightHandSide, Differenti applyBoundaryConditions( const RealType& time, DofVectorPointer& uDofs ) { const MeshPointer& mesh = this->getMesh(); if( this->cudaKernelType == "templated" ) { this->bindDofs( uDofs ); Loading Loading @@ -698,6 +697,7 @@ applyBoundaryConditions( const RealType& time, #endif userData.boundaryConditions = &this->boundaryConditionPointer.template getData< Devices::Cuda >(); Meshes::Traverser< MeshType, Cell > meshTraverser; const MeshPointer& mesh = this->getMesh(); // */ /*meshTraverser.template processBoundaryEntities< BoundaryEntitiesProcessor > ( mesh, Loading
src/Benchmarks/ODESolvers/Merson.hpp +2 −2 Original line number Diff line number Diff line Loading @@ -244,7 +244,7 @@ void Merson< Problem, SolverMonitor >::computeKFunctions( DofVectorPointer& u, RealType* _k2 = k2->getData(); RealType* _k3 = k3->getData(); RealType* _k4 = k4->getData(); RealType* _k5 = k5->getData(); //RealType* _k5 = k5->getData(); RealType* _kAux = kAux->getData(); RealType* _u = u->getData(); Loading Loading @@ -353,7 +353,6 @@ typename Problem :: RealType Merson< Problem, SolverMonitor >::computeError( con const RealType* _k3 = k3->getData(); const RealType* _k4 = k4->getData(); const RealType* _k5 = k5->getData(); RealType* _kAux = kAux->getData(); RealType eps( 0.0 ), maxEps( 0.0 ); if( std::is_same< DeviceType, Devices::Host >::value ) Loading Loading @@ -383,6 +382,7 @@ typename Problem :: RealType Merson< Problem, SolverMonitor >::computeError( con if( std::is_same< DeviceType, Devices::Cuda >::value ) { #ifdef HAVE_CUDA RealType* _kAux = kAux->getData(); dim3 cudaBlockSize( 512 ); const IndexType cudaBlocks = Cuda::getNumberOfBlocks( size, cudaBlockSize.x ); const IndexType cudaGrids = Cuda::getNumberOfGrids( cudaBlocks ); Loading
src/Examples/transport-equation/transportEquationProblem_impl.h +1 −1 Original line number Diff line number Diff line Loading @@ -188,7 +188,7 @@ getExplicitUpdate( const RealType& time, */ const MeshPointer& mesh = this->getMesh(); typedef typename MeshType::Cell Cell; int count = ::sqrt(mesh->template getEntitiesCount< Cell >()); //int count = ::sqrt(mesh->template getEntitiesCount< Cell >()); this->bindDofs( _u ); Solvers::PDE::ExplicitUpdater< Mesh, MeshFunctionType, DifferentialOperator, BoundaryCondition, RightHandSide > explicitUpdater; Pointers::SharedPointer< MeshFunctionType > u( mesh, _u ); Loading
src/TNL/Containers/Segments/ChunkedEllpackView.hpp +1 −1 Original line number Diff line number Diff line Loading @@ -622,7 +622,7 @@ segmentsReductionKernel( IndexType gridIdx, const IndexType sliceOffset = sliceInfo.pointer; const IndexType chunkSize = sliceInfo.chunkSize; const IndexType chunkIdx = sliceIdx * chunksInSlice + threadIdx.x; //const IndexType chunkIdx = sliceIdx * chunksInSlice + threadIdx.x; bool compute( true ); if( Organization == RowMajorOrder ) Loading
src/TNL/Containers/Segments/SlicedEllpackView.hpp +2 −2 Original line number Diff line number Diff line Loading @@ -253,7 +253,7 @@ forSegments( IndexType first, IndexType last, Function& f, Args... args ) const auto l = [=] __cuda_callable__ ( const IndexType segmentIdx, Args... args ) mutable { const IndexType sliceIdx = segmentIdx / SliceSize; const IndexType segmentInSliceIdx = segmentIdx % SliceSize; const IndexType segmentSize = sliceSegmentSizes_view[ sliceIdx ]; //const IndexType segmentSize = sliceSegmentSizes_view[ sliceIdx ]; const IndexType begin = sliceOffsets_view[ sliceIdx ] + segmentInSliceIdx; const IndexType end = sliceOffsets_view[ sliceIdx + 1 ]; IndexType localIdx( 0 ); Loading Loading @@ -312,7 +312,7 @@ segmentsReduction( IndexType first, IndexType last, Fetch& fetch, const Reductio auto l = [=] __cuda_callable__ ( const IndexType segmentIdx, Args... args ) mutable { const IndexType sliceIdx = segmentIdx / SliceSize; const IndexType segmentInSliceIdx = segmentIdx % SliceSize; const IndexType segmentSize = sliceSegmentSizes_view[ sliceIdx ]; //const IndexType segmentSize = sliceSegmentSizes_view[ sliceIdx ]; const IndexType begin = sliceOffsets_view[ sliceIdx ] + segmentInSliceIdx; const IndexType end = sliceOffsets_view[ sliceIdx + 1 ]; RealType aux( zero ); Loading