added option "blockdiag" in 'CreateL' to supress returning blockdiagonal matrix containin Lsubk and added imported function of Matrix in NAMESPACE