Making Druinsky-Toledo matrix nonsingular (it needs to be further modified to be well-conditioned) and adding the ability to set/query the algorithmic blocksize from python