Source code for our paper "Solving Discrete Logarithms in Smooth-Order Groups with CUDA" from SHARCS 2012. The cuda-11.5 branch will compile on CUDA 11.5 and now runs with two A100 GPUs, though is perhaps not tuned optimally.

Steven Engler 3a802c8e75 Improved the logging for dlrho. %!s(int64=8) %!d(string=hai) anos
.gitignore 3a50b8f3a9 Added a '.gitignore' file. %!s(int64=8) %!d(string=hai) anos
COPYING e2973619c9 Commit cudadl-0.8 to git %!s(int64=14) %!d(string=hai) anos
Makefile 0e11651c7d The worker's 'cuda_dev_id' argument selects a GPU. %!s(int64=8) %!d(string=hai) anos
README 9fca8d29fd cudadl-0.9 %!s(int64=14) %!d(string=hai) anos
atomic_iostream.h 3a802c8e75 Improved the logging for dlrho. %!s(int64=8) %!d(string=hai) anos
controller.cc 71bde25be1 Output number of kernel launches for each subproblem. %!s(int64=8) %!d(string=hai) anos
controller.h edfd20a685 Add new cmdline options to controller %!s(int64=14) %!d(string=hai) anos
controller_main.cc cde6e66784 Include 'unistd.h' in source files. %!s(int64=8) %!d(string=hai) anos
cudadl.h 71bde25be1 Output number of kernel launches for each subproblem. %!s(int64=8) %!d(string=hai) anos
dlrho.cc 3a802c8e75 Improved the logging for dlrho. %!s(int64=8) %!d(string=hai) anos
dpnode.cc c081555e83 Fix the MPI driver to not conflict with libevent %!s(int64=14) %!d(string=hai) anos
dpnode.h 5b62cbb2c4 Split off the three main()s in preparation for MPI wrapper %!s(int64=14) %!d(string=hai) anos
dpnode_main.cc 5b62cbb2c4 Split off the three main()s in preparation for MPI wrapper %!s(int64=14) %!d(string=hai) anos
dpstream.cu 3a802c8e75 Improved the logging for dlrho. %!s(int64=8) %!d(string=hai) anos
evutils.cc 0e999113e0 Reconstruct answer and launch new problems %!s(int64=14) %!d(string=hai) anos
evutils.h 4a0de9516a Enable libevent thread safety for the worker %!s(int64=14) %!d(string=hai) anos
gen_N.cc 9fca8d29fd cudadl-0.9 %!s(int64=14) %!d(string=hai) anos
gencios_reg_20 771b13c42c Use macros to include the inline assembly instead of inline functions %!s(int64=14) %!d(string=hai) anos
logparse.pl 9fa203afbe Add a program to parse the logfile for some interesting statistics %!s(int64=14) %!d(string=hai) anos
mpi.cc ed7cb8a165 The first argument to 'mpi' specifies number of workers. %!s(int64=8) %!d(string=hai) anos
parrhoasm.cu 3a802c8e75 Improved the logging for dlrho. %!s(int64=8) %!d(string=hai) anos
run_timing_multiprocess.sh ce3141f7e1 Added scripts for collecting timing data. %!s(int64=8) %!d(string=hai) anos
subproblem.h e5decaf5df Communicate the dpfreq to the workers correctly %!s(int64=14) %!d(string=hai) anos
summarize_data_multiprocess.py ce3141f7e1 Added scripts for collecting timing data. %!s(int64=8) %!d(string=hai) anos
worker.cc 71bde25be1 Output number of kernel launches for each subproblem. %!s(int64=8) %!d(string=hai) anos
worker.h 0e11651c7d The worker's 'cuda_dev_id' argument selects a GPU. %!s(int64=8) %!d(string=hai) anos
worker_main.cc 0e11651c7d The worker's 'cuda_dev_id' argument selects a GPU. %!s(int64=8) %!d(string=hai) anos

README

cudadl-0.9
21 Mar 2012
Ryan Henry and Ian Goldberg
{rhenry,iang}@cs.uwaterloo.ca
http://crysp.uwaterloo.ca/software/

This package contains the source code to our CUDA implementation of
van Oorschot and Wiener's parallel version of the Pollard rho discrete
log algorithm. It is intended for use on 1536-bit moduli that are
RSA numbers with smooth totient; that is, the modulus N=pq, where p and q
are 768-bit primes, and the prime factors of p-1 and q-1 are all
distinct and less then B, for a parameter B. [The value 1536 is
hardcoded as "WORDS = 24" (24*32*2 = 1536) in the Makefile; it is easy
to change this value and recompile if desired.] Note that this means
the totient of N = \phi(N) = (p-1)(q-1) has all prime factors less than
B; that is, \phi(n) is "B-smooth".

Usage:

1. Build the software. You'll need:

NTL
GMP
NVIDIA CUDA Toolkit 3.1
2 M2050 (or other compute capability level 2.0) CUDA cards
[If you have more or just 1, you'll need to modify dlrho.cc,
unfortunately.]

Hopefully just typing "make" should work. It will build gen_N and
dlrho.

2. Create the modulus N as, for example, a 1536-bit RSA number whose
totient is 2^50-smooth:

./gen_N 1536 50 > N

3. Generate a DL problem mod N and solve it:

./dlrho < N

This software is described in "Solving Discrete Logarithms in
Smooth-Order Groups with CUDA", CACR technical report 2012-02,
http://www.cacr.math.uwaterloo.ca/techreports/2012/cacr2012-02.pdf

This program is covered under version 3 of the GNU General Public
Licence; see the file COPYING for more information.

Changelog:

0.9 (21 Mar 2012)
Extend the code to handle smoothness levels (B) larger than 2^60. Now
we can handle up to 2^92. We have successfully run a test with
B = 2^80.

0.8 (23 Jan 2012)
Initial public release