void Main::stopSimulation(void)

{

  timer->stop();

  cerr << "Stopping simulation" << endl;

  statusBar()->message("Simulation paused");

  running = false;

}

void Main::togglePaused()

{

  bool s = run->isItemChecked(paused_id);

  if (s) {

    cerr << "Calling start simulation" << endl;

    startSimulation();

  } else {

    cerr << "Calling stop simulation" << endl;

    stopSimulation();

  }

}

void Main::setFluxArrowSize(int size)

{

  flux_arrow_size = size/100.;

}

void Main::enlarge()

{

  canvas.setSceneRect( QRectF( 0,0, canvas.width()*4./3., canvas.height()*4./3.) );

}

void Main::shrink()

{

  canvas.setSceneRect( QRectF( 0,0, canvas.width()*3/4, canvas.height()*3/4) );

}

void Main::scale(double factor)

{

  QMatrix m = editor->matrix();

  m.scale(factor, factor);

  editor->setMatrix( m );

}

void Main::zoomIn()

{

  QMatrix m = editor->matrix();

  m.scale( 1.1, 1.1 );

  editor->setMatrix( m );

}

void Main::zoomOut()

{

  QMatrix m = editor->matrix();

  m.scale( 0.9, 0.9 );

  editor->setMatrix( m );

}

void Main::print()

{

  if ( !printer ) printer = new QPrinter;

  if ( printer->setup(this) ) {

    //    extern Mesh mesh;

    Vector bbll,bbur;

    mesh.BoundingBox(bbll,bbur);

#ifdef QDEBUG

    qDebug() << "bbll = " << bbll << endl;

    qDebug() << "bbur = " << bbur << endl;

#endif

    double cw = (bbur.x - bbll.x);

    double ch = (bbur.y - bbll.y);

    QPainter pp(printer);

    QRect vp=pp.viewport();

#ifdef QDEBUG

    qDebug() << "Paper width = " << vp.width() << " x " << vp.height() << endl;

#endif

    // Note that Cell is also translated...

    pp.translate(-bbur.x,-bbur.y);

    if (cw>ch) {

      pp.scale(vp.width()/(2*cw*Cell::Magnification()), vp.width()/(2*cw*Cell::Magnification()));

    } else {

      pp.scale(vp.height()/(2*ch*Cell::Magnification()), vp.height()/(2*ch*Cell::Magnification()));

    }

    canvas.render(&pp, QRectF(), QRectF(0.,0.,canvas.width(),canvas.height()));

  }

}

void Main::TimeStepWrap(void)

{

  static int t;

  stringstream fname;

  TimeStep();

  t = (int)mesh.getTime();

  if ((par.export_interval > 0) && !(t%par.export_interval)){

    fname << par.datadir << "/" << par.export_fn_prefix;

    fname.fill('0');

    fname.width(6);

    fname << t << ".csv";

    this->exportCellData(QString(fname.str().c_str()));

  }

  // check number of timesteps

  if (t == par.nit) {

    emit SimulationDone();

  }

}

void Main::RestartSim(void)

{

  stopSimulation();

  if ( QMessageBox::question(

			     this,

			     tr("Restart simulation?"),

			     tr("Restart simulation.\n"

				"Are you sure?"),

			     QMessageBox::Yes | QMessageBox::No, QMessageBox::NoButton ) == QMessageBox::Yes ) {

    cerr << "Restarting simulation" << endl;

    //    extern Mesh mesh;

    mesh.Clear();

    Init();

    Plot();

    editor->FullRedraw();

  } 

  //startSimulation();

}

void Main::FitCanvasToWindow(void)

{

  double scale_factor_x = (double)editor->width()/(double)canvas.width();

  double scale_factor_y = (double)editor->height()/(double)canvas.height();

  double scale_factor = scale_factor_x > scale_factor_y ? scale_factor_x : scale_factor_y;

  QMatrix m = editor->matrix();

#ifdef QDEBUG  

  qDebug() << "editor->width() = " << editor->width() << endl;

  qDebug() << "editor->height() = " << editor->height() << endl;

  qDebug() << "scale_factor = " << scale_factor << endl;

  qDebug() << "scale_factor_x = " << scale_factor_x << endl;

  qDebug() << "scale_factor_y = " << scale_factor_y << endl;

#endif

  m.scale( scale_factor, scale_factor );

  editor->setMatrix( m );

  editor->show();

}

void Main::PauseIfRunning(void)

{

  if (running) {

    timer->stop();

  }

}

void Main::ContIfRunning(void)

{

  if (running) {

    timer->start( 0 );

  }

}

void Main::FitLeafToCanvas(void) 

{

  Vector ll,ur;

  mesh.BoundingBox(ll, ur);

  ll*=Cell::Magnification(); ur*=Cell::Magnification();

  // give the leaf some space

  Vector border = ((ur-ll)/5.);

  QRectF bb( ll.x - border.x, ll.y - border.y, ur.x-ll.x + 2*border.x, ur.y-ll.y + 2*border.y );

  // cerr << ur << ", " << ll << endl;

  // editor->fitInView(bb, Qt::KeepAspectRatio);

  editor->ensureVisible(bb);

  //editor->setTransform(viewport);

}

void Main::CleanMesh(void) 

{

  vector<double> clean_chem(Cell::NChem());

  vector<double> clean_transporters(Cell::NChem());

  for (int i=0;i<Cell::NChem();i++) {

    clean_chem[i]=par.initval[i];

    clean_transporters[i]=0.;		

  }

  mesh.CleanChemicals(clean_chem);

  mesh.CleanTransporters(clean_transporters);

  mesh.setTime(0);

  Plot();

  editor->FullRedraw();

}

