void abcSaf2::print(funkyPars *p){
if(p->numSites==0||(doSaf2==0))
return;
// fprintf(stderr,"newDim:%d doSaf2:%d\n",newDim,doSaf2);
if(doSaf2==1){
realRes2 *r=(realRes2 *) p->extras[index];
int id=0;
if(outputBanded)
printSparse(p,index,outfileSFS,outfileSFSPOS,header->name[p->refId]);
}
}
/* The gateway function */
void mexFunction( int nlhs, mxArray *plhs[],
int nrhs, const mxArray *prhs[])
{
mwIndex *irs,*jcs;
double *si,*sr;
int nnz, m,n;
nnz = mxGetNzmax(prhs[0]);
sr = mxGetPr(prhs[0]);
si = mxGetPi(prhs[0]);
irs = mxGetIr(prhs[0]);
jcs = mxGetJc(prhs[0]);
m = mxGetM(prhs[0]);
n = mxGetN(prhs[0]);
printSparse(m,n,nnz,sr,si,irs,jcs);
}
