MAINSRC = $${TARGET}.cpp

QT -= network sql xml

QT += qt3support

!win32 {

 GRAPHICS = qt #qwt

}

win32 {

 CONFIG -= console

 LIBXML2DIR = $${LIBDIR}\libxml2

 LIBICONVDIR = $${LIBDIR}\libiconv

 LIBZDIR = $${LIBDIR}\libz

 GRAPHICS = qt 

 RC_FILE = VirtualLeaf.rc

 QMAKE_CXXFLAGS += -DLIBXML_STATIC

 QMAKE_CXXFLAGS += -I$${LIBXML2DIR}\include -I$${LIBICONVDIR}\include -I$${LIBZDIR}\include

 LIBS += -L$${LIBXML2DIR}\lib -lxml2 -L$${LIBICONVDIR}\lib -L$${LIBZDIR}\lib  -lz -lm -lwsock32 -liconv

}

macx:release {

 LIBS+= -dead_strip

}

unix {

 CC = /usr/bin/gcc 

 QWTDIR = /ufs/guravage/opt/qwt-5.2.1-svn

 QMAKE_LIBDIR += $$QWTDIR/lib 

 QMAKE_CXXFLAGS += -fPIC -I/usr/include/libxml2

 QMAKE_LFLAGS += -fPIC

 LIBS += -lxml2 -lz -lm 

}

# Input

HEADERS += \

# apoplastitem.h \

 canvas.h \

 cellbase.h \

 cell.h \

 cellitem.h \

 forwardeuler.h \

       hull.h \ 

 infobar.h \

 mainbase.h \

 mainbase.h \

 matrix.h \

 mesh.h \

 modelcatalogue.h \

 Neighbor.h \

 node.h \

 nodeitem.h \

 nodeset.h \

 OptionFileDialog.h \

 output.h \

 parameter.h \

 pardialog.h \

 parse.h \

 pi.h \

 qcanvasarrow.h \

 random.h \

 rungekutta.h \

 simitembase.h \

 simplugin.h \

 sqr.h \

 tiny.h \

 transporterdialog.h \

 UniqueMessage.h \

 vector.h \

 wallbase.h \

 wall.h \

 wallitem.h \

 warning.h \

 xmlwrite.h \

 $${PARTMPL}

SOURCES += \

# apoplastitem.cpp \

 canvas.cpp \

 cellbase.cpp \

 cell.cpp \

 cellitem.cpp \

 forwardeuler.cpp \

 hull.cpp \

 mainbase.cpp \

 matrix.cpp \

 mesh.cpp \

 modelcatalogue.cpp \

 Neighbor.cpp \

 node.cpp \

 nodeitem.cpp \

 nodeset.cpp \

 OptionFileDialog.cpp \

 output.cpp \

 parameter.cpp \

 pardialog.cpp \

 parse.cpp \

 random.cpp \

 rungekutta.cpp \

 simitembase.cpp \

 transporterdialog.cpp \

 UniqueMessage.cpp \

 vector.cpp \

 wallbase.cpp \

 wall.cpp \

 wallitem.cpp \

 warning.cpp \

 xmlwrite.cpp \

 $$MAINSRC

contains( TARGET, leaf_fleming ) {

 DEFINES += FLEMING	

}

contains(GRAPHICS, qwt) {

 #macx:LIBS += -L$$QWTDIR/lib -lqwt

 #win32:LIBS += -L$$QWTDIR/lib -lqwt5

 #LIBS += -L$$QWTDIR/lib -lqwt

 INCLUDEPATH += $$QWTDIR/include

 DEFINES += HAVE_QWT

 HEADERS += data_plot.h

 SOURCES += data_plot.cpp

}

contains( GRAPHICS, qt ) {

 message( "Building Qt executable" )

 QMAKE_CXXFLAGS += -DQTGRAPHICS # -fpermissive

}

contains( GRAPHICS, xfig ) {

 message("Building Xfig executable (background runnable).")

 QMAKE_CXXFLAGS += -DXFIGGRAPHICS

}

