okay, this compiles, let's see how broken it is
works!  now check for unnecessary rounding errors, and maybe try alternative starting point
okay, I guess the compiler actually can't be trusted with stuff like mult_matrix_dxn_vect_n
permutation fix
now to restore nosnp
cluster fix