/* invert33.c - program to invert 3x3 matrices.
 * (C) Copyright 1999, by John Halleck
 * All rights reserved.
 */
/* Version of August 16th, 1999 */

#include <stdio.h>
 /* We have to print status messages  */

#include "errors.h"
 /* We have to print final status messages */

#include "matdebug.h"
 /* The matrix debugging routines would be handy */

#include "qinv.h"
 /* And we need, of course, the routines we were testing. */


static double final[3][3], given[3][3], working[3][3];

int main() {
  int status;
  printf ("non-singular 3 by 3 matrix to invert:\n");
  scanf ("%lf%lf%lf%lf%lf%lf%lf%lf%lf",
          &given[0][0],&given[0][1],&given[0][2],
          &given[1][0],&given[1][1],&given[1][2],
          &given[2][0],&given[2][1],&given[2][2]);
  matprint ("Inverting:", 3, 3, given);
  if ((status = invns (3, final, given, working))) {
     printf ("call failed...\n");
     printerror (status);
     return status;
  } 
  matprint ("Inverse:", 3, 3, final);

  return 0;
}