contains( GRAPHICS, x11 ) {

/*

 *

 *  This file is part of the Virtual Leaf.

 *

 *  VirtualLeaf is free software: you can redistribute it and/or modify

 *  it under the terms of the GNU General Public License as published by

 *  the Free Software Foundation, either version 3 of the License, or

 *  (at your option) any later version.

 *

 *  VirtualLeaf is distributed in the hope that it will be useful,

 *  but WITHOUT ANY WARRANTY; without even the implied warranty of

 *  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the

 *  GNU General Public License for more details.

 *

 *  You should have received a copy of the GNU General Public License

 *  along with the Virtual Leaf.  If not, see <http://www.gnu.org/licenses/>.

 *

 *  Copyright 2010 Roeland Merks.

 *

 */

#include <time.h>

#include <string>

#include <fstream>

#include <streambuf>

#include <QGraphicsScene>

#include <QGraphicsView>

#include <qdatetime.h>

#include <q3mainwindow.h>

#include <qstatusbar.h>

#include <qmessagebox.h>

#include <qmenubar.h>

#include <qapplication.h>

#include <qpainter.h>

#include <qprinter.h>

#include <qlabel.h>

#include <qimage.h>

#include <q3progressdialog.h>

#include <qtimer.h>

#include <qslider.h>

#include <qpixmap.h>

#include <qfile.h>

#include <qdir.h>

#include <q3filedialog.h>

#include <QGraphicsItem>

#include <QList>

#include <QDir>

#include <QDebug>

#include <set>

//Added by qt3to4:

#include <Q3ValueList>

#include <Q3PopupMenu>

#include <QMouseEvent>

#include <typeinfo>

#include <cstring>

#include <q3process.h>

#include <qlayout.h>

#include <qspinbox.h>

#include <fstream>

#include <sstream>

#include "pardialog.h"

#include "parameter.h"

#include "canvas.h"

#include "node.h"

#include "nodeset.h"

#include "nodeitem.h"

#include "cellitem.h"

#include "wallitem.h"

#include "mesh.h"

#include "xmlwrite.h"

#include "OptionFileDialog.h"

#include <cstdlib>

#include <cstdio>

#include "modelcatalogue.h"

#include <algorithm>

// Include VIB and PSB logos

#include "psb.xpm"

#include "cwi.xpm"

static const std::string _module_id("$Id$");

using namespace std;

// We use a global variable to save memory - all the brushes and pens in

// the mesh are shared.

#define FNAMESIZE 100

#define QUOTE_ME(s) QUOTE_ME_2NDLEV(s)

#define QUOTE_ME_2NDLEV(s) #s

static QColor dark_red("darkRed");

static const int imageRTTI = 984376;

extern Parameter par;

const QString Main::caption("Virtual leaf");

const QString Main::caption_with_file("Virtual leaf: %1");

FigureEditor::FigureEditor(

			   QGraphicsScene& c, Mesh &m, QWidget* parent,

			   const char* name, Qt::WindowFlags f) :

  QGraphicsView(&c,parent), mesh(m)

{

  intersection_line = 0;

  //angle_line = 0;

  setInteractive(true);

  moving = 0;  

  rotation_mode = false;

}

void FigureEditor::clear()

{

  QList<QGraphicsItem *> list = scene()->items();

  QList<QGraphicsItem *>::Iterator it = list.begin();

  for (; it != list.end(); ++it) {

  // emit this signal

  QObject::connect(   this, SIGNAL( ParsChanged() ), pardial, SLOT( Reset() ) );

}

void Main::savePars()

{

  stopSimulation();

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

  fd->setMode( Q3FileDialog::AnyFile );

  fd->setFilter( "Parameter files (*.par)");

  QString fileName;

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

    fileName = fd->selectedFile();

    ofstream parfile((const char *)fileName);

    par.Write(parfile);

  }

  startSimulation();

}

void Main::readPars()

{

  stopSimulation();

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

  fd->setMode( Q3FileDialog::ExistingFile );

  fd->setFilter( "Parameter files (*.par)");

  QString fileName;

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

    fileName = fd->selectedFile();

    par.Read((const char *)fileName);

  }

  emit ParsChanged();

}

void Main::saveStateXML()

{

  stopSimulation();

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

  fd->setMode( Q3FileDialog::AnyFile );

  QString fileName;

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

    fileName = fd->selectedFile();

    if ( QFile::exists( fileName ) &&

	 QMessageBox::question(

			       this,

			       tr("Overwrite File? -- Leaf Growth"),

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

				  " Do you want to overwrite it?")

			       .arg( fileName ),

			       tr("&Yes"), tr("&No"),

			       QString::null, 1, 1 ) ) {

      return saveStateXML();

    } else {

      mesh.XMLSave((const char *)fileName, XMLSettingsTree());

    }

  }

}

void Main::snapshot()

{

  stopSimulation();

  Q3FileDialog *fd = new Q3FileDialog( this, "Save snapshot", TRUE );

  fd->setMode( Q3FileDialog::AnyFile );

  QString fileName;

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

    fileName = fd->selectedFile();

    if ( QFile::exists( fileName ) &&

	 QMessageBox::question(

			       this,

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

      return snapshot();

    } else {

      // Save bitmaps at 1024x768

      Save((const char *)fileName, extension, 1024, 768);

    }

  }

}

void Main::readPrevStateXML()

