QuickMatch | Evaluation of parallel regular expression matching on GPU | GPU library

I need an R package or function that will allow me to match controls to cases for a large dataset, 5 million subjects. I have tried a few packages, my problems are summarized below. I only tried to match on a single covariate and I most likely will need to match on several.

Package MatchIt: The nearest neighbor, optimal, and genetic methods all just run for hours and hours. The "cem" method runs really quickly but I need to know which cases were matched/unmatched so I can do further analysis with the matched subset. Running the match.data() on the cem results only supplies the weights to be used in a regression and not the matched subset. The paired function in cem would work if I wanted one to one matching but I want to retain as many controls as possible.

matchControls() in the e1071 package: runs for a long time and them returns "not able to allocate vector of size 1352 GB"

Match() function from Matching package: Just runs and runs...

quickmatch() from the quickmatch package: It ran quickly but I am not sure I'm using the function correctly or how to extract the matched data from the "qm_matching" object returned. Below is my attempt using quickmatch on fake data.