function init_doubleWeight(csf_i, sOrb) result(doubleWeight)
! obj has the same structure as the semi-start weight, reuse them!
type(CSF_Info_t), intent(in) :: csf_i
integer, intent(in) :: sOrb
type(WeightObj_t) :: doubleWeight
character(*), parameter :: this_routine = "init_doubleWeight"
ASSERT(sOrb > 0 .and. sOrb <= nSpatOrbs)
doubleWeight%dat%F = endFx(csf_i, sOrb)
doubleWeight%dat%G = endGx(csf_i, sOrb)
doubleWeight%proc%minus => getMinus_double
doubleWeight%proc%plus => getPlus_double
doubleWeight%proc%zero => getZero_double
doubleWeight%initialized = .true.
end function init_doubleWeight