/* nag_hermitian_eigensystem(f02axc) Example Program
 *
 * Copyright 1991 Numerical Algorithms Group.
 *
 * Mark 2, 1991.
 */

#include <nag.h>
#include <stdio.h>
#include <nag_stdlib.h>
#include <nagf02.h>

#define NMAX 4
#define TDA  NMAX
#define TDV NMAX

int main(void)
{
  Integer i, j, n;
  Complex a[NMAX][TDA], v[NMAX][TDV];
  double r[NMAX];

  Vprintf("f02axc Example Program Results\n");
  Vscanf("%*[^\n]"); /* Skip heading in data file */
  Vscanf("%ld", &n);
  if (n>=1 || n<=NMAX)
    {
      for (i=0; i<n; i++)
        for (j=0; j<n; j++)
          Vscanf(" ( %lf , %lf ) ", &a[i][j].re,&a[i][j].im);
      f02axc(n, &a[0][0], (Integer)TDA, r, &v[0][0],
             (Integer)TDV, NAGERR_DEFAULT);
      Vprintf("Eigenvalues\n");
      for (i=0; i<n; i++)
        Vprintf("%9.4f", r[i]);
      Vprintf("\nEigenvectors\n");
      for (i=0; i<n; i++)
        for (j=0; j<n; j++)
          Vprintf("(%7.3f %7.3f )%s",v[i][j].re,v[i][j].im,
                  (j%4==3 || j==n-1) ? "\n" : " ");
      return EXIT_SUCCESS;
    }
  else
    {
      Vfprintf(stderr,"n is out of range: n = %4ld,\n",n);
      return EXIT_FAILURE;
    }
}