NAME

burg2 - Burg 2D convolution


SYNOPSIS

call burg2(adj,add,lag1,lag2,data,n1,n2,aa,a1,a2,residual)


INPUT PARAMETERS

adj - integer
      adjoint=1, forward=0
add - integer
      add=0 add to current values, 0=erase
lag1 - integer
      lag in first dimension
lag2 - integer
      lag in second dimension
data - real(n1,n2)
      data
n1 - integer
      data length along axis=1
n2 - integer
      data length along axis=2
aa - real(a1,a2)
      filter
a1 - integer
      filter length axis=1
a2 - integer
      filter length axis=2
residual - real(n1,n2*2)
      residual


DESCRIPTION

 Burg 2-D conv with (a1,a2) filter (monoplane annihilator if a2=2)


COMMENTS

       output residual partitioned into normal and backward parts.
       output adjoint to FILTER.
       output residual(,) aligns with data(,) at filter coef aa(lag1,lag2)


SEE ALSO

burg


LIBRARY

sepfilterf90