Download this file

2 lines (2 with data), 60 Bytes

1
2
function [y] = precondUpper(r,U,D)
y = U \ (D .* (U' \ r));