ScaLAPACK PDSYEVX Example
  =========================

  This example calls the PDSYEV simple driver routine (Parallel, Double-precision, SYmmetric,
  Eigenvalue decomposition driVer.)  When compiled and run, it produces output that can be 
  pasted directly into Matlab.



  Load the acml module:

    % module load acml

  Copy the files in the EXAMPLES/pdsyevx directory to a convenient working directory:

    % cp $SCALAPACK/EXAMPLES/pdsyevx/* .

  Clean the test directory and make a new executable:

    % make clean
    % make


  Run
  ===

  Submit the PBS job script:

    % qsub pdsyevx.pbs

  Solution data is validated. Review the contents of the pdsyev.o<pbs_jobid> file.


  Files
  =====

  Input files:

    none

  Output files:

    "stdout"         decomposition data and Matlab statements


  Details
  =======

  The "stdout" file contains decomposition data and a cluster of statements that may
  be pasted directly into Matlab for further analysis.

    Z(     1,     1)=     -0.485377931880384694D-01
    Z(     2,     1)=     -0.122222718667358635D+00
    Z(     3,     1)=     -0.282485135303397361D+00
    Z(     4,     1)=      0.950214627337748530D+00
    Z(     1,     2)=      0.912072227031435201D-01
    Z(     2,     2)=      0.426620092092790393D+00
    Z(     3,     2)=     -0.877038303257524188D+00
    Z(     4,     2)=     -0.201197301593937117D+00
    Z(     1,     3)=      0.495193043045526249D+00
    Z(     2,     3)=     -0.798642041261893088D+00
    Z(     3,     3)=     -0.298844131980119498D+00
    Z(     4,     3)=     -0.166273644422072125D+00
    Z(     1,     4)=      0.862617629821305298D+00
    Z(     2,     4)=      0.406482218545086926D+00
    Z(     3,     4)=      0.248391118466086747D+00
    Z(     4,     4)=      0.170190725350427008D+00
      N =           4
      A =  hilb(N) + diag([1:-1/N:1/N])
      W(           1 )=  0.304814036060551      ;
      W(           2 )=  0.581961202142018      ;
      W(           3 )=  0.908498110005126      ;
      W(           4 )=   2.38091712798278      ;
      backerror = A - Z * diag(W) * Z' 
      resid = A * Z - Z * diag(W)
      ortho = Z' * Z - eye(N)
      norm(backerror)
      norm(resid)
      norm(ortho)