public void TestInitialize()
{
//dx Set up a test problem
g = DenseMatrix.Identity(n);
for (int jj = 0; jj <5; jj++)
{
Vector gg = Normal(size : n);
var hh = gg.PointwiseMultiply(gg);
g.AddRowVector(hh, g);
g.AddRowVector(Normal(size : n).PointwiseMultiply( Normal(size : n)), g);
}
}
