Better Implementation of Wolfgang Schramm for Finding GCD?

Looking for ways to make this program faster, probably running Ramanujan's sum on different thread and share the data?