From e88cde26971d724869134828d3974ab23a9e4707 Mon Sep 17 00:00:00 2001 From: David Minton Date: Tue, 25 May 2021 13:45:17 -0400 Subject: [PATCH] Simplified objective function --- src/symba/symba_frag_pos.f90 | 6 +----- 1 file changed, 1 insertion(+), 5 deletions(-) diff --git a/src/symba/symba_frag_pos.f90 b/src/symba/symba_frag_pos.f90 index 9b61636b0..23f1e7eb7 100644 --- a/src/symba/symba_frag_pos.f90 +++ b/src/symba/symba_frag_pos.f90 @@ -730,11 +730,7 @@ function ke_objective_function(v_r_mag, v_r_unit, v_t_mag, v_t_unit, m_frag, vco fval = fval + 0.5_DP * m_frag(i) * dot_product(v_shift(:, i), v_shift(:, i)) end do ! The following ensures that fval = 0 is a local minimum, which is what the BFGS method is searching for - if (fval < 0.0_DP) then - fval = fval**2 - else - fval = fval**2 - end if + fval = fval**2 return