2010-10-14    <merks@cwi.nl>

	* cell.cpp (DivideWalls): accomodated for rename of Circumference -> WallCircumference

	* hull.h: added an operator< to sort Points

	* hull.cpp: added an operator< to sort Points

	* cellbase.cpp (ExactCircumference): I added a new function ExactCircumference, yielding the circumference of the cell along its wall_elements

	* VirtualLeaf.cpp: adjust info_string to accomodate for new name of function CellBase::Circumference -> CellBase::WallCircumference

	* mesh.cpp: corrected Mesh::Compactness, the boundary coordinates need to be sorted in x,y order for the convex hull algorithm (thanks Margriet!). I updated CSVExportCellData so it exports the circumferences of hull and boundary_polygon.

2010-10-14    <guravage@caterpie.sen.cwi.nl>

	* VirtualLeafpar.tmpl (export_fn_prefix): changed to 'cell.'

	* canvas.cpp (TimeStepWrap): Removed TimeStamp(). Cell data

	filename now incorporates iteration number.

	* canvas.h (MainBase): Removed TimeStamp().

2010-10-08    <guravage@caterpie.sen.cwi.nl>

title = Parameter values for The Virtual Leaf / title

label = <b>Visualization</b> / label

arrowcolor = white / string

arrowsize = 100 / double

textcolor = red / string

cellnumsize = 1 / int

nodenumsize = 1 / int

node_mag = 1.0 / double 

outlinewidth = 1.0 / double

cell_outline_color = forestgreen / string

resize_stride = 0 / int

label = / label

label = <b>Cell mechanics</b> / label

T = 1.0 / double

lambda_length = 100. / double

lambda_celllength = 0. / double

target_length = 60. / double 

cell_expansion_rate = 1. / double

cell_div_expansion_rate = 0. / double

auxin_dependent_growth = true / bool

ode_accuracy = 1e-4 / double

mc_stepsize = 0.4 / double

mc_cell_stepsize = 0.2 / double

energy_threshold = 1000. / double

bend_lambda = 0. / double

alignment_lambda = 0. / double

rel_cell_div_threshold = 2. / double

rel_perimeter_stiffness = 2 / double

collapse_node_threshold = 0.05 / double

morphogen_div_threshold = 0.2 / double

morphogen_expansion_threshold = 0.01 / double

copy_wall = true / bool

label = / label

label = <b>Auxin transport and PIN1 dynamics</b> / label

source = 0. / double

D = 0., 0.0, 0.0, 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0. / doublelist

initval = 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0. / doublelist

k1 = 1. / double

k2 = 0.3 / double

	       cerr << "Old axis is " << old_axis << ", new axis is " << norm_long_axis << endl; 

	    */

	    alignment_dh += alignment_before - alignment_after;

	    /* cerr << "alignment_dh  = " << alignment_dh << endl;

	       cerr << "cellvec = " << c.cellvec << endl;*/

	  } else {

	    // if we have no length constraint, still need to update area

	    delta_intgrl_list.push_back(DeltaIntgrl(delta_A,0,0,0,0,0));

	  }

	  old_l1=(old_p-i_min_1).Norm();

	  old_l2=(old_p-i_plus_1).Norm();

	  new_l1=(new_p-i_min_1).Norm();

	  new_l2=(new_p-i_plus_1).Norm();

	  static int count=0;

	  // Insertion of nodes (cell wall yielding)

	  if (!node.fixed) {

	      node_insertion_queue.push( Edge(cit->nb1, &node) );

	    }

	      node_insertion_queue.push( Edge(&node, cit->nb2 ) );

	    }

	    count++;

	  }

	  length_dh += 2*Node::target_length * ( w1*(old_l1 - new_l1) + 

						 w2*(old_l2 - new_l2) ) +

	    w1*(DSQR(new_l1) 

		- DSQR(old_l1)) 

	    + w2*(DSQR(new_l2) 

		  - DSQR(old_l2));

	}

	// bending energy also holds for outer boundary

	// first implementation. Can probably be done more efficiently

	// calculate circumcenter radius (gives local curvature)

