add GPLv3 license - this is the same license used by the original author on github: https://github.com/dqnykamp/sonets

- add arguments for different algorithms
  algorithms can be:
    Standard
    Triangular (i-th row and i-th column)
    Block B## based: ## indicates block size
       and can be used to optimize cache use

fix calc_alpha - conversion to Bernoulli Matrix must change W(i,j) values to 0 and 1 in memory before estimating the alpha values

avoid storing extra (Bernoulli) matrix - but do the Bernoulli conversion only when writing the matrix

initial import, assumed to be closest to the original code; obtained from /fs3/group/jonasgrp/Jose/SONET