diff --git a/src/symba/symba_frag_pos.f90 b/src/symba/symba_frag_pos.f90
index 0233aa7ed..087c345e4 100644
--- a/src/symba/symba_frag_pos.f90
+++ b/src/symba/symba_frag_pos.f90
@@ -701,7 +701,7 @@ subroutine set_fragment_radial_velocities(lerr)
       ! Arguments
       logical,                intent(out)   :: lerr
       ! Internals
-      real(DP), parameter                   :: TOL = 2e-8_DP
+      real(DP), parameter                   :: TOL = 1e-10_DP
       integer(I4B)                          :: i, j
       real(DP), dimension(:), allocatable   :: v_r_initial, v_r_sigma
       real(DP), dimension(:,:), allocatable :: v_r
diff --git a/src/util/util_minimize_bfgs.f90 b/src/util/util_minimize_bfgs.f90
index 84f63499e..5f820fec1 100644
--- a/src/util/util_minimize_bfgs.f90
+++ b/src/util/util_minimize_bfgs.f90
@@ -29,7 +29,7 @@ function util_minimize_bfgs(f, N, x0, eps, lerr) result(x1)
    ! Internals
    integer(I4B) ::  i, j, k, l, conv, num
    integer(I4B), parameter :: MAXLOOP = 1000 !! Maximum number of loops before method is determined to have failed 
-   real(DP), parameter     :: graddelta = 1e-5_DP !! Delta x for gradient calculations
+   real(DP), parameter     :: graddelta = 1e-4_DP !! Delta x for gradient calculations
    real(DP), dimension(N) :: S               !! Direction vectors 
    real(DP), dimension(N,N) :: H             !! Approximated inverse Hessian matrix 
    real(DP), dimension(N) :: grad1           !! gradient of f