#include "output.h"

#include "parse.h"

#include "xmlwrite.h"

#include "warning.h"

#include <QLocale>

#include <QDir>

#include <QStringList>

using namespace std;

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

Parameter::Parameter() {

  arrowcolor = strdup("white");

  arrowsize = 100;

  textcolor = strdup("red");

  cellnumsize = 1;

  nodenumsize = 1;

  node_mag = 1.0;

  outlinewidth = 1.0;

  cell_outline_color = strdup("forestgreen");

  resize_stride = 0;

  T = 1.0;

  lambda_length = 100.;

  lambda_celllength = 0.;

  target_length = 60.;

  cell_expansion_rate = 1.;

  cell_div_expansion_rate = 0.;

  auxin_dependent_growth = true;

  ode_accuracy = 1e-4;

  mc_stepsize = 0.4;

  mc_cell_stepsize = 0.2;

  energy_threshold = 1000.;

  bend_lambda = 0.;

  alignment_lambda = 0.;

  rel_cell_div_threshold = 2.;

  rel_perimeter_stiffness = 2;

  collapse_node_threshold = 0.05;

  morphogen_div_threshold = 0.2;

  morphogen_expansion_threshold = 0.01;

  copy_wall = true;

  source = 0.;

  D = new double[15];

  D[0] = 0.;

  D[1] = 0.0;

  D[2] = 0.0;

  D[3] = 0.;

  D[4] = 0.;

  static bool ReadP=false;

  if (ReadP) {

    //throw "Run Time Error in parameter.cpp: Please Read parameter file only once!!";

    CleanUp();

  } else

    ReadP=true;

  FILE *fp=OpenReadFile(filename);

  arrowcolor = sgetpar(fp, "arrowcolor", "white", true);

  arrowsize = fgetpar(fp, "arrowsize", 100, true);

  textcolor = sgetpar(fp, "textcolor", "red", true);

  cellnumsize = igetpar(fp, "cellnumsize", 1, true);

  nodenumsize = igetpar(fp, "nodenumsize", 1, true);

  node_mag = fgetpar(fp, "node_mag", 1.0, true);

  outlinewidth = fgetpar(fp, "outlinewidth", 1.0, true);

  cell_outline_color = sgetpar(fp, "cell_outline_color", "forestgreen", true);

  resize_stride = igetpar(fp, "resize_stride", 0, true);

  T = fgetpar(fp, "T", 1.0, true);

  lambda_length = fgetpar(fp, "lambda_length", 100., true);

  lambda_celllength = fgetpar(fp, "lambda_celllength", 0., true);

  target_length = fgetpar(fp, "target_length", 60., true);

  cell_expansion_rate = fgetpar(fp, "cell_expansion_rate", 1., true);

  cell_div_expansion_rate = fgetpar(fp, "cell_div_expansion_rate", 0., true);

  auxin_dependent_growth = bgetpar(fp, "auxin_dependent_growth", true, true);

  ode_accuracy = fgetpar(fp, "ode_accuracy", 1e-4, true);

  mc_stepsize = fgetpar(fp, "mc_stepsize", 0.4, true);

  mc_cell_stepsize = fgetpar(fp, "mc_cell_stepsize", 0.2, true);

  energy_threshold = fgetpar(fp, "energy_threshold", 1000., true);

  bend_lambda = fgetpar(fp, "bend_lambda", 0., true);

  alignment_lambda = fgetpar(fp, "alignment_lambda", 0., true);

  rel_cell_div_threshold = fgetpar(fp, "rel_cell_div_threshold", 2., true);

  rel_perimeter_stiffness = fgetpar(fp, "rel_perimeter_stiffness", 2, true);

  collapse_node_threshold = fgetpar(fp, "collapse_node_threshold", 0.05, true);

  morphogen_div_threshold = fgetpar(fp, "morphogen_div_threshold", 0.2, true);

  morphogen_expansion_threshold = fgetpar(fp, "morphogen_expansion_threshold", 0.01, true);

  copy_wall = bgetpar(fp, "copy_wall", true, true);

  source = fgetpar(fp, "source", 0., true);

  D = dgetparlist(fp, "D", 15, true);

  initval = dgetparlist(fp, "initval", 15, true);

  k1 = fgetpar(fp, "k1", 1., true);

  k2 = fgetpar(fp, "k2", 0.3, true);

  r = fgetpar(fp, "r", 1., true);

  kr = fgetpar(fp, "kr", 1., true);

    return true_str;

  else

    return false_str;

}

