anchor.diff:          endif
           w1=max(0.,w1); w2=max(0.,w2); w3=max(0.,w3);
   
  -        if (abs((w1+w2+w3)-1)>1e-4) then
  +        if (abs((w1+w2+w3)-1)>1e-3) then
               write(*,*) "Error, w1+w2+w3 != 1"
               write(*,*) w1,w2,w3
               write(*,*) "Point"
  
positive.code:         real :: internal_precision
          double precision :: slope, offset
  
          internal_precision = 1e-7
          if (present(precision)) internal_precision = precision
  
          online=.False.
  
negative.code:         real :: x0,y0,x1,y1
          double precision :: slope, x_line
  
          internal_precision = 1e-4
          if (present(precision)) internal_precision = precision
          n = size(poly,2)
  
  