Update for block coordinate descent of full BLIN model D is matrix of lagged Y (covariates) Y is observations
update_MLE_additive(A, B, D, DDT, DTD, DYT, DTY, type = "biten")