{

  // if we have already read a file, read the next file

  if (!currentFile.isEmpty() && working_dir) {

    QString next_file;

    QStringList xml_files = working_dir->entryList("*.xml");

    QString currentFile_nopath = currentFile.section( '/', -1 );

    QString currentFile_path = currentFile.section( '/', 0, -2 );

    QList<QString>::iterator f = xml_files.find( currentFile_nopath );

    if (f == xml_files.end()) {

      return;

    }

    if (f==xml_files.begin()) {

      QMessageBox mb( "Read previous leaf",

		      "No more files",

		      QMessageBox::Information,

		      QMessageBox::Ok | QMessageBox::Default,

		      QMessageBox::NoButton,

		      QMessageBox::NoButton);

      mb.exec();

      return;

    }

    next_file = *(--f);

    next_file = currentFile_path+"/"+next_file;

    readStateXML((const char *)next_file);

  }

}

int Main::readStateXML(const char *filename, bool geometry, bool pars, bool simtime)

{

  try {

    xmlNode *settings;

    mesh.XMLRead((const char *)filename, &settings, geometry, pars, simtime);

#ifdef QDEBUG

  if (!currentFile.isEmpty() && working_dir) {

    QString next_file;

    QStringList xml_files = working_dir->entryList("*.xml");

    QString currentFile_nopath = currentFile.section( '/', -1 );

    QString currentFile_path = currentFile.section( '/', 0, -2 );

    next_file = xml_files.back();

    next_file = currentFile_path+"/"+next_file;

    readStateXML((const char*)next_file);

  }

}

void Main::readFirstStateXML()

{

  // if we have already read a file, read the next file

  if (!currentFile.isEmpty() && working_dir) {

    QString next_file;

    QStringList xml_files = working_dir->entryList("*.xml");

    QString currentFile_nopath = currentFile.section( '/', -1 );

    QString currentFile_path = currentFile.section( '/', 0, -2 );

    next_file = xml_files.front();

    next_file = currentFile_path+"/"+next_file;

    readStateXML((const char*)next_file);

  }

}

void Main::readStateXML()

{

  //  extern Mesh mesh;

  stopSimulation();

#ifdef QDEBUG

  qDebug() << "Trying to open an OptionFileDialog" << endl;

#endif

  OptionFileDialog *fd = new OptionFileDialog( this, "read dialog", TRUE );

  fd->setMode( OptionFileDialog::ExistingFile );

  if (working_dir) {

    fd->setDir(*working_dir);

  }

  QString fileName;

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

    fileName = fd->selectedFile();

    if (working_dir) {

      delete working_dir;

    }

    working_dir = fd->dir();

    if (readStateXML((const char *)fileName,fd->readGeometryP(), fd->readParametersP()) )

      return readStateXML(); // user can try again

  }

}

void Main::clear()

{

  editor->clear();

}

void Main::about()

{

  static QMessageBox* about = new QMessageBox

    ( "VirtualLeaf V1.0",

      "<h3>VirtualLeaf V1.0</h3>"

      "<p>"

      "An Open Source framework for cell-based modeling of plant tissue growth and development <br>"

      "(c) 2005-2008, Roeland Merks <i>et al.</i><br>"

      "    <a href=\"http://www.psb.vib-ugent.be\">VIB Department Plant Systems Biology</a>, "

      "Ghent, Belgium <br>"

      "(c) 2008-2010, <a href=\"http://www.cwi.nl/~merks\">Roeland Merks <i>et al.</i></a> <br>"

      "    <a href=\"http://www.cwi.nl\">Centrum Wiskunde & Informatica</a> and <a href=\"http://www.ncsb.nl\">Netherlands Consortium for Systems Biology</a>, Amsterdam, Netherlands <br>"

"<br>VirtualLeaf is free software: you can redistribute it and/or modify it under the terms of the GNU General Public License as published by the Free Software Foundation, either version 3 of the License, or (at your option) any later version.<br>"

      "<br>If you use this code for your projects, please cite our paper:"

"<br>Merks, Guravage, Inze, and Beemster. An Open Source framework for cell-based modeling of plant tissue growth and development. <i>Plant Physiology,</i> submitted.<br> . <br>"

      "<br>Please share your model plugins and extensions at <a href=\"http://virtualleaf.googlecode.com\">http://virtualleaf.googlecode.com</a>.", 

      QMessageBox::Information, 1, 0, 0, this, 0, FALSE );

  about->setButtonText( 1, "Dismiss" );

  about->show();

}

void Main::gpl()

{

  static QMessageBox* gpl = new QMessageBox ( "GPL License", "", 

  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;

}

#define QDEBUG

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.CSVExportMeshData(stream);

    file.close();

  }

}

void Main::exportCellData() {

  QString fileName;

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

  stopSimulation();

  fd->setMode( Q3FileDialog::AnyFile );

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

    fileName = fd->selectedFile();

    if ( QFile::exists( fileName ) &&

			       this,

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

      return exportCellData();

    } else {

      exportCellData(fileName);

    }

  }

}