diff -N -r -c HLib-1.3orig/Library/rkmatrix.c HLib-1.3/Library/rkmatrix.c
*** HLib-1.3orig/Library/rkmatrix.c	2004-12-12 17:42:33.000000000 +0100
--- HLib-1.3/Library/rkmatrix.c	2005-10-31 16:28:11.098320384 +0100
***************
*** 1814,1828 ****
    double x1,x2,y;
    double *sigma;
    double ytol;
!   int i, size,k;
    int rows, cols;
  
    assert(r != 0x0);
    rows = r->rows;
    cols = r->cols;
    k = r->kt;
    sigma = NULL;
!   size = rows; if(cols<rows) size = cols;
  
    switch(fill) {
      default:
--- 1814,1829 ----
    double x1,x2,y;
    double *sigma;
    double ytol;
!   int i, size, k, ksigma;
    int rows, cols;
  
    assert(r != 0x0);
    rows = r->rows;
    cols = r->cols;
    k = r->kt;
+   ksigma = imin(k, imin(rows, cols));
    sigma = NULL;
!   size = imin(rows, cols);
  
    switch(fill) {
      default:
***************
*** 1851,1859 ****
        if(k > 0) {
  	sigma = getsigma_rkmatrix(r);
  
! 	for(i=0; i<k; i++) {
! 	  x1 =(0.05+0.9*i/k)*cols;
! 	  x2 =(0.05+0.9*(i+1)/k)*cols;
  	  y = 1.0;
  	  if(sigma[i]>0.0) {
  	    y = log(sigma[i])/log(10.0);
--- 1852,1860 ----
        if(k > 0) {
  	sigma = getsigma_rkmatrix(r);
  
! 	for(i=0; i<ksigma; i++) {
! 	  x1 =(0.05+0.9*i/ksigma)*cols;
! 	  x2 =(0.05+0.9*(i+1)/ksigma)*cols;
  	  y = 1.0;
  	  if(sigma[i]>0.0) {
  	    y = log(sigma[i])/log(10.0);
***************
*** 1890,1898 ****
  
      case 3:
        if(k > 0) {
! 	for(i=0; i<k; i++) {
! 	  x1 =(0.05+0.9*i/k)*cols;
! 	  x2 =(0.05+0.9*(i+1)/k)*cols;
  	  y = 1.0;
  	  if(sigma[i]>0.0) {
  	    y = log(sigma[i])/log(10);
--- 1891,1899 ----
  
      case 3:
        if(k > 0) {
! 	for(i=0; i<ksigma; i++) {
! 	  x1 =(0.05+0.9*i/ksigma)*cols;
! 	  x2 =(0.05+0.9*(i+1)/ksigma)*cols;
  	  y = 1.0;
  	  if(sigma[i]>0.0) {
  	    y = log(sigma[i])/log(10);