anchor.diff:                 ra(i), ra(i+1), phinew(i) % p, crA(pa(i)), curcr, curcr, phitmp)
                  call zinit_sfun(rnew, rnew, ra(i+1), ry(i+1), ry(i+1), phitmp, phinew(i+1)%p)
                  anorm = znormest(rnew * ry(i+1), 4, zsfun_matvec, zsfun_matvec_transp)
  -               ! anorm_t = znorm_true(rnew * ry(i+1), zsfun_matvec, zsfun_matvec_transp)
  -               ! print *, 'anorm, anorm_t = ', anorm, anorm_t
                  call zexp_mv(rnew * ry(i+1), order, -tau0, R, Stmp, eps, anorm, zsfun_matvec)
                  call zgemm('n', 'n', rnew, n(i+1) * ry(i+2), ry(i+1), ONE, &
                  Stmp, rnew, crnew(i+1) % p, ry(i+1), &
  
positive.code:     integer :: max_phi_size, max_core_size, max_R_size
      integer :: order
      real(8) :: ermax, min_res
      real(8) anorm, anorm_t
      complex(8) ZERO, ONE
      parameter( ZERO=(0.0d0,0.0d0), ONE=(1.0d0,0.0d0) )
  
  
negative.code:                ra(i), ra(i + 1), ry(i) * n(i) * ry(i + 1), ry(i) * n(i) * ry(i + 1))
             call zinit_bfun_main(phinew(i) % p, crA(pa(i):pa(i+1)-1), phinew(i+1) % p)
             anorm = znormest(ry(i) * n(i) * ry(i+1), 4, zmatvec, zmatvec_transp)
             call zexp_mv(ry(i) * n(i) * ry(i+1), order, tau0, &
                         crnew(i)%p, curcr, eps, anorm, zmatvec)
             if ( i > 1 ) then
  