function opt_func(p)

   use opt_com_mod, only: demerit, rp
   implicit none

   real(rp), dimension(:), intent(in) :: p
   real(rp) opt_func
   integer end_flag

   opt_func = demerit(p, end_flag)

end function opt_func