Line data    Source code

       1   1968757377 :   fsign = 1.d0
       2   1968757377 :   t_val = time
       3   1968757377 :   IF ( t_val .LT. 0.d0 ) THEN
       4              :     ! if t < 0 this%F(t,activeF) = -this%F(t+beta,activeF)
       5    990402494 :     t_val = t_val + this%beta
       6    990402494 :     fsign  = -1.d0 
       7              :   END IF
       8              : 
       9   1968757377 :   ratioh = t_val * this%inv_dt + 1.d0 ! t / delta t
      10   1968757377 :   index = INT(ratioh)
      11   1968757377 :   ratioh = ratioh - DBLE(index)
      12              : 
      13              :   !  interpolation first order
      14              :   ! index + 1 can be sizeHybrid + 1 => value is 0 (cf. BathOperator_init routine)
      15              :   hybrid = fsign * (this%F(index,activeF) + &
      16   1968757377 :       ratioh * ( this%F(index+1,activeF) - this%F(index,activeF) ))