condition_number()

void condition_number	(	DOUBLE *	norm,
		DOUBLE *	cond
	)

                                                 {
  int descY[9];
  COMPLEX* work;
  int lwork;
  DOUBLE* rwork;
  int lrwork;
  int info;
  int gRows,gCols;
  int ione=1;
  int izero=0;

  int locRows, locCols;
  DOUBLE _norm,_cond;

  gRows  = modgeominfo->ncntr + modgeominfo->nchief;
  gCols  = modgeominfo->ncntr;

  work  = calloc(2,sizeof(COMPLEX));
  rwork = calloc(2,sizeof(DOUBLE));

  lwork = -1;
  lrwork = -1;

  locRows  = numroc_(&gRows ,&runinfo->row_block_size,&runinfo->myrow,&izero,&runinfo->nprows);
  locCols  = numroc_(&gCols ,&runinfo->col_block_size,&runinfo->mycol,&izero,&runinfo->npcols);

  descinit_(descY   ,&gRows,&gCols,&runinfo->row_block_size,&runinfo->col_block_size,&izero,&izero,&runinfo->ctxt,&locRows,&info); 

  pzgecon_("I", &gRows, solution->matY, &ione ,&ione, descY,&_norm,&_cond, work,&lwork, rwork, &lrwork,&info);
  lwork  = (int) work[0];  free(work);  work  = calloc(lwork ,sizeof(COMPLEX));
  lrwork = (int) rwork[0]; free(rwork); rwork = calloc(lrwork,sizeof(DOUBLE));

  pzgecon_("I", &gRows, solution->matY, &ione ,&ione, descY,&_norm,&_cond, work,&lwork, rwork, &lrwork,&info);
      printf(" _NORM(Y)  = %10.8E\n",_norm);
      printf(" _RCOND(Y) = %10.8E   ; _COND(Y) = %10.8E\n",_cond,1.0/_cond);

}