anchor.diff:      enddo
       Parent(Neqns+1) = 0
   
  -    if ( IDBg1/=0 ) then
  +    if ( IDBg_sym/=0 ) then
         write (6,"(' parent')")
         write (6,"(2I6)") (i,Parent(i),i=1,Neqns)
       endif
  
positive.code:     !*Allocations
      allocate (val(NDEg*NDEg),stat=ierr)
      if ( ierr/=0 ) stop "Allocation error:val"
      !write (6,*) "nuform:stage = ", STAge
      kk = 0
      ndof2 = ndof*ndof
  
  
negative.code:   private
    public :: HECMW_SOLVE_DIRECT
  
    real(kind=kreal), parameter :: RMIn = 4.941D-300
  
    integer(kind=kint), parameter :: IDBg_ini = 0
    integer(kind=kint), parameter :: IDBg_sym = 0
    integer(kind=kint), parameter :: IDBg_num = 0
  
  contains
    !----------------------------------------------------------------------
  