main()

int main	(	int	iargc,
		char *	argv[]
	)

Definition at line 113 of file lapacktest.cc.

References DUMP, LAPACKTEST_VERSION, Options::nsize, test(), and Options::verbose.

{

  // define usage information
  char usage_text[]=
  {
    LAPACKTEST_VERSION "\n"
    "usage: lapacktest [-n n]" "\n"
    "   or: lapacktest --help|-h" "\n"
  };

  // define full help text
  char help_text[]=
  {
    "\n"
    "-v     be verbose" "\n"
    "-n N   test system of size NxN" "\n"
  };

  // define commandline options
  using namespace tfxx::cmdline;
  static Declare options[]= 
  {
    // 0: print help
    {"help",arg_no,"-"},
    // 1: verbose mode
    {"v",arg_no,"-"},
    // 2: verbose mode
    {"n",arg_yes,"3"},
    {NULL}
  };

  // no arguments? print usage...
  if (iargc<2) 
  {
    cerr << usage_text << endl;
    exit(0);
  }

  // collect options from commandline
  Commandline cmdline(iargc, argv, options);

  // help requested? print full help text...
  if (cmdline.optset(0))
  {
    cerr << usage_text << endl;
    cerr << help_text << endl;
    exit(0);
  }

  /*
  // dummy operation: print option settings
  for (int iopt=0; iopt<2; iopt++)
  {
    cout << "option: '" << options[iopt].opt_string << "'" << endl;
    if (cmdline.optset(iopt)) {  cout << "  option was set"; }
    else { cout << "option was not set"; }
    cout << endl;
    cout << "  argument (string): '" << cmdline.string_arg(iopt) << "'" << endl;
    cout << "     argument (int): '" << cmdline.int_arg(iopt) << "'" << endl;
    cout << "    argument (long): '" << cmdline.long_arg(iopt) << "'" << endl;
    cout << "   argument (float): '" << cmdline.float_arg(iopt) << "'" << endl;
    cout << "  argument (double): '" << cmdline.double_arg(iopt) << "'" << endl;
    cout << "    argument (bool): '";
    if (cmdline.bool_arg(iopt))
    { cout << "true"; } else { cout << "false"; }
    cout << "'" << endl;
  }
  while (cmdline.extra()) { cout << cmdline.next() << endl; }

  // dummy operation: print rest of command line
  while (cmdline.extra()) { cout << cmdline.next() << endl; }
  */

  Options opt;
  opt.verbose=cmdline.optset(1);
  opt.nsize=cmdline.int_arg(2);

  /*----------------------------------------------------------------------*/

  {

    if (opt.verbose)
    {
      cout << "Define system of size " << opt.nsize << "x" << opt.nsize << endl;
    }

    Tmatrix M(opt.nsize,opt.nsize);
    Tmatrix R(opt.nsize,1);

    // set up random number generator
    tfxx::numeric::RNGgaussian rng;

    {
      // set up system matrix
      aff::Iterator<Tmatrix> I(M);
      while (I.valid())
      {
        *I = rng();
        ++I;
      }
    }

    {
      // set up system vector
      aff::Iterator<Tmatrix> I(R);
      while (I.valid())
      {
        *I = rng();
        ++I;
      }
    }
    test(M,R);
  }

  /*----------------------------------------------------------------------*/

  {
    Tmatrix M(3,3), V(3,1), W(3,2);
    M(1,1)=1.;
    M(1,2)=0.4;
    M(1,3)=0.6;
    M(2,1)=-0.2;
    M(2,2)=-1.;
    M(2,3)=0.7;
    M(3,1)=-1.;
    M(3,2)=-0.3;
    M(3,3)=1.;
    DUMP(M);
    V(1)=-80.;
    V(2)=70.;
    V(3)=-35.;
    DUMP(V);
    test(M,V);
    W(1,1)=-80.;
    W(2,1)=70.;
    W(3,1)=-35.;
    W(1,2)=0.;
    W(2,2)=2.;
    W(3,2)=-3.;
    DUMP(W);
    test(M,W);
  }
}

Here is the call graph for this function: