Added fix for a bug in numpy that does not do a += a.T correctly if a is bigger than 91 x 91. I changed it to a = a + a.T.