module peftc {
  use tcai1
  use tcaf1
  integer, private  :: nx
#% _init (aa, yy)
  real, dimension (:), pointer :: aa, yy
  nx = size (yy)
  call tcai1_init (aa)
  call tcaf1_init (yy)
#% _lop (x, r)
  integer stat1, stat2
  stat1 = tcai1_lop (adj, .true., x(:nx),   r)
  stat2 = tcaf1_lop (adj, .true., x(nx+1:), r)
}