void Parameter::Write(ostream &os) const {

  if (arrowcolor) 

  os << " arrowcolor = " << arrowcolor << endl;

  os << " arrowsize = " << arrowsize << endl;

  if (textcolor) 

  os << " textcolor = " << textcolor << endl;

  os << " cellnumsize = " << cellnumsize << endl;

  os << " nodenumsize = " << nodenumsize << endl;

  os << " node_mag = " << node_mag << endl;

  os << " outlinewidth = " << outlinewidth << endl;

  if (cell_outline_color) 

  os << " cell_outline_color = " << cell_outline_color << endl;

  os << " resize_stride = " << resize_stride << endl;

  os << " T = " << T << endl;

  os << " lambda_length = " << lambda_length << endl;

  os << " lambda_celllength = " << lambda_celllength << endl;

  os << " target_length = " << target_length << endl;

  os << " cell_expansion_rate = " << cell_expansion_rate << endl;

  os << " cell_div_expansion_rate = " << cell_div_expansion_rate << endl;

  os << " auxin_dependent_growth = " << sbool(auxin_dependent_growth) << endl;

  os << " ode_accuracy = " << ode_accuracy << endl;

  os << " mc_stepsize = " << mc_stepsize << endl;

  os << " mc_cell_stepsize = " << mc_cell_stepsize << endl;

  os << " energy_threshold = " << energy_threshold << endl;

  os << " bend_lambda = " << bend_lambda << endl;

  os << " alignment_lambda = " << alignment_lambda << endl;

  os << " rel_cell_div_threshold = " << rel_cell_div_threshold << endl;

  os << " rel_perimeter_stiffness = " << rel_perimeter_stiffness << endl;

  os << " collapse_node_threshold = " << collapse_node_threshold << endl;

  os << " morphogen_div_threshold = " << morphogen_div_threshold << endl;

  os << " morphogen_expansion_threshold = " << morphogen_expansion_threshold << endl;

  os << " copy_wall = " << sbool(copy_wall) << endl;

  os << " source = " << source << endl;

  os << " D = "<< D[0] << ", " << D[1] << ", " << D[2] << ", " << D[3] << ", " << D[4] << ", " << D[5] << ", " << D[6] << ", " << D[7] << ", " << D[8] << ", " << D[9] << ", " << D[10] << ", " << D[11] << ", " << D[12] << ", " << D[13] << ", " << D[14] << endl;

  os << " initval = "<< initval[0] << ", " << initval[1] << ", " << initval[2] << ", " << initval[3] << ", " << initval[4] << ", " << initval[5] << ", " << initval[6] << ", " << initval[7] << ", " << initval[8] << ", " << initval[9] << ", " << initval[10] << ", " << initval[11] << ", " << initval[12] << ", " << initval[13] << ", " << initval[14] << endl;

  os << " k1 = " << k1 << endl;

  os << " k2 = " << k2 << endl;

  os << " r = " << r << endl;

  os << " kr = " << kr << endl;

}

{

  xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);

  xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "resize_stride" );

  ostringstream text;

    text << resize_stride;

xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());

}

