# 1 "<stdin>"
# 1 "<built-in>"
# 1 "<command-line>"
# 1 "<stdin>"
module conv
  use helix
  use compress
  implicit none
  contains
  function conv1( aa, bb) result ( ss)
    type( filter) :: aa, bb, ss
    integer :: na, nb, ns, ia, ib, i
    na = size (aa%lag)
    nb = size (bb%lag)
    ns = maxval (aa%lag) + maxval (bb%lag)
    call allocatehelix (ss, ns)
    ss%lag = (/ (i, i = 1, ns) /)
    do ia = 1, na
      i = aa%lag (ia)
      ss%flt (i) = aa%flt (ia)
    end do
    do ib = 1, nb
      i = bb%lag (ib)
      ss%flt (i) = ss%flt (i) + bb%flt (ib)
    end do
    do ia = 1, na
      do ib = 1, nb
        i = aa%lag (ia) + bb%lag (ib)
        ss%flt (i) = ss%flt (i) + aa%flt (ia) * bb%flt (ib)
      end do
    end do
    call resize (ss)
  end function
end module