module nicaf1 { # Nonstationary Internal Convolution, Adjoint is the Filter real, dimension (:), pointer :: yy integer :: n, na #% _init (yy, n, na) #% _lop (aa (n,na), rr (n)) integer a, r, y do a = 1, na { do r = na, n { y = r - a + 1 if( adj) aa (r, a) += rr (r) * yy (y) else rr (r) += aa (r, a) * yy (y) }} }