fix error in calculating moments: A = G^-1*C, not C^-1*G, which is the "drift" component we use in simulation