void update(int value){ reduce(&P_odd,even_to_odd); P_odd.normalize(); reduce(&P_even,odd_to_even); P_even.normalize(); glutPostRedisplay(); glutTimerFunc(2, tick, 0); };
void time_step(){ reduce(&P_odd,even_to_odd); P_odd.normalize(); drift(&P_odd, &N_tensor, 1000); reduce(&P_even,odd_to_even); P_even.normalize(); drift(&P_even, &N_tensor, 1000); };