//------------------------------------------------------------------------------------------------
//------------------------------------------------------------------------------------------------
void Field::selectConstituent()
{
int iconst, res;
QStringList items;
LOG_MSG("selectConstituent");
get_fieldinfo(&NX, &axis, &fraction, &nsites, &nconst, const_used, &res);
for (iconst=0; iconst<MAX_CONC+1; iconst++) {
if (iconst == constituent) continue;
if (const_used[iconst] == 1) {
items << const_name[iconst];
}
}
bool ok;
QString item = QInputDialog::getItem(this, tr("QInputDialog::getItem()"),
tr("Constituent:"), items, 0, false, &ok);
if (ok && !item.isEmpty()) {
for (iconst=0; iconst<MAX_CONC+1; iconst++) {
if (item == const_name[iconst]) {
constituent = iconst;
if (useConcPlot)
updateConcPlot();
break;
}
}
}
}
static LLVMValueRef make_fieldptr(compile_t* c, LLVMValueRef l_value,
ast_t* l_type, ast_t* right)
{
pony_assert(ast_id(l_type) == TK_NOMINAL);
pony_assert(ast_id(right) == TK_ID);
ast_t* def;
ast_t* field;
uint32_t index;
get_fieldinfo(l_type, right, &def, &field, &index);
if(ast_id(def) != TK_STRUCT)
index++;
if(ast_id(def) == TK_ACTOR)
index++;
return LLVMBuildStructGEP(c->builder, l_value, index, "");
}
//-----------------------------------------------------------------------------------------
// New version, site/cell size is fixed, the blob grows
//-----------------------------------------------------------------------------------------
void Field::displayField(int hr, int *res)
{
QGraphicsScene* scene = new QGraphicsScene(QRect(0, 0, CANVAS_WIDTH, CANVAS_WIDTH));
QBrush brush;
int i, xindex, yindex, ix, iy, w, rgbcol[3];
double xp, yp, d0, d, volume, scale, cmin, cmax, rmax;
double a, b, Wc;
int Nc;
bool growthRate;
// LOG_MSG("displayField");
*res = 0;
hour = hr;
if (slice_changed) {
get_fieldinfo(&NX, &axis, &fraction, &nsites, &nconst, const_used, res);
if (*res != 0) return;
this->data = (FIELD_DATA *)malloc(nsites*sizeof(FIELD_DATA));
get_fielddata(&axis, &fraction, &nsites, &nconst, this->data, res);
if (*res != 0) return;
slice_changed = false;
}
if (constituent == GROWTH_RATE)
growthRate = true;
else
growthRate = false;
if (axis == X_AXIS) { // Y-Z plane
xindex = 1;
yindex = 2;
} else if (axis == Y_AXIS) { // X-Z plane
xindex = 0;
yindex = 2;
} else if (axis == Z_AXIS) { // X-Y plane
xindex = 0;
yindex = 1;
}
/*
NX = size of lattice
Nc = # of sites to fill the canvas from side to side (or top to bottom) = (2/3)NX
Wc = canvas width (pixels)
w = site width = Wc/Nc
xp = a.ix + b
yp = a.iy + b
blob centre at (NX/2,NX/2) maps to canvas centre at (Wc/2,Wc/2)
=> Wc/2 = a.NX/2 + b
The width of Nc sites maps to the canvas width
=> Wc = a.Nc
=> a = Wc/Nc, b = Wc/2 - a.NX/2
*/
Nc = (2*NX)/3;
Wc = CANVAS_WIDTH;
w = Wc/Nc;
a = w;
b = Wc/2 - a*NX/2;
d0 = w*dfraction;
cmin = 1.0e10;
cmax = 0;
rmax = 0;
for (i=0; i<nsites; i++) {
rmax = MAX(rmax,data[i].dVdt);
cmin = MIN(MAX(cmin,0),data[i].conc[constituent]);
cmax = MAX(cmax,data[i].conc[constituent]);
}
brush.setStyle(Qt::SolidPattern);
brush.setColor(QColor(0,0,0));
scene->addRect(0,0,CANVAS_WIDTH,CANVAS_WIDTH,Qt::NoPen, brush);
view->setScene(scene);
view->setGeometry(QRect(0, 0, 700, 700));
if (cmax == 0) {
view->show();
return;
}
for (i=0; i<nsites; i++) {
ix = this->data[i].site[xindex];
iy = this->data[i].site[yindex];
xp = int(a*ix + b - w);
yp = int(a*iy + b - w);
chooseFieldColor(data[i].conc[constituent],cmin,cmax,use_log,rgbcol);
brush.setColor(QColor(rgbcol[0],rgbcol[1],rgbcol[2]));
scene->addRect(xp,yp,w,w,Qt::NoPen, brush);
}
for (i=0; i<nsites; i++) {
ix = this->data[i].site[xindex];
iy = this->data[i].site[yindex];
xp = int(a*ix + b - w);
yp = int(a*iy + b - w);
volume = this->data[i].volume; // = 0 if there is no cell
if (volume > 0) {
scale = pow(volume,0.3333);
d = scale*d0; // fix this - need to change d0
double f;
if (rmax > 0) {
f = data[i].dVdt/rmax;
} else {
f = 0.01;
}
chooseRateColor(f,rgbcol);
brush.setColor(QColor(rgbcol[0],rgbcol[1],rgbcol[2]));
scene->addEllipse(xp+(w-d)/2,yp+(w-d)/2,d,d,Qt::NoPen, brush);
}
}
view->show();
if (save_images) {
scene->clearSelection(); // Selections would also render to the file
scene->setSceneRect(scene->itemsBoundingRect()); // Re-shrink the scene to it's bounding contents
QImage image(scene->sceneRect().size().toSize(), QImage::Format_ARGB32); // Create the image with the exact size of the shrunk scene
image.fill(Qt::transparent); // Start all pixels transparent
QPainter painter(&image);
scene->render(&painter);
ifield++;
char filename[] = "image/field0000.png";
char numstr[5];
sprintf(numstr,"%04d",hour);
for (int i=0; i<4; i++)
filename[11+i] = numstr[i];
image.save(filename);
}
}