f
          the transfer function matrix G, stored in a similar way
          as the real parts of the zeros.

  POLESR  (output) DOUBLE PRECISION array, dimension (P*M*NPZ)
          This array contains the real parts of the poles of the
          transfer function matrix G, stored in the same way as
          the zeros. Note that only the first NP(i,j) entries are
          initialized for the (i,j) transfer function.

  POLESI  (output) DOUBLE PRECISION array, dimension (P*M*NPZ)
          This array con