Changed order of gradient for the narrow band and default value of get_upwind_gradient to 1.

src/FiniteDifference/Upwind_gradient.hpp

@@ -210,7 +210,7 @@ static bool ghost_width_is_sufficient(gridtype & grid, size_t required_width)
  * @param grid Grid, on which the gradient should be computed.
  */
 template <size_t Field_in, size_t Sign, size_t Gradient_out, typename gridtype>
-void get_upwind_gradient(gridtype & grid, const size_t order=5, const bool one_sided_BC=true)
+void get_upwind_gradient(gridtype & grid, const size_t order=1, const bool one_sided_BC=true)
 {
 	grid.template ghost_get<Field_in>();
 	grid.template ghost_get<Sign>();

src/level_set/redistancing_Sussman/NarrowBand.hpp

@@ -340,7 +340,7 @@ private:
 		copy_gridTogrid<Phi_SDF, Phi_SDF_temp>(grid_in, g_temp); // Copy Phi_SDF from the input grid to the temorary grid
 		init_sign_prop<Phi_SDF_temp, Phi_sign_temp>(g_temp); // initialize Phi_sign_temp with the sign of the
 		// input Phi_SDF
-		get_upwind_gradient<Phi_SDF_temp, Phi_sign_temp, Phi_grad_temp>(g_temp);   // Get initial gradients
+		get_upwind_gradient<Phi_SDF_temp, Phi_sign_temp, Phi_grad_temp>(g_temp, 1, true);   // Get initial gradients
 	}
 	/**@brief Checks if a value for Phi_SDF lays within the narrow band.
 	 *