BUGFIX: use lil matrices for J.
Use sparse matrices only for big (n>400) circuits by default