OppLsqr = NewOperator("Lsqr", NF, IndexUp, IndexDn)