burg2 - Burg 2D convolution
call burg2(adj,add,lag1,lag2,data,n1,n2,aa,a1,a2,residual)
adjoint=1, forward=0
add=0 add to current values, 0=erase
lag in first dimension
lag in second dimension
real(n1,n2)
data
data length along axis=1
data length along axis=2
real(a1,a2)
filter
filter length axis=1
filter length axis=2
real(n1,n2*2)
residual
Burg 2-D conv with (a1,a2) filter (monoplane annihilator if a2=2)
output residual partitioned into normal and backward parts. output adjoint to FILTER. output residual(,) aligns with data(,) at filter coef aa(lag1,lag2)
burg
sepfilterf90