Function main

SRC/matrix/main.cpp:36–239

OPS_Stream *opserrPtr = &sserr;

int main()
{

  ID data(0);

  opserr << data << endln;
  opserr << data.insert(1) << endln;
  opserr << data << endln;
  opserr << data.insert(2) << endln;
  opserr << data << endln;
  opserr << data.insert(4) << endln;
  opserr << data << endln;
  opserr << data.insert(5) << endln;
  opserr << data << endln;
  opserr << data.insert(6) << endln;
  opserr << data << endln;
  opserr << data.insert(3) << endln;
  opserr << data << endln;
  opserr << data.insert(-1) << endln;
  opserr << data << endln;

  opserr << data.insert(3) << endln;
  opserr << data << endln;
  exit(0);

  opserr.setPrecision(14);
  opserr.setFloatField(SCIENTIFIC);

  Matrix a(3,3);
  a(0,0)= 1.35691200000000e+08; a(0,1)= 1.49035621973463e+07; a(0,2) =  -3.72529029846191e-09;
  a(1,0) = 1.49035621973463e+07; a(1,1)= 1.63692388430620e+06; a(1,2)= -4.65661287307739e-10;
  a(2,0) =-3.72529029846191e-09; a(2,1)=  -4.65661287307739e-10; a(2,2)= 3.19209024960000e+07;

  opserr << a;

  Matrix b(3,3);

  a.Invert(b);
  opserr << b;

  opserr << a*b;

  Matrix I(3,3);
  for (int i=0; i<3; i++)
    I(i,i) = 1.0;

  opserr << a.Solve(I,b);
  
  opserr << b;
  opserr << a*b;
  opserr << I;

  /*
    // Vector  TEST
    cout << "\n***** Vector TEST *****  \n";
    Vector a;
    cout << "Vector a;     a -> " << a;
    Vector b(5);