How to add Gaussian noise to a sequence of bits in C++

I have generated a bit sequence, encoded it and modulated using BPSK in C++.I want to add Gaussian noise to the modulated symbols. Could anyone please help me on how this can be realized in C++?