void Main::CleanMeshChemicals(void) 

{

  vector<double> clean_chem(Cell::NChem());

  for (int i=0;i<Cell::NChem();i++) {

    clean_chem[i]=par.initval[i];

  }

  mesh.CleanChemicals(clean_chem);

  mesh.setTime(0);

  Plot();

  editor->FullRedraw();

}

void Main::CleanMeshTransporters(void) 

{

  vector<double> clean_transporters(Cell::NChem());

  for (int i=0;i<Cell::NChem();i++) {

    clean_transporters[i]=0.;

  }

  mesh.CleanTransporters(clean_transporters);

  mesh.setTime(0);

  Plot();

  editor->FullRedraw();

}

void Main::RandomizeMesh(void) 

{

  vector<double> max_chem(Cell::NChem());

  vector<double> max_transporters(Cell::NChem());

  for (int i=0;i<Cell::NChem();i++) {

    max_transporters[i]=0.;

    max_chem[i]=par.initval[i];

  }

  // Amount of PIN1 at walls

  max_transporters[1] = 0.01;

  mesh.RandomizeChemicals(max_chem, max_transporters);

  Plot();

}

void Main::XMLReadSettings(xmlNode *settings) 

{

  MainBase::XMLReadSettings(settings);

  view->setItemChecked(com_id, showcentersp);

  view->setItemChecked(mesh_id, showmeshp);

  view->setItemChecked(border_id, showbordercellp);

  view->setItemChecked(node_number_id, shownodenumbersp);

  view->setItemChecked(cell_number_id, showcellnumbersp);

  view->setItemChecked(cell_axes_id, showcellsaxesp);

  view->setItemChecked(cell_strain_id, showcellstrainp);

  view->setItemChecked(movie_frames_id, movieframesp);

  view->setItemChecked(only_boundary_id, showboundaryonlyp);

  view->setItemChecked(fluxes_id, showfluxesp);

  view->setItemChecked(hide_cells_id, hidecellsp);

  options->setItemChecked(dyn_cells_id, dynamicscellsp);

  view->setItemChecked( cell_walls_id, showwallsp);

 // view->setItemChecked( apoplasts_id, showapoplastsp);

  editor->setTransform(viewport);

}

xmlNode *Main::XMLSettingsTree(void) 

{

  showcentersp = view->isItemChecked(com_id);

  showmeshp = view->isItemChecked(mesh_id);

  showbordercellp =  view->isItemChecked(border_id);

  shownodenumbersp =  view->isItemChecked(node_number_id);

  showcellnumbersp =  view->isItemChecked(cell_number_id);

  showcellsaxesp = view->isItemChecked( cell_axes_id );

  showcellstrainp = view->isItemChecked( cell_strain_id );

  movieframesp = view->isItemChecked(movie_frames_id);;

  showboundaryonlyp =  view->isItemChecked(only_boundary_id);

  showfluxesp = view->isItemChecked(fluxes_id);

  dynamicscellsp = options->isItemChecked(dyn_cells_id);

  showwallsp = view->isItemChecked( cell_walls_id);

  //showapoplastsp = view->isItemChecked( apoplasts_id);

  hidecellsp = view->isItemChecked( hide_cells_id);

  xmlNode *settings = MainBase::XMLSettingsTree();

  QTransform viewport(editor->transform());

  xmlAddChild(settings, XMLViewportTree(viewport));

  return settings;

}

void Main::exportCellData(QString fileName) {

#ifdef QDEBUG  

  qDebug() << "exportCellData fileName: " << fileName << endl;

#endif

  QFile file(fileName);

  if ( file.open( IO_WriteOnly ) ) {

    QTextStream stream( &file );

    mesh.CSVExportCellData(stream);

    mesh.CSVExportWallData(stream);

    mesh.CSVExportMeshData(stream);

    file.close();

  }

}

void Main::exportCellData() {

  QString fileName;

  Q3FileDialog *fd = new Q3FileDialog( this, "file dialog", TRUE );

  fd->setDir(par.datadir); 

  stopSimulation();

  fd->setMode( Q3FileDialog::AnyFile );

  if ( fd->exec() == QDialog::Accepted ) {

    fileName = fd->selectedFile();

    // extract extension from filename

    QFileInfo fi(fileName);

    QString extension = fi.suffix();

    if (extension.isEmpty()) {

      extension = "csv";

      fileName += ".";

      fileName += extension;

    }

    if (extension!="csv" && extension!="CSV") {

      if (

	  QMessageBox::question(

				 this,

				 tr("Change extension? -- Cell data export"),

				 tr("VirtualLeaf can only export data to CSV format (Excel-compatible)."

				    "Do you want to change the file extension to \".csv\"?"),

				 QMessageBox::Yes, QMessageBox::No ) == QMessageBox::No 

	  ) {

	return exportCellData();

      } else {

	fileName.replace(extension,"csv");

	extension="csv";

      }

    }

    if ( QFile::exists( fileName ) &&

	 QMessageBox::question(

			       this,

			       tr("Overwrite file? -- Cell data export"),

			       tr("A file called %1 already exists."

				  " Do you want to overwrite it?").arg( fileName ),

			       QMessageBox::Yes, QMessageBox::No 

			       ) == QMessageBox::No

	 ) {

      return exportCellData();

    } else {

      exportCellData(fileName);

      QString status_message = QString("Wrote data file to %1").arg(fileName);

      statusBar()->showMessage(status_message);

    }

  }

}