double *work;
int  workSize;