Making small changes to objfun and glmm
Fixed objfunc so that value is correct but...

Hessian and gradient are incorrect for the pieces relating the var comps.
Correcting meow in objfun to be calc correctly
Cleaning up and improving code.
Continuing to clean up code, remove unneeded parts
More minor changes to glmm rpack
Implemented Matrix package's Diagonal function

Many matrices are diagonal (D.star,D.star.inv,cdouble), so to save memory I've implemented a sparse matrix package.
More modifications to objfun and glmm