anchor.diff:    call elem_symmetric_axbx( A, B, w, X )
   
     ! Print the eigenvalues and eigenvectors
  -  call elem_print_dist_mat( X )
     call elem_print_dist_mat_star_star( w )
   
     ! Shut down Elemental and MPI
  
positive.code:   call elem_register_dist_mat( X, mPad, nPad, iZero, iZero, &
      & XLoc, mPadLoc, grid )
  
    ! I do not know of a good way to flush the output from F90, as the flush
    ! command is not standard. Thus, I chose not to write to stdout from F90.
  
    ! Print the input matrices 
    call elem_print_dist_mat( A )
    call elem_print_dist_mat( B )
  
    ! Set the algorithmic blocksize to 'nb'
    call elem_set_blocksize( nb )
  
    ! For tuning purposes choose one of the following
    ! (see http://elemental.googlecode.com/hg/doc/build/html/advanced/tuning.html
    !  for more details)
    !
    !call elem_set_normal_tridiag_approach
    call elem_set_square_tridiag_approach
  
negative.code:   call elem_create_grid( grid, comm )
    call elem_grid_height( grid, r )
    call elem_grid_width( grid, c )
    call elem_grid_row( grid, row )
    call elem_grid_col( grid, col )
  
    ! Create buffers for passing into data for distributed matrices 
    call elem_length( n, row, r, mLoc )
  