{

  xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);

  xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "T" );

  ostringstream text;

    text << T;

xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());

}

{

  xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);

  xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "lambda_length" );

  ostringstream text;

    text << lambda_length;

xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());

}

{

  xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);

  xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "lambda_celllength" );

  ostringstream text;

    text << lambda_celllength;

xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());

}

{

  xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);

  xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "target_length" );

  ostringstream text;

    text << target_length;

xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());

}

{

  xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);

  xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "cell_expansion_rate" );

  ostringstream text;

    text << cell_expansion_rate;

xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());

}

{

  xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);

  xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "cell_div_expansion_rate" );

  ostringstream text;

    text << cell_div_expansion_rate;

if (!strcmp(namec, "node_mag")) {

  node_mag = standardlocale.toDouble(valc, &ok);

  if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'node_mag' from XML file.",valc); }

}

if (!strcmp(namec, "outlinewidth")) {

  outlinewidth = standardlocale.toDouble(valc, &ok);

  if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'outlinewidth' from XML file.",valc); }

}

if (!strcmp(namec, "cell_outline_color")) {

  if (cell_outline_color) { free(cell_outline_color); }

  cell_outline_color=strdup(valc);

}

if (!strcmp(namec, "resize_stride")) {

  resize_stride = standardlocale.toInt(valc, &ok);

  if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to integer while reading parameter 'resize_stride' from XML file.",valc); }

}

if (!strcmp(namec, "T")) {

  T = standardlocale.toDouble(valc, &ok);

  if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'T' from XML file.",valc); }

}

if (!strcmp(namec, "lambda_length")) {

  lambda_length = standardlocale.toDouble(valc, &ok);

  if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'lambda_length' from XML file.",valc); }

}

if (!strcmp(namec, "lambda_celllength")) {

  lambda_celllength = standardlocale.toDouble(valc, &ok);

  if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'lambda_celllength' from XML file.",valc); }

}

if (!strcmp(namec, "target_length")) {

  target_length = standardlocale.toDouble(valc, &ok);

  if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'target_length' from XML file.",valc); }

}

if (!strcmp(namec, "cell_expansion_rate")) {

  cell_expansion_rate = standardlocale.toDouble(valc, &ok);

  if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'cell_expansion_rate' from XML file.",valc); }

}

if (!strcmp(namec, "cell_div_expansion_rate")) {

  cell_div_expansion_rate = standardlocale.toDouble(valc, &ok);

  if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'cell_div_expansion_rate' from XML file.",valc); }

}

if (!strcmp(namec, "auxin_dependent_growth")) {

auxin_dependent_growth = strtobool(valc);

}

if (!strcmp(namec, "ode_accuracy")) {

  ode_accuracy = standardlocale.toDouble(valc, &ok);

  if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'ode_accuracy' from XML file.",valc); }

}

