Skip to content
This repository was archived by the owner on Aug 28, 2024. It is now read-only.

Commit

Permalink
Browse files Browse the repository at this point in the history
Use fixed size array for the objection function vector
  • Loading branch information
daminton committed May 13, 2021
1 parent 9d8bc1a commit edd784a
Showing 1 changed file with 1 addition and 2 deletions.
3 changes: 1 addition & 2 deletions src/symba/symba_frag_pos.f90
Original file line number Diff line number Diff line change
Expand Up @@ -410,8 +410,7 @@ function symba_frag_pos_ke_objective_function(v_r_mag, m_frag, v_r_unit, L_lin_t
real(DP) :: fnorm !! The objective function result: norm of the vector composed of the tangential momentum and energy
!! Minimizing this brings us closer to our objective
! Internals
integer(I4B), parameter :: Neqs = 4
real(DP), dimension(:), allocatable :: f_vec
real(DP), dimension(4) :: f_vec

f_vec(1) = sum(m_frag(:) * v_r_mag(:) * v_r_unit(1,:))
f_vec(2) = sum(m_frag(:) * v_r_mag(:) * v_r_unit(2,:))
Expand Down

0 comments on commit edd784a

Please sign in to comment.