Commit fde31a70 authored by Tomas Sobotik's avatar Tomas Sobotik
Browse files

3D parallel iterative solver almost fully fixed v2

parent 337fb813

examples/hamilton-jacobi-parallel/tnlParallelEikonalSolver3D_impl.h

+1 −1
Original line number Diff line number Diff line
@@ -1013,7 +1013,7 @@ void tnlParallelEikonalSolver<3,SchemeHost, SchemeDevice, Device, double, int>:: 
	  if(abs(fu) > 0.0)

		  sharedTau[l]=abs(cfl/fu);



	  if(u[l]*fu < 0.0 && abs(fu*sharedTau[l]) >abs(u[l])) sharedTau[l] = 0.5*abs(u[l]/fu) + this->subMesh.template getSpaceStepsProducts< 1, 0, 0 >()*this->subMesh.template getSpaceStepsProducts< 1, 0, 0 >();

	 /* if(u[l]*fu < 0.0 && abs(fu*sharedTau[l]) >abs(u[l])) sharedTau[l] = 0.9*abs(u[l]/fu)/* + this->subMesh.template getSpaceStepsProducts< 1, 0, 0 >()*this->subMesh.template getSpaceStepsProducts< 1, 0, 0 >()*/;



      if(l == 0)

      {

src/operators/godunov-eikonal/parallelGodunovEikonal3D_impl.h

+4 −4
Original line number Diff line number Diff line
@@ -377,9 +377,9 @@ Real parallelGodunovEikonalScheme< tnlGrid< 3, MeshReal, Device, MeshIndex >, Re 


	   d = 1.0 - sqrt(xf*xf + xb*xb + yf*yf + yb*yb + zf*zf + zb*zb)*ihx; /*upwind*/



//	   if(Sign(d) > 0.0 )

//		   return Sign(u[cellIndex])*d;

//	   else

	   if(Sign(d) > 0.0 )

		   return Sign(u[cellIndex])*d;

	   else

		   return signui*d;

}