diff -N -c -r HLib-1.3/Library/fem2d.c HLib-1.3a/Library/fem2d.c
*** HLib-1.3/Library/fem2d.c	2004-12-12 17:42:32.000000000 +0100
--- HLib-1.3a/Library/fem2d.c	2005-01-07 23:27:30.000000000 +0100
***************
*** 65,71 ****
    fac = new_clusterfactory(n, g2->vertices, 2);
    idx = (int *) allocmem((size_t) sizeof(int) * n);
  
!   /* Create an array "perm" describing the subset of active dofs and
       initialize x, smin and smax */
    j = 0;
    for(i=0; i<g2->vertices; i++) {
--- 65,71 ----
    fac = new_clusterfactory(n, g2->vertices, 2);
    idx = (int *) allocmem((size_t) sizeof(int) * n);
  
!   /* Create an array "idx" describing the subset of active dofs and
       initialize x, smin and smax */
    j = 0;
    for(i=0; i<g2->vertices; i++) {
***************
*** 238,247 ****
      /* Compute the element stiffness matrix by multiplying the gradients */
      for(i=0; i<3; i++)
        for(j=0; j<3; j++)
! 	elem[i][j] = fabs(det) * (g[i][0] * (a[0][0] * g[j][0] +
! 					     a[0][1] * g[j][1]) +
! 				  g[i][1] * (a[1][0] * g[j][0] +
! 					     a[1][1] * g[j][1]) +
  				  (b[0] * g[j][0] +
  				   b[1] * g[j][1]) / 6.0 +
  				  c * (i == j ? 2.0 : 1.0) / 24.0);
--- 238,247 ----
      /* Compute the element stiffness matrix by multiplying the gradients */
      for(i=0; i<3; i++)
        for(j=0; j<3; j++)
! 	elem[i][j] = fabs(det) * ((g[i][0] * (a[0][0] * g[j][0] +
! 					      a[0][1] * g[j][1]) +
! 				   g[i][1] * (a[1][0] * g[j][0] +
! 					      a[1][1] * g[j][1])) / 2.0 +
  				  (b[0] * g[j][0] +
  				   b[1] * g[j][1]) / 6.0 +
  				  c * (i == j ? 2.0 : 1.0) / 24.0);