module signoi {
  use helicon
  use polydiv
  use solver_prc_mod
  use cgstep_mod
  integer, private                      :: i
  integer                               :: niter, nd
  real                                  :: eps
  type( filter), dimension( :), pointer :: nn, ss
  real, dimension (nd), allocatable     :: dd
#% _init (nn, ss, niter, nd, eps)
i = 0
#% _lop (data, sign)
    integer stat1; i = i + 1
    call helicon_init (nn (i))
    call polydiv_init (nd, ss (i)) 
    stat1 = helicon_lop (.false., .false., data, dd)
    call solver_prc(m=sign, d=dd, Fop=helicon_lop, stepper=cgstep,
    niter=niter, Sop=polydiv_lop, nSop=nd, eps=eps)
    call cgstep_close()
#% _close
    call polydiv_close()
}