diff -N -r -c HLib-1.3orig/Library/surfacebem.c HLib-1.3/Library/surfacebem.c
*** HLib-1.3orig/Library/surfacebem.c	2005-06-15 18:16:32.000000000 +0200
--- HLib-1.3/Library/surfacebem.c	2006-02-24 15:53:10.202063472 +0100
***************
*** 1094,1100 ****
    const double *gdet = bfactory->gr->gdet;
    int qrhs = bfactory->qrhs;
    const double (*xq)[2] = (usecolbasis ? (const double (*)[2]) bfactory->xq_col[qrhs] : (const double (*)[2]) bfactory->xq_row[qrhs]);
!   const double *wq = (usecolbasis ? bfactory->wq_col[qrhs] : bfactory->wq_row[qrhs]);
    int q = (usecolbasis ? bfactory->q_col[qrhs] : bfactory->q_row[qrhs]);
    double (*basiseval)(int i, double s, double t) = (usecolbasis ? bfactory->coleval : bfactory->roweval);
    int triangles = bfactory->gr->triangles;
--- 1094,1100 ----
    const double *gdet = bfactory->gr->gdet;
    int qrhs = bfactory->qrhs;
    const double (*xq)[2] = (usecolbasis ? (const double (*)[2]) bfactory->xq_col[qrhs] : (const double (*)[2]) bfactory->xq_row[qrhs]);
!   const double *wq = (usecolbasis ? bfactory->wqd_col[qrhs] : bfactory->wqd_row[qrhs]);
    int q = (usecolbasis ? bfactory->q_col[qrhs] : bfactory->q_row[qrhs]);
    double (*basiseval)(int i, double s, double t) = (usecolbasis ? bfactory->coleval : bfactory->roweval);
    int triangles = bfactory->gr->triangles;