simplified the model matrix variance to reduce computations
Removed redundant jn property
simplfied the initialization and update of the matrices