/* nag_complex_lin_eqn_mult_rhs(f04adc) Example Program
 *
 * Copyright 1990 Numerical Algorithms Group.
 *
 * Mark 1A revised, (Oct 1990).
 */

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

int main(void)
{
#define NMAX 5
#define TDA NMAX
#define TDB 1
#define TDX 1

  Complex a[NMAX][TDA], b[NMAX][TDB], x[NMAX][TDX];
  Integer i, j, n, nrhs = 1;
  static NagError fail;

  fail.print = TRUE;
  Vprintf("f04adc Example Program Results\n");
  Vscanf("%*[^\n]");  /* Skip heading in data file */
  Vscanf("%ld", &n);
  if (n > 0 && 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);
      for (i = 0; i < n; i++)
        for (j = 0; j < TDX; j++)
          Vscanf(" ( %lf , %lf ) ", &b[i][j].re, &b[i][j].im);
      f04adc(n, nrhs, &a[0][0], (Integer)TDA, &b[0][0], (Integer)TDB,
             &x[0][0], (Integer)TDX, &fail);
      if (fail.code!=NE_NOERROR)
        return EXIT_FAILURE;
      else
        {
          Vprintf("Solution\n");
          for (i=0; i<n; i++)
            Vprintf("(%7.4f, %7.4f)\n", x[i][0].re, x[i][0].im);
        }
    }
  else
    {
      Vfprintf(stderr, "Error: n is out of range: n = %ld\n", n);
      return EXIT_FAILURE;
    }
  return EXIT_SUCCESS;
}