diff --git a/examples/hamilton-jacobi-parallel/tnlParallelEikonalSolver_impl.h b/examples/hamilton-jacobi-parallel/tnlParallelEikonalSolver_impl.h index b3e09b7294df4f7a1051fcb039d34868e497fd33..9101c846ff04138cd3259abb25d29e126673b3e0 100644 --- a/examples/hamilton-jacobi-parallel/tnlParallelEikonalSolver_impl.h +++ b/examples/hamilton-jacobi-parallel/tnlParallelEikonalSolver_impl.h @@ -800,10 +800,13 @@ void tnlParallelEikonalSolver<SchemeHost, SchemeDevice, Device, double, int>::ru { // u[l] = value; if(boundaryCondition==FROM_NORTH) - u[l] = ; + u[l] = u[i + blockDim.x*(blockDim.y-1)]; else if(boundaryCondition==FROM_SOUTH) + u[l] = u[i]; else if(boundaryCondition==FROM_EAST) + u[l] = u[blockDim.x-1 + blockDim.x*j]; else if(boundaryCondition==FROM_WEST) + u[l] = u[blockDim.x*j]; } }