We now compute the GRDA quadratic forms in 'predict' in a single matrix multiplication rather than for each test observation. Improves #1 but does not do enough
In GRDA CV function we now update estimates rather than recomputing expensive quantities for each set of tuning parameters. Closes #1