Demonstration programs to project N,J.  The programs require
python version 2.*  and the numpy library.  The pfaffian 
function in the program is from
Gonzalez-Ballestro, et al., Comp. Phys. Comm. 182, 2213 (2011).

Eq8.py     input:  input.dat, uv.*.dat
                 output:  olp.out


Computes the Eq. 8 pfaffian and overlaps on a grid in 
gauge angle phiN and cosine of  rotation angle csb.  Input includes
the even-number parity UV matrices and block/column identifier for the
odd quasiparticle operator.

Eq5.py     input:  input.dat, uv.*.dat
                 output:  olp.out

Same as above, using Eq. 5 and the odd-number parity UV matrices.  Here
the odd quasiparticle is implicit in the matrices.

nzj_project.py   input:  olp.out
                 
The values of (N,Z,J) are set within the program.  In this code, 
Z is a dummy variable.

To reproduce the number in Table I, type on the command line

$./run.sh

or 

$python Eq8.py  
$python nzj_project.py

or

$python Eq5.py
$python nzj_project.py

