implement matrix assembly by simply returning the identity matrix.  This should a sort of projected gradient descent.  It is too late at night for me to compute the actual matrix

[[Imported from SVN: r5545]]