anchor.diff:                          nspawn = sign(nspawn, parent_sign)
                       end if
   
  -                    ! 5. Set connection to the address of the spawned element
  +                    ! Set connection to the address of the spawned element
                       connection = connection_jk
   
                   else
  
positive.code:             nspawn_ki = int(pspawn_ki)
              pspawn_ki = pspawn_ki - nspawn_ki
  
              if (pspawn_ki > psuccess) nspawn_ki = nspawn_ki + 1
  
              if(nspawn_ki > 0 ) then
              ! Successful spawning on ki
  
                  ! 1.2 Generate the second random excitation 
                  call gen_excit_ptr(cdet, Pgen_jk, connection_jk, hmatel_jk)
  
                  ! Calculate P_gen for the second excitation
                  pspawn_jk = Xo_ * abs(hmatel_jk) / Pgen_jk
                  psuccess = genrand_real2()
  
                  ! Multiple offspring again...
                  nspawn_jk = int(pspawn_jk)
                  pspawn_jk = pspawn_jk - nspawn_jk
  
  
                  if (pspawn_jk > psuccess) nspawn_jk = nspawn_jk + 1
  
                  if (nspawn_jk > 0) then
  
negative.code:                     ! Combine the spawning numbers
                      nspawn = nspawn_ki * nspawn_jk
  
                      ! If H_ij is positive, then the spawned walker is of opposite
                      ! sign to the parent, otherwise the spawned walkers if of the same
                      ! sign as the parent.
                      if (hmatel_ki*hmatel_jk > 0.0_p) then
  