anchor.diff:    call elem_create_grid( grid, comm )
     call elem_grid_height( grid, r )
     call elem_grid_width( grid, c )
  -  call elem_grid_size( grid, p )
     call elem_grid_row( grid, row )
     call elem_grid_col( grid, col )
  -  call elem_grid_rank( grid, rank )
   
     ! Create buffers for passing into data for distributed matrices 
     call elem_length( n, row, r, mLoc )
  
positive.code:   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
  
negative.code: 
    ! Useful constants
    integer :: iOne = 1, iZero = 0
    integer :: n = 1000              ! problem size
    integer :: nb = 96               ! algorithmic blocksize
    integer :: comm = MPI_COMM_WORLD ! global communicator
  
  