Magma code to determine Sophie Germain primes, which are excellent candidates for Brands scheme for the Mammon project.

germain.magma

bits := 32; batchSize := 100;
maxQ := -1;
for index in [1..batchSize] do
    q := RandomPrime(bits - 1);
    if q gt maxQ and IsPrime(2*q + 1) then
        maxQ := q;
        print 2*q + 1, q;
    end if;
end for;