A copy of the source code, together with a Makefile for the programme qsearch is available here.
When run from a command line, the programme presents a menu of options for searching for quaternionic solutions to the Yang-Baxter equations:
These solutions may be used by the braid programme to calculate quaternionic invariants of virtiual knots and links, as described in [1].
[1] A. Bartholomew and R. Fenn. Quaternionic Invariants of Virtual Knots and Links (J Knot Theory Ramifications Vol. 17 No. 2 (2008) 231-251).