ids.hs

inner_decentral_sino_st ::  ProblemParams -> (Double,Double,Double) -> (Double,Double,Double)
inner_decentral_sino_st pp (simin, simed, simax) = (simin', simed', simax')
  where
    p  = fromIntegral $ _pp_p pp :: Double
    t  = fromIntegral $ _pp_t pp :: Double
    h  = fromIntegral $ _pp_h pp :: Double
    myf_val = myfunN1 simed pp
    thresh  = (p*t)/(p+h) 

    if myf_val < thresh then
      (simin',  simax') = (simed, simax)
      else
      (simin', simax') = (simin, simed)

    simed' = average simin' simax'