function [ w ] = AtAv(v, K, H, R, dt, nt, mt, rho )
%
% compute A'*A*v
%
% Matthias Heinkenschloss
% Jan 9, 2017 (revised Jan 14, 2020)
%
z  = Av(v, K, H, R, dt, nt, mt, rho );
w  = Atz(z, K, H, R, dt, nt, mt, rho );


