     sys[sysIndex]->setSolverType(2);
     sys[sysIndex]->setPrecondType(0);
     sys[sysIndex]->fullJacobian = 1;
-
 
     if(argc>1) {
       sys[sysIndex]->setNumPartitions((int)atoi(argv[1]));


  // if you don't want to visualize, then output the data
  int fileIndex = 0;
#pragma omp parallel shared(workingThread, g_lock)
  {
    int tid = omp_get_thread_num();
    while(true)

        {
          cout << "SWITCH SYSTEM " << loc_working_thread << " -> " << abs(1-loc_working_thread) << "... ";
          sys[loc_working_thread]->transferState(sys[abs(1-loc_working_thread)]);
          updateDone = false;
          cout << "SWITCH COMPLETE." << endl;
        }
      }