if (!strcmp(namec, "mc_stepsize")) {

 class Parameter {

 public: 

   Parameter();

   ~Parameter();

   void CleanUp(void);

   void Read(const char *filename);

   void Write(ostream &os) const;

   void XMLAdd(xmlNode *root) const;

   void XMLRead(xmlNode *root);

   void AssignValToPar(const char *namec, const char *valc);

   void AssignValArrayToPar(const char *namec, vector<double> valarray);

  char * arrowcolor;

  double arrowsize;

  char * textcolor;

  int cellnumsize;

  int nodenumsize;

  double node_mag;

  double outlinewidth;

  char * cell_outline_color;

  int resize_stride;

  double T;

  double lambda_length;

  double lambda_celllength;

  double target_length;

  double cell_expansion_rate;

  double cell_div_expansion_rate;

  bool auxin_dependent_growth;

  double ode_accuracy;

  double mc_stepsize;

  double mc_cell_stepsize;

  double energy_threshold;

  double bend_lambda;

  double alignment_lambda;

  double rel_cell_div_threshold;

  double rel_perimeter_stiffness;

  double collapse_node_threshold;

  double morphogen_div_threshold;

  double morphogen_expansion_threshold;

  bool copy_wall;

  double source;

  double * D;

  double * initval;

  double k1;

  double k2;

  double r;

  double kr;

#include "pardialog.h"

#include "parameter.h"

#include <cstring>

#include <qdialog.h>

#include <qlabel.h>

#include <qlineedit.h>

#include <qmessagebox.h>

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

ParameterDialog::ParameterDialog(QWidget *parent, const char *name, Qt::WindowFlags f) : QDialog(parent,name,false,f) {

    extern Parameter par;

  arrowcolor_edit = new QLineEdit( QString("%1").arg(par.arrowcolor), this, "arrowcolor_edit" );

  arrowsize_edit = new QLineEdit( QString("%1").arg(par.arrowsize), this, "arrowsize_edit" );

  textcolor_edit = new QLineEdit( QString("%1").arg(par.textcolor), this, "textcolor_edit" );

  cellnumsize_edit = new QLineEdit( QString("%1").arg(par.cellnumsize), this, "cellnumsize_edit" );

  nodenumsize_edit = new QLineEdit( QString("%1").arg(par.nodenumsize), this, "nodenumsize_edit" );

  node_mag_edit = new QLineEdit( QString("%1").arg(par.node_mag), this, "node_mag_edit" );

  outlinewidth_edit = new QLineEdit( QString("%1").arg(par.outlinewidth), this, "outlinewidth_edit" );

  cell_outline_color_edit = new QLineEdit( QString("%1").arg(par.cell_outline_color), this, "cell_outline_color_edit" );

  resize_stride_edit = new QLineEdit( QString("%1").arg(par.resize_stride), this, "resize_stride_edit" );

  T_edit = new QLineEdit( QString("%1").arg(par.T), this, "T_edit" );

  lambda_length_edit = new QLineEdit( QString("%1").arg(par.lambda_length), this, "lambda_length_edit" );

  lambda_celllength_edit = new QLineEdit( QString("%1").arg(par.lambda_celllength), this, "lambda_celllength_edit" );

  target_length_edit = new QLineEdit( QString("%1").arg(par.target_length), this, "target_length_edit" );

  cell_expansion_rate_edit = new QLineEdit( QString("%1").arg(par.cell_expansion_rate), this, "cell_expansion_rate_edit" );

  cell_div_expansion_rate_edit = new QLineEdit( QString("%1").arg(par.cell_div_expansion_rate), this, "cell_div_expansion_rate_edit" );

  auxin_dependent_growth_edit = new QLineEdit( QString("%1").arg(sbool(par.auxin_dependent_growth)), this, "auxin_dependent_growth_edit" );

  ode_accuracy_edit = new QLineEdit( QString("%1").arg(par.ode_accuracy), this, "ode_accuracy_edit" );

  mc_stepsize_edit = new QLineEdit( QString("%1").arg(par.mc_stepsize), this, "mc_stepsize_edit" );

  mc_cell_stepsize_edit = new QLineEdit( QString("%1").arg(par.mc_cell_stepsize), this, "mc_cell_stepsize_edit" );

  energy_threshold_edit = new QLineEdit( QString("%1").arg(par.energy_threshold), this, "energy_threshold_edit" );

  bend_lambda_edit = new QLineEdit( QString("%1").arg(par.bend_lambda), this, "bend_lambda_edit" );

  alignment_lambda_edit = new QLineEdit( QString("%1").arg(par.alignment_lambda), this, "alignment_lambda_edit" );

  rel_cell_div_threshold_edit = new QLineEdit( QString("%1").arg(par.rel_cell_div_threshold), this, "rel_cell_div_threshold_edit" );

  rel_perimeter_stiffness_edit = new QLineEdit( QString("%1").arg(par.rel_perimeter_stiffness), this, "rel_perimeter_stiffness_edit" );

  collapse_node_threshold_edit = new QLineEdit( QString("%1").arg(par.collapse_node_threshold), this, "collapse_node_threshold_edit" );

  morphogen_div_threshold_edit = new QLineEdit( QString("%1").arg(par.morphogen_div_threshold), this, "morphogen_div_threshold_edit" );

  morphogen_expansion_threshold_edit = new QLineEdit( QString("%1").arg(par.morphogen_expansion_threshold), this, "morphogen_expansion_threshold_edit" );

  copy_wall_edit = new QLineEdit( QString("%1").arg(sbool(par.copy_wall)), this, "copy_wall_edit" );

  source_edit = new QLineEdit( QString("%1").arg(par.source), this, "source_edit" );

  QString D_string("%1,%2,%3,%4,%5,%6,%7,%8,%9,%10,%11,%12,%13,%14,%15");

  D_string = D_string.arg(par.D[0]).arg(par.D[1]).arg(par.D[2]).arg(par.D[3]).arg(par.D[4]).arg(par.D[5]).arg(par.D[6]).arg(par.D[7]).arg(par.D[8]).arg(par.D[9]).arg(par.D[10]).arg(par.D[11]).arg(par.D[12]).arg(par.D[13]).arg(par.D[14]);

  D_edit = new QLineEdit( D_string, this, "D_edit" );

  QString initval_string("%1,%2,%3,%4,%5,%6,%7,%8,%9,%10,%11,%12,%13,%14,%15");

  initval_string = initval_string.arg(par.initval[0]).arg(par.initval[1]).arg(par.initval[2]).arg(par.initval[3]).arg(par.initval[4]).arg(par.initval[5]).arg(par.initval[6]).arg(par.initval[7]).arg(par.initval[8]).arg(par.initval[9]).arg(par.initval[10]).arg(par.initval[11]).arg(par.initval[12]).arg(par.initval[13]).arg(par.initval[14]);

  initval_edit = new QLineEdit( initval_string, this, "initval_edit" );

  grid->addWidget( new QLabel( "arrowcolor", this ),4, 0 );

  grid->addWidget( arrowcolor_edit, 4, 0+1  );

  grid->addWidget( new QLabel( "arrowsize", this ),5, 0 );

  grid->addWidget( arrowsize_edit, 5, 0+1  );

  grid->addWidget( new QLabel( "textcolor", this ),6, 0 );

  grid->addWidget( textcolor_edit, 6, 0+1  );

  grid->addWidget( new QLabel( "cellnumsize", this ),7, 0 );

  grid->addWidget( cellnumsize_edit, 7, 0+1  );

  grid->addWidget( new QLabel( "nodenumsize", this ),8, 0 );

  grid->addWidget( nodenumsize_edit, 8, 0+1  );

  grid->addWidget( new QLabel( "node_mag", this ),9, 0 );

  grid->addWidget( node_mag_edit, 9, 0+1  );

  grid->addWidget( new QLabel( "outlinewidth", this ),10, 0 );

  grid->addWidget( outlinewidth_edit, 10, 0+1  );

  grid->addWidget( new QLabel( "cell_outline_color", this ),11, 0 );

  grid->addWidget( cell_outline_color_edit, 11, 0+1  );

  grid->addWidget( new QLabel( "resize_stride", this ),12, 0 );

  grid->addWidget( resize_stride_edit, 12, 0+1  );

  grid->addWidget( new QLabel( "", this), 13, 0, 1, 2 );

  grid->addWidget( new QLabel( " <b>Cell mechanics</b>", this), 14, 0, 1, 2 );

  grid->addWidget( new QLabel( "T", this ),15, 0 );

  grid->addWidget( T_edit, 15, 0+1  );

  grid->addWidget( new QLabel( "lambda_length", this ),16, 0 );

  grid->addWidget( lambda_length_edit, 16, 0+1  );