anchor.diff:  c
   clater              if(ierr.ne.izero) return
   c
  -clater            else if(lgdefp.eq.ione.and.ibbarp.eq.ione) then
  +clater            else if(lgdefp.eq.ione.and.intord.eq.ithree) then
   clater              if(reform) call matinit_drv(
   clater     &         alnz,ja,nnz,neq,                                         ! sparse
   clater     &         x,d,iwink,wink,numnp,nwink,                              ! global
  
positive.code: c
                if(ierr.ne.izero) return
  c
  clater            else if(lgdefp.eq.ione.and.ibbarp.eq.izero) then
  clater              if(reform) call matinit_drv(
  clater     &         alnz,ja,nnz,neq,                                         ! sparse
  clater     &         x,d,iwink,wink,numnp,nwink,                              ! global
  
negative.code: c     requested, compute forces due to applied displacements and split
  c     nodes, and perform iterative solution.
  c
              if(lgdefp.eq.izero.and.intord.ne.ithree) then
                if(reform) call matinit_drv(
       &         alnz,ja,nnz,neq,                                         ! sparse
       &         x,d,iwink,wink,numnp,nwink,                              ! global
  