add wdist_dmatrix.c to wrap matrix inversion and multiply
time to check dgels cross-platform differences