function [ r, y, j] = linear_forward_model( K, q, r, x, iter )
j = K;
y = K * x;
end