#$
#$=head1 NAME
#$
#$mshelicon - convolution using mutli scale helix filters
#$
#$=head1 SYNOPSIS
#$
#$Initializer - C<call mshelicon_init(aa)>
#$
#$Operator    - C<ierr=mshelicon_lop(adj,add,xx,yy)>
#$
#$=head1 PARAMETERS
#$
#$=over 4
#$
#$=item aa - type(msfilter)  
#$
#$      helix filter to perform convolution with
#$
#$=item adj,add,xx,yy - 
#$
#$      standard operators parameters
#$
#$=back
#$
#$=head1 DESCRIPTION
#$
#$ Multi-scale onvolution, inverse to deconvolution.
#$ Requires the filter be causal with an implicit "1." at the onset.
#$
#$
#$=head1 SEE ALSO
#$
#$L<mshelix>,L<mshconest>,L<helicon>,L<nhelicon>
#$
#$=head1 LIBRARY
#$
#$B<geef90>
#$
#$=cut
#$

module mshelicon {                                  # Multi-scale convolution
  use mshelix
  use helicon
  integer           :: nx, ns
  type( msfilter)   :: aa
#%  _init (nx, ns, aa)
#%  _lop ( xx( nx),  yy( nx, ns))
  integer is, stat1
  do is = 1, ns { 
     call helicon_init( onescale( is, aa))
     stat1 = helicon_lop( adj, .true., xx, yy(:,is))
     }
}