anchor.diff:    endif
   
     vector_length = 2**order
  -!  allocate density and fill in some phantasy values
  +!  allocate density and fill in some fantasy values
     allocate (density(vector_length, nr_variables, NR_POINTS))
     density = 0.0d0
     do ipoint = 1, NR_POINTS
  
positive.code:   endif
  
    vector_length = 2**order
  !  allocate density and fill in some phantasy values
    allocate (density(vector_length, nr_variables, NR_POINTS))
    density = 0.0d0
    do ipoint = 1, NR_POINTS
  
negative.code:   endif
  
    vector_length = 2**order
  !  allocate density and fill in some fantasy values
    allocate (density(vector_length, nr_variables, NR_POINTS))
    density = 0.0d0
    do ipoint = 1, NR_POINTS
  