c-----pre-tabulation of some elements for fast-Hessian calculation
      real s(0:size)
      real al(dimD,0:size),alsq(dimD,0:size),alcub(dimD,0:size),
     &     alfour(dimD,0:size),alfive(dimD,0:size)
      real sinal(dimD,0:size),cosal(dimD,0:size),sinal_al(dimD,0:size),
     &     sinal_alsq(dimD,0:size),sinal_alcub(dimD,0:size),
     &     sinal_alfour(dimD,0:size),sinal_alfive(dimD,0:size),
     &     cosal_al(dimD,0:size),cosal_alsq(dimD,0:size),
     &     cosal_alcub(dimD,0:size),cosal_alfour(dimD,0:size)
      real ft_c(dimB,0:size),strip,striptre,stripD,stripB
      COMMON /PRETAB/s,al,alsq,alcub,alfour,alfive,sinal,
     &              cosal,sinal_al,sinal_alsq,sinal_alcub,
     &              sinal_alfour,sinal_alfive,cosal_al,
     &              cosal_alsq,cosal_alcub,cosal_alfour,ft_c,
     &              striptre,stripD