Commit 56d3cc60 authored by Tomas Sobotik's avatar Tomas Sobotik
Browse files

Parallel narrow band working, still needs a bit more optimisation.

parent 4dc39335

examples/fast-sweeping/main.h

+2 −2
Original line number Diff line number Diff line
@@ -17,9 +17,9 @@ 


#include "MainBuildConfig.h"

	//for HOST versions:

#include "tnlFastSweeping.h"

//#include "tnlFastSweeping.h"

	//for DEVICE versions:

//#include "tnlFastSweeping_CUDA.h"

#include "tnlFastSweeping_CUDA.h"

#include "fastSweepingConfig.h"

#include <solvers/tnlBuildConfigTags.h>

examples/fast-sweeping/tnlFastSweeping2D_CUDA_v4_impl.h

+83 −76

File changed.

Preview size limit exceeded, changes collapsed.

examples/fast-sweeping/tnlFastSweeping3D_CUDA_impl.h

+2 −2
Original line number Diff line number Diff line
@@ -227,8 +227,8 @@ bool tnlFastSweeping< tnlGrid< 3,MeshReal, Device, MeshIndex >, Real, Index > :: 
	Entity.refresh();

	int gid = Entity.getIndex();



	if(abs(cudaDofVector[gid]) < 1.8*h)

		cudaDofVector2[gid] = cudaDofVector[gid];

	if(abs(cudaDofVector[gid]) < 1.0*h)

		cudaDofVector2[gid] = 0.5*h;//cudaDofVector[gid];

	else

		cudaDofVector2[gid] = INT_MAX*Sign(cudaDofVector[gid]);

examples/hamilton-jacobi-parallel/tnlParallelEikonalSolver2D_impl.h

+1 −1
Original line number Diff line number Diff line
@@ -25,7 +25,7 @@ template< typename SchemeHost, typename SchemeDevice, typename Device> 
tnlParallelEikonalSolver<2,SchemeHost, SchemeDevice, Device, double, int>::tnlParallelEikonalSolver()

{

	cout << "a" << endl;

	this->device = tnlHostDevice;  /////////////// tnlCuda Device --- vypocet na GPU, tnlHostDevice   ---    vypocet na CPU

	this->device = tnlCudaDevice;  /////////////// tnlCuda Device --- vypocet na GPU, tnlHostDevice   ---    vypocet na CPU



#ifdef HAVE_CUDA

	if(this->device == tnlCudaDevice)

examples/hamilton-jacobi-parallel/tnlParallelEikonalSolver3D_impl.h

+1 −1
Original line number Diff line number Diff line
@@ -25,7 +25,7 @@ template< typename SchemeHost, typename SchemeDevice, typename Device> 
tnlParallelEikonalSolver<3,SchemeHost, SchemeDevice, Device, double, int>::tnlParallelEikonalSolver()

{

	cout << "a" << endl;

	this->device = tnlCudaDevice;  /////////////// tnlCuda Device --- vypocet na GPU, tnlHostDevice   ---    vypocet na CPU

	this->device = tnlHostDevice;  /////////////// tnlCuda Device --- vypocet na GPU, tnlHostDevice   ---    vypocet na CPU



#ifdef HAVE_CUDA

	if(this->device == tnlCudaDevice)