Commit 8ae96a98 authored by Tomas Sobotik's avatar Tomas Sobotik
Browse files

Dimension for added for Fast-sweeping

parent 9437c8d7

examples/fast-sweeping/fastSweepingConfig.h

+1 −0
Original line number Diff line number Diff line
@@ -29,6 +29,7 @@ class fastSweepingConfig 
         config.addDelimiter( "Parallel Eikonal solver settings:" );

         config.addEntry        < tnlString > ( "problem-name", "This defines particular problem.", "fast-sweeping" );

         config.addRequiredEntry        < tnlString > ( "initial-condition", "Initial condition for solver");

         config.addRequiredEntry        < tnlString > ( "dim", "Dimension of problem.");

         config.addEntry       < tnlString > ( "mesh", "Name of mesh.", "mesh.tnl" );

         config.addEntry       < tnlString > ( "exact-input", "Are the function values near the curve equal to the SDF? (yes/no)", "no" );

      }

examples/fast-sweeping/main.h

+28 −8
Original line number Diff line number Diff line
@@ -43,6 +43,23 @@ int main( int argc, char* argv[] ) 
   if( ! parseCommandLine( argc, argv, configDescription, parameters ) )

      return false;



   const tnlString& dim = parameters.getParameter< tnlString >( "dim" );



   if(dim == "2")

   {

		tnlFastSweeping<tnlGrid<2,double,tnlHost, int>, double, int> solver;

		if(!solver.init(parameters))

	   {

			cerr << "Solver failed to initialize." << endl;

			return EXIT_FAILURE;

	   }

		checkCudaDevice;

	   cout << "-------------------------------------------------------------" << endl;

	   cout << "Starting solver..." << endl;

	   solver.run();

   }

   else    if(dim == "3")

   {

		tnlFastSweeping<tnlGrid<3,double,tnlHost, int>, double, int> solver;

		if(!solver.init(parameters))

	   {
@@ -53,6 +70,9 @@ int main( int argc, char* argv[] ) 
	   cout << "-------------------------------------------------------------" << endl;

	   cout << "Starting solver..." << endl;

	   solver.run();

   }

   else

	   cerr << "Unsopported no. of dimensions: " << dim << "!" << endl;





   time(&stop);