about 50x speedup by precomputing norm2(A) (obvious)